Changeset 3449 for palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

Timestamp:

Oct 29, 2018 7:36:56 PM (5 years ago)

Author:

suehring

Message:

Branch resler -r 3439 re-integrated into current trunk: RTM 3.0, transpiration of plant canopy, output fixes in USM

File:

• palm/trunk/SOURCE/basic_constants_and_equations_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! +degc_to_k 
+! 
+! 3361 2018-10-16 20:39:37Z knoop
 ! New module (introduced with modularization of bulk cloud physics model)
 !
@@ -43 +46 @@
     IMPLICIT NONE
 
-    REAL(wp), PARAMETER ::  pi = 3.141592654_wp  !< PI
-
     REAL(wp), PARAMETER ::  c_p = 1005.0_wp                           !< heat capacity of dry air (J kg-1 K-1)
+    REAL(wp), PARAMETER ::  degc_to_k = 273.15_wp                     !< temperature (in K) of 0 deg C (K)
     REAL(wp), PARAMETER ::  g = 9.81_wp                               !< gravitational acceleration (m s-2)
     REAL(wp), PARAMETER ::  kappa = 0.4_wp                            !< von Karman constant
@@ -55 +57 @@
     REAL(wp), PARAMETER ::  molecular_weight_of_nh4no3 = 0.08004_wp   !< mol. m. ammonium sulfate (kg mol-1)
     REAL(wp), PARAMETER ::  molecular_weight_of_water = 0.01801528_wp !< mol. m. H2O (kg mol-1)
+    REAL(wp), PARAMETER ::  pi = 3.141592654_wp                       !< PI
     REAL(wp), PARAMETER ::  rho_l = 1.0E3_wp                          !< density of water (kg m-3)
     REAL(wp), PARAMETER ::  rho_nacl = 2165.0_wp                      !< density of NaCl (kg m-3)
@@ -61 +64 @@
     REAL(wp), PARAMETER ::  r_d = 287.0_wp                            !< sp. gas const. dry air (J kg-1 K-1)
     REAL(wp), PARAMETER ::  r_v = 461.51_wp                           !< sp. gas const. water vapor (J kg-1 K-1)
-    REAL(wp), PARAMETER ::  sigma_sb = 5.67E-08_wp                    !< Stefan-Boltzmann constant
+    REAL(wp), PARAMETER ::  sigma_sb = 5.67037E-08_wp                 !< Stefan-Boltzmann constant
     REAL(wp), PARAMETER ::  solar_constant = 1368.0_wp                !< solar constant at top of atmosphere
     REAL(wp), PARAMETER ::  vanthoff_nacl = 2.0_wp                    !< van't Hoff factor for NaCl
@@ -143 +146 @@
 !
 !--    Saturation vapor pressure for a specific temperature:
-       magnus_0d =  611.2_wp * EXP( 17.62_wp * ( t - 273.15_wp ) /             &
+       magnus_0d =  611.2_wp * EXP( 17.62_wp * ( t - degc_to_k ) /             &
                                                    ( t - 29.65_wp  ) )
 
@@ -163 +166 @@
 !
 !--    Saturation vapor pressure for a specific temperature:
-       magnus_1d =  611.2_wp * EXP( 17.62_wp * ( t - 273.15_wp ) /             &
+       magnus_1d =  611.2_wp * EXP( 17.62_wp * ( t - degc_to_k ) /             &
                                                ( t - 29.65_wp  ) )