Skip to content
Snippets Groups Projects
Commit 58d931f7 authored by tomrink's avatar tomrink
Browse files

fix some cf issues

parent 1a494d8e
Branches
No related tags found
No related merge requests found
......@@ -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)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment