/[lmdze]/trunk/Sources/phylmd/Interface_surf/calcul_fluxs.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Interface_surf/calcul_fluxs.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC revision 206 by guez, Tue Aug 30 12:52:46 2016 UTC
# Line 6  contains Line 6  contains
6    
7    SUBROUTINE calcul_fluxs(dtime, tsurf, p1lay, cal, beta, coef1lay, ps, &    SUBROUTINE calcul_fluxs(dtime, tsurf, p1lay, cal, beta, coef1lay, ps, &
8         qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, petAcoef, &         qsurf, radsol, dif_grnd, t1lay, q1lay, u1lay, v1lay, petAcoef, &
9         peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, fluxsens, &         peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, flux_t, &
10         dflux_s, dflux_l)         dflux_s, dflux_l)
11    
12      ! Cette routine calcule les fluxs en h et q à l'interface et une      ! Cette routine calcule les fluxs en h et q à l'interface et une
# Line 20  contains Line 20  contains
20      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2
21    
22      real, intent(IN):: dtime      real, intent(IN):: dtime
23      real, intent(IN):: tsurf(:) ! (knon) temperature de surface      real, intent(IN):: tsurf(:) ! (knon) température de surface
24      real, intent(IN):: p1lay(:) ! (knon) pression 1er niveau (milieu de couche)  
25        real, intent(IN):: p1lay(:) ! (knon)
26        ! pression première couche (milieu de couche)
27    
28      real, intent(IN):: cal(:) ! (knon) capacité calorifique du sol      real, intent(IN):: cal(:) ! (knon) capacité calorifique du sol
29      real, intent(IN):: beta(:) ! (knon) evap reelle      real, intent(IN):: beta(:) ! (knon) évaporation réelle
30      real, intent(IN):: coef1lay(:) ! (knon) coefficient d'échange      real, intent(IN):: coef1lay(:) ! (knon) coefficient d'échange
31      real, intent(IN):: ps(:) ! (knon) pression au sol      real, intent(IN):: ps(:) ! (knon) pression au sol
32      real, intent(OUT):: qsurf(:) ! (knon) humidite de l'air au dessus du sol      real, intent(OUT):: qsurf(:) ! (knon) humidité de l'air au-dessus du sol
33      real, intent(IN):: radsol(:) ! (knon) rayonnement net au sol (LW + SW)  
34        real, intent(IN):: radsol(:) ! (knon)
35        ! rayonnement net au sol (longwave + shortwave)
36    
37      real, intent(IN):: dif_grnd(:) ! (knon)      real, intent(IN):: dif_grnd(:) ! (knon)
38      ! coefficient diffusion vers le sol profond      ! coefficient de diffusion vers le sol profond
39    
40      real, intent(IN):: t1lay(:), q1lay(:), u1lay(:), v1lay(:) ! (knon)      real, intent(IN):: t1lay(:), q1lay(:), u1lay(:), v1lay(:) ! (knon)
41    
42      real, intent(IN):: petAcoef(:), peqAcoef(:) ! (knon)      real, intent(IN):: petAcoef(:), peqAcoef(:) ! (knon)
43      ! coefficients A de la résolution de la couche limite pour t et q      ! coefficients A de la résolution de la couche limite pour T et q
44    
45      real, intent(IN):: petBcoef(:), peqBcoef(:) ! (knon)      real, intent(IN):: petBcoef(:), peqBcoef(:) ! (knon)
46      ! coeff. B de la resolution de la CL pour t et q      ! coefficients B de la résolution de la couche limite pour t et q
47    
48      real, intent(OUT):: tsurf_new(:) ! (knon) température au sol      real, intent(OUT):: tsurf_new(:) ! (knon) température au sol
49      real, intent(OUT):: evap(:) ! (knon)      real, intent(OUT):: evap(:) ! (knon)
50    
51      real, intent(OUT):: fluxlat(:), fluxsens(:) ! (knon)      real, intent(OUT):: fluxlat(:), flux_t(:) ! (knon)
52      ! flux de chaleur latente et sensible      ! flux de chaleur latente et sensible
53    
54      real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon)      real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon)
# Line 53  contains Line 58  contains
58      ! Local:      ! Local:
59      integer i      integer i
60      integer knon ! nombre de points a traiter      integer knon ! nombre de points a traiter
61      real, dimension(size(ps)):: mh, oh, mq, nq, oq      real, dimension(size(ps)):: mh, oh, mq, nq, oq, dq_s_dt, coef ! (knon)
62      real, dimension(size(ps)):: dq_s_dt, coef      real qsat(size(ps)) ! mass fraction
63      real qsat(size(ps)) ! qsat en kg/kg      real sl(size(ps)) ! chaleur latente d'évaporation ou de sublimation
     real sl(size(ps)) ! chaleur latente d'evaporation ou de sublimation  
64      logical delta      logical delta
65      real zcor      real zcor
66      real, parameter:: t_grnd = 271.35, t_coup = 273.15      real, parameter:: t_grnd = 271.35, t_coup = 273.15
# Line 67  contains Line 71  contains
71           size(coef1lay), size(ps), size(qsurf), size(radsol), size(dif_grnd), &           size(coef1lay), size(ps), size(qsurf), size(radsol), size(dif_grnd), &
72           size(t1lay), size(q1lay), size(u1lay), size(v1lay), size(petAcoef), &           size(t1lay), size(q1lay), size(u1lay), size(v1lay), size(petAcoef), &
73           size(peqAcoef), size(petBcoef), size(peqBcoef), size(tsurf_new), &           size(peqAcoef), size(petBcoef), size(peqBcoef), size(tsurf_new), &
74           size(evap), size(fluxlat), size(fluxsens), size(dflux_s), &           size(evap), size(fluxlat), size(flux_t), size(dflux_s), &
75           size(dflux_l)/), "calcul_fluxs knon")           size(dflux_l)/), "calcul_fluxs knon")
76    
77      ! Traitement humidite du sol      ! Traitement de l'humidité du sol
78    
79      IF (thermcep) THEN      IF (thermcep) THEN
80         DO i = 1, knon         DO i = 1, knon
# Line 114  contains Line 118  contains
118    
119      evap = - mq - nq * tsurf_new      evap = - mq - nq * tsurf_new
120      fluxlat = - evap * sl      fluxlat = - evap * sl
121      fluxsens = mh + dflux_s * tsurf_new      flux_t = mh + dflux_s * tsurf_new
122      dflux_l = sl * nq      dflux_l = sl * nq
123    
124      ! Nouvelle valeur de l'humidité au dessus du sol :      ! Nouvelle valeur de l'humidité au dessus du sol :
Legend:
Removed from v.178  
changed lines
  Added in v.206
  ViewVC Help
Powered by ViewVC 1.1.21