Replace q_vap_saturation_from_density with q_vap_from_p_vap (#653)

* Replace q_vap_saturation_from_density with q_vap_from_p_vap

* Update Thermodynamics package version to 0.15

---------

Co-authored-by: Navid C. Constantinou <[email protected]>

Author: szy21

Project.toml

@@ -67,7 +67,7 @@ Scratch = "1"
 SeawaterPolynomials = "0.3.5"
 StaticArrays = "1"
 Statistics = "1.9"
-Thermodynamics = "0.14"
+Thermodynamics = "0.15"
 ZipFile = "0.10"
 julia = "1.10"
 

src/OceanSeaIceModels/InterfaceComputations/interface_states.jl

@@ -73,7 +73,7 @@ end
     ρₛ = convert(CT, ρₛ)
     phase = formulation.phase
     p★ = Thermodynamics.saturation_vapor_pressure(ℂₐ, Tₛ, phase)
-    q★ = Thermodynamics.q_vap_saturation_from_density(ℂₐ, Tₛ, ρₛ, p★)
+    q★ = Thermodynamics.q_vap_from_p_vap(ℂₐ, Tₛ, ρₛ, p★)
 
     # Compute saturation specific humidity according to Raoult's law
     χ_H₂O = compute_water_mole_fraction(formulation.water_mole_fraction, Sₛ)