simplified further

CliMA · Jan 29, 2024 · 9a9a019 · 9a9a019
1 parent 26cf486
commit 9a9a019
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/docs/src/IceNucleationParcel0D.md b/docs/src/IceNucleationParcel0D.md
@@ -79,7 +79,7 @@ Accounting for conservation of water, i.e. ``\frac{dq_v}{dt} + \frac{dq_l}{dt} +
 and rearranging the terms
 ```math
 \begin{equation}
-\frac{dS_l}{dt} = a_1 w S_l - \left(a_2 + a_3 S_l \right) \frac{dq_{l,vap}}{dt} - \left(a_2 + a_4  S_l\right) \frac{dq_{i,subl}}{dt} - a_5 S_l \frac{dq_{i,fus}}{dt}
+\frac{dS_l}{dt} = a_1 w S_l - \left(a_2 + a_3 \right)S_l  \frac{dq_{l,vap}}{dt} - \left(a_2 + a_4 \right) S_l \frac{dq_{i,subl}}{dt} - a_5 S_l \frac{dq_{i,fus}}{dt}
 \end{equation}
 ```
 where:
@@ -90,7 +90,7 @@ a_1 = \frac{L_v g}{c_p T^2 R_v} - \frac{g}{R_a T}
 ```
 ```math
 \begin{equation}
-a_2 = \frac{S_l}{q_v}
+a_2 = \frac{1}{q_v}
 \end{equation}
 ```
 ```math

diff --git a/parcel/parcel.jl b/parcel/parcel.jl
@@ -60,7 +60,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 = S_liq / q_vap
+    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)
@@ -144,7 +144,7 @@ function parcel_model(dY, Y, p, t)
     # Update the tendecies
     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
-    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
+    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
     dq_vap_dt = -dq_ice_dt - dq_liq_dt