diff --git a/CONTRIBUTORS.md b/CONTRIBUTORS.md index c1d9657..f3f5e3e 100644 --- a/CONTRIBUTORS.md +++ b/CONTRIBUTORS.md @@ -5,3 +5,4 @@ | james-bruten-mo | James Bruten | Met Office | 2025-12-09 | | t00sa | Sam Clarke-Green | Met Office | 2026-03-02 | | Pierre-siddall| Pierre Siddall | Met Office | 2026-03-11 | +| iboutle | Ian Boutle | Met Office | 2026-08-12 | diff --git a/src/activation.F90 b/src/activation.F90 index 3c2ec85..a0f8d85 100644 --- a/src/activation.F90 +++ b/src/activation.F90 @@ -261,14 +261,11 @@ subroutine activate(dt, cloud_mass, cloud_number, w, rho, dnumber, dmac, T, p, activated_arg = sum(dnccn_all(:)) + dactive = zero_real_wp + dmass_d = zero_real_wp if (Smax_act > smax_act_min .and. .not. l_warm) then - if (iopt_inuc < 4) then - ! For lower-order ice nucleation options, need to initialise - ! dactive and dmass_d - dactive = zero_real_wp - dmass_d = zero_real_wp - else + if (iopt_inuc > 3) then ! For higher-order ice nucleation schemes, dactive and dmass_d ! can be based on dustphys dactive = 0.01*dustphys%N(1)