- Timestamp:
- Jul 24, 2024, 4:23:34 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_fluxs_mod.F90
r5116 r5117 18 18 USE dimphy, ONLY: klon 19 19 USE indice_sol_mod 20 usesens_heat_rain_m, ONLY: sens_heat_rain20 USE sens_heat_rain_m, ONLY: sens_heat_rain 21 21 22 22 INCLUDE "clesphys.h" … … 77 77 REAL, INTENT(IN) :: fqsat ! correction factor on qsat (generally 0.98 over salty water, 1 everywhere else) 78 78 79 real, intent(in), optional:: rhoa(:) ! (knon)79 REAL, INTENT(IN), optional:: rhoa(:) ! (knon) 80 80 ! density of moist air (kg / m3) 81 81 … … 89 89 REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new, evap, fluxsens, fluxlat 90 90 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l 91 REAL, intent(out), OPTIONAL:: sens_prec_liq(:), sens_prec_sol(:) ! (knon)91 REAL, INTENT(OUT), OPTIONAL:: sens_prec_liq(:), sens_prec_sol(:) ! (knon) 92 92 REAL, DIMENSION(klon), OPTIONAL :: lat_prec_liq, lat_prec_sol 93 93 … … 142 142 dflux_s = 0. 143 143 dflux_l = 0. 144 if(PRESENT(lat_prec_liq)) lat_prec_liq = 0.145 if(PRESENT(lat_prec_sol)) lat_prec_sol = 0.144 IF (PRESENT(lat_prec_liq)) lat_prec_liq = 0. 145 IF (PRESENT(lat_prec_sol)) lat_prec_sol = 0. 146 146 147 147 ! zx_qs = qsat en kg/kg … … 258 258 259 259 ! calcul de l'enthalpie des precipitations liquides et solides 260 if(PRESENT(sens_prec_liq)) sens_prec_liq(i) &260 IF (PRESENT(sens_prec_liq)) sens_prec_liq(i) & 261 261 = - sens_heat_rain(precip_rain(i) + precip_snow(i), t1lay(i), & 262 262 q1lay(i), rhoa(i), rlvtt, tsurf_new(i), ps(i)) 263 if(PRESENT(sens_prec_sol)) sens_prec_sol(i) = 0.263 IF (PRESENT(sens_prec_sol)) sens_prec_sol(i) = 0. 264 264 ! On calcule par rapport a T=0 265 265 !! sens_prec_liq(i) = rcw * (t1lay(i) - RTT) * precip_rain(i) 266 266 !! sens_prec_sol(i) = rcs * (t1lay(i) - RTT) * precip_snow(i) 267 267 268 if(PRESENT(lat_prec_liq)) &268 IF (PRESENT(lat_prec_liq)) & 269 269 lat_prec_liq(i) = precip_rain(i) * (RLVTT - RLVTT) 270 if(PRESENT(lat_prec_sol)) &270 IF (PRESENT(lat_prec_sol)) & 271 271 lat_prec_sol(i) = precip_snow(i) * (RLSTT - RLVTT) 272 272 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.