Skip to content

Commit

Permalink
chgres_cube fix to correctly read vertical coordinate type from GRIB2…
Browse files Browse the repository at this point in the history
… data

Operational RAP GRIB2 files on the native grid (and possibly other GRIB2 data types)
are provided on isobaric coordinates, but contain some data on the 1.0 hybrid 
coordinate.  This fixes an error with chgres_cube proceeding as if all of the data 
were on hybrid coordinates. 

Addresses Issue #222
  • Loading branch information
LarissaReames-NOAA authored Nov 18, 2020
1 parent bbb662f commit 91159ad
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion sorc/chgres_cube.fd/input_data.F90
Original file line number Diff line number Diff line change
Expand Up @@ -2519,7 +2519,7 @@ subroutine read_input_atm_grib2_file(localpet)
if (.not.lret) call error_handler("OPENING GRIB2 ATM FILE.", iret)

print*,"- READ VERTICAL COORDINATE."
iret = grb2_inq(the_file,inv_file,":var0_2","_0_0:"," hybrid ")
iret = grb2_inq(the_file,inv_file,":var0_2","_0_0:",":10 hybrid level:")

if (iret <= 0) then
lvl_str = "mb:"
Expand Down

0 comments on commit 91159ad

Please sign in to comment.