NaN in tdump file

Post any defects you find in the HYSPLIT software here. The HYSPLIT Developers carefully monitor this list and will work diligently to repair any reported problems. When posting a bug report, please specify both the HYSPLIT version and operating system you are using.
Post Reply
croc
Posts: 23
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

NaN in tdump file

Post by croc »

Dear Hysplit forum,
I am getting NaN fields for lat, lon and height running hyts_std v5.3.2 for
linux.

CONTROL file content is:
22 04 17 04
1 #NUMBER OF STARTING LOCATIONS
-89.997 -24.8 10 #STARTING 1: LATITUDE LONGITUDE HEIGHT (m-agl)
-240 #TOTAL RUN TIME (backwards)
0 #VERTICAL MOTION CALCULATION METHOD
20000 #TOP OF MODEL DOMAIN (m-AGL)
11 #NUMBER OF INPUT DATA GRIDS
/data1/hysplit/meteo/era5/arl/
20220417.arl
/data1/hysplit/meteo/era5/arl/
20220416.arl
/data1/hysplit/meteo/era5/arl/
20220415.arl
/data1/hysplit/meteo/era5/arl/
20220414.arl
/data1/hysplit/meteo/era5/arl/
20220413.arl
/data1/hysplit/meteo/era5/arl/
20220412.arl
/data1/hysplit/meteo/era5/arl/
20220411.arl
/data1/hysplit/meteo/era5/arl/
20220410.arl
/data1/hysplit/meteo/era5/arl/
20220409.arl
/data1/hysplit/meteo/era5/arl/
20220408.arl
/data1/hysplit/meteo/era5/arl/
20220407.arl
./
tdump


SETUP.CFG content is:
&SETUP
khmax = 240,
tm_mixd = 1,
delt = 1,
kagl = 0,
/


Meteo input files are ERA5, produced by ERA52ARL Nov 2019 version.g.1
(downloaded 12 Aug 2024).
arldata.cfg files contain:
Model Type: ERA5
Grid Numb: 99
Vert Coord: 2
Pole Lat: 90.00
Pole Lon: 179.75
Ref Lat: 0.25
Ref Lon: 0.25
Grid Size: 0.00
Orientation: 0.00
Cone Angle: 0.00
Sync X Pt: 1.00
Sync Y Pt: 1.00
Sync Lat: -90.00
Sync Lon: 180.00
Reserved: 0.00
Numb X pt: 1440
Numb Y pt: 721
Numb Levels: 38
Level 1: .00000 11 T02M V10M U10M PRSS PBLH CAPE SHGT TPP1 DSWF LTHF DIFR
Level 2: 1000.0 7 HGTS TEMP UWND VWND WWND RELH DIFW
..................
Level 38: 1.0000 7 HGTS TEMP UWND VWND WWND RELH DIFW


tdump content is:
11 1
ERA5 22 4 7 0 0
ERA5 22 4 8 0 0
ERA5 22 4 9 0 0
ERA5 22 4 10 0 0
ERA5 22 4 11 0 0
ERA5 22 4 12 0 0
ERA5 22 4 13 0 0
ERA5 22 4 14 0 0
ERA5 22 4 15 0 0
ERA5 22 4 16 0 0
ERA5 22 4 17 0 0
1 BACKWARD OMEGA
22 4 17 4 -89.997 -24.800 10.0
2 PRESSURE MIXDEPTH
1 1 22 4 17 4 0 0 0.0 -89.997 -24.800 2799.8 677.8 24.7
1 1 22 4 17 3 0 0 -1.0 -89.993 -32.935 2799.4 677.8 24.4
1 1 22 4 17 2 0 0 -2.0 -89.992 76.521 2801.2 677.4 24.1
1 1 22 4 17 1 0 0 -3.0 -89.993 -106.602 2799.2 677.4 23.5
1 1 22 4 17 0 0 0 -4.0 -89.997 -87.655 2799.6 677.4 23.8
1 1 22 4 16 23 0 0 -5.0 NaN NaN NaN NaN NaN
1 1 22 4 16 22 0 0 -6.0 NaN NaN NaN NaN NaN
..........


Can this problem be fixed in some way?
Best regards,

Stefano
croc
Posts: 23
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

Re: NaN in tdump file

Post by croc »

Dear Hysplit forum,
I think I have a clue to locate the bug.
The problem might be associated with the following warning i get during the execution:
WARNING: advpol.f could not move off pole point

Moreover, I am getting NaN also in MESSAGE. The following is
where the NaN starts:
NOTICE main: 64319006 1 1 -274 -2970360.50 7.91646385 1.00000000
NOTICE main: 64319005 1 1 -275 -2960726.00 112.302284 0.655598700
NOTICE main: 64319004 1 1 -276 -2919918.25 3277.06274 1.00000000
NOTICE main: 64319003 1 1 -277 NaN NaN 1.00000000
NOTICE main: 64319002 1 1 -278 NaN NaN 1.00000000

If I understood correctly the file format, 278 should be the minutes
spent from the beginning of the trajectory, which corresponds to 4.6 hour,
close to the 5th hour where the NaNs make the first appearance in the
tdump file.

Since the warning message is written by the grdsiz routine,
could it be related to this
viewtopic.php?f=12&t=2785&p=7861#p7861
other problem?

Regards,

Stefano
Post Reply

Return to “Bugs”