[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