Skip to content

Commit b579233

Browse files
dmeydmey
committed
add wrf-teb integration
1 parent d96531b commit b579233

18 files changed

+3524
-27
lines changed

‎Registry/Registry.EM_COMMON‎

Lines changed: 62 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -901,6 +901,56 @@ state real CHC_SFCDIF ij misc 1 - r "C
901901
state real CMGR_SFCDIF ij misc 1 - r "CMGR_SFCDIF" "" ""
902902
state real CHGR_SFCDIF ij misc 1 - r "CHGR_SFCDIF" "" ""
903903

904+
# Town Energy Balance (TEB)
905+
state real TEB_TI_BLD ij misc 1 - r "TEB_TI_BLD" "inside building temp computed with its equation evolution" "K"
906+
state real TEB_TI_BLD_EQ ij misc 1 - r "TEB_TI_BLD_EQ" "inside building temperature" "K"
907+
state real TEB_TI_BLDWFR ij misc 1 - r "TEB_TI_BLDWFR" "inside building temperature without heating" "K"
908+
state real TEB_T_CANYON ij misc 1 - r "TEB_T_CANYON" "outdoor temperature without heating" "K"
909+
state real TEB_Q_CANYON ij misc 1 - r "TEB_Q_CANYON" "outdoor air specific humidity" "kg kg-1"
910+
state real TEB_WS_ROOF ij misc 1 - rh "TEB_WS_ROOF" "roof water content" "(kg/m2)"
911+
state real TEB_WS_ROAD ij misc 1 - rh "TEB_WS_ROAD" "road water content" "(kg/m2)"
912+
state real TEB_WSNOW_ROOF i{teb_snly}j misc 1 Z r "TEB_WSNOW_ROOF" "roof snow layers reservoir" "kg m{-2}"
913+
state real TEB_TSNOW_ROOF i{teb_snly}j misc 1 Z r "TEB_TSNOW_ROOF" "roof snow layers temperature" "K"
914+
state real TEB_RSNOW_ROOF i{teb_snly}j misc 1 Z r "TEB_RSNOW_ROOF" "roof snow layers density" "kg m{-3}"
915+
state real TEB_ASNOW_ROOF ij misc 1 - r "TEB_ASNOW_ROOF" "roof snow albedo" "dimensionless"
916+
state real TEB_ESNOW_ROOF ij misc 1 - r "TEB_ESNOW_ROOF" "roof snow emissivity" "dimensionless"
917+
state real TEB_TSSNOW_ROOF ij misc 1 - r "TEB_TSSNOW_ROOF" "roof snow surface temperature" "K"
918+
state real TEB_WSNOW_ROAD i{teb_snly}j misc 1 Z r "TEB_WSNOW_ROAD" "road snow layers reservoir" "kg m{-2}"
919+
state real TEB_TSNOW_ROAD i{teb_snly}j misc 1 Z r "TEB_TSNOW_ROAD" "road snow layers temperature" "K"
920+
state real TEB_RSNOW_ROAD i{teb_snly}j misc 1 Z r "TEB_RSNOW_ROAD" "road snow layers density" "kg m{-3}"
921+
state real TEB_ASNOW_ROAD ij misc 1 - r "TEB_ASNOW_ROAD" "road snow albedo" "dimensionless"
922+
state real TEB_ESNOW_ROAD ij misc 1 - r "TEB_ESNOW_ROAD" "road snow emissivity" "dimensionless"
923+
state real TEB_TSSNOW_ROAD ij misc 1 - r "TEB_TSSNOW_ROAD" "road snow surface temperature" "K"
924+
state real TEB_T_WIN1 ij misc 1 - r "TEB_T_WIN1" "outdoor window temperature" "K"
925+
state real TEB_T_WIN2 ij misc 1 - r "TEB_T_WIN2" "indoor window temperature" "K"
926+
state real TEB_AUX_MAX ij misc 1 - r "TEB_AUX_MAX" "auxiliar variable for autosize calcs" "W m-2"
927+
state real TEB_QI_BLD ij misc 1 - r "TEB_QI_BLD" "indoor air specific humidity" "kg kg-1"
928+
state real TEB_THER_PRODC_DAY ij misc 1 - r "TEB_THER_PRODC_DAY" "present day integrated thermal production of energy" "J/m2"
929+
state real TEB_T_FLOOR i{teb_flly}j misc 1 Z r "TEB_T_FLOOR" "floor layers temperatures" "K"
930+
state real TEB_T_MASS i{teb_flly}j misc 1 Z r "TEB_T_MASS" "internal mass layers temperatures" "K"
931+
state real TEB_T_ROAD i{teb_rdly}j misc 1 Z r "TEB_T_ROAD" "road layers temperatures" "K"
932+
state real TEB_T_ROOF i{teb_roly}j misc 1 Z r "TEB_T_ROOF" "roof layers temperatures" "K"
933+
state real TEB_T_WALL_A i{teb_waly}j misc 1 Z r "TEB_T_WALL_A" "wall layers temperatures (wall A)" "K"
934+
state real TEB_T_WALL_B i{teb_waly}j misc 1 Z r "TEB_T_WALL_B" "wall layers temperatures (wall B)" "K"
935+
state real TEB_HVAC_COOL ij misc 1 - h "TEB_HVAC_COOL" "energy consumption of the cooling system" "W m-2"
936+
state real TEB_HVAC_HEAT ij misc 1 - h "TEB_HVAC_HEAT" "energy consumption of the heating system" "W m-2"
937+
state real TEB_THER_PROD_PANEL ij misc 1 - h "TEB_THER_PROD_PANEL" "thermal energy production of solar panel on roofs" "W m-2"
938+
state real TEB_PHOT_PROD_PANEL ij misc 1 - h "TEB_PHOT_PROD_PANEL" "photovoltaic energy production of solar panel on roofs" "W m-2"
939+
state real TEB_TSK_RURAL ij misc 1 - r "TEB_TSK_RURAL" "TSK for rural fraction (TEB only)" "K"
940+
941+
# BEGIN: INPUT quantities used for evaluating WRF-TEB with TEB OFFLINE
942+
state real TEB_INPUT_TA ij misc 1 - h "TEB_INPUT_TA" " " " "
943+
state real TEB_INPUT_PS ij misc 1 - h "TEB_INPUT_PS" " " " "
944+
state real TEB_INPUT_QA ij misc 1 - h "TEB_INPUT_QA" " " " "
945+
state real TEB_INPUT_WIND ij misc 1 - h "TEB_INPUT_WIND" " " " "
946+
state real TEB_INPUT_DIR ij misc 1 - h "TEB_INPUT_DIR" " " " "
947+
state real TEB_INPUT_DIR_SW ij misc 1 - h "TEB_INPUT_DIR_SW" " " " "
948+
state real TEB_INPUT_SCA_SW ij misc 1 - h "TEB_INPUT_SCA_SW" " " " "
949+
state real TEB_INPUT_LW ij misc 1 - h "TEB_INPUT_LW" " " " "
950+
state real TEB_INPUT_RAIN ij misc 1 - h "TEB_INPUT_RAIN" " " " "
951+
state real TEB_INPUT_SNOW ij misc 1 - h "TEB_INPUT_SNOW" " " " "
952+
state real TEB_INPUT_DIR_CO2 ij misc 1 - h "TEB_INPUT_DIR_CO2" " " " "
953+
# END: INPUT quantities used for evaluating WRF-TEB with TEB OFFLINE
904954

