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

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

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

revision 286 by guez, Tue Jul 24 15:22:48 2018 UTC revision 288 by guez, Tue Jul 24 16:27:12 2018 UTC
# Line 4  module coefkz2_m Line 4  module coefkz2_m
4    
5  contains  contains
6    
7    SUBROUTINE coefkz2(nsrf, paprs, pplay, t, pcfm, pcfh)    SUBROUTINE coefkz2(nsrf, paprs, pplay, t, coefm, coefh)
8    
9      ! J'introduit un peu de diffusion sauf dans les endroits o\`u une      ! J'introduit un peu de diffusion sauf dans les endroits o\`u une
10      ! forte inversion est pr\'esente. On peut dire que la diffusion      ! forte inversion est pr\'esente. On peut dire que la diffusion
# Line 24  contains Line 24  contains
24            
25      REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K)      REAL, intent(in):: t(:, :) ! (knon, klev) temperature (K)
26    
27      REAL, intent(out):: pcfm(:, 2:) ! (knon, 2:klev) coefficient vitesse      REAL, intent(out):: coefm(:, 2:) ! (knon, 2:klev) coefficient vitesse
28    
29      REAL, intent(out):: pcfh(:, 2:) ! (knon, 2:klev)      REAL, intent(out):: coefh(:, 2:) ! (knon, 2:klev)
30      ! coefficient chaleur et humidite)      ! coefficient chaleur et humidite)
31    
32      ! Local:      ! Local:
# Line 56  contains Line 56  contains
56      ! Initialiser les sorties      ! Initialiser les sorties
57      DO k = 2, klev      DO k = 2, klev
58         DO i = 1, knon         DO i = 1, knon
59            pcfm(i, k) = 0.0            coefm(i, k) = 0.0
60            pcfh(i, k) = 0.0            coefh(i, k) = 0.0
61         ENDDO         ENDDO
62      ENDDO      ENDDO
63    
# Line 88  contains Line 88  contains
88                 ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN                 ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
89               zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &               zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
90                    /(paprs(i, 2)-paprs(i, klev+1))))**2                    /(paprs(i, 2)-paprs(i, klev+1))))**2
91               pcfm(i, k)= zl2(i)* kstable               coefm(i, k)= zl2(i)* kstable
92               pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different               coefh(i, k) = coefm(i, k) /prandtl ! h et m different
93            ENDIF            ENDIF
94         ENDDO         ENDDO
95      ENDDO      ENDDO
Legend:
Removed from v.286  
changed lines
  Added in v.288
  ViewVC Help
Powered by ViewVC 1.1.21