|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 ! ######### 00002 SUBROUTINE DEFAULT_IDEAL_FLUX(KFORCF, KFORCT, PTIMEF, PTIMET, & 00003 PSFTH, PSFTQ, PSFCO2, & 00004 HUSTARTYPE, PUSTAR, PZ0M, PALB, & 00005 PEMIS, PTSRAD) 00006 ! ######################################################################## 00007 ! 00008 !!**** *DEFAULT_IDEAL_FLUX* - routine to set default values for the configuration for ISBA FLUX assimilation scheme 00009 !! 00010 !! PURPOSE 00011 !! ------- 00012 !! 00013 !!** METHOD 00014 !! ------ 00015 !! 00016 !! EXTERNAL 00017 !! -------- 00018 !! 00019 !! 00020 !! IMPLICIT ARGUMENTS 00021 !! ------------------ 00022 !! 00023 !! REFERENCE 00024 !! --------- 00025 !! 00026 !! 00027 !! AUTHOR 00028 !! ------ 00029 !! L. Jarlan *Meteo France* 00030 !! 00031 !! MODIFICATIONS 00032 !! ------------- 00033 !! Original 02/2005 00034 !------------------------------------------------------------------------------- 00035 ! 00036 !* 0. DECLARATIONS 00037 ! ------------ 00038 ! 00039 USE MODD_SURF_PAR, ONLY : XUNDEF 00040 USE MODD_CSTS, ONLY : XTT 00041 ! 00042 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00043 USE PARKIND1 ,ONLY : JPRB 00044 ! 00045 IMPLICIT NONE 00046 ! 00047 !* 0.1 Declarations of arguments 00048 ! ------------------------- 00049 INTEGER, INTENT(OUT) :: KFORCF 00050 INTEGER, INTENT(OUT) :: KFORCT 00051 REAL, DIMENSION(:), INTENT(OUT) :: PTIMEF 00052 REAL, DIMENSION(:), INTENT(OUT) :: PTIMET 00053 REAL, DIMENSION(:), INTENT(OUT) :: PSFTH ! hourly data of heat surface flux (W/m2) 00054 REAL, DIMENSION(:), INTENT(OUT) :: PSFTQ ! hourly data of water vapor surface flux (kg/m2/s) or (W/m2) 00055 REAL, DIMENSION(:), INTENT(OUT) :: PSFCO2 ! hourly data of CO2 surface flux (kg/m2/s) 00056 CHARACTER(LEN=5), INTENT(OUT) :: HUSTARTYPE ! type of computation for friction 00057 ! 'USTAR' 00058 ! 'Z0 ' 00059 REAL, DIMENSION(:), INTENT(OUT) :: PUSTAR ! hourly data of friction (m2/s2) 00060 REAL, INTENT(OUT) :: PZ0M ! roughness length (m) 00061 REAL, INTENT(OUT) :: PALB ! albedo (-) 00062 REAL, INTENT(OUT) :: PEMIS ! emissivity (-) 00063 REAL, DIMENSION(:),INTENT(OUT) :: PTSRAD ! radiative temperature (K) 00064 ! 00065 !* 0.2 declarations of local variables 00066 ! 00067 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00068 ! 00069 !------------------------------------------------------------------------------- 00070 IF (LHOOK) CALL DR_HOOK('DEFAULT_IDEAL_FLUX',0,ZHOOK_HANDLE) 00071 ! 00072 KFORCF = 2. 00073 KFORCT = 2. 00074 ! 00075 PTIMEF(1) = 0. 00076 PTIMET(1) = 0. 00077 PTIMEF(2) = XUNDEF 00078 PTIMET(2) = XUNDEF 00079 00080 ! 00081 !---------------------------------------------------------------------------------- 00082 ! 00083 !* 1. HOURLY surface theta flux (NFORC+1 values from 00UTC to 24UTC) 00084 ! ------------------------- 00085 ! 00086 !* unit: W/m2 00087 ! 00088 PSFTH(:) = 0. 00089 ! 00090 !---------------------------------------------------------------------------------- 00091 ! 00092 !* 2. HOURLY surface vapor mixing ratio flux (NFORC+1 values from 00UTC to 24UTC) 00093 ! -------------------------------------- 00094 ! 00095 !* unit: kg/m2/s 00096 ! 00097 PSFTQ(:) = 0. 00098 ! 00099 !---------------------------------------------------------------------------------- 00100 ! 00101 !* 4. HOURLY surface CO2 flux (NFORC+1 values from 00UTC to 24UTC) 00102 ! ----------------------- 00103 ! 00104 !* unit: kg/m2/s 00105 ! 00106 PSFCO2(:) = 0. 00107 ! 00108 !---------------------------------------------------------------------------------- 00109 ! 00110 !* 5. Type of definition for friction fluxes 00111 ! -------------------------------------- 00112 ! 00113 !* HUSTARTYPE = 'Z0 ' ! friction is defined using a roughness length formulation 00114 !* = 'USTAR' ! friction is prescribed via the friction velocity u* 00115 ! 00116 HUSTARTYPE = 'Z0 ' 00117 ! 00118 !---------------------------------------------------------------------------------- 00119 ! 00120 !* 6. Roughness length (used if XUSTARTYPE = 'Z0 ') 00121 ! ---------------- 00122 ! 00123 PZ0M = 0.01 ! unit in meters 00124 ! 00125 !---------------------------------------------------------------------------------- 00126 ! 00127 !* 6. Friction (used if XUSTARTYPE = 'USTAR') 00128 ! -------- 00129 ! 00130 !* unit: m2/s2 00131 ! 00132 PUSTAR = 0. 00133 ! 00134 ! 00135 !---------------------------------------------------------------------------------- 00136 ! 00137 !* 7. HOURLY surface radiative temperature (NFORC+1 values from 00UTC to 24UTC) 00138 ! ------------------------- 00139 ! 00140 00141 PTSRAD(:) = XTT ! radiative surface temperature, (unit is K) 00142 ! 00143 !---------------------------------------------------------------------------------- 00144 ! 00145 !* 8. Radiative fields constant and uniform values 00146 ! -------------------------------------------- 00147 ! 00148 PALB = 0. ! albedo, no unit 00149 PEMIS = 1. ! emissivity, no unit 00150 ! 00151 IF (LHOOK) CALL DR_HOOK('DEFAULT_IDEAL_FLUX',1,ZHOOK_HANDLE) 00152 ! 00153 END SUBROUTINE DEFAULT_IDEAL_FLUX
1.8.0