905955
# solar location variables from radiation driver
906956
state real COSZEN ij misc 1 - rh "COSZEN" "COS of SOLAR ZENITH ANGLE" "dimensionless"
@@ -2320,6 +2370,14 @@ rconfig integer maxpatch namelist,physics 1 10
23202370
rconfig integer num_snow_layers namelist,physics 1 3 irh "num_snow_layers" "" ""
23212371
rconfig integer num_snso_layers namelist,physics 1 7 irh "num_snso_layers" "" ""
23222372

2373+
# TEB
2374+
rconfig integer teb_num_floor_layers namelist,physics 1 5 irh "teb_num_floor_layers" "" ""
2375+
rconfig integer teb_num_road_layers namelist,physics 1 5 irh "teb_num_road_layers" "" ""
2376+
rconfig integer teb_num_roof_layers namelist,physics 1 5 irh "teb_num_roof_layers" "" ""
2377+
rconfig integer teb_num_wall_layers namelist,physics 1 5 irh "teb_num_wall_layers" "" ""
2378+
rconfig integer teb_num_snow_layers namelist,physics 1 1 irh "teb_num_snow_layers" "" ""
2379+
rconfig integer teb_test_integration namelist,physics 1 0 irh "teb_test_integration" "" ""
2380+
23232381
rconfig integer num_urban_ndm derived 1 1 irh "num_urban_ndm" "maximum number of street dimensions (ndm in BEP or BEM header)" ""
23242382
rconfig integer num_urban_ng derived 1 1 irh "num_urban_ng" "number of grid levels in the ground (ng_u in BEP or BEM header)" ""
23252383
rconfig integer num_urban_nwr derived 1 1 irh "num_urban_nwr" "number of grid levels in the walls or roof (nwr_u in BEP or BEM header)" ""
@@ -2900,6 +2958,10 @@ package sfclayscheme sf_sfclay_physics==91 - -
29002958
package noahucmscheme sf_urban_physics==1 - state:trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d,sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d,a_u_bep,a_v_bep,a_t_bep,a_q_bep,a_e_bep,b_u_bep,b_v_bep,b_t_bep,b_q_bep,b_e_bep,dlg_bep,dl_u_bep,sf_bep,vl_bep,mh_urb2d,stdh_urb2d,lf_urb2d,lp_urb2d,hgt_urb2d,lb_urb2d,tgr_urb2d,cmcr_urb2d,drelr_urb2d,drelb_urb2d,drelg_urb2d,flxhumr_urb2d,flxhumb_urb2d,flxhumg_urb2d,tgrl_urb3d,smr_urb3d,cmgr_sfcdif,chgr_sfcdif,trl_urb3d,tgl_urb3d,tbl_urb3d
29012959
package bepscheme sf_urban_physics==2 - state:a_u_bep,a_v_bep,a_t_bep,a_q_bep,a_e_bep,b_u_bep,b_v_bep,b_t_bep,b_q_bep,b_e_bep,dlg_bep,dl_u_bep,sf_bep,vl_bep,trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d,sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d,hi_urb2d,lp_urb2d,hgt_urb2d,lb_urb2d,trl_urb3d,tgl_urb3d,tbl_urb3d,tsk_rural
29022960
package bep_bemscheme sf_urban_physics==3 - state:a_u_bep,a_v_bep,a_t_bep,a_q_bep,a_e_bep,b_u_bep,b_v_bep,b_t_bep,b_q_bep,b_e_bep,dlg_bep,dl_u_bep,sf_bep,vl_bep,trb_urb4d,tw1_urb4d,tw2_urb4d,tgb_urb4d,tlev_urb3d,qlev_urb3d,tw1lev_urb3d,tw2lev_urb3d,tglev_urb3d,tflev_urb3d,sf_ac_urb3d,lf_ac_urb3d,cm_ac_urb3d,sfvent_urb3d,lfvent_urb3d,sfwin1_urb3d,sfwin2_urb3d,sfw1_urb3d,sfw2_urb3d,sfr_urb3d,sfg_urb3d,hi_urb2d,lp_urb2d,hgt_urb2d,lb_urb2d,trl_urb3d,tgl_urb3d,tbl_urb3d,tsk_rural
2961+
package tebscheme sf_urban_physics==4 - state:lp_urb2d,hgt_urb2d,lb_urb2d,trl_urb3d,tgl_urb3d,tbl_urb3d,TEB_TI_BLD,TEB_T_CANYON,TEB_Q_CANYON,TEB_TI_BLD_EQ,TEB_TI_BLDWFR,TEB_WS_ROOF,TEB_WS_ROAD,TEB_WSNOW_ROOF,TEB_TSNOW_ROOF,TEB_RSNOW_ROOF,TEB_ASNOW_ROOF,TEB_TSSNOW_ROOF,TEB_ESNOW_ROOF,TEB_WSNOW_ROAD,TEB_TSNOW_ROAD,TEB_RSNOW_ROAD,TEB_ASNOW_ROAD,TEB_TSSNOW_ROAD,TEB_ESNOW_ROAD,TEB_T_WIN1,TEB_T_WIN2,TEB_AUX_MAX,TEB_QI_BLD,TEB_THER_PRODC_DAY,TEB_T_FLOOR,TEB_T_MASS,TEB_T_ROAD,TEB_T_ROOF,TEB_T_WALL_A,TEB_T_WALL_B,TEB_HVAC_COOL,TEB_HVAC_HEAT,TEB_TSK_RURAL
2962+
2963+
package teb_notest teb_test_integration==0 - -
2964+
package teb_test teb_test_integration==1 - state:TEB_INPUT_TA,TEB_INPUT_PS,TEB_INPUT_QA,TEB_INPUT_WIND,TEB_INPUT_DIR,TEB_INPUT_DIR_SW,TEB_INPUT_SCA_SW,TEB_INPUT_LW,TEB_INPUT_RAIN,TEB_INPUT_SNOW,TEB_INPUT_DIR_CO2
29032965

