Issue with ml2pl
Added by James King over 4 years ago
Hi all,
I have some data on hybrid verical levels which I would like to interpolate to pressure levels. When I use
cdo ml2pl,100000,97500,95000...
I get the error
cdo ml2pl (Warning): No 3D variable with hybrid sigma pressure coordinate found!
Any ideas would be most welcome. I can't upload the file because it's too large, but cdo sinfon returns the following information:
File format : netCDF2
-1 : Institut Source Ttype Levels Num Points Num Dtype : Parameter name
1 : unknown unknown instant 60 1 204800 1 F32 : HGT_GDS4_HYBL_S123
Grid coordinates :
1 : gaussian : points=204800 (640x320) np=0
g4_lon_3 : 0 to 359.438 by 0.562501 degrees_east circular
g4_lat_2 : 89.5701 to -89.5701 degrees_north
Vertical coordinates :
1 : generic : levels=60
lv_HYBL1 : 1 to 60 by 1 number
Time coordinate : 396 stepsJames
Replies (4)
RE: Issue with ml2pl - Added by Ralf Mueller over 4 years ago
you can upload the first timestep - this should be small enough for an upload
RE: Issue with ml2pl - Added by Ralf Mueller over 4 years ago
Hi James!
netcdf jra55_zg_1981-2013_hybrid_time1 {
dimensions:
g4_lon_3 = 640 ;
g4_lat_2 = 320 ;
lv_HYBL1 = 60 ;
initial_time0_hours = UNLIMITED ; // (1 currently)
variables:
float g4_lon_3(g4_lon_3) ;
g4_lon_3:standard_name = "longitude" ;
g4_lon_3:long_name = "longitude" ;
g4_lon_3:units = "degrees_east" ;
g4_lon_3:axis = "X" ;
float g4_lat_2(g4_lat_2) ;
g4_lat_2:standard_name = "latitude" ;
g4_lat_2:long_name = "latitude" ;
g4_lat_2:units = "degrees_north" ;
g4_lat_2:axis = "Y" ;
double lv_HYBL1(lv_HYBL1) ;
lv_HYBL1:long_name = "hybrid level" ;
lv_HYBL1:units = "number" ;
lv_HYBL1:axis = "Z" ;
double initial_time0_hours(initial_time0_hours) ;
initial_time0_hours:standard_name = "time" ;
initial_time0_hours:long_name = "initial time" ;
initial_time0_hours:units = "hours since 1800-1-1 00:00:00" ;
initial_time0_hours:calendar = "standard" ;
initial_time0_hours:axis = "T" ;
float HGT_GDS4_HYBL_S123(initial_time0_hours, lv_HYBL1, g4_lat_2, g4_lon_3) ;
HGT_GDS4_HYBL_S123:long_name = "Geopotential height" ;
HGT_GDS4_HYBL_S123:units = "gpm" ;
HGT_GDS4_HYBL_S123:grid_type = "gaussian" ;
HGT_GDS4_HYBL_S123:_FillValue = 1.e+20f ;
HGT_GDS4_HYBL_S123:missing_value = 1.e+20f ;
HGT_GDS4_HYBL_S123:N = 124, 112, 124, 120, 124, 120, 124, 124, 120, 124, 120, 124 ;
HGT_GDS4_HYBL_S123:statistical_process_duration = "instantaneous (beginning at reference time at intervals of 6 hours)" ;
HGT_GDS4_HYBL_S123:statistical_process_descriptor = "average of N uninitialized analyses" ;
HGT_GDS4_HYBL_S123:forecast_time_units = "hours" ;
HGT_GDS4_HYBL_S123:forecast_time = 0 ;
HGT_GDS4_HYBL_S123:parameter_number = 7 ;
HGT_GDS4_HYBL_S123:parameter_table_version = 200 ;
HGT_GDS4_HYBL_S123:gds_grid_type = 4 ;
HGT_GDS4_HYBL_S123:level_indicator = 109 ;
HGT_GDS4_HYBL_S123:center = "Japanese Meteorological Agency - Tokyo (RSMC)" ;
HGT_GDS4_HYBL_S123:sub_center = "241" ;
// global attributes:
:CDI = "Climate Data Interface version ?? (http://mpimet.mpg.de/cdi)" ;
:Conventions = "None" ;
:history = "Tue Sep 24 16:07:52 2019: cdo -selmon,1 -selyear,1981 jra55_zg_1981-2013_hybrid.nc jra55_zg_1981-2013_hybrid_time1.nc\n",
"Tue Sep 24 14:21:33 2019: cdo selname,HGT_GDS4_HYBL_S123 jra55_zg_1981-2013_hybrid.nc jra55_zg_1981-2013_hybrid_1.nc\n",
"Tue Sep 24 13:53:59 2019: cdo selyear,1981/2013 jra55_zg_1958-2013_hybrid.nc jra55_zg_1981-2013_hybrid.nc\n",
"Tue Sep 24 13:39:10 2019: cdo mergetime /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.195801_195812.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.195901_195912.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196001_196012.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196101_196112.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196201_196212.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196301_196312.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196401_196412.engelstaedter383169.nc /soge-home/data/analysis/jra-55/TL319L60/monthly/hgt/nc/anl_mdl.007_hgt.reg_tl319.196501_196512.engelstaedter383169.nc" ;
:creation_date = "Thu Aug 8 04:58:34 MDT 2019" ;
:NCL_Version = "6.6.2" ;
:system = "Linux casper10 3.10.0-693.21.1.el7.x86_64 #1 SMP Wed Mar 7 19:03:37 UTC 2018 x86_64 x86_64 x86_64 GNU/Linux" ;
:grib_source = "anl_mdl.007_hgt.reg_tl319.195801_195812.engelstaedter383169.grb" ;
:title = "NCL: convert-GRIB-to-netCDF" ;
:CDO = "Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo)" ;
data:
I cannot see any information about the vertical coordinate except that is has 60 levels. no coefficients, no formular. CDO cannot make anything with it, so it declares the z-axis as generic. You might check the CF-convention on how to declare it in your case.
BTW: CDO can handle GRIB right away - no need to convert it and potentially loose information.
hth
ralf
RE: Issue with ml2pl - Added by James King over 4 years ago
Hi Ralf,
Thnaks for that. Maybe I'm chasing a lost cause here.