Skip to content

noresm3_0_025_cam6_4_121: Re-introduce the Bennartz mg diagnostics #250

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
gold2718 merged 3 commits into from
Jan 28, 2026
Merged

noresm3_0_025_cam6_4_121: Re-introduce the Bennartz mg diagnostics #250

Changes from 1 commit
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
d3b7326
re-introduce the Bennartz mg diagnostics CCN_B ACTNL_B and FCTL_B
Ovewh Dec 19, 2025
0e47f71
only do Bennartz calculation if those history fields are active
Ovewh Jan 5, 2026
b69bba6
Refactor the calculation of ACTNL_B FCTL_B and CCN_B to conditionally…
Ovewh Jan 28, 2026
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
47 changes: 47 additions & 0 deletions src/physics/cam7/micro_pumas_cam.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1152,6 +1152,11 @@ subroutine micro_pumas_cam_init(pbuf2d)
call addfld ('FCTL', horiz_only, 'A', 'fraction', 'Fractional occurrence of cloud top liquid', sampled_on_subcycle=.true.)
call addfld ('FCTI', horiz_only, 'A', 'fraction', 'Fractional occurrence of cloud top ice', sampled_on_subcycle=.true.)

! Bennartz cloud top diagnostics for comparing to Bernartz CDNC concentrations
call addfld ('ACTNL_B', horiz_only, 'A', 'm-3', 'Average Cloud Top droplet number (Bennartz)', sampled_on_subcycle=.true.)
call addfld ('FCTL_B', horiz_only, 'A', 'fraction','Fractional occurrence of cloud top liquid (Bennartz)', sampled_on_subcycle=.true.)
call addfld ('CCN_B', horiz_only, 'A', 'm-3', 'Average Cloud Top liquid CCN (Bennartz)', sampled_on_subcycle=.true.)

! New frequency arrays for mixed phase and supercooled liquid (only and mixed) for (a) Cloud Top and (b) everywhere..
call addfld ('FREQM', (/ 'lev' /), 'A', 'fraction', 'Fractional occurrence of mixed phase', sampled_on_subcycle=.true.)
call addfld ('FREQSL', (/ 'lev' /), 'A', 'fraction', 'Fractional occurrence of only supercooled liquid', sampled_on_subcycle=.true.)
Expand Down Expand Up @@ -1287,6 +1292,17 @@ subroutine micro_pumas_cam_init(pbuf2d)
call add_default ('FREQS ', 1, ' ')
call add_default ('FREQL ', 1, ' ')
call add_default ('FREQI ', 1, ' ')
call add_default ('ACTNL ', 1, ' ')
call add_default ('ACTREL ', 1, ' ')
call add_default ('FCTL ', 1, ' ')
call add_default ('ACTNI ', 1, ' ')
call add_default ('ACTREI ', 1, ' ')
call add_default ('FCTI ', 1, ' ')
!These are for comparing to Bennartz
call add_default ('FCTL_B ', 1, ' ')
call add_default ('ACTNL_B ', 1, ' ')
call add_default ('CCN_B ', 1, ' ')

do m = 1, ncnst
call cnst_get_ind(cnst_names(m), mm)
call add_default(cnst_name(mm), 1, ' ')
Expand Down Expand Up @@ -1813,6 +1829,10 @@ subroutine micro_pumas_cam_tend(state, ptend, dtime, pbuf)

real(r8) :: ftem_grid(pcols,pver)

real(r8) :: fctl_b(pcols) !frequency of occurrence for Bennartz
real(r8) :: ctnl_b(pcols) !cdnc [/m3] for Bennartz
real(r8) :: ccn_b(pcols) !ccm [/m3] defined as for cdnc for Bennartz

! Variables for precip efficiency calculation
real(r8) :: minlwp ! LWP threshold

Expand Down Expand Up @@ -3492,6 +3512,10 @@ subroutine micro_pumas_cam_tend(state, ptend, dtime, pbuf)
fctsl_grid = 0._r8
fctslm_grid= 0._r8

fctl_b = 0._r8
ctnl_b = 0._r8
ccn_b = 0._r8

do i = 1, ngrdcol
do k = top_lev, pver
if ( liqcldf_grid(i,k) > 0.01_r8 .and. icwmrst_grid(i,k) > 1.e-7_r8 ) then
Expand Down Expand Up @@ -3538,6 +3562,25 @@ subroutine micro_pumas_cam_tend(state, ptend, dtime, pbuf)
if (qcsevap_idx > 0 ) qcsevapout_grid_ptr = qcsevapout_grid
if (qisevap_idx > 0 ) qisevapout_grid_ptr = qisevapout_grid

!Calculate values for comparing with Bennartz 2017
do i = 1, ncol
do k = top_lev, pver
!Criterions for Bennartz (2017) to use values from a column
!1) 268 < T < 300 [K]
!2) liquid cloud fraction > 10 %
if ( liqcldf_grid(i,k) > 0.1_r8 &
.and. state_loc%t(i,k) > 268.0_r8 &
.and. state_loc%t(i,k) < 300.0_r8 ) then
!Save cloud fraction and in-cloud number conc
ctnl_b(i) = icwnc_grid(i,k) * liqcldf_grid(i,k)
fctl_b(i) = liqcldf_grid(i,k)
ccn_b(i) = ncal(i,k) * liqcldf_grid(i,k)
exit !==> Go out to i=1,ncol-loop
end if
!--IH
end do
end do

! --------------------------------------------- !
! General outfield calls for microphysics !
! --------------------------------------------- !
Expand Down Expand Up @@ -3828,6 +3871,10 @@ subroutine micro_pumas_cam_tend(state, ptend, dtime, pbuf)
call outfld('FCTSL', fctsl_grid, pcols, lchnk)
call outfld('FCTSLM', fctslm_grid, pcols, lchnk)

call outfld( 'FCTL_B' , fctl_b, pcols, lchnk )
call outfld( 'ACTNL_B' , ctnl_b, pcols, lchnk )
call outfld( 'CCN_B' , ccn_b, pcols, lchnk )

if (micro_mg_version > 2) then
call outfld('PRACGO', pracgo_grid, pcols, lchnk)
call outfld('PSACRO', psacro_grid, pcols, lchnk)
Expand Down