Skip to content

urban updates: urb solar panel, green roof, building drag coeff; local climate zone #1385

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
andreazonato wants to merge 4 commits into from
Closed

urban updates: urb solar panel, green roof, building drag coeff; local climate zone #1385

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
a33b5f3
implementation of LCZ, Green Roof, Photovoltaic panels, new building …
Jan 26, 2021
dc4a843
Solved bugs on use_wudapt_lcz
Jan 28, 2021
9754795
Solved - issue on Registry.EM_COMMON
Jan 29, 2021
8cac4ef
Modification of nz_um=18,nbu_max=15
andreazonato Feb 25, 2021
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
30 changes: 24 additions & 6 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -894,6 +894,21 @@ 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
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 @@ -957,9 +972,9 @@ state real WSTAR_YSU ij misc 1 - - "WSTAR
state real DELTA_YSU ij misc 1 - - "DELTA_YSU" "entrainment layer depth from ysupbl" "m"

# MYJ PBL variables; GBM PBL: EXCH_H, EXCH_M
state real EXCH_H ikj misc 1 Z r "EXCH_H" "SCALAR EXCHANGE COEFFICIENTS " "m2 s-1"
state real EXCH_M ikj misc 1 Z r "EXCH_M" "EXCHANGE COEFFICIENTS " "m2 s-1"
state real CT ij misc 1 - r "CT" "COUNTERGRADIENT TERM" "K"
state real EXCH_H ikj misc 1 Z r "EXCH_H" "SCALAR EXCHANGE COEFFICIENTS " "m2 s-1"
state real EXCH_M ikj misc 1 Z r "EXCH_M" "EXCHANGE COEFFICIENTS " "m2 s-1"
state real CT ij misc 1 - r "CT" "COUNTERGRADIENT TERM" "K"
state real THZ0 ij misc 1 - r "THZ0" "POTENTIAL TEMPERATURE AT ZNT" "K"
state real Z0 ij misc 1 - r "Z0" "Background ROUGHNESS LENGTH" "m"
state real QZ0 ij misc 1 - r "QZ0" "SPECIFIC HUMIDITY AT ZNT" "kg kg-1"
Expand Down Expand Up @@ -2328,6 +2343,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" "number of grid levels in the green roof (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,7 +2354,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_zgrd" ""
rconfig integer num_urban_hi namelist,physics 1 15 irh "num_urban_hi" "" ""
rconfig integer num_months namelist,physics 1 12 irh "num_months" "" ""
rconfig integer sf_surface_mosaic namelist,physics 1 0 rh "sf_surface_mosaic" "1= mosaic, 0=no mosaic method, add by danli" ""
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,9 +2918,10 @@ 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 nolsmscheme sf_surface_physics==0
Comment thread
andreazonato marked this conversation as resolved.
Outdated
package slabscheme sf_surface_physics==1 - -
package lsmscheme sf_surface_physics==2 - state:flx4,fvb,fbur,fgsn,smcrel,xlaidyn
package ruclsmscheme sf_surface_physics==3 - state:smfr3d,keepfr3dflag,soilt1,rhosnf,snowfallac,precipfr,acrunoff
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 derived 1 1 irh "num_urban_ngr" "number of grid levels in the green roof (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_zgrd" ""

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
16 changes: 16 additions & 0 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,6 +681,7 @@ 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
Expand All @@ -699,6 +701,20 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
& ,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
13 changes: 0 additions & 13 deletions dyn_em/module_initialize_real.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2983,19 +2983,6 @@ SUBROUTINE init_domain_rk ( grid &
IF ( ( config_flags%sf_urban_physics == 1 ) .OR. ( config_flags%sf_urban_physics == 2 ) .OR. ( config_flags%sf_urban_physics == 3 ) ) THEN
DO j = jts , MIN(jde-1,jte)
DO i = its , MIN(ide-1,ite)
IF ( MMINLU == 'NLCD40' .OR. MMINLU == 'MODIFIED_IGBP_MODIS_NOAH') THEN
IF ( grid%FRC_URB2D(i,j) .GE. 0.5 .AND. &
(grid%ivgtyp(i,j).NE.13 .AND. grid%ivgtyp(i,j).NE.24 .AND. grid%ivgtyp(i,j).NE.25 .AND. grid%ivgtyp(i,j).NE.26)) grid%ivgtyp(i,j)=13
ELSE IF ( MMINLU == "USGS" ) THEN
IF ( grid%FRC_URB2D(i,j) .GE. 0.5 .AND. &
grid%ivgtyp(i,j).NE.1 ) grid%ivgtyp(i,j)=1
ENDIF

IF ( grid%FRC_URB2D(i,j) == 0. ) THEN
IF ( (MMINLU == 'NLCD40' .OR. MMINLU == 'MODIFIED_IGBP_MODIS_NOAH') .AND. &
(grid%ivgtyp(i,j)==24 .OR. grid%ivgtyp(i,j)==25 .OR. grid%ivgtyp(i,j)==26 .OR. grid%ivgtyp(i,j)==13) ) grid%FRC_URB2D(i,j) = 0.9
IF ( MMINLU == 'USGS' .AND. grid%ivgtyp(i,j)==1 ) grid%FRC_URB2D(i,j) = 0.9
ENDIF
IF ( skip_middle_points_t ( ids , ide , jds , jde , i , j , em_width , hold_ups ) ) CYCLE
grid%LP_URB2D(i,j) = grid%URB_PARAM(i,91,j)
grid%LB_URB2D(i,j) = grid%URB_PARAM(i,95,j)
Expand Down
Loading