/[lmdze]/trunk/phylmd/Interface_surf/pbl_surface.f90
ViewVC logotype

Diff of /trunk/phylmd/Interface_surf/pbl_surface.f90

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

revision 332 by guez, Tue Aug 13 09:19:22 2019 UTC revision 341 by guez, Mon Oct 21 06:11:44 2019 UTC
# Line 4  module pbl_surface_m Line 4  module pbl_surface_m
4    
5  contains  contains
6    
7    SUBROUTINE pbl_surface(pctsrf, t, q, u, v, julien, mu0, ftsol, cdmmax, &    SUBROUTINE pbl_surface(pctsrf, t_seri, q_seri, u, v, julien, mu0, ftsol, &
8         cdhmax, ftsoil, qsol, paprs, play, fsnow, fqsurf, falbe, fluxlat, &         cdmmax, cdhmax, ftsoil, qsol, paprs, play, fsnow, fqsurf, falbe, &
9         rain_fall, snow_fall, frugs, agesno, rugoro, d_t, d_q, d_u, d_v, &         fluxlat, rain_fall, snow_fall, frugs, agesno, rugoro, d_t, d_q, d_u, &
10         flux_t, flux_q, flux_u, flux_v, cdragh, cdragm, q2, dflux_t, dflux_q, &         d_v, flux_t, flux_q, flux_u, flux_v, cdragh, cdragm, q2, dflux_t, &
11         coefh, t2m, q2m, u10m_srf, v10m_srf, pblh, capcl, oliqcl, cteicl, pblt, &         dflux_q, coefh, t2m, q2m, u10m_srf, v10m_srf, pblh, capcl, oliqcl, &
12         therm, plcl, fqcalving, ffonte, run_off_lic_0, albsol, sollw, solsw, &         cteicl, pblt, therm, plcl, fqcalving, ffonte, run_off_lic_0, albsol, &
13         tsol)         sollw, solsw, tsol)
14    
15      ! From phylmd/clmain.F, version 1.6, 2005/11/16 14:47:19      ! From phylmd/clmain.F, version 1.6, 2005/11/16 14:47:19
16      ! Author: Z. X. Li (LMD/CNRS)      ! Author: Z. X. Li (LMD/CNRS)
# Line 42  contains Line 42  contains
42      REAL, INTENT(inout):: pctsrf(:, :) ! (klon, nbsrf)      REAL, INTENT(inout):: pctsrf(:, :) ! (klon, nbsrf)
43      ! pourcentages de surface de chaque maille      ! pourcentages de surface de chaque maille
44    
45      REAL, INTENT(IN):: t(klon, klev) ! temperature (K)      REAL, INTENT(IN):: t_seri(:, :) ! (klon, klev) air temperature, in K
46      REAL, INTENT(IN):: q(klon, klev) ! vapeur d'eau (kg / kg)      REAL, INTENT(IN):: q_seri(:, :) ! (klon, klev) mass fraction of water vapor
47      REAL, INTENT(IN):: u(klon, klev), v(klon, klev) ! vitesse      REAL, INTENT(IN):: u(klon, klev), v(klon, klev) ! vitesse
48      INTEGER, INTENT(IN):: julien ! jour de l'annee en cours      INTEGER, INTENT(IN):: julien ! jour de l'annee en cours
49      REAL, intent(in):: mu0(klon) ! cosinus de l'angle solaire zenithal          REAL, intent(in):: mu0(klon) ! cosinus de l'angle solaire zenithal    
# Line 82  contains Line 82  contains
82      REAL, INTENT(IN):: rugoro(klon)      REAL, INTENT(IN):: rugoro(klon)
83    
84      REAL, intent(out):: d_t(:, :), d_q(:, :) ! (klon, klev)      REAL, intent(out):: d_t(:, :), d_q(:, :) ! (klon, klev)
85      ! changement pour t et q      ! changement pour t_seri et q_seri
86    
87      REAL, intent(out):: d_u(klon, klev), d_v(klon, klev)      REAL, intent(out):: d_u(klon, klev), d_v(klon, klev)
88      ! changement pour "u" et "v"      ! changement pour "u" et "v"
# Line 121  contains Line 121  contains
121      REAL capcl(klon, nbsrf)      REAL capcl(klon, nbsrf)
122      REAL oliqcl(klon, nbsrf)      REAL oliqcl(klon, nbsrf)
123      REAL cteicl(klon, nbsrf)      REAL cteicl(klon, nbsrf)
124      REAL, INTENT(inout):: pblt(klon, nbsrf) ! T au nveau HCL      REAL, INTENT(inout):: pblt(klon, nbsrf) ! temp\'erature au nveau HCL
125      REAL therm(klon, nbsrf)      REAL therm(klon, nbsrf)
126      REAL plcl(klon, nbsrf)      REAL plcl(klon, nbsrf)
127    
# Line 159  contains Line 159  contains
159      REAL rugmer(klon)      REAL rugmer(klon)
160      REAL ytsoil(klon, nsoilmx)      REAL ytsoil(klon, nsoilmx)
161      REAL yts(klon), ypctsrf(klon), yz0_new(klon)      REAL yts(klon), ypctsrf(klon), yz0_new(klon)
162      real yrugos(klon) ! longueur de rugosite (en m)      real yrugos(klon) ! longueur de rugosit\'e, en m
163      REAL yalb(klon)      REAL yalb(klon)
164      REAL snow(klon) ! column-density of mass of snow at the surface, in kg m-2      REAL snow(klon) ! column-density of mass of snow at the surface, in kg m-2
165      real yqsurf(klon), yagesno(klon)      real yqsurf(klon), yagesno(klon)
# Line 200  contains Line 200  contains
200      REAL ypblt(klon)      REAL ypblt(klon)
201      REAL ytherm(klon)      REAL ytherm(klon)
202      REAL u1(klon), v1(klon)      REAL u1(klon), v1(klon)
203      REAL tair1(klon), qair1(klon), tairsol(klon)      REAL tair1(klon)
     REAL psfce(klon), patm(klon)  
     REAL zgeo1(klon)  
