Project

General

Profile

unexpected hybrid coordinate formula modification

Added by Andres Luhamaa over 7 years ago

Hi
I tried to interpolate NWP model file from LCC to PlateCarree and it seems to do what I want:

cdo remapbil,targetgrid infile outfile

but as a side effect, cdo has also changed representation of vertical coordinates and IMO incorrectly.

In the original file, there are dimensions "hybrid1 = 1" and "p0 = 1", which have been transformed to "lev = 1" and "bnds = 2", so that variables that were given at "hybrid1" are now given at "lev". The problem is that the transformation is not done correctly. In the original file, the hybrid1 is described as:

double hybrid1(hybrid1) ;
hybrid1:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
hybrid1:formula = "p(n,k,j,i) = ap(k) + b(k)*ps(n,j,i)" ;
hybrid1:formula_terms = "ap: ap1 b: b1 ps: surface_air_pressure p0: p01" ;

where the formula_terms refers to actually existing variables, but in the interpolated file there is:

double lev(lev) ;
lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
lev:long_name = "hybrid sigma pressure coordinate" ;
lev:formula = "p = ap + b*ps" ;
lev:formula_terms = "ap: ap b: b ps: " ;

and there are no such variables as "ap" or "b" in the file and my final application actually tries to interpret this formula and obviously fails.

Is there a way to turn this unexpected conversion of vertical coordinate off or fix it somewhere?