Problem using CFSR in grb2 format
Added by Alessio Spassiani almost 11 years ago
Hi,
I am attempting to remap CFSR data in grb2 format using the command:
cdo remap,r144x73,remapweights_cfsr.nc cfsr.$year$month$day.grb2 remap.cfsr.$year$month$day.grb2
however cdo doesn't map all the variables, specifically the v-component of wind (VGRDprs) which seems odd since it is able to remap the u-component of wind (UGRDprs).
My data file is too large to attached, but it was downloaded from the CISL Research Data Archive. Below is the control file for the data I am working with (before remapping).
dset ^./cfsr.%y4%m2%d2.grb2 index ^./cfsr.20080703.grb2.idx undef 9.999E+20 title ./cfsr.20080703.grb2 * produced by g2ctl v0.0.8.8 * command line options: -0 ./cfsr.%y4%m2%d2.grb2 * griddef=1:0:(720 x 361):grid_template=0:winds(N/S): lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.500000 lon 0.000000 to 359.500000 by 0.500000 #points=259920:winds(N/S) dtype grib2 options template ydef 361 linear -90.000000 0.5 xdef 720 linear 0.000000 0.500000 tdef 3449 linear 00Z03jul2008 6hr * PROFILE hPa zdef 37 levels 100000 97500 95000 92500 90000 87500 85000 82500 80000 77500 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 22500 20000 17500 15000 12500 10000 7000 5000 3000 2000 1000 700 500 300 200 100 options pascals vars 6 HGTprs 37,100 0,3,5 ** (1000 975 950 925 900.. 7 5 3 2 1) Geopotential Height [gpm] PRESmsl 0,101,0 0,3,0 ** mean sea level Pressure [Pa] PRMSLmsl 0,101,0 0,3,1 ** mean sea level Pressure Reduced to MSL [Pa] TMPprs 37,100 0,0,0 ** (1000 975 950 925 900.. 7 5 3 2 1) Temperature [K] UGRDprs 37,100 0,2,2 ** (1000 975 950 925 900.. 7 5 3 2 1) U-Component of Wind [m/s] VGRDprs 37,100 0,2,3 ** (1000 975 950 925 900.. 7 5 3 2 1) V-Component of Wind [m/s] ENDVARS
It appears cdo doesn't not recognize the v-component of the wind. When I use the “sinfon” command I get the following:
cdo sinfon cfsr.20080703.grb2 File format: GRIB2 JPEG -1 : Institut Source Ttype Levels Num Gridsize Num Dtype : Parameter name 1 : NCEP unknown instant 1 1 259920 1 P16z : msl 2 : NCEP unknown instant 37 2 259920 1 P16z : gh 3 : NCEP unknown instant 37 2 259920 1 P10z : t 4 : NCEP unknown instant 37 2 259920 1 P14z : u 5 : NCEP unknown instant 1 1 259920 1 P17z : prmsl Grid coordinates : 1 : lonlat > size : dim = 259920 nx = 720 ny = 361 lon : first = 0 last = 359.5 inc = 0.5 degrees_east circular lat : first = 90 last = -90 inc = -0.5 degrees_north Vertical coordinates : 1 : meansea level : 0 2 : pressure Pa : 100 200 300 500 700 1000 2000 3000 5000 7000 10000 12500 15000 17500 20000 22500 25000 30000 35000 40000 45000 50000 55000 60000 65000 70000 75000 77500 80000 82500 85000 87500 90000 92500 95000 97500 100000 Time coordinate : unlimited steps RefTime = 2008-07-03 00:00:00 Units = hours Calendar = proleptic_gregorian YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss 2008-07-03 00:00:00 2008-07-03 06:00:00 2008-07-03 12:00:00 2008-07-03 18:00:00 cdo sinfon: Processed 5 variables over 4 timesteps ( 0.41s )
Is there a reason why CDO will recognize the u-component of the wind but not the v-component. Is there a possible solution?
Below is the current version of CDO I am using.
Climate Data Operators version 1.6.3rc3 (http://code.zmaw.de/projects/cdo) Compiler: gcc -std=gnu99 -g -O2 -pthread version: gcc (GCC) 4.4.7 20120313 (Red Hat 4.4.7-3) Compiled: by localuser on tiger.meteo.mcgill.ca (x86_64-unknown-linux-gnu) Jan 23 2014 12:12:26 Features: PTHREADS NC4 Z JASPER UDUNITS2 PROJ.4 XML2 CURL Libraries: proj/4.7 xml2/2.7.6 curl/7.19.7 Filetypes: srv ext ieg grb grb2 nc nc2 nc4 nc4c CDI library version : 1.6.3 of Jan 23 2014 12:12:16 CGRIBEX library version : 1.6.3 of Jan 8 2014 19:55:18 GRIB_API library version : 1.9.16 netCDF library version : 4.1.1 of Jan 30 2014 05:31:22 $ HDF5 library version : 1.8.5 SERVICE library version : 1.3.1 of Jan 23 2014 12:11:57 EXTRA library version : 1.3.1 of Jan 23 2014 12:11:51 IEG library version : 1.3.1 of Jan 23 2014 12:11:55 FILE library version : 1.8.2 of Jan 23 2014 12:11:51
Thanks in advance for the help.
Sincerely,
Alessio
Replies (1)
RE: Problem using CFSR in grb2 format - Added by Uwe Schulzweida almost 11 years ago
Hi Alessio,
CDO does not support the GRIB2 Multi-field feature. Multi-field is a GRIB2 record with more than one field. u and v are often stored in a Multi-field record. CDO reads only the first field of a Multi-field record.
Cheers,
Uwe