scale_const.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_const
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: const_setup

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(RP), public            :: const_pi      = 3.14159265358979_rp 
  real(RP), public            :: const_d2r
  real(RP), public            :: const_eps     = 1.e-16_rp           
  real(RP), public            :: const_eps1    = 0.99999999999999_rp 
  real(RP), public            :: const_huge    = 1.e+30_rp           
  !$acc declare create(CONST_PI, CONST_EPS, CONST_EPS1, CONST_HUGE)

  integer,  public, parameter :: const_undef2  = -32768
  real(SP), public, parameter :: const_undef4  = -9.9999e30          
  real(DP), public, parameter :: const_undef8  = -9.9999d30          
  real(RP), public            :: const_undef

  ! adopted constants
  real(RP), public            :: const_radius  = 6.37122e+6_rp       
  real(RP), public            :: const_ohm     = 7.2920e-5_rp        
  real(RP), public            :: const_grav    = 9.80665_rp          

  ! physical constants
  real(RP), public, parameter :: const_stb     = 5.670373e-8_rp      
  real(RP), public, parameter :: const_karman  = 0.4_rp              
  real(RP), public, parameter :: const_r       = 8.31436_rp          

  ! dry air constants
  real(RP), public            :: const_mdry    =  28.966_rp          
  real(RP), public            :: const_rdry    =  287.04_rp          
  real(RP), public            :: const_cpdry   = 1004.64_rp          
  real(RP), public            :: const_cvdry
  real(RP), public            :: const_laps    = 6.5e-3_rp           
  real(RP), public            :: const_lapsdry

  ! water constants
  real(RP), public            :: const_mvap    =  18.016_rp          
  real(RP), public, parameter :: const_rvap    =  461.50_rp          
  real(RP), public, parameter :: const_cpvap   = 1846.00_rp          
  real(RP), public            :: const_cvvap
  real(RP), public, parameter :: const_cl      = 4218.0_rp           
  real(RP), public, parameter :: const_ci      = 2106.0_rp           

  real(RP), public            :: const_epsvap
  real(RP), public            :: const_epstvap

  real(RP), public, parameter :: const_emelt   = 3.4e5_rp
  real(RP), public, parameter :: const_tmelt   = 273.15_rp

  real(RP), public, parameter :: const_lhv0    = 2.501e+6_rp         
  real(RP), public            :: const_lhv00
  real(RP), public, parameter :: const_lhs0    = 2.834e+6_rp         
  real(RP), public            :: const_lhs00
  real(RP), public            :: const_lhf0
  real(RP), public            :: const_lhf00
  real(RP), public, parameter :: const_psat0   =  610.78_rp          
  real(RP), public, parameter :: const_dwatr   = 1000.0_rp           
  real(RP), public, parameter :: const_dice    =  916.8_rp           

  real(RP), public            :: const_sound

  real(RP), public            :: const_pstd    = 101325.0_rp         
  real(RP), public            :: const_pre00   = 100000.0_rp         
  real(RP), public            :: const_tstd    = 288.15_rp           
  real(RP), public, parameter :: const_tem00   = 273.15_rp           
  real(RP), public, parameter :: const_ppm     = 1.e-6_rp            

  integer,  public            :: const_i_lw    = 1
  integer,  public            :: const_i_sw    = 2

  character(len=H_SHORT), public :: const_thermodyn_type = 'EXACT'

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  logical, private :: initialized = .false.
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine const_setup
    use scale_prc, only: &
       prc_abort
    implicit none

    real(RP) :: CONST_SmallPlanetFactor = 1.0_rp 

    namelist / param_const / &
       const_radius,           &
       const_ohm,              &
       const_grav,             &
       const_rdry,             &
       const_cpdry,            &
       const_laps,             &
       const_pstd,             &
       const_pre00,            &
       const_tstd,             &
       const_thermodyn_type,   &
       const_smallplanetfactor

    integer :: ierr
    !---------------------------------------------------------------------------

    if ( initialized ) return
    initialized = .true.

    log_newline
    log_info("CONST_setup",*) 'Setup'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_const,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("CONST_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("CONST_setup",*) 'Not appropriate names in namelist PARAM_CONST. Check!'
       call prc_abort
    endif
    log_nml(param_const)

    if    ( rp == sp ) then
       const_undef = real(const_undef4,kind=rp)
    elseif( rp == dp ) then
       const_undef = real(const_undef8,kind=rp)
    else
       log_error("CONST_setup",*) 'unsupported precision: ', rp
       call prc_abort
    endif

    const_pi      = 4.0_rp * atan( 1.0_rp )
    const_d2r     = const_pi / 180.0_rp
    const_eps     =          epsilon(0.0_rp)
    const_eps1    = 1.0_rp - epsilon(0.0_rp)
    const_huge    =             huge(0.0_rp)

    const_radius  = const_radius / const_smallplanetfactor
    const_ohm     = const_ohm    * const_smallplanetfactor

    const_cvdry   = const_cpdry - const_rdry
    const_lapsdry = const_grav / const_cpdry

    const_cvvap   = const_cpvap - const_rvap
    const_epsvap  = const_rdry / const_rvap
    const_epstvap = 1.0_rp / const_epsvap - 1.0_rp

    const_lhf0    = const_lhs0 - const_lhv0

    const_lhv00   = const_lhv0 - ( const_cpvap - const_cl ) * const_tem00
    const_lhs00   = const_lhs0 - ( const_cpvap - const_ci ) * const_tem00
    const_lhf00   = const_lhf0 - ( const_cl    - const_ci ) * const_tem00

    const_sound = sqrt( const_cpdry * const_rdry / ( const_cpdry - const_rdry ) * const_tem00 )

    log_newline
    log_info("CONST_setup",*) 'List of constants '
    log_info_cont(*) 'PI                                                : PI      = ', const_pi
    log_info_cont(*) 'Small number                                      : EPS     = ', const_eps
    log_info_cont(*) 'Small number (1-EPS)                              : EPS1    = ', const_eps1
    log_info_cont(*) 'Huge  number                                      : HUGE    = ', const_huge
    log_info_cont(*) 'undefined number (INT2)                           : UNDEF2  = ', const_undef2
    log_info_cont(*) 'undefined number (REAL,general use)               : UNDEF   = ', const_undef
    log_info_cont(*) 'undefined number (REAL4)                          : UNDEF4  = ', const_undef4
    log_info_cont(*) 'undefined number (REAL8)                          : UNDEF8  = ', const_undef8

    log_info_cont(*) 'radius of the planet                          [m] : RADIUS  = ', const_radius
    log_info_cont(*) 'angular velocity of the planet              [1/s] : OHM     = ', const_ohm
    log_info_cont(*) 'standard acceleration of gravity           [m/s2] : GRAV    = ', const_grav

    log_info_cont(*) 'Stefan-Boltzman constant                [W/m2/K4] : STB     = ', const_stb
    log_info_cont(*) 'von Karman constant                               : KARMAN  = ', const_karman
    log_info_cont(*) 'universal gas constant                  [J/mol/K] : R       = ', const_r

    log_info_cont(*) 'mass weight (dry air)                     [g/mol] : Mdry    = ', const_mdry
    log_info_cont(*) 'specific gas constant (dry air)          [J/kg/K] : Rdry    = ', const_rdry
    log_info_cont(*) 'specific heat (dry air, const. pressure) [J/kg/K] : CPdry   = ', const_cpdry
    log_info_cont(*) 'specific heat (dry air, const. volume)   [J/kg/K] : Cvdry   = ', const_cvdry
    log_info_cont(*) 'lapse rate of ISA                           [K/m] : LAPS    = ', const_laps
    log_info_cont(*) 'dry adiabatic lapse rate                    [K/m] : LAPSdry = ', const_lapsdry

    log_info_cont(*) 'mass weight (water vapor)                 [g/mol] : Mvap    = ', const_mvap
    log_info_cont(*) 'specific gas constant (water vapor)      [J/kg/K] : Rvap    = ', const_rvap
    log_info_cont(*) 'specific heat (vapor, const. pressure)   [J/kg/K] : CPvap   = ', const_cpvap
    log_info_cont(*) 'specific heat (vapor, const. volume)     [J/kg/K] : CVvap   = ', const_cvvap
    log_info_cont(*) 'specific heat (liquid water)             [J/kg/K] : CL      = ', const_cl
    log_info_cont(*) 'specific heat (ice)                      [J/kg/K] : CI      = ', const_ci
    log_info_cont(*) 'Rdry / Rvap                                       : EPSvap  = ', const_epsvap
    log_info_cont(*) '1 / EPSvap - 1                                    : EPSTvap = ', const_epstvap

    log_info_cont(*) 'latent heat of vaporizaion at 0C           [J/kg] : LHV0    = ', const_lhv0
    log_info_cont(*) 'latent heat of sublimation at 0C           [J/kg] : LHS0    = ', const_lhs0
    log_info_cont(*) 'latent heat of fusion      at 0C           [J/kg] : LHF0    = ', const_lhf0
    log_info_cont(*) 'latent heat of vaporizaion at 0K           [J/kg] : LHV00   = ', const_lhv00
    log_info_cont(*) 'latent heat of sublimation at 0K           [J/kg] : LHS00   = ', const_lhs00
    log_info_cont(*) 'latent heat of fusion      at 0K           [J/kg] : LHF00   = ', const_lhf00
    log_info_cont(*) 'Thermodynamics calculation type : ', trim(const_thermodyn_type)
    log_info_cont(*) 'saturate pressure of water vapor at 0C       [Pa] : PSAT0   = ', const_psat0
    log_info_cont(*) 'density of water                          [kg/m3] : DWATR   = ', const_dwatr
    log_info_cont(*) 'density of ice                            [kg/m3] : DICE    = ', const_dice

    log_info_cont(*) 'speed of sound (dry air at 0C)              [m/s] : SOUND   = ', const_sound
    log_info_cont(*) 'standard pressure                            [Pa] : Pstd    = ', const_pstd
    log_info_cont(*) 'pressure reference                           [Pa] : PRE00   = ', const_pre00
    log_info_cont(*) 'standard temperature (15C)                    [K] : Tstd    = ', const_tstd
    log_info_cont(*) 'temperature reference (0C)                    [K] : TEM00   = ', const_tem00

    return
  end subroutine const_setup

end module scale_const