Skip to content

WRF-urban updates for the next WRF release #1309

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.

Sign up for GitHub

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

Closed
cenlinhe wants to merge 9 commits into from
Closed

WRF-urban updates for the next WRF release #1309

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
d93bfb8
Include new urban codes for green roof, solar panel, and Local climat…
cenlinhe Oct 30, 2020
b5e99d4
Bug fix for NMM buids
cenlinhe Nov 3, 2020
0b552c1
Bug Fix for urban calculation (line 4370 in module_sf_bep_bem.F)
cenlinhe Nov 3, 2020
e2b4665
Remove unnecessary files e.g., compile.log,...
cenlinhe Nov 3, 2020
b8b2805
Remove compile.log & compile_nmm.log
cenlinhe Nov 3, 2020
3d747e8
Bug Fix for negative qr value in roof temperature calculation in modu…
cenlinhe Nov 5, 2020
0c389c0
Add new registry urban variable to the bep_bem package; The original …
cenlinhe Nov 5, 2020
1079049
Add urban category namelist flag for using new LCZ and error message …
cenlinhe Nov 17, 2020
39ab23a
Bug fix for calling wrf error message in the LCZ inconsistency checki…
cenlinhe Nov 18, 2020
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
21 changes: 19 additions & 2 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -885,7 +885,7 @@ state real TGLEV_URB3D i{umap8}j misc 1 Z r "T
state real TFLEV_URB3D i{umap9}j misc 1 Z r "TFLEV_URB3D" "FLOOR TEMPERATURE" "K"
state real SF_AC_URB3D ij misc 1 - r "SF_AC_URB3D" "SENSIBLE HEAT FLUX FROM THE AIR COND." "W m{-2}"
state real LF_AC_URB3D ij misc 1 - r "LF_AC_URB3D" "LATENT HEAT FLUX FROM THE AIR COND." "W m{-2}"
state real CM_AC_URB3D ij misc 1 - r "CM_AC_URB3D" "CONSUMPTION OF THE AIR COND." "W m{-2}"
state real CM_AC_URB3D ij misc 1 - rh "CM_AC_URB3D" "CONSUMPTION OF THE AIR COND." "W m{-2}"
state real SFVENT_URB3D ij misc 1 - r "SFVENT_URB3D" "SENSIBLE HEAT FLUX FROM URBAN VENTILATION" "W m{-2}"
state real LFVENT_URB3D ij misc 1 - r "LFVENT_URB3D" "LATENT HEAT FLUX FROM URBAN VENTILATION" "W m{-2}"
state real SFWIN1_URB3D i{umap7}j misc 1 Z r "SFWIN1_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC WINDOW" "W m{-2}"
Expand All @@ -894,6 +894,20 @@ state real SFW1_URB3D i{umap4}j misc 1 Z r "S
state real SFW2_URB3D i{umap4}j misc 1 Z r "SFW2_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real SFR_URB3D i{umap5}j misc 1 Z r "SFR_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real SFG_URB3D i{umap0}j misc 1 Z r "SFG_URB3D" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real EP_PV_URB3D ij misc 1 - rh "EP_PV_URB3D" "ELEC. PRODUCTION OF ROOFTOP PV PANELS" "W m{-2}" !PVP
Copy link
Contributor

@davegill davegill Nov 5, 2020

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

All of these variables need to be in a package

Copy link
Contributor Author

@cenlinhe cenlinhe Nov 5, 2020

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks! will add to a package in an updated version.

