From 58d931f7f61a0628dde156be01d7ea002cf541b9 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Wed, 6 Jul 2022 20:56:21 -0500
Subject: [PATCH] fix some cf issues
---
modules/util/util.py | 6 +++++-
1 file changed, 5 insertions(+), 1 deletion(-)
diff --git a/modules/util/util.py b/modules/util/util.py
index 52c4c368..a3b94576 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -663,12 +663,13 @@ rho_water = 1.
rho_ice = 0.917
-def compute_lwc_iwc(iphase, reff, opd):
+def compute_lwc_iwc(iphase, reff, opd, geo_dz):
xy_shape = iphase.shape
iphase = iphase.flatten()
reff = reff.flatten()
opd = opd.flatten()
+ geo_dz = geo_dz.flatten()
lwp_dcomp = np.zeros(reff.shape[0])
iwp_dcomp = np.zeros(reff.shape[0])
@@ -681,6 +682,9 @@ def compute_lwc_iwc(iphase, reff, opd):
iwp_dcomp[ice] = 0.667 * opd[ice] * rho_ice * reff[ice]
lwp_dcomp[no_ice] = 0.55 * opd[no_ice] * rho_water * reff[no_ice]
+ iwp_dcomp /= geo_dz
+ lwp_dcomp /= geo_dz
+
lwp_dcomp = lwp_dcomp.reshape(xy_shape)
iwp_dcomp = iwp_dcomp.reshape(xy_shape)
--
GitLab