Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Code change for RRFS application #4

Closed
wants to merge 14 commits into from
Closed
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 22 additions & 2 deletions sorc/chgres_cube.fd/input_data.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand All @@ -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."
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down