[MITgcm-support] Initial Conditions on Staggered Grid

Martin Losch Martin.Losch at awi.de
Wed Jul 23 06:32:34 EDT 2014


Hi Arjumand,

The model expects nx,ny grid points also for the u and v points. It will fill nx+1 or ny+1 by itself. So you should only save u_sta(1:256,1:256,:) and v_sta(1:256,1:256,:)

Martin

On Jul 23, 2014, at 10:29 AM, Muhammad Arjumand Masood <masood at g.harvard.edu> wrote:

> Hi,
> 
> I am having some trouble getting the gcm to the read in my U,V initial conditions when I specify them according to a staggered grid. To be specific, U is read in incorrectly, whereas V is read in as I intend. 
> 
> The components of flow are staggered in space in the form of an Arakawa C grid. (http://mitgcm.org/sealion/online_documents/node45.html). When I specify the U,V initial conditions, in an un-staggered Arakawa A-grid from, the variables are read in basically correctly. Since the gcm is eventually outputting U,V in the Arakawa C grid form, shouldn’t I be able to specify exactly how I want the initial U,V to be read in? 
> 
> Maybe this is an issue in the way I am writing my .bin file for the initial conditions because I don’t know exactly how the .bin files are read in by the gcm. I create the initial conditions in matlab and get: 
> 
> >> size(u_sta)
>    257   256    50
> 
> >> size(v_sta)
>    256   257    50
> 
> I’ve checked what u_sta and v_sta look like in matlab, and they are both fine. 
> 
> My command for creating the bin files is simply: 
> 
> fid=fopen('MyU_sta.bin','w','b'); fwrite(fid,u_sta,'real*4'); fclose(fid);
> fid=fopen('MyV_sta.bin','w','b'); fwrite(fid,v_sta,'real*4'); fclose(fid);
> 
> Any help would be appreciated!
> 
> Thanks,
> Arjumand
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list