29042966
package nolsmscheme sf_surface_physics==0 - -
29052967
package slabscheme sf_surface_physics==1 - -

‎Registry/registry.dimspec‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -53,6 +53,12 @@ dimspec snly 2 namelist=num_snow_layers z snow_layers
5353
dimspec l 2 namelist=num_soil_layers z soil_layers
5454
dimspec snsl 2 namelist=num_snso_layers z snso_layers
5555

56+
dimspec teb_flly 2 namelist=teb_num_floor_layers z teb_floor_layers
57+
dimspec teb_rdly 2 namelist=teb_num_road_layers z teb_road_layers
58+
dimspec teb_roly 2 namelist=teb_num_roof_layers z teb_roof_layers
59+
dimspec teb_waly 2 namelist=teb_num_wall_layers z teb_wall_layers
60+
dimspec teb_snly 2 namelist=teb_num_snow_layers z teb_snow_layers
61+
5662
dimspec umap0 2 namelist=num_urban_ndm z num_urban_ndm
5763
dimspec umap1 2 namelist=urban_map_zrd z urban_map_zrd
5864
dimspec umap2 2 namelist=urban_map_zwd z urban_map_zwd

‎arch/Config.pl‎

Lines changed: 32 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -15,6 +15,8 @@
1515
$sw_jasperinc_path="";
1616
$sw_esmflib_path="";
1717
$sw_esmfinc_path="";
18+
$sw_teblib_path="";
19+
$sw_tebinc_path="";
1820
$sw_ldflags="";
1921
$sw_compileflags="";
2022
$sw_opt_level="";
@@ -321,6 +323,17 @@
321323
$sw_esmf_ldflag = "yes" ;
322324
}
323325

