scienceconstants.F90

Go to the documentation of this file.

module science_constants_mod
  use datadefn_mod, only : default_precision
  use optionsdatabase_mod, only : options_get_string, options_get_real
  use state_mod, only : model_state_type
  implicit none

#ifndef TEST_MODE
  private
#endif

  real(kind=DEFAULT_PRECISION) :: smallp=1.0e-14, von_karman_constant, z0, z0th, alphah, betam, betah, &
       gammam, gammah, pi, surface_vapour_mixing_ratio, cp , rlvap, rlvap_over_cp, r, r_over_cp, g,&
       convective_limit, ratio_mol_wts, rlargep

  real(kind=DEFAULT_PRECISION) :: seconds_in_a_day=86400.0
  public smallp, von_karman_constant, z0, z0th, alphah, betam, betah, gammam, gammah, pi, cp, &
       rlvap, rlvap_over_cp, r, r_over_cp, g, convective_limit, ratio_mol_wts, rlargep, initialise_science_constants, &
       seconds_in_a_day
contains

  subroutine initialise_science_constants(current_state)
    type(model_state_type), target, intent(inout) :: current_state

    von_karman_constant=options_get_real(current_state%options_database, "von_karman_constant")
    z0=options_get_real(current_state%options_database, "z0")
    z0th=options_get_real(current_state%options_database, "z0th")
    alphah=options_get_real(current_state%options_database, "alphah")
    betam=options_get_real(current_state%options_database, "betam")
    betah=options_get_real(current_state%options_database, "betah")
    gammam=options_get_real(current_state%options_database, "gammam")
    gammah=options_get_real(current_state%options_database, "gammah")
    pi=options_get_real(current_state%options_database, "pi")
    cp=options_get_real(current_state%options_database, "cp")
    rlvap=options_get_real(current_state%options_database, "rlvap")
    r=options_get_real(current_state%options_database, "r")
    g=options_get_real(current_state%options_database, "G")
    convective_limit=options_get_real(current_state%options_database, "convective_limit")
    ratio_mol_wts=options_get_real(current_state%options_database, "ratio_mol_wts")
    rlargep=options_get_real(current_state%options_database, "rlargep")

    rlvap_over_cp=rlvap/cp
    r_over_cp=r/cp
  end subroutine initialise_science_constants  
end module science_constants_mod