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

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

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

revision 288 by guez, Tue Jul 24 16:27:12 2018 UTC revision 289 by guez, Tue Jul 24 16:40:15 2018 UTC
# Line 48  contains Line 48  contains
48           size(v, 1), size(q, 1), size(t, 1), size(ts), size(cdragm), &           size(v, 1), size(q, 1), size(t, 1), size(ts), size(cdragm), &
49           size(zgeop, 1), size(coefm, 1), size(coefh, 1), size(q2, 1)], &           size(zgeop, 1), size(coefm, 1), size(coefh, 1), size(q2, 1)], &
50           "coef_diff_turb knon")           "coef_diff_turb knon")
       
     CALL coefkz(nsrf, paprs, pplay, ts, u, v, t, q, zgeop, coefm, coefh)  
   
     IF (iflag_pbl == 1) THEN  
        CALL coefkz2(nsrf, paprs, pplay, t, coefm0, coefh0)  
        coefm = max(coefm, coefm0)  
        coefh = max(coefh, coefh0)  
     END IF  
   
     IF (ok_kzmin) THEN  
        ! Calcul d'une diffusion minimale pour les conditions tres stables  
        CALL coefkzmin(paprs, pplay, u, v, t, q, cdragm, coefh0)  
        coefm = max(coefm, coefh0)  
        coefh = max(coefh, coefh0)  
     END IF  
51    
52      IF (iflag_pbl >= 6) THEN      IF (iflag_pbl >= 6) THEN
53         ! Mellor et Yamada adapt\'e \`a Mars, Richard Fournier et         ! Mellor et Yamada adapt\'e \`a Mars, Richard Fournier et
# Line 84  contains Line 69  contains
69    
70         CALL yamada4(dtime, zlev, zlay, u, v, teta, q2, coefm, coefh, &         CALL yamada4(dtime, zlev, zlay, u, v, teta, q2, coefm, coefh, &
71              ustarhb(u(:, 1), v(:, 1), cdragm))              ustarhb(u(:, 1), v(:, 1), cdragm))
72        else
73           CALL coefkz(nsrf, paprs, pplay, ts, u, v, t, q, zgeop, coefm, coefh)
74    
75           IF (iflag_pbl == 1) THEN
76              CALL coefkz2(nsrf, paprs, pplay, t, coefm0, coefh0)
77              coefm = max(coefm, coefm0)
78              coefh = max(coefh, coefh0)
79           END IF
80    
81           IF (ok_kzmin) THEN
82              ! Calcul d'une diffusion minimale pour les conditions tres stables
83              CALL coefkzmin(paprs, pplay, u, v, t, q, cdragm, coefh0)
84              coefm = max(coefm, coefh0)
85              coefh = max(coefh, coefh0)
86           END IF
87      END IF      END IF
88    
89    end subroutine coef_diff_turb    end subroutine coef_diff_turb
Legend:
Removed from v.288  
changed lines
  Added in v.289
  ViewVC Help
Powered by ViewVC 1.1.21