xtrct_grid issue

[sherryj]

Hello All,

I attempted to use 'xtrct_grid' tool to obtain a subgrid of raw HRRR file and have written a shell script (

download_hrrr_test.txt

(4.86 KiB) Downloaded 370 times

) 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

[sonny.zinn]

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".

[sherryj]

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

[sonny.zinn]

The input for directory must end with a slash (/). Please check your script if it includes a trailing slash.

[sherryj]

Thank you! After adding a slash (/) at the end, it works!

Best,
Sherry