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

[Jody Klymak]

Hi Martin,

On Feb 7, 2013, at  2:15 AM, Martin Losch <Martin.Losch at awi.de> wrote:
> 
> 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.

OK, that looks reasonable, though I may have to do both 3-D and 2-D: I'll want to save occasional 3-D snapshots....

> 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.

Maybe this is the best approach - certainly easier to program.  My problem is that once I have things on disk they tend to stay there ;-)

Thanks,  Jody

> 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
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

--
Jody Klymak    
http://web.uvic.ca/~jklymak/