GFS 0p25 File Error

Topics about the HYSPLIT trajectory model.
StonesInAPond
Posts: 1
Joined: May 12th, 2020, 12:05 pm
Registered HYSPLIT User: Yes

GFS 0p25 File Error

Post by StonesInAPond »

Hi Everyone,

I am new to this board so have patience with me in case something feels a bit off with my post.

Anyway, I am using the GFS archive to conduct HYSPLIT modelling. I have been using the archive to model backwards atmospheric trajectories in hourly timesteps. I have been receiving errors using the 2019-06-03 GFS 0p25 archive file (3rd of June 2019). My error code has been “ERROR metset: meteorological data time interval varies. Check input meteo file.” With it also stating “Changed from 180 min to 4500 min, At day/hr 3 6”. Would anyone be able to advice if there is anything wrong with this file? I have tried checking the logs of missing data but without any success.

Cheers,
Stones
barbara.stunder
Posts: 451
Joined: November 9th, 2012, 4:23 pm
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by barbara.stunder »

Thank you for bringing this to our attention. Yes, there is a problem with that file. I deleted it from the archive and noted that in the file readme_gfs0p25_missings.txt (ftp://ftp.arl.noaa.gov/archives/gfs0p25.v1). Unfortunately we cannot fix the file, so as an alternative, use /archives/gdas0p5/20190603_gdas0p5 instead.
croc
Posts: 17
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by croc »

Hi,
I am sorry to add a new message to an old thread but my problem
is strictly related to the suggested method.
I am trying to run backtrajectories starting 19 06 02 22, therefore
I am using meteo files:
20190603_gfs0p25
20190602_gfs0p25
20190601_gfs0p25
20190531_gfs0p25
.....

nevertheless, I am getting the error:

HYSPLIT - Initialization
HYSPLIT version: hysplit.v5.2.0
Last Changed Date: 2022-01-16
NOTICE: using namelist file SETUP.CFG

WARNING metset: Old format meteo data grid 12 1
At line 126 of file metold.f
Fortran runtime error: Bad value during integer read

I made many other test but with no luck. Could you please
suggest me what I might be mistaken.
Regards,

Stefano

NOTE: to complete this message, this specific problem is due,
as sonny.zinn suggested (kudo), to an incomplete meteo file.
Yet, I cannot calculate even very short trajetories starting between
2 Jun 2019 22:00 and 3 Jun 2019 00:00 as I describe in my next post.
Last edited by croc on December 1st, 2022, 6:29 am, edited 1 time in total.
sonny.zinn
Posts: 354
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by sonny.zinn »

Since 20190603_gfs0p25 does not exist, we suppose that the missing file is substituted with 20190603_gdas0p5. Our cursory checks on files from May 1, 2019 to June 2, 2019 did not find an issue. How far back does your backward trajectories go?

We will create the listing.md5.txt file containing MD5 checksums in the gfs0p25.v1 archive directory. The checksums may be used to see if one or more files you have are different from our files. Once the checksum file becomes available, we will add a post to this thread. It may take a few days to create the checksum file.
croc
Posts: 17
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by croc »

Many thanks for your answer. You are right: I sym linked
20190603_gfs0p25 to 20190603_gdas0p5. Sorry for the mistake.

I further reduced the CONTROL file to reproduce the problem:

test1:
----------------------------------------------------------------------
- cat CONTROL
19 06 02 22
1 #NUMBER OF STARTING LOCATIONS
50.00 12.00 50 #STARTING 1: LATITUDE LONGITUDE HEIGHT (m-agl)
-10 #TOTAL RUN TIME (backwards)
0 #VERTICAL MOTION CALCULATION METHOD
11000 #TOP OF MODEL DOMAIN (m-AGL)
2 #NUMBER OF INPUT DATA GRIDS
./
20190603_gdas0p5
./
20190602_gfs0p25
./
tdump

- hyts_std
HYSPLIT - Initialization
HYSPLIT version: hysplit.v5.2.0
Last Changed Date: 2022-01-16
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 1 1 GFSQ
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 2 1 GFSG
Calculation Started ... please be patient
*ERROR* metpos: start point not within (x,y,t) any data file
- start time before start of meteorology data
- start time after end of meteorology data
STOP 900
----------------------------------------------------------------------


test2:
----------------------------------------------------------------------
- cat CONTROL
19 06 02 20
1 #NUMBER OF STARTING LOCATIONS
50.00 12.00 50 #STARTING 1: LATITUDE LONGITUDE HEIGHT (m-agl)
-10 #TOTAL RUN TIME (backwards)
0 #VERTICAL MOTION CALCULATION METHOD
11000 #TOP OF MODEL DOMAIN (m-AGL)
2 #NUMBER OF INPUT DATA GRIDS
./
20190603_gdas0p5
./
20190602_gfs0p25
./
tdump

- hyts_std
HYSPLIT - Initialization
HYSPLIT version: hysplit.v5.2.0
Last Changed Date: 2022-01-16
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 1 1 GFSQ
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 2 1 GFSG
Calculation Started ... please be patient
Percent complete: 10.0
Percent complete: 20.0
Percent complete: 30.0
Percent complete: 40.0
Percent complete: 50.0
Percent complete: 60.0
Percent complete: 70.0
Percent complete: 80.0
Percent complete: 90.0
Percent complete: 100.0
Complete Hysplit
----------------------------------------------------------------------


test3:
----------------------------------------------------------------------
- cat CONTROL
19 06 03 1
1 #NUMBER OF STARTING LOCATIONS
50.00 12.00 50 #STARTING 1: LATITUDE LONGITUDE HEIGHT (m-agl)
-10 #TOTAL RUN TIME (backwards)
0 #VERTICAL MOTION CALCULATION METHOD
11000 #TOP OF MODEL DOMAIN (m-AGL)
2 #NUMBER OF INPUT DATA GRIDS
./
20190603_gdas0p5
./
20190602_gfs0p25
./
tdump

- hyts_std
HYSPLIT - Initialization
HYSPLIT version: hysplit.v5.2.0
Last Changed Date: 2022-01-16
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 1 1 GFSQ
*WARNING*: metini - vertical velocity missing, switching to divergence
Input data file: 2 1 GFSG
Calculation Started ... please be patient
Percent complete: 10.0
Percent complete: 20.0
Percent complete: 30.0
Percent complete: 40.0
Percent complete: 50.0
Percent complete: 60.0
Percent complete: 70.0
Percent complete: 80.0
Percent complete: 90.0
Percent complete: 100.0
Complete Hysplit
----------------------------------------------------------------------


It looks to me like something related to the mixed use of GFS and
GDAS meteo files.
Regards,
sonny.zinn
Posts: 354
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by sonny.zinn »

The listing.md5.txt file is added to the gfs0p25.v1 archive. MD5 checksums are listed in the file along with their corresponding filenames. You might want to compare file checksums of your copy of gfs0p25.v1 files with those in listing.md5.txt. For Linux computers, the md5sum program can be used to compute MD5 checksums. Windows users may use CertUtil. The checksums should match.
croc
Posts: 17
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by croc »

Thanks to sonny.zinn for calculating md5sum of all the gfs0p25.v1
archive. This will help me to avoid using corrupted/partially
downloaded files.
Now I can confirm the md5 checksum of one of the two
files (20190602_gfs0p25) used to reproduce the problem
shown in test1-test3 of my last post (the other is from gdas0p5
archive but I don't expect to be corrupted, at least on the
basis of its size).

Do you have any clue about the metpos error?
Regards,
sonny.zinn
Posts: 354
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by sonny.zinn »

Please see the response on December 6, 2022 below which supersedes this post (December 6, 2022)


Interpolation across meteorological files of different grids and/or resolution does not appear to be working. If the CONTROL file is changed to use GDAS0p5 files, as shown below, no runtime error occurs. We will open a ticket about the interpolation issue.

Code: Select all

19 06 03 1
1 #NUMBER OF STARTING LOCATIONS
50.00 12.00 50 #STARTING 1: LATITUDE LONGITUDE HEIGHT (m-agl)
-10 #TOTAL RUN TIME (backwards)
0 #VERTICAL MOTION CALCULATION METHOD
11000 #TOP OF MODEL DOMAIN (m-AGL)
2 #NUMBER OF INPUT DATA GRIDS
./
20190603_gdas0p5
./
20190602_gdas0p5
./
tdump
croc
Posts: 17
Joined: January 22nd, 2015, 10:25 am
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by croc »

Thank you sonny.zinn for opening the ticket. Could you please
post here its progresses when available or point us to the related thread,
if public?

Regards,
sonny.zinn
Posts: 354
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: GFS 0p25 File Error

Post by sonny.zinn »

Use of an additional meteorological ameliorates the interpolation issue. Your CONTROL file posted on December 1, 2022 in the above used two meteorological files, namely, 20190603_gdas0p5 and 20190602_gfs0p25. Adding 20190602_gdas0p5 to the CONTROL file resolves the runtime issue: see the CONTROL file below. HYSPLIT will use the finer grid, 20190602_gfs0p25, for computation of trajectory on June 2 and will use 20190602_gdas0p5 when interpolating meteorological variables from June 3rd to 2nd.

Code: Select all

19 06 02 22
1
50.00 12.00 50
-10
0
11000
3
/pub/archives/gdas0p5/
20190603_gdas0p5
/pub/archives/gdas0p5/
20190602_gdas0p5
/pub/archives/gfs0p25.v1/
20190602_gfs0p25
./
tdump
Post Reply

Return to “Trajectory Model”