state real T_PV_URB3D i{umap5}j misc 1 - r "T_PV_URB3D" "PHOTOVOLTAIC PANELS TEMPERATURE " "K" !PVP
state real TRV_URB4D i{umap10}j misc 1 Z r "TRV_URB4D" "GREEN ROOF LAYER TEMPERATURE" "K"
state real QR_URB4D i{umap10}j misc 1 Z r "QR_URB4D" "GREEN ROOF LAYER MOISTURE" "dimensionless"
state real QGR_URB3D ij misc 1 Z rh "QGR_URB4D" "GREEN ROOF LAYER MOISTURE OUTPUT" "dimensionless"
state real TGR_URB3D ij misc 1 Z rh "TGR_URB4D" "GREEN ROOF LAYER TEMPERATURE OUTPUT" "dimensionless"
state real DRAIN_URB4D i{umap5}j misc 1 Z r "DRAIN_URB4D" "GREEN ROOF DRAINAGE" "mm"
state real DRAINGR_URB3D ij misc 1 - rh "DRAINGR_URB3D" "ACCUMULATED GREEN ROOF DRAINAGE" "mm"
state real SFRV_URB3D i{umap5}j misc 1 Z r "SFR_URB3D" "SENSIBLE HEAT FLUX FROM GREEN ROOF" "W m{-2}"
state real LFRV_URB3D i{umap5}j misc 1 Z r "LFR_URB3D" "LATENT HEAT FLUX FROM GREEN ROOF" "W m{-2}"
state real DGR_URB3D i{umap5}j misc 1 Z r "DGR_URB3D" "ROOF LAYER DEPTH WATER RETENTION" "mm"
state real DG_URB3D i{umap0}j misc 1 Z r "DG_URB4D" "ROOF LAYER DEPTH WATER RETENTION" "mm"
state real LFR_URB3D i{umap5}j misc 1 Z r "LFR_URB3D" "LATENT HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real LFG_URB3D i{umap0}j misc 1 Z r "LFG_URB3D" "LATENT HEAT FLUX FROM URBAN SFC" "W m{-2}"
state real CMR_SFCDIF ij misc 1 - r "CMR_SFCDIF" "" ""
state real CHR_SFCDIF ij misc 1 - r "CHR_SFCDIF" "" ""
state real CMC_SFCDIF ij misc 1 - r "CMC_SFCDIF" "" ""
Expand Down Expand Up @@ -2328,6 +2342,7 @@ rconfig integer num_urban_ngb derived 1 1
rconfig integer num_urban_nf derived 1 1 irh "num_urban_nf" "number of grid levels in the floors (nf_u in BEM header)" ""
rconfig integer num_urban_nz derived 1 1 irh "num_urban_nz" "maximum number of vertical levels in the urban grid (nz_um in BEP or BEM header)" ""
rconfig integer num_urban_nbui derived 1 1 irh "num_urban_nbui" "maximum number of types of buildings in an urban class (nbui_max in BEM header)" ""
rconfig integer num_urban_ngr derived 1 1 irh "num_urban_ngr" "maximum number of types of buildings in an urban class (ngr_u in BEP or BEM header)" ""

rconfig integer urban_map_zrd derived 1 1 rh "urban_map_zrd" "urban mapping 1: ind_zrd" ""
rconfig integer urban_map_zwd derived 1 1 rh "urban_map_zwd" "urban mapping 2: ind_zwd" ""
Expand All @@ -2338,6 +2353,7 @@ rconfig integer urban_map_bd derived 1 1
rconfig integer urban_map_wd derived 1 1 rh "urban_map_wd" "urban mapping 7: ind_wd" ""
rconfig integer urban_map_gbd derived 1 1 rh "urban_map_gbd" "urban mapping 8: ind_gbd" ""
rconfig integer urban_map_fbd derived 1 1 rh "urban_map_fbd" "urban mapping 9: ind_fbd" ""
rconfig integer urban_map_zgrd derived 1 1 rh "urban_map_zgrd" "urban mapping 10: ind_grd" ""
Copy link
Contributor

@davegill davegill Nov 2, 2020

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

the variables urban_map_{zrd, zwd, gd, zd, zdf, bd, wd, gbd, fdb} are all defined in the NMM Registry. To be conventional, urban_map_zgrd should also be there.


rconfig integer num_urban_hi namelist,physics 1 15 irh "num_urban_hi" "" ""
rconfig integer num_months namelist,physics 1 12 irh "num_months" "" ""
Expand All @@ -2358,6 +2374,7 @@ rconfig integer imomentum namelist,physics 1 0
rconfig integer ishallow namelist,physics 1 0 rh "ishallow" "shallow convection in G3 scheme" ""
rconfig real convtrans_avglen_m namelist,physics 1 30 rh "convtrans_avglen_m" "averaging time for convective transport output variables (minutes)" ""
rconfig integer num_land_cat namelist,physics 1 21 - "num_land_cat" "" ""
rconfig integer use_wudapt_lcz namelist,physics 1 0 - "use_wudapt_lcz" "" ""
rconfig integer num_soil_cat namelist,physics 1 16 - "num_soil_cat" "" ""
rconfig integer mp_zero_out namelist,physics 1 0 - "mp_zero_out" "microphysics fields set to zero 0=no action taken, 1=all fields but Qv, 2=all fields including Qv" "flag"
rconfig real mp_zero_out_thresh namelist,physics 1 1.e-8 - "mp_zero_out_thresh" "minimum threshold for non-Qv moist fields, below are set to zero" "kg/kg"
Expand Down Expand Up @@ -2901,7 +2918,7 @@ package sfclayscheme sf_sfclay_physics==91 - -

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
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
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
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,ep_pv_urb3d,t_pv_urb3d,trv_urb4d,qr_urb4d,qgr_urb3d,tgr_urb3d,drain_urb4d,draingr_urb3d,sfrv_urb3d,lfrv_urb3d,dgr_urb3d,dg_urb3d,lfr_urb3d,lfg_urb3d

