Skip to content

Hygroscopicity MS11 ratio option #31

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
Pierre Siddall (Pierre-siddall) merged 15 commits into from
May 20, 2026
Merged

Hygroscopicity MS11 ratio option #31

Show file tree
Hide file tree
Changes from 11 commits
Commits
Show all changes
15 commits
Select commit Hold shift + click to select a range
daf5462
Added plumbing to allow ratio to be transferred from namelists in UM …
JoeCartonKelly-MO Mar 31, 2026
335fb3e
more plumbing plus a check to make sure ratio adds up to one and defa…
JoeCartonKelly-MO Mar 31, 2026
37d746f
Introduced ratio for hygroscopicity routine, allowing growth to be di…
JoeCartonKelly-MO Mar 31, 2026
f1e20da
Merge branch 'MetOffice:main' into hygro_ratios_option
JoeCartonKelly-MO Mar 31, 2026
3a351f6
Updated styling to match UMDP better and improve readability
JoeCartonKelly-MO Apr 7, 2026
3af1d9f
umdp3_fixer.py run and files tidied
JoeCartonKelly-MO Apr 8, 2026
2fdb5c4
Added myself to contributors list
JoeCartonKelly-MO Apr 15, 2026
f99165d
Corrected bug to use md not mdcopy and removed bc from routine
JoeCartonKelly-MO Apr 22, 2026
84ca6a5
updated description of ratio to make clear it only impacts wet part o…
JoeCartonKelly-MO Apr 28, 2026
f6495b1
updated to make mdcopy use consistent - won't impact output as all va…
JoeCartonKelly-MO Apr 28, 2026
e457333
Fixed bug where water content routine wouldnt run for seasalt/organic…
JoeCartonKelly-MO May 12, 2026
4a08b69
Simplified things following Dan's advice
JoeCartonKelly-MO May 12, 2026
771b35b
fixed small bug
JoeCartonKelly-MO May 12, 2026
1ec35b5
umdp3 checker changes
JoeCartonKelly-MO May 12, 2026
c804c50
Merge branch 'main' into hygro_ratios_option
Pierre-siddall May 20, 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
1 change: 1 addition & 0 deletions CONTRIBUTORS.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,3 +6,4 @@
| andrewcoughtrie | Andrew Coughtrie | Met Office | 2026-02-10 |
| jennyhickson | Jenny Hickson | Met Office | 2026-03-02 |
| Pierre-siddall | Pierre Siddall | Met Office | 2026-03-13 |
| JoeCartonKelly-MO | Joseph Carton-Kelly | Met Office | 2026-04-15 |
9 changes: 9 additions & 0 deletions src/control/core/interface/ukca_config_specification_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -443,6 +443,10 @@ MODULE ukca_config_specification_mod
REAL :: acc_cor_scav_scaling ! Scaling factor for scavenging
! parameters for the accumulation and
! coarse modes
REAL :: solinsol_hygro_ratio(4) ! SOL/INSOL hygroscopicity ratios
! cp_su, cp_cl, cp_bc, cp_oc
! This ratio only affects the wet
! part of the aerosol

! -- GLOMAP deposition configuration options --
LOGICAL :: l_ddepaer ! True for aerosol dry deposition
Expand Down Expand Up @@ -996,6 +1000,7 @@ SUBROUTINE init_ukca_configuration()
glomap_config%l_dust_mp_ageing = .FALSE.
glomap_config%dry_depvel_acc_scaling = rmdi
glomap_config%acc_cor_scav_scaling = rmdi
glomap_config%solinsol_hygro_ratio(:) = rmdi

