diff --git a/AirSeaFluxCode.py b/AirSeaFluxCode.py index 00aadf084afb354c584102cca84b0a0f922a4814..fce5b2ade17e5f7143478f000ddaa8ff7b2399c2 100644 --- a/AirSeaFluxCode.py +++ b/AirSeaFluxCode.py @@ -524,7 +524,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, wind[ind] = np.copy(spd[ind]) u10n[ind] = wind[ind]-usr[ind]/kappa*(np.log(h_in[0, ind]/10) - psim[ind]) - u10n[u10n < 0] = np.nan + u10n = np.where(u10n < 0, np.nan, u10n) itera[ind] = np.ones(1)*it sensible = -rho*cp*usr*tsr latent = -rho*lv*usr*qsr