Results 39 comments of cuill

@calquigs In your hgrid.gr3, no need to add an extra node to make it closed, which caused inconsistent between the number of nodes read from header information and the actual...

@SorooshMani-NOAA This error happens when the requested data is not available yet on aws.

@SorooshMani-NOAA I would suggest you use gfs2 and hrrr3 because there are issues with original gfs/hrrr data. Below are plots for hrrr and hrrr3. As you can see, hrrr_orig puts...

@SorooshMani-NOAA With your script, rnday is 3.5 days. HRRR has 49 forecast records. hrrr3 considered both hindcast and forecast scenarios: For hindcast, rnday is the real run period and set...

Here is the script I just used to generate hrrr_2022093006.nc now = datetime.now() last_cycle = np.datetime64(pd.DatetimeIndex([now-timedelta(hours=2)]).floor('6H').values[0], 'h').tolist() print(last_cycle) rnday = 1 record = 2 hgrid = Hgrid.open('../../static/hgrid.gr3', crs='epsg:4326') pscr='/sciclone/pscr/lcui01/HRRR/' hrrr...

With your script, it generates **start** and **today** as follows: last_cycle is 2022-09-30 12:00:00 start is 2022-09-29 00:00:00 rnday is 3.5 So you can set rnday = 2, because only...

No, I would say the remaining 0.5 days will be covered by sflux1, 1 day covered by hrrr_2022092900.nc, and 2 days covered by hrrr_2022093000.nc

Here is the code to generate sflux1 from gfs2: **now = datetime.now() last_cycle = np.datetime64(pd.DatetimeIndex([now-timedelta(hours=2)]).floor('6H').values[0], 'h').tolist() start = (last_cycle - timedelta(days=1)).replace(hour=0) rnday = 2 record = 3 hgrid = Hgrid.open('../../static/hgrid.gr3',...

@SorooshMani-NOAA Thanks for your suggestion. I'll add it to the to-do list but I don't have ETA at the moment. I have other priorities to do. You are welcome to...