! -- GLOMAP deposition configuration options --
glomap_config%l_ddepaer = .FALSE.
Expand Down Expand Up @@ -1122,6 +1127,7 @@ SUBROUTINE ukca_get_config( &
hno3_uptake_coeff, &
sigwmin, &
sigma_updraught_scaling, &
solinsol_hygro_ratio, &
l_cal360, &
l_ukca_chem_aero, &
l_ukca_mode, &
Expand Down Expand Up @@ -1324,6 +1330,7 @@ SUBROUTINE ukca_get_config( &
REAL, OPTIONAL, INTENT(OUT) :: ph_fit_coeff_b
REAL, OPTIONAL, INTENT(OUT) :: ph_fit_intercept
REAL, OPTIONAL, INTENT(OUT) :: hno3_uptake_coeff
REAL, OPTIONAL, INTENT(OUT) :: solinsol_hygro_ratio(4)

LOGICAL, OPTIONAL, INTENT(OUT) :: l_cal360
LOGICAL, OPTIONAL, INTENT(OUT) :: l_ukca_chem_aero
Expand Down Expand Up @@ -1748,6 +1755,8 @@ SUBROUTINE ukca_get_config( &
IF (PRESENT(mode_incld_so2_rfrac)) &
mode_incld_so2_rfrac = glomap_config%mode_incld_so2_rfrac
IF (PRESENT(l_aero_rainout)) l_aero_rainout = glomap_config%l_aero_rainout
IF (PRESENT(solinsol_hygro_ratio)) solinsol_hygro_ratio(:) = &
glomap_config%solinsol_hygro_ratio(:)
IF (PRESENT(l_cv_rainout)) l_cv_rainout = glomap_config%l_cv_rainout
IF (PRESENT(i_mode_nucscav)) i_mode_nucscav = glomap_config%i_mode_nucscav
IF (PRESENT(l_impc_scav)) l_impc_scav = glomap_config%l_impc_scav
Expand Down
23 changes: 23 additions & 0 deletions src/control/core/top_level/ukca_setup_mod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -129,6 +129,7 @@ SUBROUTINE ukca_setup(error_code, &
ph_fit_intercept, &
sigwmin, &
sigma_updraught_scaling, &
solinsol_hygro_ratio, &
const_rmol, &
const_tfs, &
const_rho_water, &
Expand Down Expand Up @@ -461,6 +462,7 @@ SUBROUTINE ukca_setup(error_code, &
REAL, OPTIONAL, INTENT(IN) :: sigwmin
REAL, OPTIONAL, INTENT(IN) :: hno3_uptake_coeff
REAL, OPTIONAL, INTENT(IN) :: sigma_updraught_scaling
REAL, OPTIONAL, INTENT(IN) :: solinsol_hygro_ratio(4)
REAL, OPTIONAL, INTENT(IN) :: const_rmol
REAL, OPTIONAL, INTENT(IN) :: const_tfs
REAL, OPTIONAL, INTENT(IN) :: const_rho_water
Expand Down Expand Up @@ -600,6 +602,8 @@ SUBROUTINE ukca_setup(error_code, &
LOGICAL :: l_strat_scheme_selected ! True if a Stratospheric scheme is
! selected for chemistry

REAL :: sum_solinsol_hygro_ratio ! Sum of solinsol hygroscopicity ratios (=1)

INTEGER (KIND=jpim), PARAMETER :: zhook_in = 0 ! DrHook tracing entry
INTEGER (KIND=jpim), PARAMETER :: zhook_out = 1 ! DrHook tracing exit
REAL (KIND=jprb) :: zhook_handle ! DrHook tracing
Expand Down Expand Up @@ -1313,6 +1317,9 @@ SUBROUTINE ukca_setup(error_code, &

END IF

IF (PRESENT(solinsol_hygro_ratio)) &
glomap_config%solinsol_hygro_ratio(:) = solinsol_hygro_ratio(:)

! -- GLOMAP deposition configuration options --

! Dry deposition
Expand Down Expand Up @@ -1721,6 +1728,22 @@ SUBROUTINE ukca_setup(error_code, &

END IF

! Ensure the sum of solinsol_hygro_ratio is 1
IF (ukca_config%l_ukca_mode .AND. &
glomap_config%i_mode_setup == 11) THEN
sum_solinsol_hygro_ratio = 0.0
DO i = 1,4
sum_solinsol_hygro_ratio = sum_solinsol_hygro_ratio + &
glomap_config%solinsol_hygro_ratio(i)
END DO
IF (PRESENT(solinsol_hygro_ratio)) THEN
glomap_config%solinsol_hygro_ratio(:) = &
glomap_config%solinsol_hygro_ratio(:) / sum_solinsol_hygro_ratio
ELSE
glomap_config%solinsol_hygro_ratio(:) = [1.0, 0.0, 0.0, 0.0]
END IF
END IF

! Initialise chemical definition arrays
CALL ukca_chem1_init()

Expand Down
66 changes: 47 additions & 19 deletions src/science/core/aerosols/glomap/ukca_volume_mode.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -156,8 +156,8 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
USE ukca_config_constants_mod, ONLY: rho_water, avogadro, rho_so4

USE ukca_mode_setup, ONLY: nmodes, nanion, ncation, &
cp_su, cp_oc, cp_cl, &
cp_so, cp_nh4, cp_no3, cp_nn
cp_su, cp_oc, cp_cl, cp_so, &
cp_nh4, cp_no3, cp_nn

USE yomhook, ONLY: lhook, dr_hook
USE parkind1, ONLY: jprb, jpim
Expand Down Expand Up @@ -218,6 +218,7 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
INTEGER :: iimode
INTEGER :: iicp
INTEGER :: ierr(nbox)
INTEGER :: dim_1_comp, dim_2_comp ! dimensions to allocate for component_copy
LOGICAL (KIND=log_small) :: mask(nbox)
LOGICAL (KIND=log_small) :: mask_sol(nbox)
LOGICAL (KIND=log_small) :: mask_nosol(nbox)
Expand All @@ -240,6 +241,8 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
REAL :: mmwrhow
REAL :: f_ao
REAL :: cl(nbox,-nanion:ncation) !ION CONCS (MOL/CC OF AIR)
REAL :: mdcopy(nbox,nmodes,glomap_variables_local%ncp) !Redistributes ions for solinsol
LOGICAL, ALLOCATABLE :: component_copy(:,:)

REAL :: tmp1(nbox,nmodes)
!
Expand Down Expand Up @@ -280,6 +283,14 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
soluble => glomap_variables_local%soluble
x => glomap_variables_local%x

! Create component_copy which can be used for SOL/INSOL wet aerosol
! growth. It copies the values of component but can be updated later
! in the routine
dim_1_comp = SIZE(component,1)
dim_2_comp = SIZE(component,2)
ALLOCATE(component_copy(dim_1_comp,dim_2_comp))
component_copy(:,:) = component(:,:)

!at this point in the code, the value of RP does not matter
rp(:)=100.0e-9 ! dummy value
CALL ukca_vapour(nbox,t,pmid,s,rp,wts,rhosol_strat)
Expand All @@ -297,6 +308,7 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &

denom(:)=0.0 ! define on all points to avoid error
denom2(:)=0.0
mdcopy(:,:,:) = md(:,:,:) ! assigning md input values to mdcopy

! Correct relative humidities to lie within the range of 10-90%
corrh(:)=rh(:)
Expand Down Expand Up @@ -348,21 +360,35 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
cl(:,i)=0.0 ! set all concentrations to zero initially
END DO

IF (component(imode,cp_su)) THEN ! assume all H2SO4 --> SO4
cl(:,-2)=md(:,imode,cp_su)/avogadro ! [SO4] in moles/cc (air)
! If SOLINSOL then redistribute sulphate md into other species
! hygro_ratio(3) is black carbon which is important to include
! in the ratio but isn't used in this routine
IF (glomap_config%i_mode_setup == 11) THEN
mdcopy(:,imode,cp_su) = glomap_config%solinsol_hygro_ratio(1)* &
(md(:,imode,cp_su))
mdcopy(:,imode,cp_cl) = glomap_config%solinsol_hygro_ratio(2)* &
(md(:,imode,cp_su))
mdcopy(:,imode,cp_oc) = glomap_config%solinsol_hygro_ratio(4)* &
(md(:,imode,cp_su))
component_copy(imode, cp_cl) = .TRUE.
component_copy(imode, cp_oc) = .TRUE.
END IF

IF (component_copy(imode,cp_su)) THEN ! assume all H2SO4 --> SO4
cl(:,-2)=mdcopy(:,imode,cp_su)/avogadro ! [SO4] in moles/cc (air)
END IF

IF (component(imode,cp_so)) THEN
cl(:,-2)=cl(:,-2)+(fhyg_aom/avogadro)*(md(:,imode,cp_so)/f_ao)
IF (component_copy(imode,cp_so)) THEN
Comment thread
JoeCartonKelly-MO marked this conversation as resolved.
Outdated
cl(:,-2)=cl(:,-2)+(fhyg_aom/avogadro)*(mdcopy(:,imode,cp_so)/f_ao)
! .. Increment concentration of SO4 ions to represent the
! .. presence of hygroscopic aged organic aerosol mass in CP_SO.
! .. Assume it has uptake behaviour at fraction FHYG_AOM of SO4.
! .. Need to divide by F_AO because MD of CP_SO is in "moles" of POM
! .. whereas CL needs to be in moles of aged organic (MM=0.15kg/mol).
END IF

IF (component(imode,cp_oc)) THEN
cl(:,-2)=cl(:,-2)+(fhyg_aom/avogadro)*(md(:,imode,cp_oc)/f_ao)
IF (component_copy(imode,cp_oc)) THEN
Comment thread
JoeCartonKelly-MO marked this conversation as resolved.
Outdated
cl(:,-2)=cl(:,-2)+(fhyg_aom/avogadro)*(mdcopy(:,imode,cp_oc)/f_ao)
! .. Increment concentration of SO4 ions to represent the
! .. presence of hygroscopic aged organic aerosol mass in CP_OC.
! .. Assume it has uptake behaviour at fraction FHYG_AOM of SO4.
Expand All @@ -376,26 +402,26 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
!
END IF

IF (component(imode,cp_cl)) THEN ! assume complete dissociation
cl(:,3)=md(:,imode,cp_cl)/avogadro ! [Na] in moles per cc (air)
cl(:,-4)=md(:,imode,cp_cl)/avogadro ! [Cl] in moles per cc (air)
IF (component_copy(imode,cp_cl)) THEN ! assume complete dissociation
cl(:,3)=mdcopy(:,imode,cp_cl)/avogadro ! [Na] in moles per cc (air)
cl(:,-4)=mdcopy(:,imode,cp_cl)/avogadro ! [Cl] in moles per cc (air)
Comment thread
JoeCartonKelly-MO marked this conversation as resolved.
END IF

! Nitrate scheme
IF (UBOUND(component,DIM=2) >= cp_no3) THEN
IF (UBOUND(component_copy,DIM=2) >= cp_no3) THEN

IF (component(imode,cp_nh4)) THEN ! assume complete dissociation
cl(:,2)=md(:,imode,cp_nh4)/avogadro ! [NH4] in moles per cc (air)
IF (component_copy(imode,cp_nh4)) THEN ! assume complete dissociation
cl(:,2)=mdcopy(:,imode,cp_nh4)/avogadro ! [NH4] in moles per cc (air)
END IF

IF (component(imode,cp_no3)) THEN ! assume complete dissociation
cl(:,-3)=md(:,imode,cp_no3)/avogadro ! [NO3] in moles per cc (air)
IF (component_copy(imode,cp_no3)) THEN ! assume complete dissociation
cl(:,-3)=mdcopy(:,imode,cp_no3)/avogadro ! [NO3] in moles per cc (air)
END IF

IF (component(imode,cp_nn)) THEN ! assume complete dissociation
cl(:,3)=cl(:,3)+md(:,imode,cp_nn)/avogadro
IF (component_copy(imode,cp_nn)) THEN ! assume complete dissociation
cl(:,3)=cl(:,3)+mdcopy(:,imode,cp_nn)/avogadro
! [Na] in moles per cc (air)
cl(:,-3)=cl(:,-3)+md(:,imode,cp_nn)/avogadro
cl(:,-3)=cl(:,-3)+mdcopy(:,imode,cp_nn)/avogadro
! [NO3] in moles per cc (air)
END IF

Expand Down Expand Up @@ -880,6 +906,8 @@ SUBROUTINE ukca_volume_mode( glomap_variables_local, nbox, nd, md, mdt, &
END DO
END IF

DEALLOCATE(component_copy)

IF (lhook) CALL dr_hook(ModuleName//':'//RoutineName,zhook_out,zhook_handle)
RETURN
END SUBROUTINE ukca_volume_mode
Expand Down