Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

(RESOLVED) ERA5 from CDS - g2print does not properly read timesteps ?

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

K

Klemens

New member
Dear all,

I try to run WRF with ERA5 data on model levels downloaded from CDS.
The download script for one variable is below:
Code: 
  c.retrieve("reanalysis-era5-complete",
  {
  "class": "ea",
  "expver": "1",
  "stream": "oper",
  "type": "an",
  "param"   : parameter_str+".128",
  "levtype": "ml",
  "levelist": "1/to/137",
  "date": date,
  "time": "00/to/23/by/1",
  "format": "grib"
   },
   output_name)

Data should be in grib2 with hourly data for one day.
But when I run g2print.exe on it I get only data for the first timestep:

Code: 
tauruslogin6:~/WRF3.8.1/WPS/ungrib> ./g2print.exe /scratch/p_convection/ERA5/ERA5_ml_129_20060707_rec.grb
 ungrib - grib edition num           2
 reading from grib file =
 /scratch/p_convection/ERA5/ERA5_ml_129_20060707_rec.grb

      ECMWF
 -------------------------------------------------------------------------------
 --------
  rec Prod Cat Param  Lvl    Lvl      Lvl     Prod    Name            Time
     Fcst
  num Disc     num    code   one      two     Templ
     hour
 -------------------------------------------------------------------------------
 --------
 Gaussian Grid: Dx,Dy,lat,lon  0.1875000       480.0000       89.85654
  0.0000000E+00
   1   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   2   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   3   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   4   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   5   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   6   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   7   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   8   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
   9   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  10   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  11   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  12   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  13   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  14   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  15   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  16   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  17   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  18   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  19   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  20   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  21   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  22   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  23   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00
  24   0    3   4     105       1       0       0     GP       2006-07-07_00:00:00   00


   Successful completion of g2print

When I convert the file with CDOs and make an ncdump, I get the correct number of timesteps:

Code: 
// global attributes:
                :CDI = "Climate Data Interface version 1.9.4 (http://mpimet.mpg.de/cdi)" ;
                :Conventions = "CF-1.6" ;
                :history = "Fri Mar 15 10:27:25 2019: cdo -f nc copy ERA5_ml_129_20060707_rec.grb ERA5_ml_129_20060707_rec.nc" ;
                :institution = "European Centre for Medium-Range Weather Forecasts" ;
                :CDO = "Climate Data Operators version 1.9.4 (http://mpimet.mpg.de/cdo)" ;
data:

 time = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19,
    20, 21, 22, 23 ;

 lon = 0, 0.1875, 0.375, 0.5625, 0.75, 0.9375, 1.125, 1.3125, 1.5, 1.6875,
    1.875, 2.0625, 2.25, 2.4375, 2.625, 2.8125, 3, 3.1875, 3.375, 3.5625,
    3.75, 3.9375, 4.125, 4.3125, 4.5, 4.6875, 4.875, 5.0625, 5.25, 5.4375,
    5.625, 5.8125, 6, 6.1875, 6.375, 6.5625, 6.75, 6.9375, 7.125, 7.3125,
    7.5, 7.6875, 7.875, 8.0625, 8.25, 8.4375, 8.625, 8.8125, 9, 9.1875,
    9.375, 9.5625, 9.75, 9.9375, 10.125, 10.3125, 10.5, 10.6875, 10.875,
    11.0625, 11.25, 11.4375, 11.625, 11.8125, 12, 12.1875, 12.375, 12.5625,
    12.75, 12.9375, 13.125, 13.3125, 13.5, 13.6875, 13.875, 14.0625, 14.25,
    14.4375, 14.625, 14.8125, 15, 15.1875, 15.375, 15.5625, 15.75, 15.9375,
    16.125, 16.3125, 16.5, 16.6875, 16.875, 17.0625, 17.25, 17.4375, 17.625,

What's going wrong ?

Many thanks for your help in advance !

Best regards,
Klemens
 
Klemens,
This topic has previously been covered on this forum. Please take a look here:
http://forum.mmm.ucar.edu/phpBB3/viewtopic.php?p=1127#p1127
 
I have read this topic and had contact with Andrea. In her case the problem was due to the definition of the region in the download script.
The only other information potentially related to my problem from this thread is the violation of the grib1 standards but since I am expecting grib2 data from CDS, I thought it might not be relevant ?

Klemens
 
Hi,
Thanks for posting the blog you found. Hopefully this will be helpful to others who may run across the same problem.
 
You must log in or register to reply here.
Top