[MITgcm-support] Advice on writing 2-D vertical sections from 3-D model

[Martin Losch]

Hi Jody

as far as I know there is no generic code for writing slabs in a 
convenient (and efficient) way. [There is a very inefficient way using 
the diagStats machinery, but that will be very slow, basically each of 
your grid points in your xz-slab is a separate region.]

For snap-shots of U,V,W, and T, I would call WRITE_REC_XZ_RL in 
model/src/write_state.F and comment out writing the 3D fields, and then 
you can use the regular dumpFreq to control the save intervals.
I think you can use the existing code as templates.

When I am trying to subsample large output (e.g. just a vertical slice 
or a couple of grid points), I usually write the full 3D arrays 
temporarily to scratch-disk and do a post processing step in a separate 
program that reads the 3D file(s), cuts out and saves what I want and 
then deletes the 3D file(s). That works fairly well in "chain-jobs" 
where I have to restart the model anyway.

Martin

On 02/06/2013 11:15 PM, Jody Klymak wrote:
>
> Hi All,
>
> Whats the easiest way to save 2-D vertical sections of data from 3-D simulations?
>
> I'm running some 3-D simulations, and I really don't care about saving all the information in the y-direction, but would like to save xz-slices of U,V,W and T (snapshots will be fine).  I didn't see that diagnostics would do that for me out of the box.  Any other suggestions?  I can see that there is WRITE_REC_XZ_RL, so that'll do what I want - is there a good place to call it and specify the save intervals?
>
> Thanks,  Jody
>
> --
> Jody Klymak
> http://web.uvic.ca/~jklymak/
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>

-- 
Martin Losch
Alfred Wegener Institute for Polar and Marine Research
Postfach 120161, 27515 Bremerhaven, Germany;
Tel./Fax: ++49(0471)4831-1872/1797