Skip to content

cam6_4_141: Aerosol hygroscopicity update #1462

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
cacraigucar merged 4 commits into from
Dec 30, 2025
Merged

cam6_4_141: Aerosol hygroscopicity update #1462

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
a20a8f5
use hygro from aero_props (phys props files)
fvitt Dec 16, 2025
9892570
remove sol_factb limiters
fvitt Dec 23, 2025
a229835
remove the sol_factb_interstitial namelist setting
fvitt Dec 23, 2025
f1d029c
Merge tag 'cam6_4_140' into aerosol_hygroscopicity
cacraigucar Dec 29, 2025
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: 0 additions & 1 deletion bld/build-namelist
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3660,7 +3660,6 @@ if ($clubb_sgs =~ /$TRUE/io) {
# Tuning for wet scavenging of modal aerosols
if ($chem =~ /_mam/) {
add_default($nl, 'sol_facti_cloud_borne');
add_default($nl, 'sol_factb_interstitial');
add_default($nl, 'sol_factic_interstitial');
}
if ($carma =~ /trop_strat/) {
Expand Down
28 changes: 3 additions & 25 deletions src/chemistry/aerosol/aerosol_state_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -927,7 +927,7 @@ function sol_factb_interstitial(self, bin_ndx, ncol, nlev, aero_props) result(so
real(r8) :: totmmr(ncol,nlev)
real(r8) :: solmmr(ncol,nlev)
integer :: ispc
character(len=aero_name_len) :: spectype
real(r8) :: spechygro

sol_factb(:,:) = 0.0_r8

Expand All @@ -936,40 +936,18 @@ function sol_factb_interstitial(self, bin_ndx, ncol, nlev, aero_props) result(so

do ispc = 1, aero_props%nspecies(bin_ndx)

call aero_props%species_type(bin_ndx, ispc, spectype)
call aero_props%get(bin_ndx, ispc, hygro=spechygro)
call self%get_ambient_mmr(ispc, bin_ndx, aer_mmr)

totmmr(:ncol,:) = totmmr(:ncol,:) + aer_mmr(:ncol,:)

if (trim(spectype) == 'sulfate') then
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*0.5_r8
end if
if (trim(spectype) == 'p-organic') then
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*0.2_r8
end if
if (trim(spectype) == 's-organic') then
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*0.2_r8
end if
if (trim(spectype) == 'dust') then
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*0.1_r8
end if
if (trim(spectype) == 'seasalt') then
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*0.8_r8
end if
solmmr(:ncol,:) = solmmr(:ncol,:) + aer_mmr(:ncol,:)*spechygro

end do !nspec

where ( totmmr > 0._r8 )
sol_factb = solmmr/totmmr
end where

where ( sol_factb > 0.8_r8 )
sol_factb = 0.8_r8
end where
where ( sol_factb < 0.1_r8 )
sol_factb = 0.1_r8
end where

end function sol_factb_interstitial


Expand Down