Mod_PHY____dat.f90

Go to the documentation of this file.

     module mod_phy____dat

!------------------------------------------------------------------------------+
!                                                         Sat  8-Jun-2013  MAR |
!    module Mod_PHY____dat  contains specific constants for                    |
!           MAR PHYsical    Parameterizations                                  |
!                                                                              |
!                     These constants may be slightly different in HOST Model  |
!                                 and may be modified in  PHY________INI       |
!                                                                              |
!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
!           Last Modification by H. Gallee,               Sat  8-Jun-2013      |
!                                                                              |
!------------------------------------------------------------------------------+



! General Variables
! =================

      use mod_real


      IMPLICIT NONE


! Physical Constants (may be changed in the interface with the Host Model)
! ==================

      real(kind=real8), SAVE  ::  zer0   =    0.00         !                                !    0.00
      real(kind=real8), SAVE  ::  half   =    0.50         !                                !    0.50
      real(kind=real8), SAVE  ::  un_1   =    1.00         !                                !    1.00
      real(kind=real8), SAVE  ::  pinmbr =    3.1416       ! pi                             !    3.1416
      real(kind=real8), SAVE  ::  dg2rad = 1745.33d-5      ! pi / 180                       !           
      real(kind=real8), SAVE  ::  epsq   =    1.e-6        ! Arbirary Small Value, H2O      !    1.e-6
      real(kind=real8), SAVE  ::  eps1   =    1.e-1        ! Arbirary Small Value           !    1.e-1
      real(kind=real8), SAVE  ::  eps6   =    1.e-6        ! Arbirary Small Value           !    1.e-6
      real(kind=real8), SAVE  ::  epsn   =    1.e-9        ! Arbirary Small Value           !    1.e-9
      real(kind=real8), SAVE  ::  epsp   =    1.e-12       ! Arbirary Small Value           !    1.e-12
      real(kind=real8), SAVE  ::  r_1by3 =    0.333333     ! 1 / 3                          !
      real(kind=real8), SAVE  ::  r_5by3 =    1.666666     ! 5 / 3                          !
      real(kind=real8), SAVE  ::  r_1000 =    1.e+3        !                                !    1.e+3
      real(kind=real8), SAVE  ::  ea_max =   50.           ! MAX allowed exponential Argum. ! computed by HOST
      real(kind=real8), SAVE  ::  ea_min =  -50.           ! MIN allowed exponential Argum. ! computed by HOST
      real(kind=real8), SAVE  ::  a_molv =    1.35d-5      ! Air Molecular Viscosity        !    1.35d-5 m2/s 
      real(kind=real8), SAVE  ::  rhoice =  920.00d+0      ! Ice         Specific Mass      !  920.00d+0 kg/m3
      real(kind=real8), SAVE  ::  bsnoro =  255.00d+0      ! Blown Snow  Specific Mass      !  255.00d+0 kg/m3
      real(kind=real8), SAVE  ::  lhvh2o = 2508.00d+3      ! Latent Heat of Vapor. of Rain  ! 2500.00d+3 J/kg
      real(kind=real8), SAVE  ::  lhfh2o =    3.34d+5      ! Latent Heat of Fusion of Snow  !    3.34d+5 J/kg
      real(kind=real8), SAVE  ::  lhsh2o = 2833.60d+3      ! Latent Heat of Sublim.of Snow  ! 2833.60d+3 J/kg
      real(kind=real8), SAVE  ::  cpdair = 1004.708845     ! dry air specific heat at cst p ! 1004.00    J/kg/K
      real(kind=real8), SAVE  ::  r_dair =  287.05967      ! dry air perfect gas law  cst   !  287.      J/kg/K
      real(kind=real8), SAVE  ::  rcp    =    0.285857     ! R / Cp                         ! 287./1004. -
      real(kind=real8), SAVE  ::  p0_kap =    3.730037     !                                ! 100 kPa ** (R/Cp)
      real(kind=real8), SAVE  ::  lv_cpd                   ! LhvH2O      /  CpdAir          !
      real(kind=real8), SAVE  ::  ls_cpd                   ! LhsH2O      /  CpdAir          !
      real(kind=real8), SAVE  ::  lc_cpd                   ! LhfH2O      /  CpdAir          !
      real(kind=real8), SAVE  ::  hc_wat = 4186.00d+0      ! Water Heat Capacity            ! 4186.00d+0 J/kg/K
      real(kind=real8), SAVE  ::  rhowat = 1000.00d+0      ! Water Specific Mass            ! 1000.00d+0 kg/m3
      real(kind=real8), SAVE  ::  tf_sno =  273.16         ! Snow  Melting  Point           !  273.16    K
      real(kind=real8), SAVE  ::  tf_sea =  271.2          ! Sea   Melting  Point           !  271.2     K
      real(kind=real8), SAVE  ::  stefbo =    5.67d-8      ! Stefan-Boltzman Constant       !    5.67d-8 W/m2/K4
      real(kind=real8), SAVE  ::  grav_f =    9.81         !     Gravitational  Force       !    9.81    m/s2
!     real(kind=real8), SAVE  ::  Grav_I =    0.101937     ! 1 /(Gravitational  Force)      ! 1 /9.81    s2/m
      real(kind=real8), SAVE  ::  grav_i                   ! 1 /(Gravitational  Force)      ! 1 /9.81    s2/m
      real(kind=real8), SAVE  ::  gravf2                   !    (Gravitational  Force) ** 2 !    9.81    m2/s4
      real(kind=real8), SAVE  ::  vonkrm =    0.4          ! von Karman Constant            !    0.4
      real(kind=real8), SAVE  ::  a_stab =    5.8          ! Stability  Coefficient Moment  !    5.8
      real(kind=real8), SAVE  ::  ahstab =    5.4          ! Stability  Coefficient Heat    !    5.4
      real(kind=real8), SAVE  ::  asstab =    4.0          ! Stability  Coefficient Blown * !    4.0
      real(kind=real8), SAVE  ::  r_stab =    3.0          ! Turbul.Diffusivit.Ratio K*/Km  !        

      real(kind=real8), SAVE  ::  earthr = 6371.229e3      ! Earth   Radius                 !            m
      real(kind=real8), SAVE  ::  diraxx                   ! x-Axis  Direction              !   90       degrees is the most natural choice
      real(kind=real8), SAVE  ::  sh_max                   ! Highest Domain  Grid  Point    !        
      real(kind=real8), SAVE  ::  dzamin                   ! Thinest Atmosph.Layer Thickness!        

      character(len=3), dimension(0:12)  ::  labmon  !
      data labmon /'---','Jan','Feb','Mar','Apr','May','Jun'           &
     &                  ,'Jul','Aug','Sep','Oct','Nov','Dec'/

      integer, SAVE         , dimension(0:12)  ::  njyear  ! Nb of Days since Begin of Year !
      data njyear / 0   ,   0 ,  31 ,  59 ,  90 , 120 , 151            &
     &                  , 181 , 212 , 243 , 273 , 304 , 334/

      integer, SAVE         , dimension(0:12)  ::  njleap  ! Nb of added Days for Leap Year !
      data njleap / 0   ,   0 ,   0 ,   1 ,   1 ,   1 ,   1            &
     &                  ,   1 ,   1 ,   1 ,   1 ,   1 ,   1/


!     CAUTION: values in the 3rd column are purely indicative

      end module mod_phy____dat