"Error W3FI00: 4" using grib2arl for ERA-Interim data

andries · Post by **andries** » July 21st, 2014, 11:17 am

Hello,

I intend to use hysplit (in windows) for calculating back trajectories with use of ERA-Interim data (reduced Gaussian grids N128, global data). When I want to convert the ERA-Interim data with grib2arl, I receive the following error:

108
108
while executing
"exec $exec_dir/grib2arl.exe $arg1 $arg2 $arg3 $arg4 $arg9 $arg5 $arg6 $arg7 $arg8 >STDOUT"
(procedure "xtrct_data" line 60)
invoked from within
"xtrct_data "
invoked from within
".eraarch.bot.save invoke"
("uplevel" body line 1)
invoked from within
"uplevel #0 [list $w invoke]"
(procedure "tk::ButtonUp" line 24)
invoked from within
"tk::ButtonUp .eraarch.bot.save"
(command bound to event)

As suggested in a previous post with a similar problem viewtopic.php?f=14&t=241, I tried to use the following command in DOS:

grib2arl -ioutput-3d.grib -y60 -x80 -g3 -p0 -k25 -soutput-2d.grib -cei_an_geopotsurf_gg128_200911.grib

with "output-3d.grib", "output-2d.grib", and "ei_an_geopotsurf_gg128_200911.grib" containing the merged GH,T,u,v,w,RH variables on pressurelevels, the merged 2m-temp, u10m,v10m at the surface level, and the surface geopotential as the invariant, respectively. After running this command, there appears "108" at the command line, and the file "MESSAGE" contains the following text:

Command Line Arguments:
Cntr - 60 -80
Grid - 3
Size - 100 100 25
SfcP - 0
Rain - 6
Init - 1
---------------------------------------------------------
Opened file: output-3d.grib
-----------------------------------------------------
SUBROUTINE analyze ...
Input buffer allocation: 262252
Error W3FI00: 4

Running grib2arl with other parameter settings (for -y, -x, -k, -p) results in the same error.

I also tried to use grib2arl without the "i", "s", "c" in front of the files:

grib2arl -output-3d.grib -y60 -x80 -g3 -p0 -k25 -output-2d.grib -ei_an_geopotsurf_gg128_200911.grib

this results in another error message. At the command line appears:

forrtl: severe (157): Program Exception - access violation

and in the MESSAGE file:

Command Line Arguments:
Cntr - 60 80
Grid - 3
Size - 100 100 25
SfcP - 0
Rain - 6
Init - 1
---------------------------------------------------------
SUBROUTINE setgrid ...
Configured global lat/lon output grid
Output array allocation: 0 0

And in the CFG_ARL file:

Model Type: G
Grid Numb: 99
Vert Coord: 0
Pole Lat: 0.00
Pole Lon: 0.00
Ref Lat: 0.00
Ref Lon: 0.00
Grid Size: 0.00
Orientation: 0.00
Cone Angle: 0.00
Sync X Pt: 1.00
Sync Y Pt: 1.00
Sync Lat: 0.00
Sync Lon: 0.00
Reserved: 0.00
Numb X pt: 0
Numb Y pt: 0
Numb Levels: 1

I would be very grateful for any help. And of course I would be happy to provide more information if needed to find a solution.
Thank you very much,
Andries

Post by **ariel.stein** » July 21st, 2014, 11:56 am

There is probably some problem with your data files.
Can you send me a short extract of the 3d, 2d and invariant files?

andries · Post by **andries** » July 21st, 2014, 1:07 pm

Thank you for the fast response. I used wgrib -V for extracting some information of the following files and copy-pasted the first few lines below. I hope it is useful. please let me know if you need any other information.

output-3D.grib (GH, T,u,v,w,RH on pressurelevels), rec 0 and 1, and rec 38 and 39.

wgrib -V output-3d.grib | more
rec 1:0:date 2009110100 Z kpds5=129 kpds6=100 kpds7=1 levels=(0,1) grid=255 1 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 430940 472280 num bits 16 BDS_Ref 430940 DecScale 0 BinScale 0

rec 2:262256:date 2009110100 Z kpds5=129 kpds6=100 kpds7=2 levels=(0,2) grid=255 2 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 385585 419242 num bits 16 BDS_Ref 385585 DecScale 0 BinScale 0

