Changeset 3361 for palm/trunk/SOURCE/biometeorology_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/biometeorology_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -45 +45 @@
     USE arrays_3d,                                                             &
         ONLY:  hyp, p, pt, q, u, v, w
+
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  rd_d_rv
 
     USE kinds  !< to set precision of INTEGER and REAL arrays according to PALM
@@ -231 +234 @@
                    CALL calculate_pet_static(                                      &
                            pt(k,j,i) * (hyp(k) / 100000.0_wp )**0.286_wp + t_zero, &  !< Air temperature (°C)
-                           q(k,j,i) * hyp(k) / ( q(k,j,i) + 0.622_wp ) / 100._wp,  &  !< Vapor pressure (hPa)
+                           q(k,j,i) * hyp(k) / ( q(k,j,i) + rd_d_rv ) / 100._wp,   &  !< Vapor pressure (hPa)
                            SQRT( MAX( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +  &
                                       ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +  &
@@ -380 +383 @@
                     CALL calculate_pet_static(                                      &
                        pt(k,j,i) * (hyp(k) / 100000.0_wp )**0.286_wp + t_zero, &  !< Air temperature (°C)
-                       q(k,j,i) * hyp(k) / ( q(k,j,i) + 0.622_wp ) / 100.0_wp, &  !< Vapor pressure (hPa)
+                       q(k,j,i) * hyp(k) / ( q(k,j,i) + rd_d_rv ) / 100.0_wp,  &  !< Vapor pressure (hPa)
                        SQRT( MAX( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +  &
                                   ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +  &