[MITgcm-support] Regional output - diagstats

[Martin Losch]

Hi Holger,

the diagstats can do "statistics" (mean, min/max, std) of sub-regions, but not really 2 or 3D fields of these sub-regions. I suppose you want to do the latter. I don't think that there is code for doing that. For that it's probably easier to actually modify the low level I/O routines for your specific requirements (although I have to admit that MDS_WRITE_FIELD (in mdsio_write_field.F) is not my favorite routine.

M.

On Sep 13, 2012, at 12:51 AM, Holger Brix wrote:

> Hi,
> 
> I am trying to figure out a way to save output with high temporal resolution for a small sub-region of my global setup. Looking at the code it seems to me that diagstats is the place to do this. I was wondering if this is right and if someone could share an example for how to do this.
> 
> Thanks a lot,
> Holger
> -- 
> __________________________________________
> 
> Holger Brix
> Dept. of Atmospheric and Oceanic Sciences
> University of California Los Angeles
> 5845 Slichter Hall
> Los Angeles, CA 90095-1567, USA
> phone UCLA: (310) 825 4526
> phone JPL : (828) 393 1018
> fax (UCLA): (310) 206 3051
> email: hbrix at ucla.edu
> http://www.atmos.ucla.edu/~hbrix
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support