urban updates: solar panel, green roof, building drag coefficient#1427
urban updates: solar panel, green roof, building drag coefficient#1427davegill merged 9 commits intowrf-model:developfrom
Conversation
…drag coefficient, diffuse and direct radiation and building materials permeability in BEP+BEM /sf_urban_physics=3)
Conflicts: Registry/Registry.EM_COMMON dyn_em/module_initialize_real.F phys/module_physics_init.F phys/module_sf_bep.F phys/module_sf_bep_bem.F phys/module_sf_noahdrv.F phys/module_sf_noahmpdrv.F phys/module_sf_noahmplsm.F phys/module_sf_urban.F phys/module_surface_driver.F run/URBPARM.TBL run/URBPARM_LCZ.TBL Changes to be committed: modified: Registry/Registry.EM_COMMON modified: Registry/Registry.NMM modified: Registry/registry.dimspec modified: dyn_em/module_first_rk_step_part1.F modified: dyn_em/start_em.F modified: dyn_nmm/module_PHYSICS_CALLS.F modified: dyn_nmm/start_domain_nmm.F modified: phys/module_physics_init.F modified: phys/module_sf_bem.F modified: phys/module_sf_bep_bem.F modified: phys/module_sf_noahdrv.F modified: phys/module_sf_noahmpdrv.F modified: phys/module_sf_urban.F modified: phys/module_surface_driver.F modified: run/URBPARM.TBL modified: run/URBPARM_LCZ.TBL modified: share/module_check_a_mundo.F
|
jenkins |
|
@andreazonato @cenlinhe @dudhia @weiwangncar |
Yes, I think that this PR is what we want. It passes the standard jenkins tests. However, we need to see that the results are physically reasonable when the new features are activated. We would like to cut off the release branch from the develop branch this week, so testing is a priority. |
|
I will try to do tests for this PR code tonight and tomorrow morning. |
|
Dear all,
regarding the test of te physics, about it we have a paper submitted to
ACP, and still on review.
Wanna me to share this with you as "proof" of the code working?
Il giorno gio 4 mar 2021 alle ore 17:15 dudhia <[email protected]>
ha scritto:
… ***@***.**** commented on this pull request.
------------------------------
In phys/module_sf_bem.F
<#1427 (comment)>:
> do irof=1,nrof
trof(irof)=trof1D(irof)
- end do
+
added blank line
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1427 (review)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ALWNNTTX4KKECGLMRRDMUUTTB6WZJANCNFSM4YSJUKEA>
.
|
|
@davegill @andreazonato Hi Andrea, I think you can put the reference here and ask Dave to include it in the release note of the PR message above. Although it is under review, you can still put the authors, paper title, journal name, and "under review". |
phys/module_sf_bem.F
Outdated
| real em_av !averadge emissivity | ||
| real twal_int(6) !twalint | ||
| real em_av !averadge emissivity | ||
| real twal_int(6) !twalint |
There was a problem hiding this comment.
Jimy,
I lined up these also. There were tabs vs spaces used here. I have followed the existing convention, and they should line up nicely now.
There was a problem hiding this comment.
Jimy,
Here is what the source code looks like in an editor on my Mac:
integer iw,jw
real b_wall(6)
real b_wind(6)
real emwal_av !averadge emissivity of the wall
real emwin_av !averadge emissivity of the window
real em_av !averadge emissivity
real twal_int(6) !twalint
real twin(4) !twinint
The github representation is playing with tabs and spaces. My editor will not allow me to introduce a tab character in the source code.
|
I have completed my code pass through and left all my notes in it.
…On Thu, Mar 4, 2021 at 10:04 AM Dave Gill ***@***.***> wrote:
***@***.**** commented on this pull request.
------------------------------
In phys/module_surface_driver.F
<#1427 (comment)>:
> @@ -2966,7 +3001,7 @@ SUBROUTINE surface_driver( &
IF(SF_URBAN_PHYSICS.eq.1) THEN
DO j=j_start(ij),j_end(ij) !urban
DO i=i_start(ij),i_end(ij) !urban
- IF( IVGTYP(I,J) == ISURBAN .or. IVGTYP(I,J) == LCZ_1 .or. IVGTYP(I,J) == LCZ_2 .or. &
+ IF(IVGTYP(I,J) == ISURBAN .or. IVGTYP(I,J) == LCZ_1 .or. IVGTYP(I,J) == LCZ_2 .or. &
More efforts to clean up the bad indenting in the original code, in the
areas that we specifically put in during the conflict resolution
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1427 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/AEIZ77EOGTUK2FHHBJPPXK3TB64RLANCNFSM4YSJUKEA>
.
|
|
@davegill @dudhia @weiwangncar I have done a few tests using this PR code and all urban physics options work well. I checked key land surface variables (e.g., temperature, heat flux, radiation, etc.) and the results look reasonable: /glade/scratch/cenlinhe/WRF_noahmp_test/WRF_LCZ_urban1427/run/urban*/ |
|
This is the under-review paper Andrea mentioned which discussed this urban code update here: A. Zonato, A. Martilli, E. Gutierrez, F. Chen, C. He, M. Barlage, D. Zardi, and L. Giovannini (2021): Exploring the effects of rooftop mitigation strategies on urban temperatures and energy consumption, Atmospheric Chemistry and Physics, under review. |
|
@cenlinhe @dudhia @weiwangncar
Cenlin, Wei and Jimy, |
|
I have now resolved all the ones except where I think the white space
change made it worse and where there was no specific response comment.
…On Thu, Mar 4, 2021 at 1:10 PM Dave Gill ***@***.***> wrote:
@cenlinhe <https://github.com/cenlinhe> @dudhia
<https://github.com/dudhia> @weiwangncar <https://github.com/weiwangncar>
I have done a few tests using this PR code and all urban physics options
work well. I checked key land surface variables (e.g., temperature, heat
flux, radiation, etc.) and the results look reasonable
Cenlin,
Thanks for jumping on this so quickly. I think that we are now good to go.
Wei and Jimy,
The comments so far about incorrect spaces (etc) were on purpose to align
the columns and fix indenting in the original code. This is ready for a
review.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1427 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/AEIZ77H6N2Z46PBYVVORMLDTB7SLDANCNFSM4YSJUKEA>
.
|
modified: phys/module_sf_bem.F modified: phys/module_sf_bep_bem.F
|
Thanks. I had one comment left in the BEP-BEM module. |
|
Actually I had several in that module related to do..enddo alignments.
…On Thu, Mar 4, 2021 at 2:34 PM Dave Gill ***@***.***> wrote:
***@***.**** commented on this pull request.
------------------------------
In phys/module_sf_bem.F
<#1427 (comment)>:
> @@ -1330,8 +1585,8 @@ subroutine algebra_long(emwal,emwin,sigma,twalint,twinint,&
real b_wind(6)
real emwal_av !averadge emissivity of the wall
real emwin_av !averadge emissivity of the window
- real em_av !averadge emissivity
- real twal_int(6) !twalint
+ real em_av !averadge emissivity
+ real twal_int(6) !twalint
Jimy,
Here is what the source code looks like in an editor on my Mac:
integer iw,jw
real b_wall(6)
real b_wind(6)
real emwal_av !averadge emissivity of the wall
real emwin_av !averadge emissivity of the window
real em_av !averadge emissivity
real twal_int(6) !twalint
real twin(4) !twinint
The github representation is playing with tabs and spaces. My editor will
not allow me to introduce a tab character in the source code.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1427 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/AEIZ77C6JQBJHTISHHPLMI3TB74HZANCNFSM4YSJUKEA>
.
|
|
@dudhia |
|
I mean bep_bem
…On Thu, Mar 4, 2021 at 2:38 PM Jimy Dudhia ***@***.***> wrote:
Actually I had several in that module related to do..enddo alignments.
On Thu, Mar 4, 2021 at 2:34 PM Dave Gill ***@***.***> wrote:
> ***@***.**** commented on this pull request.
> ------------------------------
>
> In phys/module_sf_bem.F
> <#1427 (comment)>:
>
> > @@ -1330,8 +1585,8 @@ subroutine algebra_long(emwal,emwin,sigma,twalint,twinint,&
> real b_wind(6)
> real emwal_av !averadge emissivity of the wall
> real emwin_av !averadge emissivity of the window
> - real em_av !averadge emissivity
> - real twal_int(6) !twalint
> + real em_av !averadge emissivity
> + real twal_int(6) !twalint
>
> Jimy,
> Here is what the source code looks like in an editor on my Mac:
>
> integer iw,jw
> real b_wall(6)
> real b_wind(6)
> real emwal_av !averadge emissivity of the wall
> real emwin_av !averadge emissivity of the window
> real em_av !averadge emissivity
> real twal_int(6) !twalint
> real twin(4) !twinint
>
> The github representation is playing with tabs and spaces. My editor will
> not allow me to introduce a tab character in the source code.
>
> —
> You are receiving this because you were mentioned.
> Reply to this email directly, view it on GitHub
> <#1427 (comment)>, or
> unsubscribe
> <https://github.com/notifications/unsubscribe-auth/AEIZ77C6JQBJHTISHHPLMI3TB74HZANCNFSM4YSJUKEA>
> .
>
|
|
The lines in bep_bem left with my comments are 796, 844, 878(?), 1673
…On Thu, Mar 4, 2021 at 2:39 PM Dave Gill ***@***.***> wrote:
@dudhia <https://github.com/dudhia>
Several of the suggestions for indenting that you had that actually put
the code back to where it was supposed to be, and removed that line from
any diff. There were a couple of these in BEM, and the newer indenting
suggestion that you pointed out in BEP_BEM.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
<#1427 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/AEIZ77GIFQ6UECQ4T22SRELTB74ZZANCNFSM4YSJUKEA>
.
|
…f-model#1427) KEYWORDS: urban, green roof, solar panel, new building drag coefficient SOURCE: Andrea Zonato (University of Trento, Italy, [email protected]) DESCRIPTION OF CHANGES: This update includes new capabilities to account for urban solar panel, green roof, and new building drag coefficient. These three new capabilities are for BEP+BEM. The new building drag coefficient is based on Santiago and Martilli (2010) and Guti ́errez et al. (2015). More details are available in the WRF User's Guide. LIST OF MODIFIED FILES: M Registry/Registry.EM_COMMON M Registry/Registry.NMM M Registry/registry.dimspec M dyn_em/module_first_rk_step_part1.F M dyn_em/start_em.F M dyn_nmm/module_PHYSICS_CALLS.F M dyn_nmm/start_domain_nmm.F M phys/module_physics_init.F M phys/module_sf_bem.F M phys/module_sf_bep_bem.F M phys/module_sf_noahdrv.F M phys/module_sf_noahmpdrv.F M phys/module_sf_urban.F M phys/module_surface_driver.F M run/URBPARM.TBL M run/URBPARM_LCZ.TBL M share/module_check_a_mundo.F TESTS CONDUCTED: 1. The model compilation has been tested successfully (tested in Cheyenne by both the developer and Cenlin He, NCAR) 2. The model simulation has been tested successfully (tested in Cheyenne by the developer) 3. The model simulation has been tested successfully (tested in WRF Docker by Cenlin He, NCAR) 4. The jenkins testing is all PASS The original PR was split into two pieces. Here is a figure from the original PR (left) and the combination of the two PRs (right). We want to see nearly the exact same features. This is a figure of 2-meter temperature around the Denver area.  RELEASE NOTE: WRF-urban updates for green roof, solar panel, and new building drag coefficient for BEP+BEM. The paper describing the green roof update is under development. The new building drag coefficient is based on Santiago and Martilli (2010) and Guti ́errez et al. (2015). Santiago, J. L. and Martilli, A. (2010). A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier-Stokes Microscale Simulations. Boundary- Layer Meteorology, 137(3):417–439. Guti ́errez, E., Martilli, A., Santiago, J. L., and Gonz ́alez, J. E. (2015). A Mechanical Drag Coefficient Formulation and Urban Canopy Parameter Assimilation Technique for Complex Urban Environments. Boundary-Layer Meteorology, 157(2):333–341. A. Zonato, A. Martilli, E. Gutierrez, F. Chen, C. He, M. Barlage, D. Zardi, and L. Giovannini (2021): Exploring the effects of rooftop mitigation strategies on urban temperatures and energy consumption, Atmospheric Chemistry and Physics, under review.
KEYWORDS: urban, green roof, solar panel, new building drag coefficient
SOURCE: Andrea Zonato (University of Trento, Italy, [email protected])
DESCRIPTION OF CHANGES:
This update includes new capabilities to account for urban solar panel, green roof, and new building drag coefficient. These three new capabilities are for BEP+BEM. The new building drag coefficient is based on Santiago and Martilli (2010) and Guti ́errez et al. (2015). More details are available in the WRF User's Guide.
LIST OF MODIFIED FILES:
M Registry/Registry.EM_COMMON
M Registry/Registry.NMM
M Registry/registry.dimspec
M dyn_em/module_first_rk_step_part1.F
M dyn_em/start_em.F
M dyn_nmm/module_PHYSICS_CALLS.F
M dyn_nmm/start_domain_nmm.F
M phys/module_physics_init.F
M phys/module_sf_bem.F
M phys/module_sf_bep_bem.F
M phys/module_sf_noahdrv.F
M phys/module_sf_noahmpdrv.F
M phys/module_sf_urban.F
M phys/module_surface_driver.F
M run/URBPARM.TBL
M run/URBPARM_LCZ.TBL
M share/module_check_a_mundo.F
TESTS CONDUCTED:
The original PR was split into two pieces. Here is a figure from the original PR (left) and the combination of the two PRs
(right). We want to see nearly the exact same features. This is a figure of 2-meter temperature around the Denver area.
RELEASE NOTE: WRF-urban updates for green roof, solar panel, and new building drag coefficient for BEP+BEM. The paper describing the green roof update is under development. The new building drag coefficient is based on Santiago and Martilli (2010) and Guti ́errez et al. (2015).
Santiago, J. L. and Martilli, A. (2010). A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier-Stokes Microscale Simulations. Boundary- Layer Meteorology, 137(3):417–439.
Guti ́errez, E., Martilli, A., Santiago, J. L., and Gonz ́alez, J. E. (2015). A Mechanical Drag Coefficient Formulation and Urban Canopy Parameter Assimilation Technique for Complex Urban Environments. Boundary-Layer Meteorology, 157(2):333–341.
A. Zonato, A. Martilli, E. Gutierrez, F. Chen, C. He, M. Barlage, D. Zardi, and L. Giovannini (2021): Exploring the effects of rooftop mitigation strategies on urban temperatures and energy consumption, Atmospheric Chemistry and Physics, under review.