Page 1 of 1

Error in converting GRIB2 data using api2arl

Posted: June 4th, 2019, 3:23 pm
by cherio
I am trying to convert grib2 data into arl format and having trouble.
I appreciate it if you gave any advice in moving forward.


api2arl version:api2arl_v4

Meteorological data
  • Japan Meteorological Agency Meso-Scale Model (GPV-MSM)
  • Cycle:8/day (00,03,06,09,12,15,18,21UTC)
  • Timestep:3-hourly (for pressure surface)
  • Format:GRIB2
api2arl.cfg: created by api2arl

Code: Select all

&SETUP
 numatm = 12,
 atmgrb = 'gh','pres','t','pt','u','v','w','z','r','qv','q','etad',
 atmcat =     3 ,   3 ,    0 ,   0 ,   2 ,    2 ,    2 ,   2 ,   1 ,   1 ,    1,   2,
 atmnum =     5 ,   0 ,    0 ,   2 ,   2 ,    3 ,    8 ,   9 ,   1 ,   2 ,    2,   8,
 atmcnv =   1.0 ,0.01 ,  1.0 , 1.0 , 1.0 ,  1.0 , 0.01 , 1.0 , 1.0 , 1.0 ,  1.0, 0.01
 atmarl = 'HGTS','PRES','TEMP','THET','UWND','VWND','WWND','DZDT','RELH','SPHU','SPHU','WWND',
 numsfc = 21,
 sfcgrb = 'orog','sp','tp','tcc','hpbl','uflx','vflx','dswr','shtf','lhtf','2t','r','10u','10v','ewss','nsss','z','sshf','ssr','ncpc','acpc',
 sfccat =   3,   3,   1,  6,  19,  2, 17,  4,  0,  0,  0,  1,  2,  2,  2,  2,  3,  0,  4,   1,   1
 sfcnum =   5,   0,   8,  1,   3, 17, 18,  0, 11, 10,  0,  1,  2,  3,  2,  3,  5, 11,  0,   9,  10
 sfccnv = 1.0,0.01,.001,1.0, 1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,0.1,1.0,1.0,.001,.001
 sfcarl = 'SHGT','PRSS','TPPT','TCLD','PBLH','UMOF','VMOF','DSWF','SHTF','LHTF','T02M','RH2M','U10M','V10M','UMOF','VMOF','SHGT','SHTF','DSWF','TPPN','TPPC',
/
arldata.cfg: created by api2arl (added missing RELH by hand)

Code: Select all

Model Type:         rjtd
Grid Numb:            99
Vert Coord:            2
Pole Lat:                47.60
Pole Lon:               150.00
Ref Lat:                  0.10
Ref Lon:                  0.12
Grid Size:                0.00
Orientation:              0.00
Cone Angle:               0.00
Sync X Pt:                1.00
Sync Y Pt:                1.00
Sync Lat:                22.40
Sync Lon:               120.00
Reserved:                 0.00
Numb X pt:           241
Numb Y pt:           253
Numb Levels:          17
Level  0:           .00000  0
Level  1:           1000.0  6 HGTS TEMP UWND VWND WWND RELH
Level  2:           975.00  6 HGTS TEMP UWND VWND WWND RELH
Level  3:           950.00  6 HGTS TEMP UWND VWND WWND RELH
Level  4:           925.00  6 HGTS TEMP UWND VWND WWND RELH
Level  5:           900.00  6 HGTS TEMP UWND VWND WWND RELH
Level  6:           850.00  6 HGTS TEMP UWND VWND WWND RELH
Level  7:           800.00  6 HGTS TEMP UWND VWND WWND RELH
Level  8:           700.00  6 HGTS TEMP UWND VWND WWND RELH
Level  9:           600.00  6 HGTS TEMP UWND VWND WWND RELH
Level 10:           500.00  6 HGTS TEMP UWND VWND WWND RELH
Level 11:           400.00  6 HGTS TEMP UWND VWND WWND RELH
Level 12:           300.00  6 HGTS TEMP UWND VWND WWND RELH
Level 13:           250.00  6 HGTS TEMP UWND VWND WWND RELH
Level 14:           200.00  6 HGTS TEMP UWND VWND WWND RELH
Level 15:           150.00  6 HGTS TEMP UWND VWND WWND RELH
Level 16:           100.00  6 HGTS TEMP UWND VWND WWND RELH
Error I got

Code: Select all

 Existing decoding configuration:api2arl.cfg
 Existing encoding configuration:arldata.cfg
  NOTICE pakset:
  Number of index records =            1
  Number of records /time =           97
  NOTICE pakini: start initialization rec =            1
  NOTICE pakini:   end initialization rec =           97
 *ERROR* pakrec - at index:            1
   Argument list times    :           11           3          13           6
   Do not match initial   :           11           3          13           3
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_UNDERFLO
W_FLAG IEEE_DENORMAL
STOP 900
I cannot figure out the path forward and I really appreciate your advice.

Re: Error in converting GRIB2 data using api2arl

Posted: June 5th, 2019, 10:37 am
by Fantine
api2arl_v4 is an example to create a converter for GRID2 data. Some customization are always needed for different datasets. https://www.ready.noaa.gov/documents/Tu ... cnvrt.html

The cfg file lists the variables from the met data (atmgrb and sfcgrb) matching to those in HYSPLIT format (atmarl and sfcarl). The one you created currently is for ECMWF variables. You will need to come up with a similar table for the GPV-MSM variables and modify the code in order to read in the met data.