|
SURFEX v7.3
General documentation of Surfex
|
|
SURFEX v7.3
General documentation of Surfex
|
00001 ! ############################################################################### 00002 SUBROUTINE COUPLING_WATFLUX_n(HPROGRAM, HCOUPLING, & 00003 PTSTEP, KYEAR, KMONTH, KDAY, PTIME, KI, KSV, KSW, PTSUN, PZENITH, PZENITH2, & 00004 PAZIM, PZREF, PUREF, PZS, PU, PV, PQA, PTA, PRHOA, PSV, PCO2, HSV, & 00005 PRAIN, PSNOW, PLW, PDIR_SW, PSCA_SW, PSW_BANDS, PPS, PPA, & 00006 PSFTQ, PSFTH, PSFTS, PSFCO2, PSFU, PSFV, & 00007 PTRAD, PDIR_ALB, PSCA_ALB, PEMIS, & 00008 PPEW_A_COEF, PPEW_B_COEF, & 00009 PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF, & 00010 HTEST ) 00011 ! ############################################################################### 00012 ! 00013 !!**** *COUPLING_WATFLUX_n * - Driver for WATER_FLUX scheme for lakes 00014 !! 00015 !! PURPOSE 00016 !! ------- 00017 ! 00018 !!** METHOD 00019 !! ------ 00020 !! 00021 !! REFERENCE 00022 !! --------- 00023 !! 00024 !! 00025 !! AUTHOR 00026 !! ------ 00027 !! V. Masson 00028 !! 00029 !! MODIFICATIONS 00030 !! ------------- 00031 !! Original 01/2004 00032 !! B. Decharme 03/2009 TS_WATER could change during a run => ALB and EMIS 00033 !! V. Masson 05/2009 Implicitation of momentum fluxes 00034 !! B. Decharme 09/2009 Radiative properties at time t+1 in order to close 00035 ! the energy budget between surfex and the atmosphere 00036 !! B. Decharme 01/2010 Add XTT 00037 !! B. Decharme 09/2012 New wind implicitation 00038 !!---------------------------------------------------------------------------- 00039 ! 00040 ! 00041 USE MODD_CSTS, ONLY : XRD, XCPD, XP00, XTT, XDAY, XTTS 00042 USE MODD_SURF_PAR, ONLY : XUNDEF 00043 USE MODD_SURF_ATM, ONLY : LCPL_ESM, CIMPLICIT_WIND 00044 USE MODD_WATER_PAR 00045 ! 00046 USE MODD_WATFLUX_n, ONLY : CWAT_ALB,XTS, XZ0, XDIR_ALB, XSCA_ALB, XEMIS, TTIME, & 00047 LINTERPOL_TS, XTICE 00048 USE MODD_CH_WATFLUX_n, ONLY : CSV, CCH_DRY_DEP, XDEP, NBEQ, NSV_CHSBEG, NSV_CHSEND,& 00049 NSV_DSTBEG, NSV_DSTEND, NAEREQ, NDSTEQ, NSLTEQ, & 00050 NSV_AERBEG, NSV_AEREND, NSV_SLTBEG, NSV_SLTEND 00051 ! 00052 USE MODI_WATER_FLUX 00053 USE MODI_ADD_FORECAST_TO_DATE_SURF 00054 USE MODI_DIAG_INLINE_WATFLUX_n 00055 USE MODI_CH_AER_DEP 00056 USE MODI_CH_DEP_WATER 00057 USE MODI_DSLT_DEP 00058 USE MODI_UPDATE_RAD_SEAWAT 00059 USE MODI_INTERPOL_TS_WATER_MTH 00060 ! 00061 USE MODE_DSLT_SURF 00062 USE MODD_SLT_n, ONLY: XEMISRADIUS_SLT, XEMISSIG_SLT 00063 USE MODD_DST_n, ONLY: XEMISRADIUS_DST, XEMISSIG_DST 00064 USE MODD_DST_SURF 00065 USE MODD_SLT_SURF 00066 ! 00067 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK 00068 USE PARKIND1 ,ONLY : JPRB 00069 ! 00070 USE MODI_ABOR1_SFX 00071 ! 00072 USE MODI_COUPLING_ICEFLUX_n 00073 ! 00074 IMPLICIT NONE 00075 ! 00076 !* 0.1 declarations of arguments 00077 ! 00078 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! program calling surf. schemes 00079 CHARACTER(LEN=1), INTENT(IN) :: HCOUPLING ! type of coupling 00080 ! 'E' : explicit 00081 ! 'I' : implicit 00082 INTEGER, INTENT(IN) :: KYEAR ! current year (UTC) 00083 INTEGER, INTENT(IN) :: KMONTH ! current month (UTC) 00084 INTEGER, INTENT(IN) :: KDAY ! current day (UTC) 00085 REAL, INTENT(IN) :: PTIME ! current time since midnight (UTC, s) 00086 INTEGER, INTENT(IN) :: KI ! number of points 00087 INTEGER, INTENT(IN) :: KSV ! number of scalars 00088 INTEGER, INTENT(IN) :: KSW ! number of short-wave spectral bands 00089 REAL, DIMENSION(KI), INTENT(IN) :: PTSUN ! solar time (s from midnight) 00090 REAL, INTENT(IN) :: PTSTEP ! atmospheric time-step (s) 00091 REAL, DIMENSION(KI), INTENT(IN) :: PZREF ! height of T,q forcing (m) 00092 REAL, DIMENSION(KI), INTENT(IN) :: PUREF ! height of wind forcing (m) 00093 ! 00094 REAL, DIMENSION(KI), INTENT(IN) :: PTA ! air temperature forcing (K) 00095 REAL, DIMENSION(KI), INTENT(IN) :: PQA ! air humidity forcing (kg/m3) 00096 REAL, DIMENSION(KI), INTENT(IN) :: PRHOA ! air density (kg/m3) 00097 REAL, DIMENSION(KI,KSV),INTENT(IN) :: PSV ! scalar variables 00098 ! ! chemistry: first char. in HSV: '#' (molecule/m3) 00099 ! ! 00100 CHARACTER(LEN=6), DIMENSION(KSV),INTENT(IN):: HSV ! name of all scalar variables 00101 REAL, DIMENSION(KI), INTENT(IN) :: PU ! zonal wind (m/s) 00102 REAL, DIMENSION(KI), INTENT(IN) :: PV ! meridian wind (m/s) 00103 REAL, DIMENSION(KI,KSW),INTENT(IN) :: PDIR_SW ! direct solar radiation (on horizontal surf.) 00104 ! ! (W/m2) 00105 REAL, DIMENSION(KI,KSW),INTENT(IN) :: PSCA_SW ! diffuse solar radiation (on horizontal surf.) 00106 ! ! (W/m2) 00107 REAL, DIMENSION(KSW),INTENT(IN) :: PSW_BANDS ! mean wavelength of each shortwave band (m) 00108 REAL, DIMENSION(KI), INTENT(IN) :: PZENITH ! zenithal angle at t (radian from the vertical) 00109 REAL, DIMENSION(KI), INTENT(IN) :: PZENITH2 ! zenithal angle at t+1 (radian from the vertical) 00110 REAL, DIMENSION(KI), INTENT(IN) :: PAZIM ! azimuthal angle (radian from North, clockwise) 00111 REAL, DIMENSION(KI), INTENT(IN) :: PLW ! longwave radiation (on horizontal surf.) 00112 ! ! (W/m2) 00113 REAL, DIMENSION(KI), INTENT(IN) :: PPS ! pressure at atmospheric model surface (Pa) 00114 REAL, DIMENSION(KI), INTENT(IN) :: PPA ! pressure at forcing level (Pa) 00115 REAL, DIMENSION(KI), INTENT(IN) :: PZS ! atmospheric model orography (m) 00116 REAL, DIMENSION(KI), INTENT(IN) :: PCO2 ! CO2 concentration in the air (kg/m3) 00117 REAL, DIMENSION(KI), INTENT(IN) :: PSNOW ! snow precipitation (kg/m2/s) 00118 REAL, DIMENSION(KI), INTENT(IN) :: PRAIN ! liquid precipitation (kg/m2/s) 00119 ! 00120 ! 00121 REAL, DIMENSION(KI), INTENT(OUT) :: PSFTH ! flux of heat (W/m2) 00122 REAL, DIMENSION(KI), INTENT(OUT) :: PSFTQ ! flux of water vapor (kg/m2/s) 00123 REAL, DIMENSION(KI), INTENT(OUT) :: PSFU ! zonal momentum flux (Pa) 00124 REAL, DIMENSION(KI), INTENT(OUT) :: PSFV ! meridian momentum flux (Pa) 00125 REAL, DIMENSION(KI), INTENT(OUT) :: PSFCO2 ! flux of CO2 (kg/m2/s) 00126 REAL, DIMENSION(KI,KSV),INTENT(OUT):: PSFTS ! flux of scalar var. (kg/m2/s) 00127 ! 00128 REAL, DIMENSION(KI), INTENT(OUT) :: PTRAD ! radiative temperature (K) 00129 REAL, DIMENSION(KI,KSW),INTENT(OUT):: PDIR_ALB! direct albedo for each spectral band (-) 00130 REAL, DIMENSION(KI,KSW),INTENT(OUT):: PSCA_ALB! diffuse albedo for each spectral band (-) 00131 REAL, DIMENSION(KI), INTENT(OUT) :: PEMIS ! emissivity (-) 00132 ! 00133 REAL, DIMENSION(KI), INTENT(IN) :: PPEW_A_COEF! implicit coefficients (m2s/kg) 00134 REAL, DIMENSION(KI), INTENT(IN) :: PPEW_B_COEF! needed if HCOUPLING='I' (m/s) 00135 REAL, DIMENSION(KI), INTENT(IN) :: PPET_A_COEF 00136 REAL, DIMENSION(KI), INTENT(IN) :: PPEQ_A_COEF 00137 REAL, DIMENSION(KI), INTENT(IN) :: PPET_B_COEF 00138 REAL, DIMENSION(KI), INTENT(IN) :: PPEQ_B_COEF 00139 CHARACTER(LEN=2), INTENT(IN) :: HTEST ! must be equal to 'OK' 00140 ! 00141 !* 0.2 declarations of local variables 00142 ! 00143 REAL, DIMENSION(KI) :: ZEXNA ! Exner function at forcing level 00144 REAL, DIMENSION(KI) :: ZEXNS ! Exner function at surface level 00145 REAL, DIMENSION(KI) :: ZWIND ! Wind 00146 REAL, DIMENSION(KI) :: ZCD ! Drag coefficient 00147 REAL, DIMENSION(KI) :: ZCDN ! Neutral Drag coefficient 00148 REAL, DIMENSION(KI) :: ZCH ! Heat transfer coefficient 00149 REAL, DIMENSION(KI) :: ZRI ! Richardson number 00150 REAL, DIMENSION(KI) :: ZHU ! Near surface relative humidity 00151 REAL, DIMENSION(KI) :: ZRESA_WATER! aerodynamical resistance 00152 REAL, DIMENSION(KI) :: ZUSTAR ! friction velocity (m/s) 00153 REAL, DIMENSION(KI) :: ZUSTAR2 ! square of friction velocity (m2/s2) 00154 REAL, DIMENSION(KI) :: ZZ0H ! heat roughness length over sea 00155 REAL, DIMENSION(KI) :: ZQSAT ! humidity at saturation 00156 REAL, DIMENSION(KI) :: ZQA ! specific humidity (kg/kg) 00157 REAL, DIMENSION(KI) :: ZEMIS ! Emissivity at time t 00158 REAL, DIMENSION(KI) :: ZTRAD ! Radiative temperature at time t 00159 REAL, DIMENSION(KI) :: ZSFTH_ICE ! Sea ice flux of heat 00160 REAL, DIMENSION(KI) :: ZSFTQ_ICE ! Sea ice flux of ice sublimation 00161 ! 00162 REAL, DIMENSION(KI,KSW) :: ZDIR_ALB ! Direct albedo at time t 00163 REAL, DIMENSION(KI,KSW) :: ZSCA_ALB ! Diffuse albedo at time t 00164 ! 00165 REAL :: ZCONVERTFACM0_SLT, ZCONVERTFACM0_DST 00166 REAL :: ZCONVERTFACM3_SLT, ZCONVERTFACM3_DST 00167 REAL :: ZCONVERTFACM6_SLT, ZCONVERTFACM6_DST 00168 ! 00169 INTEGER :: ISWB ! number of shortwave spectral bands 00170 INTEGER :: JSWB ! loop counter on shortwave spectral bands 00171 REAL(KIND=JPRB) :: ZHOOK_HANDLE 00172 ! 00173 !------------------------------------------------------------------------------------- 00174 ! Preliminaries: 00175 !------------------------------------------------------------------------------------- 00176 IF (LHOOK) CALL DR_HOOK('COUPLING_WATFLUX_N',0,ZHOOK_HANDLE) 00177 IF (HTEST/='OK') THEN 00178 CALL ABOR1_SFX('COUPLING_WATFLUXN: FATAL ERROR DURING ARGUMENT TRANSFER') 00179 END IF 00180 !------------------------------------------------------------------------------------- 00181 ! 00182 ZEXNA (:) = XUNDEF 00183 ZEXNS (:) = XUNDEF 00184 ZWIND (:) = XUNDEF 00185 ZCD (:) = XUNDEF 00186 ZCDN (:) = XUNDEF 00187 ZCH (:) = XUNDEF 00188 ZRI (:) = XUNDEF 00189 ZHU (:) = XUNDEF 00190 ZRESA_WATER(:) = XUNDEF 00191 ZUSTAR (:) = XUNDEF 00192 ZUSTAR2 (:) = XUNDEF 00193 ZZ0H (:) = XUNDEF 00194 ZQSAT (:) = XUNDEF 00195 ZEMIS (:) = XUNDEF 00196 ZTRAD (:) = XUNDEF 00197 ZDIR_ALB (:,:) = XUNDEF 00198 ZSCA_ALB (:,:) = XUNDEF 00199 ! 00200 IF(LCPL_ESM)THEN 00201 ZSFTQ_ICE(:) = XUNDEF 00202 ZSFTH_ICE(:) = XUNDEF 00203 ENDIF 00204 ! 00205 !------------------------------------------------------------------------------------- 00206 ! 00207 ZEXNS(:) = (PPS(:)/XP00)**(XRD/XCPD) 00208 ZEXNA(:) = (PPA(:)/XP00)**(XRD/XCPD) 00209 ! 00210 ZWIND(:) = SQRT(PU(:)**2+PV(:)**2) 00211 ! 00212 PSFTS(:,:) = 0. 00213 ! 00214 ZHU = 1. 00215 ! 00216 ZQA(:) = PQA(:) / PRHOA(:) 00217 ! 00218 ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 00219 ! Time evolution 00220 ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 00221 ! 00222 TTIME%TIME = TTIME%TIME + PTSTEP 00223 CALL ADD_FORECAST_TO_DATE_SURF(TTIME%TDATE%YEAR,TTIME%TDATE%MONTH,TTIME%TDATE%DAY,TTIME%TIME) 00224 ! 00225 !-------------------------------------------------------------------------------------- 00226 ! Fluxes over water according to Charnock formulae 00227 !-------------------------------------------------------------------------------------- 00228 ! 00229 ! 00230 CALL WATER_FLUX(XZ0, PTA, ZEXNA, PRHOA, XTS, ZEXNS, ZQA, PRAIN, & 00231 PSNOW, XTT, ZWIND, PZREF, PUREF, PPS, ZQSAT, PSFTH, PSFTQ,& 00232 ZUSTAR, ZCD, ZCDN, ZCH, ZRI, ZRESA_WATER, ZZ0H ) 00233 ! 00234 !------------------------------------------------------------------------------------- 00235 ! Outputs: 00236 !------------------------------------------------------------------------------------- 00237 ! 00238 ! Momentum fluxes 00239 ! 00240 IF(CIMPLICIT_WIND=='OLD')THEN 00241 ! old implicitation 00242 ZUSTAR2(:) = (ZCD(:)*ZWIND(:)*PPEW_B_COEF(:)) / & 00243 (1.0-PRHOA(:)*ZCD(:)*ZWIND(:)*PPEW_A_COEF(:)) 00244 ELSE 00245 ! new implicitation 00246 ZUSTAR2(:) = (ZCD(:)*ZWIND(:)*(2.*PPEW_B_COEF(:)-ZWIND(:))) /& 00247 (1.0-2.0*PRHOA(:)*ZCD(:)*ZWIND(:)*PPEW_A_COEF(:)) 00248 ! 00249 ZWIND(:) = PRHOA(:)*PPEW_A_COEF(:)*ZUSTAR2(:) + PPEW_B_COEF(:) 00250 ZWIND(:) = MAX(ZWIND(:),0.) 00251 ! 00252 WHERE(PPEW_A_COEF(:)/= 0.) 00253 ZUSTAR2(:) = MAX( ( ZWIND(:) - PPEW_B_COEF(:) ) / (PRHOA(:)*PPEW_A_COEF(:)), 0.) 00254 ENDWHERE 00255 ! 00256 ENDIF 00257 ! 00258 PSFU = 0. 00259 PSFV = 0. 00260 WHERE (ZWIND(:)>0.) 00261 PSFU(:) = - PRHOA(:) * ZUSTAR2(:) * PU(:) / ZWIND(:) 00262 PSFV(:) = - PRHOA(:) * ZUSTAR2(:) * PV(:) / ZWIND(:) 00263 END WHERE 00264 ! 00265 ! CO2 flux 00266 ! 00267 PSFCO2(:) = 0.0 ! Assumes no CO2 emission over water bodies 00268 ! 00269 !------------------------------------------------------------------------------------- 00270 !radiative properties at t 00271 !------------------------------------------------------------------------------------- 00272 ! 00273 ISWB = SIZE(PSW_BANDS) 00274 ! 00275 DO JSWB=1,ISWB 00276 ZDIR_ALB(:,JSWB) = XDIR_ALB(:) 00277 ZSCA_ALB(:,JSWB) = XSCA_ALB(:) 00278 END DO 00279 ! 00280 ZEMIS = XEMIS 00281 ZTRAD = XTS 00282 ! 00283 !------------------------------------------------------------------------------------- 00284 !Specific fields for GELATO when using earth system model 00285 !(intermediate step before explicit sea and ice fluxes comutation) 00286 !------------------------------------------------------------------------------------- 00287 ! 00288 IF(LCPL_ESM)THEN 00289 CALL COUPLING_ICEFLUX_n(KI, PTA, ZEXNA, PRHOA, XTICE, ZEXNS, & 00290 ZQA, PRAIN, PSNOW, ZWIND, PZREF, PUREF, & 00291 PPS, XTS, XTTS, ZSFTH_ICE, ZSFTQ_ICE ) 00292 ENDIF 00293 ! 00294 !------------------------------------------------------------------------------------- 00295 ! Scalar fluxes: 00296 !------------------------------------------------------------------------------------- 00297 ! 00298 ! 00299 IF (NBEQ>0) THEN 00300 IF (CCH_DRY_DEP == "WES89") THEN 00301 CALL CH_DEP_WATER (ZRESA_WATER, ZUSTAR, PTA, ZTRAD, & 00302 PSV(:,NSV_CHSBEG:NSV_CHSEND), & 00303 CSV(NSV_CHSBEG:NSV_CHSEND), & 00304 XDEP(:,1:NBEQ) ) 00305 00306 PSFTS(:,NSV_CHSBEG:NSV_CHSEND) = - PSV(:,NSV_CHSBEG:NSV_CHSEND) & 00307 * XDEP(:,1:NBEQ) 00308 IF (NAEREQ > 0 ) THEN 00309 CALL CH_AER_DEP(PSV(:,NSV_AERBEG:NSV_AEREND),& 00310 PSFTS(:,NSV_AERBEG:NSV_AEREND),& 00311 ZUSTAR,ZRESA_WATER,PTA,PRHOA) 00312 END IF 00313 00314 ELSE 00315 PSFTS(:,NSV_CHSBEG:NSV_CHSEND) =0. 00316 IF(NSV_AERBEG.LT.NSV_AEREND) PSFTS(:,NSV_AERBEG:NSV_AEREND) =0. 00317 ENDIF 00318 ENDIF 00319 ! 00320 IF (NDSTEQ>0) THEN 00321 CALL DSLT_DEP(PSV(:,NSV_DSTBEG:NSV_DSTEND), PSFTS(:,NSV_DSTBEG:NSV_DSTEND), & 00322 ZUSTAR, ZRESA_WATER, PTA, PRHOA, XEMISSIG_DST, XEMISRADIUS_DST, & 00323 JPMODE_DST, XDENSITY_DST, XMOLARWEIGHT_DST, ZCONVERTFACM0_DST, & 00324 ZCONVERTFACM6_DST, ZCONVERTFACM3_DST, LVARSIG_DST, LRGFIX_DST, & 00325 CVERMOD ) 00326 00327 CALL MASSFLUX2MOMENTFLUX( & 00328 PSFTS(:,NSV_DSTBEG:NSV_DSTEND), & !I/O ![kg/m2/sec] In: flux of only mass, out: flux of moments 00329 PRHOA, & !I [kg/m3] air density 00330 XEMISRADIUS_DST, &!I [um] emitted radius for the modes (max 3) 00331 XEMISSIG_DST, &!I [-] emitted sigma for the different modes (max 3) 00332 NDSTMDE, & 00333 ZCONVERTFACM0_DST, & 00334 ZCONVERTFACM6_DST, & 00335 ZCONVERTFACM3_DST, & 00336 LVARSIG_DST, LRGFIX_DST ) 00337 ENDIF 00338 ! 00339 IF (NSLTEQ>0) THEN 00340 CALL DSLT_DEP(PSV(:,NSV_SLTBEG:NSV_SLTEND), PSFTS(:,NSV_SLTBEG:NSV_SLTEND), & 00341 ZUSTAR, ZRESA_WATER, PTA, PRHOA, XEMISSIG_SLT, XEMISRADIUS_SLT, & 00342 JPMODE_SLT, XDENSITY_SLT, XMOLARWEIGHT_SLT, ZCONVERTFACM0_SLT, & 00343 ZCONVERTFACM6_SLT, ZCONVERTFACM3_SLT, LVARSIG_SLT, LRGFIX_SLT, & 00344 CVERMOD ) 00345 00346 CALL MASSFLUX2MOMENTFLUX( & 00347 PSFTS(:,NSV_SLTBEG:NSV_SLTEND), & !I/O ![kg/m2/sec] In: flux of only mass, out: flux of moments 00348 PRHOA, & !I [kg/m3] air density 00349 XEMISRADIUS_SLT, &!I [um] emitted radius for the modes (max 3) 00350 XEMISSIG_SLT, &!I [-] emitted sigma for the different modes (max 3) 00351 NSLTMDE, & 00352 ZCONVERTFACM0_SLT, & 00353 ZCONVERTFACM6_SLT, & 00354 ZCONVERTFACM3_SLT, & 00355 LVARSIG_SLT, LRGFIX_SLT ) 00356 ENDIF 00357 ! 00358 !------------------------------------------------------------------------------------- 00359 ! Inline diagnostics 00360 !------------------------------------------------------------------------------------- 00361 ! 00362 CALL DIAG_INLINE_WATFLUX_n(PTSTEP,PTA, XTS, ZQA, PPA, PPS, PRHOA, PU, PV, PZREF, & 00363 PUREF, ZCD, ZCDN, ZCH, ZRI, ZHU, XZ0, ZZ0H, ZQSAT, & 00364 PSFTH, PSFTQ, PSFU, PSFV, PDIR_SW, PSCA_SW, PLW, & 00365 ZDIR_ALB, ZSCA_ALB, ZEMIS, ZTRAD, PRAIN, PSNOW, & 00366 XTICE, ZSFTH_ICE, ZSFTQ_ICE ) 00367 ! 00368 !------------------------------------------------------------------------------------- 00369 !Radiative properties at time t+1 (see by the atmosphere) in order to close 00370 !the energy budget between surfex and the atmosphere 00371 !------------------------------------------------------------------------------------- 00372 ! 00373 IF (LINTERPOL_TS.AND.MOD(TTIME%TIME,XDAY) == 0.) THEN 00374 CALL INTERPOL_TS_WATER_MTH(TTIME%TDATE%YEAR,TTIME%TDATE%MONTH,TTIME%TDATE%DAY,XTS) 00375 ENDIF 00376 ! 00377 CALL UPDATE_RAD_SEAWAT(CWAT_ALB,XTS,PZENITH2,XTT,XEMIS,XDIR_ALB, & 00378 XSCA_ALB,PDIR_ALB,PSCA_ALB,PEMIS,PTRAD ) 00379 ! 00380 !======================================================================================= 00381 ! 00382 IF (LHOOK) CALL DR_HOOK('COUPLING_WATFLUX_N',1,ZHOOK_HANDLE) 00383 ! 00384 END SUBROUTINE COUPLING_WATFLUX_n
1.8.0
