🔧

parcel/parcel.jl

@@ -67,7 +67,7 @@ function parcel_model(dY, Y, p, t)
 
     # Adiabatic parcel coefficients
     a1 = L_vap * grav / cp_a / T^2 / R_v - grav / R_a / T
-    a2 = R_v * p_a / e_sl / R_a
+    a2 = 1 / q_vap
     a3 = L_vap^2 / R_v / T^2 / cp_a
     a4 = L_vap * L_subl / R_v / T^2 / cp_a
     a5 = L_vap * L_fus / R_v / cp_a / (T^2)
@@ -204,9 +204,9 @@ function parcel_model(dY, Y, p, t)
     end
 
     # Update the tendecies
-    dq_vap_dt = -dq_ice_dt - dq_liq_dt
     dq_ice_dt = dq_ice_dt_vap_to_ice + dq_ice_dt_liq_to_ice
     dq_liq_dt = dq_liq_dt_vap_to_liq - dq_ice_dt_liq_to_ice
+    dq_vap_dt = -dq_ice_dt - dq_liq_dt
     ds_l_dt = a1 * w * s_liq - (a2 + a3) * s_liq * dq_liq_dt_vap_to_liq - (a2 + a4) * s_liq * dq_ice_dt_vap_to_ice - a5 * s_liq * dq_ice_dt_liq_to_ice
     dp_a_dt = -p_a * grav / R_a / T * w
     dT_dt = -grav / cp_a * w + L_vap / cp_a * dq_liq_dt_vap_to_liq + L_fus / cp_a * dq_ice_dt_liq_to_ice + L_subl / cp_a * dq_ice_dt_vap_to_ice