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

Timestamp:

Jul 11, 2019 1:57:56 PM (4 years ago)

Author:

Giersch

Message:

Pressure and density profile calculations revised using basic functions, comments improved, function for ideal gas law revised

File:

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

palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Comment of barometric formula improved, function for ideal gas law revised 
+! 
+! 4084 2019-07-10 17:09:11Z knoop
 ! Changed precomputed fractions to be variable based
 ! 
@@ -236 +239 @@
 !
 !--    compute density according to ideal gas law:
-       ideal_gas_law_rho_pt_0d = p / (r_d * (1.0_wp / exner_function_invers(p)) * t)
+       ideal_gas_law_rho_pt_0d = p / (r_d * exner_function(p) * t)
 
     END FUNCTION ideal_gas_law_rho_pt_0d
@@ -255 +258 @@
 !
 !--    compute density according to ideal gas law:
-       ideal_gas_law_rho_pt_1d = p / (r_d * (1.0_wp / exner_function_invers(p)) * t)
+       ideal_gas_law_rho_pt_1d = p / (r_d * exner_function(p) * t)
 
     END FUNCTION ideal_gas_law_rho_pt_1d
@@ -334 +337 @@
 ! Description:
 ! ------------
-!> Compute the barometric formula for scalar arguments.
+!> Compute the barometric formula for scalar arguments. The calculation is 
+!> based on the assumption of a polytropic atmosphere and neutral 
+!> stratification, where the temperature lapse rate is g/cp. 
 !------------------------------------------------------------------------------!
     FUNCTION barometric_formula_0d( z, t_0, p_0)
@@ -354 +359 @@
 ! Description:
 ! ------------
-!> Compute the barometric formula for 1-D array arguments.
+!> Compute the barometric formula for 1-D array arguments. The calculation is 
+!> based on the assumption of a polytropic atmosphere and neutral 
+!> stratification, where the temperature lapse rate is g/cp.
 !------------------------------------------------------------------------------!
     FUNCTION barometric_formula_1d( z, t_0, p_0)