package nolsmscheme sf_surface_physics==0 - -
package slabscheme sf_surface_physics==1 - -
Expand Down
3 changes: 3 additions & 0 deletions Registry/Registry.NMM
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1463,6 +1463,7 @@ rconfig integer num_urban_ngb namelist,physics 1 1
rconfig integer num_urban_nf namelist,physics 1 1 irh "num_urban_nf" "number of grid levels in the floors (nf_u in BEM header)" ""
rconfig integer num_urban_nz namelist,physics 1 1 irh "num_urban_nz" "maximum number of vertical levels in the urban grid (nz_um in BEP or BEM header)" ""
rconfig integer num_urban_nbui namelist,physics 1 1 irh "num_urban_nbui" "maximum number of types of buildings in an urban class (nbui_max in BEM header)" ""
rconfig integer num_urban_ngr namelist,physics 1 1 irh "num_urban_ngr" "maximum number of types of buildings in an urban class (ngr_u in BEP or BEM header)" ""

rconfig integer urban_map_zrd derived 1 1 rh "urban_map_zrd" "urban mapping 1: ind_zrd" ""
rconfig integer urban_map_zwd derived 1 1 rh "urban_map_zwd" "urban mapping 2: ind_zwd" ""
Expand All @@ -1473,6 +1474,7 @@ rconfig integer urban_map_bd derived 1 1
rconfig integer urban_map_wd derived 1 1 rh "urban_map_wd" "urban mapping 7: ind_wd" ""
rconfig integer urban_map_gbd derived 1 1 rh "urban_map_gbd" "urban mapping 8: ind_gbd" ""
rconfig integer urban_map_fbd derived 1 1 rh "urban_map_fbd" "urban mapping 9: ind_fbd" ""
rconfig integer urban_map_zgrd derived 1 1 rh "urban_map_zgrd" "urban mapping 10: ind_grd" ""

rconfig integer sf_surface_mosaic namelist,physics 1 0 rh "sf_surface_mosaic" "1= mosaic, 0=no mosaic method, add by danli" ""
rconfig integer mosaic_cat namelist,physics 1 3 rh "mosaic_cat" "works when sf_surface_mosaic=1; it is the number of mosaic tiles" ""
Expand All @@ -1487,6 +1489,7 @@ rconfig integer maxens3 namelist,physics 1 16
rconfig integer ensdim namelist,physics 1 144 irh "ensdim" "" ""
rconfig integer chem_opt namelist,physics max_domains 0 rh "chem_opt" "" ""
rconfig integer num_land_cat namelist,physics 1 24 - "num_land_cat" "" ""
rconfig integer use_wudapt_lcz namelist,physics 1 0 - "use_wudapt_lcz" "" ""
rconfig integer num_soil_cat namelist,physics 1 16 - "num_soil_cat" "" ""
rconfig integer topo_wind namelist,physics max_domains 0 - "topo_wind" "2: Use Mass sfc drag scheme, 1: improve effects topography over surface wind, 0:not" ""
rconfig integer mp_zero_out namelist,physics 1 0 - "mp_zero_out" "microphysics fields set to zero 0=no action taken, 1=all fields but Qv, 2=all fields including Qv" "flag"
Expand Down
1 change: 1 addition & 0 deletions Registry/registry.dimspec
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -63,6 +63,7 @@ dimspec umap6 2 namelist=urban_map_bd z urban_map_bd
dimspec umap7 2 namelist=urban_map_wd z urban_map_wd
dimspec umap8 2 namelist=urban_map_gbd z urban_map_gbd
dimspec umap9 2 namelist=urban_map_fbd z urban_map_fbd
dimspec umap10 2 namelist=urban_map_zgrd z urban_map_zgrd

