/[lmdze]/trunk/Sources/phylmd/clmain.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/clmain.f

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

revision 228 by guez, Fri Nov 3 12:38:47 2017 UTC revision 229 by guez, Mon Nov 6 17:20:45 2017 UTC
# Line 95  contains Line 95  contains
95      ! flux de vapeur d'eau (kg / m2 / s) à la surface      ! flux de vapeur d'eau (kg / m2 / s) à la surface
96    
97      REAL, intent(out):: flux_u(klon, nbsrf), flux_v(klon, nbsrf)      REAL, intent(out):: flux_u(klon, nbsrf), flux_v(klon, nbsrf)
98      ! tension du vent à la surface, en Pa      ! tension du vent (flux turbulent de vent) à la surface, en Pa
99    
100      REAL, INTENT(out):: cdragh(klon), cdragm(klon)      REAL, INTENT(out):: cdragh(klon), cdragm(klon)
101      real q2(klon, klev + 1, nbsrf)      real q2(klon, klev + 1, nbsrf)
# Line 369  contains Line 369  contains
369               END DO               END DO
370    
371               ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), coefm(:knon, 1))               ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), coefm(:knon, 1))
   
              ! iflag_pbl peut \^etre utilis\'e comme longueur de m\'elange  
   
372               CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), &               CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), &
373                    yu(:knon, :), yv(:knon, :), yteta(:knon, :), &                    yu(:knon, :), yv(:knon, :), yteta(:knon, :), &
374                    coefm(:knon, 1), yq2(:knon, :), ykmm(:knon, :), &                    coefm(:knon, 1), yq2(:knon, :), ykmm(:knon, :), &
375                    ykmn(:knon, :), ykmq(:knon, :), ustar(:knon), iflag_pbl)                    ykmn(:knon, :), ykmq(:knon, :), ustar(:knon))
   
376               coefm(:knon, 2:) = ykmm(:knon, 2:klev)               coefm(:knon, 2:) = ykmm(:knon, 2:klev)
377               coefh(:knon, 2:) = ykmn(:knon, 2:klev)               coefh(:knon, 2:) = ykmn(:knon, 2:klev)
378            END IF            END IF
379    
380            ! calculer la diffusion des vitesses "u" et "v"            CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), &
381            CALL clvent(knon, dtime, yu(:knon, 1), yv(:knon, 1), &                 yt(:knon, :), yu(:knon, :), ypaprs(:knon, :), &
382                 coefm(:knon, :), yt, yu, ypaprs, ypplay, ydelp, y_d_u, &                 ypplay(:knon, :), ydelp(:knon, :), y_d_u(:knon, :), &
383                 y_flux_u(:knon))                 y_flux_u(:knon))
384            CALL clvent(knon, dtime, yu(:knon, 1), yv(:knon, 1), &            CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), &
385                 coefm(:knon, :), yt, yv, ypaprs, ypplay, ydelp, y_d_v, &                 yt(:knon, :), yv(:knon, :), ypaprs(:knon, :), &
386                   ypplay(:knon, :), ydelp(:knon, :), y_d_v(:knon, :), &
387                 y_flux_v(:knon))                 y_flux_v(:knon))
388    
389            ! calculer la diffusion de "q" et de "h"            ! calculer la diffusion de "q" et de "h"
Legend:
Removed from v.228  
changed lines
  Added in v.229
  ViewVC Help
Powered by ViewVC 1.1.21