Ignore:
Timestamp:
Jul 11, 2019 1:57:56 PM (2 years ago)
Author:
Giersch
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/basic_constants_and_equations_mod.f90

    r4084 r4088  
    2525! -----------------
    2626! $Id$
     27! Comment of barometric formula improved, function for ideal gas law revised
     28!
     29! 4084 2019-07-10 17:09:11Z knoop
    2730! Changed precomputed fractions to be variable based
    2831!
     
    236239!
    237240!--    compute density according to ideal gas law:
    238        ideal_gas_law_rho_pt_0d = p / (r_d * (1.0_wp / exner_function_invers(p)) * t)
     241       ideal_gas_law_rho_pt_0d = p / (r_d * exner_function(p) * t)
    239242
    240243    END FUNCTION ideal_gas_law_rho_pt_0d
     
    255258!
    256259!--    compute density according to ideal gas law:
    257        ideal_gas_law_rho_pt_1d = p / (r_d * (1.0_wp / exner_function_invers(p)) * t)
     260       ideal_gas_law_rho_pt_1d = p / (r_d * exner_function(p) * t)
    258261
    259262    END FUNCTION ideal_gas_law_rho_pt_1d
     
    334337! Description:
    335338! ------------
    336 !> Compute the barometric formula for scalar arguments.
     339!> Compute the barometric formula for scalar arguments. The calculation is
     340!> based on the assumption of a polytropic atmosphere and neutral
     341!> stratification, where the temperature lapse rate is g/cp.
    337342!------------------------------------------------------------------------------!
    338343    FUNCTION barometric_formula_0d( z, t_0, p_0)
     
    354359! Description:
    355360! ------------
    356 !> Compute the barometric formula for 1-D array arguments.
     361!> Compute the barometric formula for 1-D array arguments. The calculation is
     362!> based on the assumption of a polytropic atmosphere and neutral
     363!> stratification, where the temperature lapse rate is g/cp.
    357364!------------------------------------------------------------------------------!
    358365    FUNCTION barometric_formula_1d( z, t_0, p_0)
Note: See TracChangeset for help on using the changeset viewer.