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, qlcrit, mprog_min=1.e-8
 
  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, qlcrit, mprog_min
  
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
 
    qlcrit=options_get_real(current_state%options_database, "qlcrit")
    ! AH - qlcrit is used in diagnostics, mprog_min is used in radar reflect
    !      these should be combined so that mprog_min is used (consistent with
    !      the UM and makes more sense)
    !mprog_min=1.e-8
    
  end subroutine initialise_science_constants  
end module science_constants_mod