/[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 237 by guez, Thu Nov 9 13:26:00 2017 UTC revision 238 by guez, Thu Nov 9 14:11:39 2017 UTC
# Line 172  contains Line 172  contains
172      REAL ycoefm0(klon, klev), ycoefh0(klon, klev)      REAL ycoefm0(klon, klev), ycoefh0(klon, klev)
173      REAL yzlay(klon, klev), zlev(klon, klev + 1), yteta(klon, klev)      REAL yzlay(klon, klev), zlev(klon, klev + 1), yteta(klon, klev)
174      REAL ykmm(klon, klev + 1), ykmn(klon, klev + 1)      REAL ykmm(klon, klev + 1), ykmn(klon, klev + 1)
     REAL ykmq(klon, klev + 1)  
175      REAL yq2(klon, klev + 1)      REAL yq2(klon, klev + 1)
176      REAL delp(klon, klev)      REAL delp(klon, klev)
177      INTEGER i, k, nsrf      INTEGER i, k, nsrf
# Line 323  contains Line 322  contains
322               ycoefh0(:knon, 1) = 0.               ycoefh0(:knon, 1) = 0.
323               coefm(:knon, :) = max(coefm(:knon, :), ycoefm0(:knon, 2:))               coefm(:knon, :) = max(coefm(:knon, :), ycoefm0(:knon, 2:))
324               coefh(:knon, :) = max(coefh(:knon, :), ycoefh0(:knon, 2:))               coefh(:knon, :) = max(coefh(:knon, :), ycoefh0(:knon, 2:))
325               ycdragm(:knon) = max(ycdragm(:knon), ycoefm0(:knon, 1))               ycdragm(:knon) = max(ycdragm(:knon), 0.)
326               ycdragh(:knon) = max(ycdragh(:knon), ycoefh0(:knon, 1))               ycdragh(:knon) = max(ycdragh(:knon), 0.)
327            END IF            END IF
328    
329            ! on met un seuil pour ycdragm et ycdragh            ! on met un seuil pour ycdragm et ycdragh
# Line 379  contains Line 378  contains
378    
379               ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), ycdragm(:knon))               ustar(:knon) = ustarhb(yu(:knon, 1), yv(:knon, 1), ycdragm(:knon))
380               CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), &               CALL yamada4(dtime, rg, zlev(:knon, :), yzlay(:knon, :), &
381                    yu(:knon, :), yv(:knon, :), yteta(:knon, :), &                    yu(:knon, :), yv(:knon, :), yteta(:knon, :), yq2(:knon, :), &
382                    ycdragm(:knon), yq2(:knon, :), ykmm(:knon, :), &                    ykmm(:knon, :), ykmn(:knon, :), ustar(:knon))
383                    ykmn(:knon, :), ykmq(:knon, :), ustar(:knon))               coefm(:knon, :) = ykmm(:knon, 2:klev)
384               coefm(:knon, 2:) = ykmm(:knon, 2:klev)               coefh(:knon, :) = ykmn(:knon, 2:klev)
              coefh(:knon, 2:) = ykmn(:knon, 2:klev)  
385            END IF            END IF
386    
387            CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), &            CALL clvent(dtime, yu(:knon, 1), yv(:knon, 1), coefm(:knon, :), &
Legend:
Removed from v.237  
changed lines
  Added in v.238
  ViewVC Help
Powered by ViewVC 1.1.21