204      REAL rugo1(klon)      REAL rugo1(klon)
205      REAL zgeop(klon, klev)      REAL zgeop(klon, klev)
206    
# Line 233  contains Line 231  contains
231      cdragm = 0.      cdragm = 0.
232      dflux_t = 0.      dflux_t = 0.
233      dflux_q = 0.      dflux_q = 0.
     yrugos = 0.  
234      ypaprs = 0.      ypaprs = 0.
235      ypplay = 0.      ypplay = 0.
236      ydelp = 0.      ydelp = 0.
# Line 311  contains Line 308  contains
308                  ydelp(j, k) = delp(i, k)                  ydelp(j, k) = delp(i, k)
309                  yu(j, k) = u(i, k)                  yu(j, k) = u(i, k)
310                  yv(j, k) = v(i, k)                  yv(j, k) = v(i, k)
311                  yt(j, k) = t(i, k)                  yt(j, k) = t_seri(i, k)
312                  yq(j, k) = q(i, k)                  yq(j, k) = q_seri(i, k)
313               END DO               END DO
314            END DO            END DO
315    
# Line 337  contains Line 334  contains
334               ycdragh(:knon) = max(ycdragh(:knon), 0.)               ycdragh(:knon) = max(ycdragh(:knon), 0.)
335            end IF            end IF
336    
           ! on met un seuil pour ycdragm et ycdragh  
337            IF (nsrf == is_oce) THEN            IF (nsrf == is_oce) THEN
338                 ! On met un seuil pour ycdragm et ycdragh :
339               ycdragm(:knon) = min(ycdragm(:knon), cdmmax)               ycdragm(:knon) = min(ycdragm(:knon), cdmmax)
340               ycdragh(:knon) = min(ycdragh(:knon), cdhmax)               ycdragh(:knon) = min(ycdragh(:knon), cdhmax)
341            END IF            END IF
# Line 385  contains Line 382  contains
382    
383            DO k = 1, klev            DO k = 1, klev
384               DO j = 1, knon               DO j = 1, knon
                 i = ni(j)  
