[MITgcm-devel] diagstats_mnc_out.F

[Ed Hill]

On Wed, 2005-07-13 at 20:43 -0400, Jean-Michel Campin wrote:
> 
> I have a few remarks about the mnc version of the diagStat output file:

Hi Jean-Michel,

For items 0--2 you've checked in fixes so they're all set.

> 3) When I do ncdump, I would prefer to see, e.g. ETAN_lv_vol
>  with 2 dimensions (Time, Region) instead of 3 (Time, region, Zd000001).
> <- I don't know how to do this.

OK, thats not too hard to do -- I'll take care of it.

> 4) vertical coordinate would be useful.
> <- I saw the "DIAGST_MNC_NEEDSWORK"

Yes, as soon as we have two vectors for the gridalt grids that are
equivalent to rC, rF then I can finish it.  This means implementing the
one-location vertical sum of dpphys which I don't completely understand.
I'll talk to Andrea.

> 5) I think it would be nice to put the right time (computed from
>  myTime, frequency & phase) to correspond to the midle of the
>  time period (if averaged) or the time of the snap-shot.
> <- not urgent

Yes.

> 6) would better not put any "units" for _vol field, since
>    "udiag" is not the right units for this variable
>   (But it's hard to know exactly what they are).
>   And the minimum, maximum, average & standard-deviation
>   are missing in the description.

OK, I'll look into it.

> 7) I need only 1 output file (and I get only 1 ascii file) but
>  get plenty of .nc files (1 per tile), all identical except bi,bj
>  in the header.
> <- This is anoying and needs to be fixed.
>    And since there is no "tile" information, we can drop the
>    tile number in the file-name suffix.

Yes, this is easy to fix.

Ed

-- 
Edward H. Hill III, PhD
office:  MIT Dept. of EAPS;  Rm 54-1424;  77 Massachusetts Ave.
             Cambridge, MA 02139-4307
emails:  eh3 at mit.edu                ed at eh3.com
URLs:    http://web.mit.edu/eh3/    http://eh3.com/
phone:   617-253-0098
fax:     617-253-4464