Hello All,
I attempted to use 'xtrct_grid' tool to obtain a subgrid of raw HRRR file and have written a shell script ( ) to automate these procedures. However, when I execute the script, the subgridding often generate files with smaller sizes compared to the correct outputs. For example, the following message is the processing log that I received when I subgrid the HRRR file "20220101_12-17_hrrr" based on the script:
Enter meteorological directory ...
/data/HRRR/
Enter meteorological file name ...
20220101_12-17_hrrr
Lower left corner: 21.1381 -122.7190
Upper right corner: 47.8419 -60.9162
Number of data levels (inc sfc): 36
Enter new lower left corner (lat,lon):
26.43,-91.06
Enter new upper right corner (lat,lon):
36.76,-74.54
Enter desired number of levels (incl sfc): 36
36
NOTICE pakset:
Number of index records = 1
Number of records /time = 299
22 1 112 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 1
NOTICE pakini: end initialization rec = 299
22 1 113 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 300
NOTICE pakini: end initialization rec = 598
22 1 114 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 599
NOTICE pakini: end initialization rec = 897
22 1 115 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 898
NOTICE pakini: end initialization rec = 1196
As you may have already noticed from the above message, the process of HRRR file missed the following procedures:
22 1 116 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 1197
NOTICE pakini: end initialization rec = 1495
22 1 117 1 0AAINDX 0 0.0000000E+00 0.0000000E+00
NOTICE pakini: start initialization rec = 1496
NOTICE pakini: end initialization rec = 1794
If I re-run the same script multiple times, it occasionally generates the correct outputs.
Additionally, when I manually run the xtrct_grid tool without using the script, it never gives me incorrect output. Does anyone know why the xtrct_grid is not stable when executed through a script?
Thank you!
Regards,
Sherry
xtrct_grid issue
-
- Posts: 361
- Joined: May 8th, 2019, 1:31 pm
- Registered HYSPLIT User: Yes
Re: xtrct_grid issue
Instead of running xtrct_grid via the expect command, you might want to directly run xtrct_grid as follows:
${EXECDIR}/xtrct_grid <<EOF
${OUTDIR}/
${hrrr1}
${lower}
${upper}
${levels}
EOF
Note there is no space between "<<" and "EOF".
${EXECDIR}/xtrct_grid <<EOF
${OUTDIR}/
${hrrr1}
${lower}
${upper}
${levels}
EOF
Note there is no space between "<<" and "EOF".
Re: xtrct_grid issue
Thank you for the suggestions! I updated the code, but get the following error message:
subregion of 20220101_00-05_hrrr
Enter meteorological directory ...
Enter meteorological file name ...
Unable to find file: 20220101_00-05_hrrr
On local directory : /data/HRRR/Raw
Check input CONTROL file for correct values
STOP 900
stat: cannot statx '/data/HRRR/extract.bin': No such file or directory
I looks like the xtrct_grid do not load "meteorological directory" and "file name" successfully using the updated EOF code.
Best,
Sherry
subregion of 20220101_00-05_hrrr
Enter meteorological directory ...
Enter meteorological file name ...
Unable to find file: 20220101_00-05_hrrr
On local directory : /data/HRRR/Raw
Check input CONTROL file for correct values
STOP 900
stat: cannot statx '/data/HRRR/extract.bin': No such file or directory
I looks like the xtrct_grid do not load "meteorological directory" and "file name" successfully using the updated EOF code.
Best,
Sherry
-
- Posts: 361
- Joined: May 8th, 2019, 1:31 pm
- Registered HYSPLIT User: Yes
Re: xtrct_grid issue
The input for directory must end with a slash (/). Please check your script if it includes a trailing slash.
Re: xtrct_grid issue
Thank you! After adding a slash (/) at the end, it works!
Best,
Sherry
Best,
Sherry