Changeset 3361 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Oct 16, 2018 8:39:37 PM (6 years ago)

Author:

knoop

Message:

Introduced global constant rd_d_rv=0.622

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -486 +486 @@
 
     USE basic_constants_and_equations_mod,                                     &
-        ONLY:  c_p, g, lv_d_cp, l_v, magnus, rho_l, r_d, r_v
+        ONLY:  c_p, g, lv_d_cp, l_v, magnus, rho_l, r_d, r_v, rd_d_rv
 
     USE calc_mean_profile_mod,                                                 &
@@ -1999 +1999 @@
 !
 !--       Calculate vapour pressure
-          e  = surf%qv1(m) * surface_pressure / ( surf%qv1(m) + 0.622_wp )
+          e  = surf%qv1(m) * surface_pressure / ( surf%qv1(m) + rd_d_rv )
           f3 = EXP ( - surf%g_d(m) * (e_s - e) )
        ELSE
@@ -2040 +2040 @@
 !
 !--    Calculate saturation water vapor mixing ratio
-       q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
+       q_s = rd_d_rv * e_s / ( surface_pressure - e_s )
 !
 !--    In case of dewfall, set evapotranspiration to zero
@@ -2079 +2079 @@
                        / ( surf_t_surface%var_1d(m) - 29.65_wp)**2 )
 
-       dq_s_dt = 0.622_wp * e_s_dt / ( surface_pressure - e_s_dt )
+       dq_s_dt = rd_d_rv * e_s_dt / ( surface_pressure - e_s_dt )
 !
 !--    Calculate net radiation radiation without longwave outgoing flux because
@@ -2350 +2350 @@
 !
 !--       Calculate mixing ratio at saturation
-          q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
+          q_s = rd_d_rv * e_s / ( surface_pressure - e_s )
 
           resistance = surf%r_a(m) / ( surf%r_a(m) + surf%r_s(m) + 1E-5_wp )