dimspec mocat 2 namelist=mosaic_cat z mosaic categories # danli
dimspec mocat2 2 namelist=mosaic_cat_soil z mosaic categories * soil layers # danli
Expand Down
42 changes: 29 additions & 13 deletions dyn_em/module_first_rk_step_part1.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -668,6 +668,7 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,RN_URB2D=grid%rn_urb2d , TS_URB2D=grid%ts_urb2d & !H urban
& ,FRC_URB2D=grid%frc_urb2d & !H urban
& ,UTYPE_URB2D=grid%utype_urb2d & !H urban
& ,SWDDIR=grid%swddir,SWDDIF=grid%swddif & !multi-layer urban _gl
! Optional urban for BEP scheme
& ,SF_URBAN_PHYSICS=config_flags%sf_urban_physics &
& ,num_urban_ndm = config_flags%num_urban_ndm & !multi-layer urban
Expand All @@ -680,25 +681,40 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,urban_map_wd = config_flags%urban_map_wd & !multi-layer urban
& ,urban_map_gbd = config_flags%urban_map_gbd & !multi-layer urban
& ,urban_map_fbd = config_flags%urban_map_fbd & !multi-layer urban
& ,urban_map_zgrd = config_flags%urban_map_zgrd & !multi-layer urban
& ,NUM_URBAN_HI=config_flags%num_urban_hi & !multi-layer urban
& ,TSK_RURAL=grid%tsk_rural & !multi-layer urban
& ,TRB_URB4D=grid%trb_urb4d,TW1_URB4D=grid%tw1_urb4d & !multi-layer urban
& ,TW2_URB4D=grid%tw2_urb4d,TGB_URB4D=grid%tgb_urb4d & !multi-layer urban
& ,TLEV_URB3D=grid%tlev_urb3d & !multi-layer urban
& ,QLEV_URB3D=grid%qlev_urb3d & !multi-layer urban
& ,TW1LEV_URB3D=grid%tw1lev_urb3d & !multi-layer urban
& ,TW2LEV_URB3D=grid%tw2lev_urb3d & !multi-layer urban
& ,TGLEV_URB3D=grid%tglev_urb3d & !multi-layer urban
& ,TFLEV_URB3D=grid%tflev_urb3d & !multi-layer urban
& ,SF_AC_URB3D=grid%sf_ac_urb3d & !multi-layer urban
& ,LF_AC_URB3D=grid%lf_ac_urb3d & !multi-layer urban
& ,CM_AC_URB3D=grid%cm_ac_urb3d & !multi-layer urban
& ,SFVENT_URB3D=grid%sfvent_urb3d & !multi-layer urban
& ,LFVENT_URB3D=grid%lfvent_urb3d & !multi-layer urban
& ,SFWIN1_URB3D=grid%sfwin1_urb3d & !multi-layer urban
& ,SFWIN2_URB3D=grid%sfwin2_urb3d & !multi-layer urban
& ,TLEV_URB3D=grid%tlev_urb3d & !multi-layer urban
& ,QLEV_URB3D=grid%qlev_urb3d & !multi-layer urban
& ,TW1LEV_URB3D=grid%tw1lev_urb3d & !multi-layer urban
& ,TW2LEV_URB3D=grid%tw2lev_urb3d & !multi-layer urban
& ,TGLEV_URB3D=grid%tglev_urb3d & !multi-layer urban
& ,TFLEV_URB3D=grid%tflev_urb3d & !multi-layer urban
& ,SF_AC_URB3D=grid%sf_ac_urb3d & !multi-layer urban
& ,LF_AC_URB3D=grid%lf_ac_urb3d & !multi-layer urban
& ,CM_AC_URB3D=grid%cm_ac_urb3d & !multi-layer urban
& ,SFVENT_URB3D=grid%sfvent_urb3d & !multi-layer urban
& ,LFVENT_URB3D=grid%lfvent_urb3d & !multi-layer urban
& ,SFWIN1_URB3D=grid%sfwin1_urb3d & !multi-layer urban
& ,SFWIN2_URB3D=grid%sfwin2_urb3d & !multi-layer urban
& ,SFW1_URB3D=grid%sfw1_urb3d,SFW2_URB3D=grid%sfw2_urb3d & !multi-layer urban
& ,SFR_URB3D=grid%sfr_urb3d,SFG_URB3D=grid%sfg_urb3d & !multi-layer urban
& ,EP_PV_URB3D=grid%ep_pv_urb3d & !GRZ
& ,T_PV_URB3D=grid%t_pv_urb3d & !GRZ
& ,TRV_URB4D=grid%trv_urb4d & !GRZ
& ,QR_URB4D=grid%qr_urb4d & !GRZ
& ,QGR_URB3D=grid%qgr_urb3d & !GRZ
& ,TGR_URB3D=grid%tgr_urb3d & !GRZ
& ,DRAIN_URB4D=grid%drain_urb4d & !GRZ
& ,DRAINGR_URB3D=grid%draingr_urb3d & !GRZ
& ,SFRV_URB3D=grid%sfrv_urb3d & !GRZ
& ,LFRV_URB3D=grid%lfrv_urb3d & !GRZ
& ,DGR_URB3D=grid%dgr_urb3d & !GRZ
& ,DG_URB3D=grid%dg_urb3d & !GRZ
& ,LFR_URB3D=grid%lfr_urb3d & !GRZ
& ,LFG_URB3D=grid%lfg_urb3d & !GRZ
& ,LP_URB2D=grid%lp_urb2d,HI_URB2D=grid%hi_urb2d & !multi-layer urban
& ,LB_URB2D=grid%lb_urb2d,HGT_URB2D=grid%hgt_urb2d & !multi-layer urban
& ,MH_URB2D=grid%mh_urb2d,STDH_URB2D=grid%stdh_urb2d & !SLUCM
Expand Down
Loading