From b60182fe9a2d798bbd3aa5b394fae78c84d2f784 Mon Sep 17 00:00:00 2001 From: Andrew Scherer Date: Thu, 20 Mar 2025 18:06:31 -0700 Subject: [PATCH] adds salinity dependence to the saturation specific humidity calculation --- pycoare/util.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/pycoare/util.py b/pycoare/util.py index c38fbc5..3ce31d4 100644 --- a/pycoare/util.py +++ b/pycoare/util.py @@ -1,5 +1,5 @@ import numpy as np -from numpy.typing import NDArray, ArrayLike +from numpy.typing import ArrayLike, NDArray def grv(lat: ArrayLike) -> NDArray[np.float64]: @@ -63,7 +63,7 @@ def qsat(t: ArrayLike, p: ArrayLike) -> NDArray[np.float64]: return es -def qsea(t: ArrayLike, p: ArrayLike) -> NDArray[np.float64]: +def qsea(t: ArrayLike, p: ArrayLike, s: ArrayLike = 35) -> NDArray[np.float64]: """Returns saturation specific humidity at sea surface from temperature and pressure. :param t: temperature (degC) @@ -74,7 +74,7 @@ def qsea(t: ArrayLike, p: ArrayLike) -> NDArray[np.float64]: :rtype: NDArray[np.float64] """ ex = qsat(t, p) # returns ex as ndarray float - es = 0.98 * ex + es = (1 - 0.02 * s / 35) * ex qs = 622 * es / (p - 0.378 * es) return qs