|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 SUBROUTINE INIT_SURFCONSPHY 00002 ! 00003 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00004 USE PARKIND1 ,ONLY : JPRB 00005 ! 00006 ! 00007 ! 00008 IMPLICIT NONE 00009 ! 00010 * INITIALIZES THE CONSTANTS FOR THE COMMONS OF THE FLXSURF3 ROUTINE FROM 00011 * CANADIAN METEOROLOGICAL CENTER 00012 * 00013 #include "surfcon.h" 00014 #include "consphy.h" 00015 00016 * CONSTANTS RELATED TO VERTICAL DIFFUSION 00017 00018 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00019 00020 IF (LHOOK) CALL DR_HOOK('INIT_SURFCONSPHY',0,ZHOOK_HANDLE) 00021 AS = 12. 00022 ASX = 5. 00023 CI = 40. 00024 BS = 1.0 00025 BETA = 1.0 00026 FACTN = 1.2 00027 HMIN = 30. 00028 ANGMAX= 0.85 00029 RAC3 = SQRT(3.) 00030 00031 * PHYSICAL CONSTANTS 00032 00033 CPD =.100546e+4 ! J K-1 kg-1 ! specific heat of dry air 00034 CPV =.186946e+4 ! J K-1 kg-1 ! specific heat of water vapour 00035 RGASD =.28705e+3 ! J K-1 kg-1 ! gas constant for dry air 00036 RGASV =.46151e+3 ! J K-1 kg-1 ! gas constant for water vapour 00037 TRPL =.27316e+3 ! K ! triple point of water 00038 TCDK =.27315e+3 ! ! conversion from kelvin to celsius 00039 RAUW =.1e+4 ! ! density of liquid H2O 00040 EPS1 =.62194800221014 ! ! RGASD/RGASV 00041 EPS2 =.3780199778986 ! ! 1 - EPS1 00042 DELTA =.6077686814144 ! ! 1/EPS1 - 1 00043 CAPPA =.28549121795 ! ! RGASD/CPD 00044 TGL =.27316e+3 ! K ! ice temperature in the atmosphere 00045 CONSOL =.1367e+4 ! W m-2 ! solar constant 00046 GRAV =.980616e+1 ! M s-2 ! gravitational acceleration 00047 RAYT =.637122e+7 ! M ! mean radius of the earth 00048 STEFAN =.566948e-7 ! J m-2 s-1 K-4 ! Stefan-Boltzmann constant 00049 PI =.314159265359e+1 ! ! PI constant = ACOS(-1) 00050 OMEGA =.7292e-4 ! s-1 ! angular speed of rotation of the earth 00051 KNAMS =.514791 ! ! conversion from knots to m/s 00052 STLO =.6628486583943e-3 ! K s2 m-2 ! Schuman-Newell Lapse Rate 00053 KARMAN =.35 ! ! Von Karman constant 00054 RIC =.2 ! ! Critical Richardson number 00055 CHLC =.2501e+7 ! J kg-1 ! latent heat of condensation 00056 CHLF =.334e+6 ! J kg-1 ! latent heat of fusion 00057 T1S =.27316e+3 ! K ! constant used to calculate L/Cp in fcn HTVOCP 00058 T2S =.25816e+3 ! K ! constant used to calculate L/Cp in fcn HTVOCP 00059 AW =.3135012829948e+4 ! ! constant used to calculate L/Cp in fcn HTVOCP 00060 BW =.2367075766316e+1 ! ! constant used to calculate L/Cp in fcn HTVOCP 00061 AI =.2864887713087e+4 ! ! constant used to calculate L/Cp in fcn HTVOCP 00062 BI =.166093131502 ! ! constant used to calculate L/Cp in fcn HTVOCP 00063 SLP =.6666666666667e-1 ! ! constant used to calculate L/Cp in fcn HTVOCP 00064 IF (LHOOK) CALL DR_HOOK('INIT_SURFCONSPHY',1,ZHOOK_HANDLE) 00065 00066 END SUBROUTINE INIT_SURFCONSPHY
1.8.0
