!copyright (C) 2001  MSC-RPN COMM  %%%RPNPHY%%%
*
      REAL CPD, CPV, RGASD, RGASV, TRPL, TCDK, RAUW, EPS1, EPS2
      REAL DELTA, CAPPA, TGL, CONSOL, GRAV, RAYT, STEFAN, PI
      REAL OMEGA
      REAL KNAMS, STLO, KARMAN, RIC, CHLC, CHLF
      REAL T1S, T2S, AW, BW, AI, BI, SLP
      LOGICAL INIT
*
      COMMON/CTESPHY/ INIT, CPD, CPV, RGASD, RGASV, TRPL, TCDK, RAUW,
     $                EPS1, EPS2, DELTA, CAPPA, TGL, CONSOL,
     $                GRAV, RAYT, STEFAN, PI, OMEGA,
     $                KNAMS, STLO, KARMAN, RIC, CHLC, CHLF,
     $                T1S, T2S, AW, BW, AI, BI, SLP
*
#if defined(DOC)
*COMMON BLOCK /CTESPHY/
*          contains all the constants read in from file of constants
*          for physics code in routine INCTPHY.
*
* INIT     logical variable to indicate if the file "CONSTANTES" were
*          read before
* CPD      =.100546e+4 J K-1 kg-1; specific heat of dry air
* CPV      =.186946e+4 J K-1 kg-1; specific heat of water vapour
* RGASD    =.28705e+3 J K-1 kg-1; gas constant for dry air
* RGASV    =.46151e+3 J K-1 kg-1; gas constant for water vapour
* TRPL     =.27316e+3 K; triple point of water
* TCDK     =.27315e+3; conversion from kelvin to celsius
* RAUW     =.1e+4; density of liquid H2O
* EPS1     =.62194800221014 ; RGASD/RGASV
* EPS2     =.3780199778986 ; 1 - EPS1
* DELTA    =.6077686814144 ; 1/EPS1 - 1
* CAPPA    =.28549121795; RGASD/CPD
* TGL      =.27316e+3 K; ice temperature in the atmosphere
* CONSOL   =.1367e+4 W m-2; solar constant
* GRAV     =.980616e+1 M s-2; gravitational acceleration
* RAYT     =.637122e+7 M; mean radius of the earth
* STEFAN   =.566948e-7 J m-2 s-1 K-4; Stefan-Boltzmann constant
* PI       =.314159265359e+1; PI constant = ACOS(-1)
* OMEGA    =.7292e-4s-1; angular speed of rotation of the earth
* KNAMS    =.514791; conversion from knots to m/s
* STLO     =.6628486583943e-3 K s2 m-2; Schuman-Newell Lapse Rate
* KARMAN   =.35; Von Karman constant
* RIC      =.2; Critical Richardson number
* CHLC     =.2501e+7 J kg-1; latent heat of condensation
* CHLF     =.334e+6 J kg-1; latent heat of fusion
* T1S      =.27316e+3 K; constant used to calculate L/Cp in fcn
*          HTVOCP
* T2S      =.25816e+3 K; constant used to calculate L/Cp in fcn
*          HTVOCP
* AW       =.3135012829948e+4; constant used to calculate L/Cp in fcn
*          HTVOCP
* BW       =.2367075766316e+1; constant used to calculate L/Cp in fcn
*          HTVOCP
* AI       =.2864887713087e+4; constant used to calculate L/Cp in fcn
*          HTVOCP
* BI       =.166093131502; constant used to calculate L/Cp in fcn
*          HTVOCP
* SLP      =.6666666666667e-1; constant used to calculate L/Cp in fcn
*          HTVOCP
*
#endif