[MITgcm-devel] Problems with Z Interface-level output

[Baylor Fox-Kemper]

Hi Ed et al,
  There's still a problem with output at non-centered Z levels.  Note 
the problems below:

1) Unlike Xp1 and Yp1, the Z coordinate is named Zl
2) The coordinate variable Zl is not filled with any values, only index 
number.

Thanks,
    -Baylor

public1:/archive/bfk/MITgcm.runs> ncdump -h 
fr.4.0/fr.4.0_0001/state.0000.all.nc
netcdf state.0000.all {
dimensions:
         T = UNLIMITED ; // (101 currently)
         X = 150 ;
         Xp1 = 151 ;
         Y = 51 ;
         Yp1 = 52 ;
         Z = 38 ;
         Zl = 38 ;
variables:
         int T(T) ;
                 T:long_name = "iteration_count" ;
         double model_time(T) ;
                 model_time:long_name = "Model Time" ;
                 model_time:units = "s" ;
         double X(X) ;
         double Xp1(Xp1) ;
         double Y(Y) ;
         double Yp1(Yp1) ;
         double Z(Z) ;
         double U(T, Z, Y, Xp1) ;
                 U:units = "m/s" ;
                 U:coordinates = "XU YU RC iter" ;
         double V(T, Z, Yp1, X) ;
                 V:units = "m/s" ;
                 V:coordinates = "XV YV RC iter" ;
         double W(T, Zl, Y, X) ;
                 W:units = "m/s" ;
                 W:coordinates = "XC YC RC iter" ;
         double Temp(T, Z, Y, X) ;
                 Temp:units = "degC" ;
                 Temp:long_name = "potential_temperature" ;
                 Temp:coordinates = "XC YC RC iter" ;
         double S(T, Z, Y, X) ;
                 S:long_name = "salinity" ;
                 S:coordinates = "XC YC RC iter" ;
         double Eta(T, Y, X) ;
                 Eta:long_name = "free-surface_r-anomaly" ;
                 Eta:units = "m" ;
                 Eta:coordinates = "XC YC iter" ;
}