385                  y_d_t(j, k) = y_d_t(j, k) * ypctsrf(j)                  y_d_t(j, k) = y_d_t(j, k) * ypctsrf(j)
386                  y_d_q(j, k) = y_d_q(j, k) * ypctsrf(j)                  y_d_q(j, k) = y_d_q(j, k) * ypctsrf(j)
387                  y_d_u(j, k) = y_d_u(j, k) * ypctsrf(j)                  y_d_u(j, k) = y_d_u(j, k) * ypctsrf(j)
# Line 448  contains Line 444  contains
444            forall (k = 2:klev) coefh(ni(:knon), k) &            forall (k = 2:klev) coefh(ni(:knon), k) &
445                 = coefh(ni(:knon), k) + ycoefh(:knon, k) * ypctsrf(:knon)                 = coefh(ni(:knon), k) + ycoefh(:knon, k) * ypctsrf(:knon)
446    
447            ! diagnostic t, q a 2m et u, v a 10m            ! Diagnostic temp\'erature, q \`a 2 m et u, v \`a 10 m:
448    
449            DO j = 1, knon            u1(:knon) = yu(:knon, 1) + y_d_u(:knon, 1)
450               i = ni(j)            v1(:knon) = yv(:knon, 1) + y_d_v(:knon, 1)
451               u1(j) = yu(j, 1) + y_d_u(j, 1)            tair1(:knon) = yt(:knon, 1) + y_d_t(:knon, 1)
452               v1(j) = yv(j, 1) + y_d_v(j, 1)  
453               tair1(j) = yt(j, 1) + y_d_t(j, 1)            IF (nsrf == is_oce) THEN
454               qair1(j) = yq(j, 1) + y_d_q(j, 1)               rugo1(:knon) = frugs(ni(:knon), is_oce)
455               zgeo1(j) = rd * tair1(j) / (0.5 * (ypaprs(j, 1) + ypplay(j, &            else
456                    1))) * (ypaprs(j, 1)-ypplay(j, 1))               rugo1(:knon) = yrugos(:knon)
457               tairsol(j) = yts(j) + y_d_ts(j)            END IF
              rugo1(j) = yrugos(j)  
              IF (nsrf == is_oce) THEN  
                 rugo1(j) = frugs(i, nsrf)  
              END IF  
              psfce(j) = ypaprs(j, 1)  
              patm(j) = ypplay(j, 1)  
           END DO  
458    
459            CALL stdlevvar(nsrf, u1(:knon), v1(:knon), tair1(:knon), qair1, &            CALL stdlevvar(nsrf, u1(:knon), v1(:knon), tair1(:knon), &
460                 zgeo1, tairsol, yqsurf(:knon), rugo1, psfce, patm, yt2m, yq2m, &                 yq(:knon, 1) + y_d_q(:knon, 1), rd * tair1(:knon) &
461                 yt10m, yq10m, wind10m(:knon), ustar(:knon))                 / (0.5 * (ypaprs(:knon, 1) + ypplay(:knon, 1))) &
462                   * (ypaprs(:knon, 1) - ypplay(:knon, 1)), &
463                   yts(:knon) + y_d_ts(:knon), yqsurf(:knon), rugo1, &
464                   ypaprs(:knon, 1), ypplay(:knon, 1), yt2m, yq2m, yt10m, yq10m, &
465                   wind10m(:knon), ustar(:knon))
466    
467            DO j = 1, knon            DO j = 1, knon
468               i = ni(j)               i = ni(j)
Legend:
Removed from v.332  
changed lines
  Added in v.341
  ViewVC Help
Powered by ViewVC 1.1.21