Skip to content
Snippets Groups Projects
Commit 7cfd9be4 authored by tomrink's avatar tomrink
Browse files

minor...

parent f8275f6a
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 15 additions and 0 deletions
modules/util/util.py
+ 15
− 0
View file @ 7cfd9be4
...@@ -662,6 +662,21 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain ...@@ -662,6 +662,21 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
rho_water = 1. rho_water = 1.
rho_ice = 0.917 rho_ice = 0.917
# real(kind=real4), parameter:: Rho_Water = 1.0 !g / m ^ 3
# real(kind=real4), parameter:: Rho_Ice = 0.917 !g / m ^ 3
#
# !--- compute
# cloud
# water
# path
# if (Iphase == 0) then
# Cwp_Dcomp(Elem_Idx, Line_Idx) = 0.55 * Tau * Reff * Rho_Water
# Lwp_Dcomp(Elem_Idx, Line_Idx) = 0.55 * Tau * Reff * Rho_Water
# else
# Cwp_Dcomp(Elem_Idx, Line_Idx) = 0.667 * Tau * Reff * Rho_Ice
# Iwp_Dcomp(Elem_Idx, Line_Idx) = 0.667 * Tau * Reff * Rho_Ice
# endif
def compute_lwc_iwc(iphase, reff, opd, geo_dz): def compute_lwc_iwc(iphase, reff, opd, geo_dz):
xy_shape = iphase.shape xy_shape = iphase.shape
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment