Project

General

Profile

Problem using CFSR in grb2 format

Added by Alessio Spassiani almost 10 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 10 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

    (1-1/1)