LMDZ
calcul_fluxs_mod.F90
Go to the documentation of this file.
1 !
2 MODULE calcul_fluxs_mod
3 
4 
5 CONTAINS
6  SUBROUTINE calcul_fluxs( knon, nisurf, dtime, &
7  tsurf, p1lay, cal, beta, cdragh, cdragq, ps, &
8  precip_rain, precip_snow, snow, qsurf, &
9  radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, gustiness, &
10  fqsat, petacoef, peqacoef, petbcoef, peqbcoef, &
11  tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
12 
13  USE dimphy, ONLY : klon
14  USE indice_sol_mod
15 
16  include "clesphys.h"
17 
18 ! Cette routine calcule les fluxs en h et q a l'interface et eventuellement
19 ! une temperature de surface (au cas ou ok_veget = false)
20 !
21 ! L. Fairhead 4/2000
22 !
23 ! input:
24 ! knon nombre de points a traiter
25 ! nisurf surface a traiter
26 ! tsurf temperature de surface
27 ! p1lay pression 1er niveau (milieu de couche)
28 ! cal capacite calorifique du sol
29 ! beta evap reelle
30 ! cdragh coefficient d'echange temperature
31 ! cdragq coefficient d'echange evaporation
32 ! ps pression au sol
33 ! precip_rain precipitations liquides
34 ! precip_snow precipitations solides
35 ! snow champs hauteur de neige
36 ! runoff runoff en cas de trop plein
37 ! petAcoef coeff. A de la resolution de la CL pour t
38 ! peqAcoef coeff. A de la resolution de la CL pour q
39 ! petBcoef coeff. B de la resolution de la CL pour t
40 ! peqBcoef coeff. B de la resolution de la CL pour q
41 ! radsol rayonnement net aus sol (LW + SW)
42 ! dif_grnd coeff. diffusion vers le sol profond
43 !
44 ! output:
45 ! tsurf_new temperature au sol
46 ! qsurf humidite de l'air au dessus du sol
47 ! fluxsens flux de chaleur sensible
48 ! fluxlat flux de chaleur latente
49 ! dflux_s derivee du flux de chaleur sensible / Ts
50 ! dflux_l derivee du flux de chaleur latente / Ts
51 !
52 
53  include "YOETHF.h"
54  include "FCTTRE.h"
55  include "YOMCST.h"
56 
57 ! Parametres d'entree
58 !****************************************************************************************
59  INTEGER, INTENT(IN) :: knon, nisurf
60  REAL , INTENT(IN) :: dtime
61  REAL, DIMENSION(klon), INTENT(IN) :: petAcoef, peqAcoef
62  REAL, DIMENSION(klon), INTENT(IN) :: petBcoef, peqBcoef
63  REAL, DIMENSION(klon), INTENT(IN) :: ps, q1lay
64  REAL, DIMENSION(klon), INTENT(IN) :: tsurf, p1lay, cal, beta, cdragh,cdragq
65  REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow ! pas utiles
66  REAL, DIMENSION(klon), INTENT(IN) :: radsol, dif_grnd
67  REAL, DIMENSION(klon), INTENT(IN) :: t1lay, u1lay, v1lay,gustiness
68  REAL, INTENT(IN) :: fqsat ! correction factor on qsat (generally 0.98 over salty water, 1 everywhere else)
69 
70 ! Parametres entree-sorties
71 !****************************************************************************************
72  REAL, DIMENSION(klon), INTENT(INOUT) :: snow ! snow pas utile
73 
74 ! Parametres sorties
75 !****************************************************************************************
76  REAL, DIMENSION(klon), INTENT(OUT) :: qsurf
77  REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new, evap, fluxsens, fluxlat
78  REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l
79 
80 ! Variables locales
81 !****************************************************************************************
82  INTEGER :: i
83  REAL, DIMENSION(klon) :: zx_mh, zx_nh, zx_oh
84  REAL, DIMENSION(klon) :: zx_mq, zx_nq, zx_oq
85  REAL, DIMENSION(klon) :: zx_pkh, zx_dq_s_dt, zx_qsat
86  REAL, DIMENSION(klon) :: zx_sl, zx_coefh, zx_coefq, zx_wind
87  REAL, DIMENSION(klon) :: d_ts
88  REAL :: zdelta, zcvm5, zx_qs, zcor, zx_dq_s_dh
89  REAL :: qsat_new, q1_new
90  REAL, PARAMETER :: t_grnd = 271.35, t_coup = 273.15
91  REAL, PARAMETER :: max_eau_sol = 150.0
92  CHARACTER (len = 20) :: modname = 'calcul_fluxs'
93  LOGICAL :: fonte_neige
94  LOGICAL, SAVE :: check = .false.
95  !$OMP THREADPRIVATE(check)
96 
97 ! End definition
98 !****************************************************************************************
99 
100  IF (check) WRITE(*,*)'Entree ', modname,' surface = ',nisurf
101 
102  IF (check) THEN
103  WRITE(*,*)' radsol (min, max)', &
104  minval(radsol(1:knon)), maxval(radsol(1:knon))
105  ENDIF
106 
107 ! Traitement neige et humidite du sol
108 !****************************************************************************************
109 !
110 !!$ WRITE(*,*)'test calcul_flux, surface ', nisurf
111 !!PB test
112 !!$ if (nisurf == is_oce) then
113 !!$ snow = 0.
114 !!$ qsol = max_eau_sol
115 !!$ else
116 !!$ where (precip_snow > 0.) snow = snow + (precip_snow * dtime)
117 !!$ where (snow > epsilon(snow)) snow = max(0.0, snow - (evap * dtime))
118 !!$! snow = max(0.0, snow + (precip_snow - evap) * dtime)
119 !!$ where (precip_rain > 0.) qsol = qsol + (precip_rain - evap) * dtime
120 !!$ endif
121 !!$ IF (nisurf /= is_ter) qsol = max_eau_sol
122 
123 
124 !
125 ! Initialisation
126 !****************************************************************************************
127  evap = 0.
128  fluxsens=0.
129  fluxlat=0.
130  dflux_s = 0.
131  dflux_l = 0.
132 !
133 ! zx_qs = qsat en kg/kg
134 !****************************************************************************************
135  DO i = 1, knon
136  zx_pkh(i) = (ps(i)/ps(i))**rkappa
137  IF (thermcep) THEN
138  zdelta=max(0.,sign(1.,rtt-tsurf(i)))
139  zcvm5 = r5les*rlvtt*(1.-zdelta) + r5ies*rlstt*zdelta
140  zcvm5 = zcvm5 / rcpd / (1.0+rvtmp2*q1lay(i))
141  zx_qs= r2es * foeew(tsurf(i),zdelta)/ps(i)
142  zx_qs=min(0.5,zx_qs)
143  zcor=1./(1.-retv*zx_qs)
144  zx_qs=zx_qs*zcor
145  zx_dq_s_dh = foede(tsurf(i),zdelta,zcvm5,zx_qs,zcor) &
146  /rlvtt / zx_pkh(i)
147  ELSE
148  IF (tsurf(i).LT.t_coup) THEN
149  zx_qs = qsats(tsurf(i)) / ps(i)
150  zx_dq_s_dh = dqsats(tsurf(i),zx_qs)/rlvtt &
151  / zx_pkh(i)
152  ELSE
153  zx_qs = qsatl(tsurf(i)) / ps(i)
154  zx_dq_s_dh = dqsatl(tsurf(i),zx_qs)/rlvtt &
155  / zx_pkh(i)
156  ENDIF
157  ENDIF
158  zx_dq_s_dt(i) = rcpd * zx_pkh(i) * zx_dq_s_dh
159  zx_qsat(i) = zx_qs
160  zx_wind(i)=min_wind_speed+sqrt(gustiness(i)+u1lay(i)**2+v1lay(i)**2)
161  zx_coefh(i) = cdragh(i) * zx_wind(i) * p1lay(i)/(rd*t1lay(i))
162  zx_coefq(i) = cdragq(i) * zx_wind(i) * p1lay(i)/(rd*t1lay(i))
163 ! zx_wind(i)=min_wind_speed+SQRT(gustiness(i)+u1lay(i)**2+v1lay(i)**2) &
164 ! * p1lay(i)/(RD*t1lay(i))
165 ! zx_coefh(i) = cdragh(i) * zx_wind(i)
166 ! zx_coefq(i) = cdragq(i) * zx_wind(i)
167  ENDDO
168 
169 
170 ! === Calcul de la temperature de surface ===
171 ! zx_sl = chaleur latente d'evaporation ou de sublimation
172 !****************************************************************************************
173 
174  DO i = 1, knon
175  zx_sl(i) = rlvtt
176  IF (tsurf(i) .LT. rtt) zx_sl(i) = rlstt
177  ENDDO
178 
179 
180  DO i = 1, knon
181 ! Q
182  zx_oq(i) = 1. - (beta(i) * zx_coefq(i) * peqbcoef(i) * dtime)
183  zx_mq(i) = beta(i) * zx_coefq(i) * &
184  (peqacoef(i) - &
185 ! conv num avec precedente version
186  fqsat * zx_qsat(i) + fqsat * zx_dq_s_dt(i) * tsurf(i)) &
187 ! fqsat * ( zx_qsat(i) - zx_dq_s_dt(i) * tsurf(i)) ) &
188  / zx_oq(i)
189  zx_nq(i) = beta(i) * zx_coefq(i) * (- fqsat * zx_dq_s_dt(i)) &
190  / zx_oq(i)
191 
192 ! H
193  zx_oh(i) = 1. - (zx_coefh(i) * petbcoef(i) * dtime)
194  zx_mh(i) = zx_coefh(i) * petacoef(i) / zx_oh(i)
195  zx_nh(i) = - (zx_coefh(i) * rcpd * zx_pkh(i))/ zx_oh(i)
196 
197 ! Tsurface
198  tsurf_new(i) = (tsurf(i) + cal(i)/(rcpd * zx_pkh(i)) * dtime * &
199  (radsol(i) + zx_mh(i) + zx_sl(i) * zx_mq(i)) &
200  + dif_grnd(i) * t_grnd * dtime)/ &
201  ( 1. - dtime * cal(i)/(rcpd * zx_pkh(i)) * ( &
202  zx_nh(i) + zx_sl(i) * zx_nq(i)) &
203  + dtime * dif_grnd(i))
204 
205 !
206 ! Y'a-t-il fonte de neige?
207 !
208 ! fonte_neige = (nisurf /= is_oce) .AND. &
209 ! & (snow(i) > epsfra .OR. nisurf == is_sic .OR. nisurf == is_lic) &
210 ! & .AND. (tsurf_new(i) >= RTT)
211 ! if (fonte_neige) tsurf_new(i) = RTT
212  d_ts(i) = tsurf_new(i) - tsurf(i)
213 ! zx_h_ts(i) = tsurf_new(i) * RCPD * zx_pkh(i)
214 ! zx_q_0(i) = zx_qsat(i) + zx_dq_s_dt(i) * d_ts(i)
215 
216 !== flux_q est le flux de vapeur d'eau: kg/(m**2 s) positive vers bas
217 !== flux_t est le flux de cpt (energie sensible): j/(m**2 s)
218  evap(i) = - zx_mq(i) - zx_nq(i) * tsurf_new(i)
219  fluxlat(i) = - evap(i) * zx_sl(i)
220  fluxsens(i) = zx_mh(i) + zx_nh(i) * tsurf_new(i)
221 
222 ! Derives des flux dF/dTs (W m-2 K-1):
223  dflux_s(i) = zx_nh(i)
224  dflux_l(i) = (zx_sl(i) * zx_nq(i))
225 
226 ! Nouvelle valeure de l'humidite au dessus du sol
227  qsat_new=zx_qsat(i) + zx_dq_s_dt(i) * d_ts(i)
228  q1_new = peqacoef(i) - peqbcoef(i)*evap(i)*dtime
229  qsurf(i)=q1_new*(1.-beta(i)) + beta(i)*qsat_new
230 !
231 ! en cas de fonte de neige
232 !
233 ! if (fonte_neige) then
234 ! bilan_f = radsol(i) + fluxsens(i) - (zx_sl(i) * evap (i)) - &
235 ! & dif_grnd(i) * (tsurf_new(i) - t_grnd) - &
236 ! & RCPD * (zx_pkh(i))/cal(i)/dtime * (tsurf_new(i) - tsurf(i))
237 ! bilan_f = max(0., bilan_f)
238 ! fq_fonte = bilan_f / zx_sl(i)
239 ! snow(i) = max(0., snow(i) - fq_fonte * dtime)
240 ! qsol(i) = qsol(i) + (fq_fonte * dtime)
241 ! endif
242 !!$ if (nisurf == is_ter) &
243 !!$ & run_off(i) = run_off(i) + max(qsol(i) - max_eau_sol, 0.0)
244 !!$ qsol(i) = min(qsol(i), max_eau_sol)
245  ENDDO
246 !
247 !****************************************************************************************
248 !
249  END SUBROUTINE calcul_fluxs
250 !
251 !****************************************************************************************
252 !
253  SUBROUTINE calcul_flux_wind(knon, dtime, &
254  u0, v0, u1, v1, gustiness, cdrag_m, &
255  acoefu, acoefv, bcoefu, bcoefv, &
256  p1lay, t1lay, &
257  flux_u1, flux_v1)
258 
259  USE dimphy
260  include "YOMCST.h"
261  include "clesphys.h"
262 
263 ! Input arguments
264 !****************************************************************************************
265  INTEGER, INTENT(IN) :: knon
266  REAL, INTENT(IN) :: dtime
267  REAL, DIMENSION(klon), INTENT(IN) :: u0, v0 ! u and v at niveau 0
268  REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness ! u and v at niveau 1
269  REAL, DIMENSION(klon), INTENT(IN) :: cdrag_m ! cdrag pour momentum
270  REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV
271  REAL, DIMENSION(klon), INTENT(IN) :: p1lay ! pression 1er niveau (milieu de couche)
272  REAL, DIMENSION(klon), INTENT(IN) :: t1lay ! temperature
273 ! Output arguments
274 !****************************************************************************************
275  REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1
276  REAL, DIMENSION(klon), INTENT(OUT) :: flux_v1
277 
278 ! Local variables
279 !****************************************************************************************
280  INTEGER :: i
281  REAL :: mod_wind, buf
282 
283 !****************************************************************************************
284 ! Calculate the surface flux
285 !
286 !****************************************************************************************
287  DO i=1,knon
288  mod_wind = min_wind_speed + sqrt(gustiness(i)+(u1(i) - u0(i))**2 + (v1(i)-v0(i))**2)
289  buf = cdrag_m(i) * mod_wind * p1lay(i)/(rd*t1lay(i))
290  flux_u1(i) = (acoefu(i) - u0(i)) / (1/buf - bcoefu(i)*dtime )
291  flux_v1(i) = (acoefv(i) - v0(i)) / (1/buf - bcoefv(i)*dtime )
292  END DO
293 
294  END SUBROUTINE calcul_flux_wind
295 !
296 !****************************************************************************************
297 !
298 END MODULE calcul_fluxs_mod
calcul_fluxs_mod::calcul_flux_wind
subroutine calcul_flux_wind(knon, dtime, u0, v0, u1, v1, gustiness, cdrag_m, AcoefU, AcoefV, BcoefU, BcoefV, p1lay, t1lay, flux_u1, flux_v1)
Definition: calcul_fluxs_mod.F90:258
dimphy::klon
integer, save klon
Definition: dimphy.F90:3
beta
!$Id!Parameters for nlm real spfac!IM cf epmax real ptcrit real omtrain real dttrig real beta
Definition: cv30param.h:5
min_wind_speed
!$Id ok_orolf LOGICAL ok_limitvrai LOGICAL ok_all_xml INTEGER iflag_ener_conserv REAL solaire RCFC12 RCFC12_act CFC12_ppt!IM ajout CFMIP2 CMIP5 LOGICAL ok_4xCO2atm RCFC12_per CFC12_ppt_per!OM correction du bilan d eau global!OM Correction sur precip KE REAL cvl_corr!OM Fonte calotte dans bilan eau LOGICAL ok_lic_melt!IM simulateur ISCCP INTEGER overlap!IM seuils cdrh REAL cdhmax!IM param stabilite s terres et en dehors REAL f_ri_cd_min!IM MAFo pmagic evap0!Frottement au f_cdrag_oce REAL min_wind_speed
Definition: clesphys.h:46
false
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL false
Definition: calcul_STDlev.h:26
calcul_fluxs_mod
Definition: calcul_fluxs_mod.F90:2
calcul_fluxs_mod::calcul_fluxs
subroutine calcul_fluxs(knon, nisurf, dtime, tsurf, p1lay, cal, beta, cdragh, cdragq, ps, precip_rain, precip_snow, snow, qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, gustiness, fqsat, petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
Definition: calcul_fluxs_mod.F90:12
dimphy
Definition: dimphy.F90:1
indice_sol_mod
Definition: indice_sol_mod.F90:1