[MITgcm-devel] something to report

[Manfredi Manizza]

Hi,

I am reporting a feture I found comparing the output
of the ECCO lat-lon-cap global set-up provided by Patrick to me.

The dimensions of the files are different depending on the format.
If in netcdf :

S =

          attributes: [1x1 struct]
     iters_from_file: [12x1 double]
         si_AREAtave: [90x1082x12 double]

If in binary :


> more pickup_seaice.0000008772.meta
  nDims = [   2 ];
  dimList = [
     90,    1,   90,
   1170,    1, 1170
  ];
  dataprec = [ 'float64' ];
  nrecords = [     6 ];
  timeStepNumber = [       8772 ];
  timeInterval = [  3.157920000000E+07 ];
  nFlds = [    6 ];
  fldList = {
  'siTICE  ' 'siAREA  ' 'siHEFF  ' 'siHSNOW ' 'siUICE  ' 'siVICE  '
  };


The grid quantities that I have , when loaded, have the dimensions
of 90 by 1170.

In fact i can use the functions of gcmfaces only using the binary 
files...

I am not sure if there is something I am missing or if
there is something wrong in the code. I just felt I needed to
send this out.

Thanks in advance for your attention.

Cheers
Manfredi