326+
# A separately-installed TEB library is required to build the WRF-TEB physics option.
327+
# User must set environment variable TEB_PATH to the build folder of TEB
328+
# to enable this feature prior to running configure.
329+
if ( $ENV{WRF_TEB} eq "1" && $ENV{TEB_PATH} )
330+
{
331+
printf "Configuring to use TEB library for WRF-TEB support...\n" ;
332+
printf(" \$TEB_PATH = %s\n",$ENV{TEB_PATH});
333+
$sw_teblib_path = $ENV{TEB_PATH};
334+
$sw_tebinc_path = "$ENV{TEB_PATH}/mod";
335+
}
336+
324337
# parse the configure.wrf file
325338

326339
$validresponse = 0 ;
@@ -556,6 +569,16 @@
556569
$_ =~ s:CONFIGURE_TERRAIN_AND_LANDUSE:-DLANDREAD_STUB=1:g;
557570
}
558571

572+
if ( $sw_teblib_path && $sw_tebinc_path )
573+
{
574+
$_ =~ s:CONFIGURE_TEB_FLAG:-DWRF_TEB:g ;
575+
$_ =~ s:CONFIGURE_TEB_INC:-I$sw_tebinc_path:g ;
576+
}
577+
else
578+
{ $_ =~ s:CONFIGURE_TEB_FLAG::g ;
579+
$_ =~ s:CONFIGURE_TEB_INC::g ;
580+
}
581+
559582
# ESMF substitutions in configure.defaults
560583
if ( $sw_esmflib_path && $sw_esmfinc_path )
561584
{
@@ -869,6 +892,15 @@
869892
$_ =~ s:CONFIGURE_GRIB2_LIB::g ;
870893
}
871894

