ini_csts.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
      SUBROUTINE ini_csts 
!     ##################
!
!!****  *INI_CSTS * - routine to initialize the module MODD_CST
!!
!!    PURPOSE
!!    -------
!       The purpose of this routine is to initialize  the physical constants
!     stored in  module MODD_CST.
!      
!
!!**  METHOD
!!    ------
!!      The physical constants are set to their numerical values 
!!     
!!
!!    EXTERNAL
!!    --------
!!      FMLOOK : to retrieve logical unit number associated to a file
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      Module MODD_CST     : contains physical constants
!!
!!    REFERENCE
!!    ---------
!!      Book2 of the documentation (module MODD_CST, routine INI_CSTS)
!!      
!!
!!    AUTHOR
!!    ------
!!      V. Ducrocq       * Meteo France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    18/05/94 
!!      J. Stein    02/01/95  add the volumic mass of liquid water
!!      J.-P. Pinty 13/12/95  add the water vapor pressure over solid ice
!!      J. Stein    29/06/97  add XTH00
!!      V. Masson   05/10/98  add XRHOLI
!!      C. Mari     31/10/00  add NDAYSEC
!!      V. Masson   01/03/03  add XCONDI
!!      A. Voldoire 01/12/09  add XTTSI, XICEC, XTTS for ESM
!!      J. Escobar  28/03/2014 for pb with emissivity/aerosol reset XSURF_TINY=1.0e-80 in real8 case
!!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_csts
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_ini_cturbs
!
USE modi_ini_ocean_csts
!
USE modi_ini_surf_csts
!
IMPLICIT NONE
!  
!-------------------------------------------------------------------------------
!
!*       1.     FUNDAMENTAL CONSTANTS
!               ---------------------
!

REAL(KIND=JPRB) :: zhook_handle

IF (lhook) CALL dr_hook('INI_CSTS',0,zhook_handle)

#ifdef SFX_MNH
#ifdef MNH_MPI_DOUBLE_PRECISION
xsurf_tiny    = 1.0e-80
#else
xsurf_tiny    = tiny(xsurf_tiny    )
#endif
#else
xsurf_tiny    = 1.0e-80
#endif
xsurf_tiny_12 = sqrt(xsurf_tiny    )
xsurf_epsilon = epsilon(xsurf_epsilon ) * 10.0

xpi         = 2.*asin(1.)
xkarman     = 0.4
xboltz      = 1.380658e-23
xlightspeed = 299792458.
xplanck     = 6.6260755e-34
xavogadro   = 6.0221367e+23
!
!-------------------------------------------------------------------------------
!
!*       2.     ASTRONOMICAL CONSTANTS
!               ----------------------
!
xday   = 86400.
xsiyea = 365.25*xday*2.*xpi/ 6.283076
xsiday = xday/(1.+xday/xsiyea)
xomega = 2.*xpi/xsiday
ndaysec = 24*3600 ! Number of seconds in a day
!
!-------------------------------------------------------------------------------!
!
!
!*       3.     TERRESTRIAL GEOIDE CONSTANTS
!               ----------------------------
!
xradius = 6371229.
xg      = 9.80665
!
!-------------------------------------------------------------------------------
!
!*       4.     REFERENCE PRESSURE
!               -------------------
!
xp00 = 1.e5
xth00 = 300.
!-------------------------------------------------------------------------------
!
!*       5.     RADIATION CONSTANTS
!               -------------------
!
!JUAN OVERFLOW XSTEFAN = 2.* XPI**5 * XBOLTZ**4 / (15.* XLIGHTSPEED**2 * XPLANCK**3)
xstefan = ( 2.* xpi**5 / 15. ) * ( (xboltz / xplanck)* xboltz ) * (xboltz/(xlightspeed*xplanck))**2 
xi0     = 1370.
!
!-------------------------------------------------------------------------------
!
!*       6.     THERMODYNAMIC CONSTANTS
!               -----------------------
!
xmd    = 28.9644e-3
xmv    = 18.0153e-3
xrd    = xavogadro * xboltz / xmd
xrv    = xavogadro * xboltz / xmv
xcpd   = 7.* xrd /2.
xcpv   = 4.* xrv
xrholw = 1000.
xrholi = 917.
xcondi = 2.22
xcl    = 4.218e+3
xci    = 2.106e+3
xtt    = 273.16
xttsi  = xtt - 1.8
xicec  = 0.5
xtts   = xtt*(1-xicec) + xttsi*xicec
xlvtt  = 2.5008e+6
xlstt  = 2.8345e+6
xlmtt  = xlstt - xlvtt
xestt  = 611.14
xgamw  = (xcl - xcpv) / xrv
xbetaw = (xlvtt/xrv) + (xgamw * xtt)
xalpw  = log(xestt) + (xbetaw /xtt) + (xgamw *log(xtt))
xgami  = (xci - xcpv) / xrv
xbetai = (xlstt/xrv) + (xgami * xtt)
xalpi  = log(xestt) + (xbetai /xtt) + (xgami *log(xtt))
!
!-------------------------------------------------------------------------------
!
!*       7.     TURBULENCE CONSTANTS
!               --------------------
!
 CALL ini_cturbs
!-------------------------------------------------------------------------------
!
!*       8.     OCEAN CONSTANTS
!               ---------------
!
 CALL ini_ocean_csts
!
!*       9.     SURFACE CONSTANTS
!               -----------------
!
 CALL ini_surf_csts
IF (lhook) CALL dr_hook('INI_CSTS',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE ini_csts