[MITgcm-support] Output 2D slice from 3D simulation

[Martin Losch]

Xiaozhou,

as far as I know, the diagnostics package can write 3D fields and 2D horizontal slices of 3D fields. I don’t think that vertical slices are possible.

If you need 2D slices you’d probably have to write your own code for this. There are subroutines in pkg/rw/write_rec.F that do vertical slices, but you need to call them in the right place. Here: pkg/obcs/obcs_output.F there are examples how to use them.

Martin

On 21. May 2018, at 23:12, Ruan Xiaozhou <saberruan at hotmail.com> wrote:
> 
> Dear all,
> 
> I’m posting this again in case someone has a solution.. I have been trying to output 2D slice (YZ) data from a 3D simulation (because the 3D field is too large to be outputted at high frequency) but did not find proper locations in the diagnostics package. I’d really appreciate it if you could point me to the right directions. 
> 
> Best,
> Xiaozhou
> --
> Xiaozhou Ruan
> Division of Geological and Planetary Sciences
> California Institute of Technology
> MC 131-24
> Pasadena, CA  91125
> 
> phone: +1 626-395-8715
> email: xiaozhour at caltech.edu
> web:   http://www.its.caltech.edu/~xruan
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support