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

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

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

revision 248 by guez, Fri Jan 5 16:40:13 2018 UTC revision 249 by guez, Fri Jan 5 17:15:05 2018 UTC
# Line 327  contains Line 327  contains
327                 yq(:knon, 1), zgeop(:knon, 1), yts(:knon), yqsurf(:knon), &                 yq(:knon, 1), zgeop(:knon, 1), yts(:knon), yqsurf(:knon), &
328                 yrugos(:knon), ycdragm(:knon), ycdragh(:knon))                 yrugos(:knon), ycdragm(:knon), ycdragh(:knon))
329    
330              IF (iflag_pbl == 1) THEN
331                 ycdragm(:knon) = max(ycdragm(:knon), 0.)
332                 ycdragh(:knon) = max(ycdragh(:knon), 0.)
333              end IF
334              
335              ! on met un seuil pour ycdragm et ycdragh
336              IF (nsrf == is_oce) THEN
337                 ycdragm(:knon) = min(ycdragm(:knon), cdmmax)
338                 ycdragh(:knon) = min(ycdragh(:knon), cdhmax)
339              END IF
340    
341            CALL coefkz(nsrf, ypaprs(:knon, :), ypplay(:knon, :), ksta, &            CALL coefkz(nsrf, ypaprs(:knon, :), ypplay(:knon, :), ksta, &
342                 ksta_ter, yts(:knon), yu(:knon, :), yv(:knon, :), yt(:knon, :), &                 ksta_ter, yts(:knon), yu(:knon, :), yv(:knon, :), yt(:knon, :), &
343                 yq(:knon, :), zgeop(:knon, :), ycoefm(:knon, :), &                 yq(:knon, :), zgeop(:knon, :), ycoefm(:knon, :), &
# Line 337  contains Line 348  contains
348                    ycoefh0(:knon, :))                    ycoefh0(:knon, :))
349               ycoefm(:knon, :) = max(ycoefm(:knon, :), ycoefm0(:knon, :))               ycoefm(:knon, :) = max(ycoefm(:knon, :), ycoefm0(:knon, :))
350               ycoefh(:knon, :) = max(ycoefh(:knon, :), ycoefh0(:knon, :))               ycoefh(:knon, :) = max(ycoefh(:knon, :), ycoefh0(:knon, :))
              ycdragm(:knon) = max(ycdragm(:knon), 0.)  
              ycdragh(:knon) = max(ycdragh(:knon), 0.)  
           END IF  
   
           ! on met un seuil pour ycdragm et ycdragh  
           IF (nsrf == is_oce) THEN  
              ycdragm(:knon) = min(ycdragm(:knon), cdmmax)  
              ycdragh(:knon) = min(ycdragh(:knon), cdhmax)  
351            END IF            END IF
352    
353            IF (ok_kzmin) THEN            IF (ok_kzmin) THEN
Legend:
Removed from v.248  
changed lines
  Added in v.249
  ViewVC Help
Powered by ViewVC 1.1.21