[MITgcm-support] 2d and 3d-diagnostics in one netcdf file?

[Martin Losch]

Hi all,

I am embarrassed to ask because it may have been discussed already  
and I just cannot find it:
1. Time-averaging with the diagnostics package: the documentation  
says that frequency > 0 produces averages, frequency < 0 produces  
snapshots. What does averagingFeq do?

2. Is it possible in the diagnostics package to have a 2D and 3D  
field in the same output stream, e.g.,

  frequency(5) = 43200.,
filename(5)  = 'state_diags',
fields(1,5)  = 'ETAN','THETA','SALT','UVEL','VVEL','WVEL',

? With these parameters I seem to be getting the level1 of ETAN  
(which is correct), but then I also get only the first levels for  
THETA, SALT, etc. here's the header of the netcdf file, note, that  
both Zmd000001 and Zld000001 = 1. Is that to be expected?

> > ncdump -h state_diags.0000072000.t001.nc
> netcdf state_diags.0000072000.t001 {
> dimensions:
>         T = UNLIMITED ; // (5 currently)
>         Zmd000001 = 1 ;
>         X = 180 ;
>         Y = 108 ;
>         Zd000001 = 1 ;
>         Xp1 = 181 ;
>         Yp1 = 109 ;
>         Zld000001 = 1 ;
> variables:
>         double T(T) ;
>                 T:long_name = "model_time" ;
>                 T:units = "s" ;
>         int iter(T) ;
>                 iter:long_name = "iteration_count" ;
>         double diag_levels(Zmd000001) ;
>                 diag_levels:description = "Idicies of vertical  
> levels within the source arrays" ;
>         double X(X) ;
>         double Y(Y) ;
>         float ETAN(T, Zd000001, Y, X) ;
>                 ETAN:description = "Surface Height Anomaly" ;
>                 ETAN:units = "m" ;
>         float THETA(T, Zmd000001, Y, X) ;
>                 THETA:description = "Potential Temperature" ;
>                 THETA:units = "degC" ;
>         float SALT(T, Zmd000001, Y, X) ;
etc.

Thanks,

Martin