diff --git a/sorc/chgres_cube.fd/input_data.F90 b/sorc/chgres_cube.fd/input_data.F90 index 497ea442f..7107ae141 100644 --- a/sorc/chgres_cube.fd/input_data.F90 +++ b/sorc/chgres_cube.fd/input_data.F90 @@ -2607,7 +2607,7 @@ subroutine read_input_atm_grib2_file(localpet) if (localpet == 0) print*,"- FILE CONTAINS SPFH." endif - if (localpet == 0) print*,"- FIND ICMR, SCLIWC, OR CICE IN FILE" + if (localpet == 0) print*,"- FIND ICMR, SCLIWC, CICE, OR CIMIXR IN FILE" iret = grb2_inq(the_file,inv_file,trac_names_grib_1(4),trac_names_grib_2(4),lvl_str_space) if (iret <= 0) then @@ -2619,7 +2619,13 @@ subroutine read_input_atm_grib2_file(localpet) if (iret <= 0) then iret = grb2_inq(the_file,inv_file, ':var0_2','_6_0:',lvl_str_space) if (iret <= 0 ) then - call handle_grib_error(vname, slevs(1),method,value,varnum,rc,var=dummy2d) + iret = grb2_inq(the_file,inv_file, ':var0_2','_1_82:',lvl_str_space) + if (iret <= 0 ) then + call handle_grib_error(vname, slevs(1),method,value,varnum,rc,var=dummy2d) + else + trac_names_grib_2(4) = '_1:82:' + if (localpet ==0) print*,"- FILE CONTAINS CIMIXR." + endif else trac_names_grib_2(4) = '_6_0' if (localpet == 0) print*,"- FILE CONTAINS CICE." @@ -5037,6 +5043,10 @@ subroutine read_input_sfc_grib2_file(localpet) vname=":VEG:" rc= grb2_inq(the_file, inv_file, vname,slev,'n=1106:',data2=dummy2d) + if (rc <= 0) then + rc= grb2_inq(the_file, inv_file, ':var2_0_231:',slev,data2=dummy2d) + endif + if (rc <= 0) then rc= grb2_inq(the_file, inv_file, vname,slev,'n=1102:',data2=dummy2d) if (rc <= 0) then @@ -5064,6 +5074,11 @@ subroutine read_input_sfc_grib2_file(localpet) vname=":VEG:" rc= grb2_inq(the_file, inv_file, vname,slev,'n=1107:',data2=dummy2d) + + if (rc <=0) then + rc= grb2_inq(the_file, inv_file, ':var2_0_232:',slev,data2=dummy2d) + endif + if (rc <=0) then rc= grb2_inq(the_file, inv_file, vname,slev,'n=1103:',data2=dummy2d) if (rc <=0) then @@ -5092,6 +5107,11 @@ subroutine read_input_sfc_grib2_file(localpet) loc=varnum) vname=":var0_7_198:" rc= grb2_inq(the_file, inv_file, vname,slev,':n=1108:',data2=dummy2d) + + if (rc <=0) then + rc= grb2_inq(the_file, inv_file, ':LAI:',':surface:', data2=dummy2d) + endif + if (rc <=0) then rc= grb2_inq(the_file, inv_file, vname,slev,':n=1104:',data2=dummy2d) if (rc <=0) then