895+
if ( $sw_teblib_path && $sw_tebinc_path )
896+
{
897+
$_ =~ s:CONFIGURE_TEB_LIB:-L$sw_teblib_path -lteb:g ;
898+
}
899+
else
900+
{
901+
$_ =~ s:CONFIGURE_TEB_LIB::g ;
902+
}
903+
872904
if ( $sw_gpfs_path ne "" )
873905
{ if (/^GPFS.*=/)
874906
{ $_ =~ s/\r|\n//g;

‎arch/postamble‎

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -19,6 +19,7 @@ ARCHFLAGS = $(COREDEFS) -DIWORDSIZE=$(IWORDSIZE) -DDWORDSIZE=$(DWORDSIZ
1919
CONFIGURE_4DVAR_FLAG \
2020
CONFIGURE_TERRAIN_AND_LANDUSE \
2121
CONFIGURE_WAVELET_FLAG \
22+
CONFIGURE_TEB_FLAG \
2223
CONFIGURE_NESTOPT \
2324
-DUSE_ALLOCATABLES \
2425
-Dwrfmodel \
@@ -53,7 +54,9 @@ INCLUDE_MODULES = $(MODULE_SRCH_FLAG) \
5354
-I$(WRF_SRC_ROOT_DIR)/wrftladj \
5455
-I$(WRF_SRC_ROOT_DIR)/chem -I$(WRF_SRC_ROOT_DIR)/inc \
5556
-I$(NETCDFPATH)/include \
56-
CONFIGURE_RTTOV_INC
57+
CONFIGURE_RTTOV_INC \
58+
CONFIGURE_TEB_INC
59+
5760
REGISTRY = Registry
5861
CC_TOOLS_CFLAGS = CONFIGURE_NMM_CORE
5962

‎arch/preamble‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -116,7 +116,7 @@ LIBWRFLIB = libwrflib.a
116116
#NOWIN $(WRF_SRC_ROOT_DIR)/frame/pack_utils.o
117117

118118
#NOWIN LIB_EXTERNAL = \
119-
#NOWIN CONFIGURE_NETCDF_LIB_PATH CONFIGURE_PNETCDF_LIB_PATH CONFIGURE_GRIB2_LIB CONFIGURE_ATMOCN_LIB CONFIGURE_HDF5_LIB_PATH
119+
#NOWIN CONFIGURE_NETCDF_LIB_PATH CONFIGURE_PNETCDF_LIB_PATH CONFIGURE_GRIB2_LIB CONFIGURE_ATMOCN_LIB CONFIGURE_HDF5_LIB_PATH CONFIGURE_TEB_LIB
120120

121121

122122
#### Architecture specific settings ####

‎dyn_em/module_first_rk_step_part1.F‎

Lines changed: 44 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -712,6 +712,50 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
712712
& ,SF_BEP=grid%sf_bep,VL_BEP=grid%vl_bep &
713713
& ,A_E_BEP=grid%a_e_bep,B_E_BEP=grid%b_e_bep,DLG_BEP=grid%dlg_bep &
714714
& ,DL_U_BEP=grid%dl_u_bep &
715+
!
716+
! BEGIN: Optional urban for TEB
717+
& ,teb_test_integration=grid%teb_test_integration &
718+
& ,teb_num_snow_layers=grid%teb_num_snow_layers &
719+
& ,teb_num_road_layers=grid%teb_num_road_layers &
720+
& ,teb_num_roof_layers=grid%teb_num_roof_layers &
721+
& ,teb_num_wall_layers=grid%teb_num_wall_layers &
722+
& ,teb_num_floor_layers=grid%teb_num_floor_layers &
723+
& ,TEB_TI_BLD=grid%TEB_TI_BLD &
724+
& ,TEB_T_CANYON=grid%TEB_T_CANYON, TEB_Q_CANYON=grid%TEB_Q_CANYON &
725+
& ,TEB_WS_ROOF=grid%TEB_WS_ROOF, TEB_WS_ROAD=grid%TEB_WS_ROAD &
726+
& ,TEB_WSNOW_ROOF=grid%TEB_WSNOW_ROOF, TEB_TSNOW_ROOF=grid%TEB_TSNOW_ROOF &
727+
& ,TEB_RSNOW_ROOF=grid%TEB_RSNOW_ROOF, TEB_ASNOW_ROOF=grid%TEB_ASNOW_ROOF &
728+
& ,TEB_TSSNOW_ROOF=grid%TEB_TSSNOW_ROOF, TEB_ESNOW_ROOF=grid%TEB_ESNOW_ROOF &
729+
& ,TEB_WSNOW_ROAD=grid%TEB_WSNOW_ROAD, TEB_TSNOW_ROAD=grid%TEB_TSNOW_ROAD &
730+
& ,TEB_RSNOW_ROAD=grid%TEB_RSNOW_ROAD, TEB_ASNOW_ROAD=grid%TEB_ASNOW_ROAD &
731+
& ,TEB_TSSNOW_ROAD=grid%TEB_TSSNOW_ROAD, TEB_ESNOW_ROAD=grid%TEB_ESNOW_ROAD &
732+
& ,TEB_T_WIN1=grid%TEB_T_WIN1, TEB_T_WIN2=grid%TEB_T_WIN2 &
733+
& ,TEB_AUX_MAX=grid%TEB_AUX_MAX, TEB_THER_PRODC_DAY=grid%TEB_THER_PRODC_DAY &
734+
& ,TEB_QI_BLD=grid%TEB_QI_BLD &
735+
& ,TEB_T_FLOOR=grid%TEB_T_FLOOR, TEB_T_MASS=grid%TEB_T_MASS &
736+
& ,TEB_T_ROAD=grid%TEB_T_ROAD, TEB_T_ROOF=grid%TEB_T_ROOF &
737+
& ,TEB_T_WALL_A=grid%TEB_T_WALL_A, TEB_T_WALL_B=grid%TEB_T_WALL_B &
738+
& ,TEB_HVAC_COOL=grid%TEB_HVAC_COOL, TEB_HVAC_HEAT=grid%TEB_HVAC_HEAT &
739+
& ,TEB_THER_PROD_PANEL=grid%TEB_THER_PROD_PANEL &
740+
& ,TEB_PHOT_PROD_PANEL=grid%TEB_PHOT_PROD_PANEL &
741+
& ,TEB_TSK_RURAL=grid%TEB_TSK_RURAL &
742+
!
743+
! BEGIN: INPUT quantities used for evaluating WRF-TEB with TEB OFFLINE
744+
& ,TEB_INPUT_TA=grid%TEB_INPUT_TA &
745+
& ,TEB_INPUT_PS=grid%TEB_INPUT_PS &
746+
& ,TEB_INPUT_QA=grid%TEB_INPUT_QA &
747+
& ,TEB_INPUT_WIND=grid%TEB_INPUT_WIND &
748+
& ,TEB_INPUT_DIR=grid%TEB_INPUT_DIR &
749+
& ,TEB_INPUT_DIR_SW=grid%TEB_INPUT_DIR_SW &
750+
& ,TEB_INPUT_SCA_SW=grid%TEB_INPUT_SCA_SW &
751+
& ,TEB_INPUT_LW=grid%TEB_INPUT_LW &
752+
& ,TEB_INPUT_RAIN=grid%TEB_INPUT_RAIN &
753+
& ,TEB_INPUT_SNOW=grid%TEB_INPUT_SNOW &
754+
& ,TEB_INPUT_DIR_CO2=grid%TEB_INPUT_DIR_CO2 &
755+
! END: INPUT quantities used for evaluating WRF-TEB with TEB OFFLINE
756+
!
757+
! END: Optional urban for TEB
758+
!
715759
& ,CMR_SFCDIF=grid%cmr_sfcdif, CHR_SFCDIF=grid%chr_sfcdif & !I/O urban
716760
& ,CMC_SFCDIF=grid%cmc_sfcdif, CHC_SFCDIF=grid%chc_sfcdif & !I/O urban
717761
& ,CMGR_SFCDIF=grid%cmgr_sfcdif, CHGR_SFCDIF=grid%chgr_sfcdif & !I/O urban

‎dyn_em/module_initialize_real.F‎

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -2957,7 +2957,8 @@ SUBROUTINE init_domain_rk ( grid &
29572957

29582958
! Split NUDAPT Urban Parameters
29592959

2960-
IF ( ( config_flags%sf_urban_physics == 1 ) .OR. ( config_flags%sf_urban_physics == 2 ) .OR. ( config_flags%sf_urban_physics == 3 ) ) THEN
2960+
IF ( ( config_flags%sf_urban_physics == 1 ) .OR. ( config_flags%sf_urban_physics == 2 ) .OR. ( config_flags%sf_urban_physics == 3 ) .OR. &
2961+
( config_flags%sf_urban_physics == 4 ) ) THEN
29612962
DO j = jts , MIN(jde-1,jte)
29622963
DO i = its , MIN(ide-1,ite)
29632964
IF ( MMINLU == 'NLCD40' .OR. MMINLU == 'MODIFIED_IGBP_MODIS_NOAH') THEN

0 commit comments

Comments
 (0)