.............................................................................................................

rec 38:9703472:date 2009110100 W kpds5=135 kpds6=100 kpds7=1 levels=(0,1) grid=255 1 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
W=Vertical velocity [Pa s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -0.00203862 0.00130974 num bits 16 BDS_Ref -0.00203862 DecScale 0 BinScale -24

rec 39:9965728:date 2009110100 W kpds5=135 kpds6=100 kpds7=2 levels=(0,2) grid=255 2 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
W=Vertical velocity [Pa s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -0.0033009 0.00412727 num bits 16 BDS_Ref -0.0033009 DecScale 0 BinScale -23

output-2D.grib (2m-temp, 10m u and v):

wgrib -V ../2D/output-2d.grib | more
rec 1:0:date 2009110100 10U kpds5=165 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
10U=10 metre U wind component [m s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -26.5468 25.1232 num bits 16 BDS_Ref -26.5468 DecScale 0 BinScale -10

rec 2:262256:date 2009110100 10V kpds5=166 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
10V=10 metre V wind component [m s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -22.181 18.5553 num bits 16 BDS_Ref -22.181 DecScale 0 BinScale -10

rec 3:524512:date 2009110100 2T kpds5=167 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
2T=2 metre temperature [K]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 223.119 312.492 num bits 16 BDS_Ref 223.119 DecScale 0 BinScale -9

rec 4:786768:date 2009110106 10U kpds5=165 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
10U=10 metre U wind component [m s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -21.9989 23.1173 num bits 16 BDS_Ref -21.9989 DecScale 0 BinScale -10

rec 5:1049024:date 2009110106 10V kpds5=166 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
10V=10 metre V wind component [m s**-1]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -19.9139 17.6095 num bits 16 BDS_Ref -19.9139 DecScale 0 BinScale -10

rec 6:1311280:date 2009110106 2T kpds5=167 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
2T=2 metre temperature [K]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 229.269 315.098 num bits 16 BDS_Ref 229.269 DecScale 0 BinScale -9

Invariant (surface geopotential)

rec 1:0:date 2009110100 Z kpds5=129 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -1140.27 53026.7 num bits 16 BDS_Ref -1140.27 DecScale 0 BinScale 0

rec 2:262320:date 2009110106 Z kpds5=129 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data -1140.27 53026.7 num bits 16 BDS_Ref -1140.27 DecScale 0 BinScale 0

here I also put the information of one of the files (geopotential height on plevels) that I used to produce the output-3D file (that contains also the other variables). When I use this geopotential file instead of the output-3D file I recieve the same error.

rec 1:0:date 2009110100 Z kpds5=129 kpds6=100 kpds7=1 levels=(0,1) grid=255 1 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 430940 472280 num bits 16 BDS_Ref 430940 DecScale 0 BinScale 0

rec 2:262320:date 2009110100 Z kpds5=129 kpds6=100 kpds7=2 levels=(0,2) grid=255 2 mb ECext=1 class=14 type=analysis stream=1025 Version=0001 anl:
Z=Geopotential [m**2 s**-2]
timerange 0 P1 0 P2 0 TimeU 1 nx 512 ny 256 GDS grid 4 num_in_ave 0 missing 0
center 98 subcenter 0 process 128 Table 128 scan: WE:NS winds(N/S)
gaussian: lat 89.463000 to -89.463000
long 0.000000 to 359.297000 by 0.703000, (512 x 256) scan 0 mode 128 bdsgrid 1
min/max data 385585 419242 num bits 16 BDS_Ref 385585 DecScale 0 BinScale 0

Post by **ariel.stein** » July 21st, 2014, 1:36 pm

Can you please send me the files?

andries · Post by **andries** » July 22nd, 2014, 6:36 am

sure, I uploaded the files with use of "copy". Please, find them under the following link:
https://copy.com/sQ9eTYvw7c8Q

Post by **ariel.stein** » July 23rd, 2014, 4:01 pm

I can reproduce the error in my computer. It seems your are extracting a grid that is not currently accepted by the code. Can you used another grid definition?

andries · Post by **andries** » July 24th, 2014, 5:09 am

What grid definition/resolution would you suggest? The ERA-Interim data files I used and sent under the link are in the default format: gaussian grid, N128 (128 gridpoints from pole to equator), in total 256 and 512 gridpoints in latitudinal and longitudinal direction, respectively.

Post by **ariel.stein** » July 24th, 2014, 9:37 am

Where are you getting the data from?
I am getting my grib data from http://data-portal.ecmwf.int/data/d/interim_daily/
The projection for that database is latitude longitude.

andries · Post by **andries** » July 25th, 2014, 10:52 am

We obtained the data via the MARS archive in a gaussian grid. I just found a solution that worked: to convert the data with use of cdo from gaussian to latlon projection with the following command:

cdo remapbil,r720x360 {ifile} {ofile}

With the ERA-Interim data in latlon projection the function grib2arl works, both in the windows and unix version of Hysplit. Thank you very much for your help in solving the issue.

singingsea · Post by **singingsea** » August 14th, 2014, 3:54 pm

I tried to get ERA data from http://data-portal.ecmwf.int/data/d/interim_daily/ to avoid the projection issue.

But I end up with some different warning messages

forrtl: severe (64): input conversion error, unit 20, file F:\hysplit4\working\CFG_ARL

Image PC Routine Line Source
grib2arl.exe 00000001400AAFF8 Unknown Unknown Unknown
grib2arl.exe 00000001400A60D9 Unknown Unknown Unknown
grib2arl.exe 000000014006B67D Unknown Unknown Unknown
grib2arl.exe 0000000140052807 Unknown Unknown Unknown
grib2arl.exe 00000001400520D1 Unknown Unknown Unknown
grib2arl.exe 000000014004CBBE Unknown Unknown Unknown
grib2arl.exe 0000000140049F05 Unknown Unknown Unknown
grib2arl.exe 0000000140017316 Unknown Unknown Unknown
grib2arl.exe 000000014000EE9A Unknown Unknown Unknown
grib2arl.exe 00000001400015E7 Unknown Unknown Unknown
grib2arl.exe 00000001400AD1DC Unknown Unknown Unknown
grib2arl.exe 0000000140094C1F Unknown Unknown Unknown
kernel32.dll 0000000076ED59ED Unknown Unknown Unknown
ntdll.dll 000000007724C541 Unknown Unknown Unknown
forrtl: severe (64): input conversion error, unit 20, file F:\hysplit4\working\CFG_ARL

Image PC Routine Line Source
grib2arl.exe 00000001400AAFF8 Unknown Unknown Unknown
grib2arl.exe 00000001400A60D9 Unknown Unknown Unknown
grib2arl.exe 000000014006B67D Unknown Unknown Unknown
grib2arl.exe 0000000140052807 Unknown Unknown Unknown
grib2arl.exe 00000001400520D1 Unknown Unknown Unknown
grib2arl.exe 000000014004CBBE Unknown Unknown Unknown
grib2arl.exe 0000000140049F05 Unknown Unknown Unknown
grib2arl.exe 0000000140017316 Unknown Unknown Unknown
grib2arl.exe 000000014000EE9A Unknown Unknown Unknown
grib2arl.exe 00000001400015E7 Unknown Unknown Unknown
grib2arl.exe 00000001400AD1DC Unknown Unknown Unknown
grib2arl.exe 0000000140094C1F Unknown Unknown Unknown
kernel32.dll 0000000076ED59ED Unknown Unknown Unknown
ntdll.dll 000000007724C541 Unknown Unknown Unknown
while executing
"exec $exec_dir/grib2arl.exe $arg1 $arg2 $arg3 $arg4 $arg9 $arg5 $arg6 $arg7 $arg8 >STDOUT"
(procedure "xtrct_data" line 60)
invoked from within
"xtrct_data "
invoked from within
".eraarch.bot.save invoke"
("uplevel" body line 1)
invoked from within
"uplevel #0 [list $w invoke]"
(procedure "tk::ButtonUp" line 24)
invoked from within
"tk::ButtonUp .eraarch.bot.save"
(command bound to event)

Could anyone give some help? Thanks!

HYSPLIT Forum: hysplitbbs.arl.noaa.gov

"Error W3FI00: 4" using grib2arl for ERA-Interim data

"Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data

Re: "Error W3FI00: 4" using grib2arl for ERA-Interim data