/[lmdze]/trunk/phylmd/pbl_surface.f
ViewVC logotype

Diff of /trunk/phylmd/pbl_surface.f

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

revision 242 by guez, Mon Nov 13 12:12:41 2017 UTC revision 243 by guez, Tue Nov 14 14:38:36 2017 UTC
# Line 310  contains Line 310  contains
310               END DO               END DO
311            END DO            END DO
312    
           ! calculer Cdrag et les coefficients d'echange  
313            CALL coefkz(nsrf, ypaprs, ypplay, ksta, ksta_ter, yts(:knon), &            CALL coefkz(nsrf, ypaprs, ypplay, ksta, ksta_ter, yts(:knon), &
314                 yrugos, yu, yv, yt, yq, yqsurf(:knon), coefm(:knon, :), &                 yrugos, yu, yv, yt, yq, yqsurf(:knon), coefm(:knon, :), &
315                 coefh(:knon, :), ycdragm(:knon), ycdragh(:knon))                 coefh(:knon, :), ycdragm(:knon), ycdragh(:knon))
# Line 416  contains Line 415  contains
415               y_dflux_q(j) = y_dflux_q(j) * ypct(j)               y_dflux_q(j) = y_dflux_q(j) * ypct(j)
416            END DO            END DO
417    
           DO k = 2, klev  
              DO j = 1, knon  
                 i = ni(j)  
                 coefh(j, k) = coefh(j, k) * ypct(j)  
                 coefm(j, k) = coefm(j, k) * ypct(j)  
              END DO  
           END DO  
           DO j = 1, knon  
              i = ni(j)  
              ycdragh(j) = ycdragh(j) * ypct(j)  
              ycdragm(j) = ycdragm(j) * ypct(j)  
           END DO  
418            DO k = 1, klev            DO k = 1, klev
419               DO j = 1, knon               DO j = 1, knon
420                  i = ni(j)                  i = ni(j)
# Line 464  contains Line 451  contains
451               agesno(i, nsrf) = yagesno(j)               agesno(i, nsrf) = yagesno(j)
452               fqcalving(i, nsrf) = y_fqcalving(j)               fqcalving(i, nsrf) = y_fqcalving(j)
453               ffonte(i, nsrf) = y_ffonte(j)               ffonte(i, nsrf) = y_ffonte(j)
454               cdragh(i) = cdragh(i) + ycdragh(j)               cdragh(i) = cdragh(i) + ycdragh(j) * ypct(j)
455               cdragm(i) = cdragm(i) + ycdragm(j)               cdragm(i) = cdragm(i) + ycdragm(j) * ypct(j)
456               dflux_t(i) = dflux_t(i) + y_dflux_t(j)               dflux_t(i) = dflux_t(i) + y_dflux_t(j)
457               dflux_q(i) = dflux_q(i) + y_dflux_q(j)               dflux_q(i) = dflux_q(i) + y_dflux_q(j)
458            END DO            END DO
# Line 491  contains Line 478  contains
478               END DO               END DO
479            END DO            END DO
480    
481            ycoefh(ni(:knon), :) = ycoefh(ni(:knon), :) + coefh(:knon, :)            forall (k = 2:klev) ycoefh(ni(:knon), k) &
482                   = ycoefh(ni(:knon), k) + coefh(:knon, k) * ypct(:knon)
483    
484            ! diagnostic t, q a 2m et u, v a 10m            ! diagnostic t, q a 2m et u, v a 10m
485    
Legend:
Removed from v.242  
changed lines
  Added in v.243
  ViewVC Help
Powered by ViewVC 1.1.21