errors in parcel

parcel/parcel.jl

@@ -86,7 +86,7 @@ function parcel_model(dY, Y, p, t)
         dqi_dt_new_depo = dN_act_dt_depo * 4 / 3 * FT(π) * r_nuc^3 * ρ_ice / ρ_air
     end
     if "MohlerRate_Deposition" in ice_nucleation_modes
-        dN_act_dt_depo = MohlerDepositionRate(aerosol, ip.deposition, S_i, T, (dY[1] * ξ), N_aer)
+        dN_act_dt_depo = CMI_het.MohlerDepositionRate(aerosol, ip.deposition, S_i, T, (dY[1] * ξ), N_aer)
         dqi_dt_new_depo = dN_act_dt_depo * 4 / 3 * FT(π) * r_nuc^3 * ρ_ice / ρ_air
     end
 
@@ -268,6 +268,7 @@ function run_parcel(IC, t_0, t_end, p)
     #! format: off
     FT = eltype(IC)
     (; const_dt, ice_nucleation_modes, growth_modes, droplet_size_distribution) = p
+    timestepper = ODE.Tsit5()
 
     println(" ")
     println("Running parcel model with: ")
@@ -278,15 +279,12 @@ function run_parcel(IC, t_0, t_end, p)
     end
     if "MohlerRate_Deposition" in ice_nucleation_modes
         print("Deposition nucleation based on rate eqn in Mohler 2006 ")
-        timestepper = ODE.Tsit5()
     end
     if "ImmersionFreezing" in ice_nucleation_modes
         print("Immersion freezing ")
-        timestepper = ODE.Tsit5()
     end
     if "HomogeneousFreezing" in ice_nucleation_modes
         print("Homogeneous freezing ")
-        timestepper = ODE.Tsit5()
     end
     print("\n")
     print("Growth modes: ")