Skip to content

MYNN PBL clouds in photolysis option 4 (TUV) #1622

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

Merged
davegill merged 7 commits into from
Jan 19, 2022
Merged

MYNN PBL clouds in photolysis option 4 (TUV) #1622

Show file tree
Hide file tree
Changes from 6 commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
4b89261
Adds capability to include sub-grid boundary layer clouds from the MY…
jordanschnell Jan 5, 2022
4c4a82b
Add check for proper namelist values (icloud = 1 and MYNNPBL )
jordanschnell Jan 5, 2022
22bbdf5
check for TUV
jordanschnell Jan 5, 2022
8c6936e
Add ifdefs around a chem nml option in check_a_mundo
davegill Jan 5, 2022
52f2cae
Stop being helpful, don't check MYNN and icloud_bl as fatal
davegill Jan 6, 2022
6474aab
TUV 18:13
davegill Jan 6, 2022
661d7ac
fix, add bl_cld/QC_BL to other QC variables.
jordanschnell Jan 6, 2022
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
1 change: 1 addition & 0 deletions Registry/registry.chem
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3816,6 +3816,7 @@ rconfig real lambda_cutoff namelist,chem max_domains
# for cloud handling in TUV photolysis scheme
rconfig integer cld_od_opt namelist,chem max_domains 1 rh "cld_od_opt" "" ""
rconfig integer pht_cldfrc_opt namelist,chem max_domains 1 rh "pht_cldfrc_opt" "" ""
rconfig logical phot_blcld namelist,chem max_domains .false. rh "phot_blcld" "" ""
rconfig integer vertmix_onoff namelist,chem max_domains 1 rh "vertmix_onoff" "" ""
rconfig integer chem_in_opt namelist,chem max_domains 0 rh "chem_in_opt" "" ""
rconfig integer phot_opt namelist,chem max_domains 0 rh "phot_opt" "" ""
Expand Down
1 change: 1 addition & 0 deletions chem/chem_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -960,6 +960,7 @@ end SUBROUTINE sum_pm_driver
chem,rho,dz8w,grid%xlat,grid%xlong, &
zmid,z_at_w, &
grid%qc_cu,grid%qi_cu, &
grid%qc_bl,grid%qi_bl, &
grid%ph_macr,grid%ph_o31d,grid%ph_o33p,grid%ph_no2, &
grid%ph_clno2, &
grid%ph_no3o2, &
Expand Down
26 changes: 22 additions & 4 deletions chem/module_phot_tuv.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,6 +75,7 @@ subroutine tuv_driver( &
gmt, julday, t_phy, moist, aerwrf, &
p8w, t8w, p_phy, chem, rho_phy, &
dz8w, xlat, xlong, z, z_at_w, gd_cloud, gd_cloud2, &
bl_cld,bl_cld2, &
tauaer1,tauaer2,tauaer3,tauaer4, &
gaer1,gaer2,gaer3,gaer4, &
waer1,waer2,waer3,waer4, &
Expand Down Expand Up @@ -117,7 +118,7 @@ subroutine tuv_driver( &
pm2_5_dry,pm2_5_water
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
INTENT(INOUT ) :: &
gd_cloud,gd_cloud2
gd_cloud,gd_cloud2,bl_cld,bl_cld2
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
INTENT(IN ) :: tauaer1, tauaer2, tauaer3, tauaer4, &
waer1, waer2, waer3, waer4, &
Expand Down Expand Up @@ -436,13 +437,23 @@ subroutine tuv_driver( &
moist(i,kts:kte,j,p_qc), moist(i,kts:kte,j,p_qi), &
moist(i,kts:kte,j,p_qs), kts, kte )
elseif( config_flags%pht_cldfrc_opt == 2 ) then
if( config_flags%cu_diag == 1 ) then
call cldfrac_fractional( cldfrac(kts:kte), &
if( config_flags%cu_diag == 1 ) then
if ( config_flags%phot_blcld .and. config_flags%icloud_bl == 1 ) then
call cldfrac_fractional( cldfrac(kts:kte), &
moist(i,kts:kte,j,p_qv), &
moist(i,kts:kte,j,p_qc) + gd_cloud(i,kts:kte,j), &
moist(i,kts:kte,j,p_qi) + gd_cloud2(i,kts:kte,j) + &
bl_cld(i,kts:kte,j) + bl_cld2(i,kts:kte,j), &
Comment thread
davegill marked this conversation as resolved.
Outdated
moist(i,kts:kte,j,p_qs), &
p_phy(i,kts:kte,j), t_phy(i,kts:kte,j), kts, kte )
else
call cldfrac_fractional( cldfrac(kts:kte), &
moist(i,kts:kte,j,p_qv), &
moist(i,kts:kte,j,p_qc) + gd_cloud(i,kts:kte,j), &
moist(i,kts:kte,j,p_qi) + gd_cloud2(i,kts:kte,j), &
moist(i,kts:kte,j,p_qs), &
p_phy(i,kts:kte,j), t_phy(i,kts:kte,j), kts, kte )
endif
else
call cldfrac_fractional( cldfrac(kts:kte), &
moist(i,kts:kte,j,p_qv), moist(i,kts:kte,j,p_qc), &
Expand Down Expand Up @@ -513,7 +524,14 @@ subroutine tuv_driver( &
qll(:) = 0.
qll(kts:kte) = moist(i,kts:kte,j,p_qc) + moist(i,kts:kte,j,p_qi)
if( config_flags%cu_diag == 1 ) then
qll(kts:kte) = qll(kts:kte) + gd_cloud(i,kts:kte,j) + gd_cloud2(i,kts:kte,j)
if ( config_flags%phot_blcld .and. config_flags%icloud_bl == 1 ) then
qll(kts:kte) = qll(kts:kte) + gd_cloud(i,kts:kte,j) + &
gd_cloud2(i,kts:kte,j) + &
bl_cld(i,kts:kte,j) + &
bl_cld2(i,kts:kte,j)
else
qll(kts:kte) = qll(kts:kte) + gd_cloud(i,kts:kte,j) + gd_cloud2(i,kts:kte,j)
endif
endif
qll(kts:kte) = 1.e3*rhoa(kts:kte)*qll(kts:kte)
where( qll(kts:kte) < 1.e-5 )
Expand Down
4 changes: 3 additions & 1 deletion chem/photolysis_driver.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,7 @@ SUBROUTINE photolysis_driver (id,curr_secs,ktau,dtstep, &
dt_cld,af_dir,af_dn,af_up,par,erythema, &
gmt,julday,t_phy,moist,aerwrf,p8w,t8w,p_phy, &
chem,rho_phy,dz8w,xlat,xlong,z,z_at_w,gd_cloud,gd_cloud2, &
bl_cld,bl_cld2, &
ph_macr,ph_o31d,ph_o33p,ph_no2, &
ph_clno2, &
ph_no3o2, &
Expand Down Expand Up @@ -99,7 +100,7 @@ SUBROUTINE photolysis_driver (id,curr_secs,ktau,dtstep, &
REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
OPTIONAL, &
INTENT(INOUT ) :: &
gd_cloud,gd_cloud2
gd_cloud,gd_cloud2,bl_cld,bl_cld2
!
! arrays that hold the aerosol optical properties
!
Expand Down Expand Up @@ -225,6 +226,7 @@ SUBROUTINE photolysis_driver (id,curr_secs,ktau,dtstep, &
gmt, julday, t_phy, moist, aerwrf, &
p8w, t8w, p_phy, chem, rho_phy, &
dz8w, xlat, xlong, z, z_at_w, gd_cloud, gd_cloud2, &
bl_cld,bl_cld2, &
tauaer1,tauaer2,tauaer3,tauaer4, &
gaer1,gaer2,gaer3,gaer4, &
waer1,waer2,waer3,waer4, &
Expand Down
25 changes: 25 additions & 0 deletions share/module_check_a_mundo.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1591,6 +1591,31 @@ END FUNCTION bep_bem_ngr_u
CALL wrf_debug ( 1, wrf_err_message )
END IF

#if (WRF_CHEM == 1)
!-----------------------------------------------------------------------
! Make sure phot_blcld is only used when icloud_bl==1 and MYNN is chosen.
!-----------------------------------------------------------------------

oops = 0
DO i = 1, model_config_rec % max_dom
IF ( .NOT. model_config_rec % grid_allowed(i) ) CYCLE
IF ( model_config_rec%phot_blcld(i) ) THEN
IF ( ( model_config_rec%icloud_bl .eq. 1 ) .AND. &
( ( model_config_rec%bl_pbl_physics(i) .EQ. MYNNPBLSCHEME2 ) .OR. &
( model_config_rec%bl_pbl_physics(i) .EQ. MYNNPBLSCHEME3 ) ) ) THEN
!CORRECTLY CONFIGURED
ELSE
oops = oops + 1
END IF
END IF
ENDDO ! Loop over domains
IF ( oops .GT. 0 ) THEN
wrf_err_message = '--- ERROR: Need MYNN PBL and icloud_bl = 1 for phot_blcld = .true.'
CALL wrf_message(wrf_err_message)
count_fatal_error = count_fatal_error + 1
END IF
#endif

!-----------------------------------------------------------------------
! We need to know if any of the cumulus schemes are active. This
! allows the model to allocate space.
Expand Down