|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 SUBROUTINE INI_CTURBS 00002 !! 00003 !!**** *INI_CTURBS* - routine to initialize the turbulence scheme 00004 !! constants. 00005 !! 00006 !! PURPOSE 00007 !! ------- 00008 ! The purpose of this routine is to initialize the turbulence 00009 ! scheme constants that are stored in module MODD_CTURB 00010 ! 00011 !! METHOD 00012 !! ------ 00013 !! The constants are set to their numerical values 00014 !! 00015 !! AUTHOR 00016 !! ------ 00017 !! V. Masson * Meteo-France * 00018 !! 00019 !! MODIFICATIONS 00020 !! ------------- 00021 !! -------------------------------------------------------------------------- 00022 ! 00023 !* 0. DECLARATIONS 00024 ! ------------ 00025 ! 00026 USE MODD_CSTS, ONLY : XKARMAN 00027 USE MODD_CANOPY_TURB, ONLY : XTKEMIN, XCEP, XCED, XALPSBL, XA0, XCMFS, XASBL,& 00028 XCSHF, XCTP 00029 ! 00030 ! 00031 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00032 USE PARKIND1 ,ONLY : JPRB 00033 ! 00034 IMPLICIT NONE 00035 ! 00036 ! --------------------------------------------------------------------------- 00037 ! 00038 ! 1. SETTING THE NUMERICAL VALUES 00039 ! ---------------------------- 00040 ! 00041 ! 1.1 Constant for dissipation of Tke 00042 ! 00043 00044 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00045 00046 IF (LHOOK) CALL DR_HOOK('INI_CTURBS',0,ZHOOK_HANDLE) 00047 XCED = 0.845 00048 ! Redelsperger-Sommeria (1981) = 0.70 00049 ! Schmidt-Schumann (1989) = 0.845 00050 ! Cheng-Canuto-Howard (2002) = 0.845 00051 ! 00052 ! 00053 ! 1.2 Constant for wind pressure-correlations 00054 ! 00055 XCEP = 2.11 00056 ! Redelsperger-Sommeria (1981) = 4. 00057 ! Schmidt-Schumann (1989) = 3.5 00058 ! Cheng-Canuto-Howard (2002) = 2.11 00059 ! 00060 ! 00061 ! 1.3 Constant a0 for wind pressure-correlations 00062 ! 00063 XA0 = 0.6 00064 ! Redelsperger-Sommeria (1981) = 0.6 00065 ! Schmidt-Schumann (1989) = 0.55 00066 ! Cheng-Canuto-Howard (2002) = 0.6 00067 ! 00068 ! 00069 ! 1.4 Constant for temperature and vapor pressure-correlations 00070 ! 00071 XCTP = 4.65 00072 ! Redelsperger-Sommeria (1981) = 4. 00073 ! Schmidt-Schumann (1989) = 3.25 00074 ! Cheng-Canuto-Howard (2002) = 4.65 00075 ! 00076 ! 00077 ! 1.5 Value related to the TKE universal function within SBL 00078 ! 00079 XALPSBL = 4.63 00080 ! Redelsperger et al 2001 = 4.63 00081 ! Wyngaard et al. 1974 = 3.75 00082 ! Stull 1988 = 4.75 00083 ! 00084 ! 00085 ! 00086 ! 2. Derivated constants 00087 ! ------------------- 00088 ! 00089 ! 2.1 Constant in fluxes equations 00090 ! 00091 XCMFS= 2./3./XCEP*(1.-XA0) !Constant for the momentum flux due to shear (RS) 00092 ! 00093 ! Redelsperger-Sommeria (1981) ......... 0.066 00094 ! Schmidt-Schumann (1989) ......... 0.086 00095 ! Cheng-Canuto-Howard (2002) ......... 0.126 00096 ! 00097 XCSHF= 2./3./XCTP !Constant for the sensible heat flux(RS) 00098 ! 00099 ! Redelsperger-Sommeria (1981) ......... 0.167 00100 ! Schmidt-Schumann (1989) ......... 0.204 00101 ! Cheng-Canuto-Howard (2002) ......... 0.143 00102 ! 00103 ! 00104 ! 2.2 Value related to the TKE universal function within SBL 00105 ! 00106 ! 00107 XASBL = 0.5*( XALPSBL**(3./2.)*XKARMAN*XCED + XKARMAN/SQRT(XALPSBL)/XCMFS ) 00108 ! Redelsperger et al 2001 00109 ! 00110 ! 00111 ! 00112 ! 3. MINIMUM VALUES 00113 ! -------------- 00114 ! 00115 XTKEMIN=1.E-4 00116 IF (LHOOK) CALL DR_HOOK('INI_CTURBS',1,ZHOOK_HANDLE) 00117 ! 00118 00119 END SUBROUTINE INI_CTURBS
1.8.0
