/[lmdze]/trunk/Sources/phylmd/vdif_kcay.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/vdif_kcay.f

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

trunk/phylmd/vdif_kcay.f revision 105 by guez, Thu Sep 4 10:40:24 2014 UTC trunk/Sources/phylmd/vdif_kcay.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 4  module vdif_kcay_m Line 4  module vdif_kcay_m
4    
5  contains  contains
6    
7    SUBROUTINE vdif_kcay(ngrid, dt, g, rconst, plev, temp, zlev, zlay, &    SUBROUTINE vdif_kcay(ngrid, dt, g, plev, zlev, zlay, u, v, teta, cd, q2, &
8         u, v, teta, cd, q2, q2diag, km, kn, ustar, l_mix)         q2diag, km, kn, ustar, l_mix)
9    
10      ! From LMDZ4/libf/phylmd/vdif_kcay.F, version 1.1 2004/06/22 11:45:36      ! From LMDZ4/libf/phylmd/vdif_kcay.F, version 1.1 2004/06/22 11:45:36
11    
12      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
13        use yamada_m, only: yamada
14    
15        INTEGER ngrid
16      ! dt : pas de temps      ! dt : pas de temps
17      ! g : g      ! g : g
18      ! zlev : altitude a chaque niveau (interface inferieure de la couche      ! zlev : altitude a chaque niveau (interface inferieure de la couche
# Line 31  contains Line 33  contains
33    
34      REAL, intent(in):: dt      REAL, intent(in):: dt
35      real, intent(in):: g      real, intent(in):: g
36      real rconst      real plev(klon, klev+1)
     real plev(klon, klev+1), temp(klon, klev)  
37      real ustar(klon), snstable      real ustar(klon), snstable
38      REAL zlev(klon, klev+1)      REAL zlev(klon, klev+1)
39      REAL zlay(klon, klev)      REAL zlay(klon, klev)
# Line 57  contains Line 58  contains
58      ! q : echelle de vitesse au bas de chaque couche      ! q : echelle de vitesse au bas de chaque couche
59      ! (valeur a la fin du pas de temps)      ! (valeur a la fin du pas de temps)
60    
61      INTEGER nlay, nlev, ngrid      INTEGER nlay, nlev
62      REAL unsdz(klon, klev)      REAL unsdz(klon, klev)
63      REAL unsdzdec(klon, klev+1)      REAL unsdzdec(klon, klev+1)
64      REAL q(klon, klev+1)      REAL q(klon, klev+1)
# Line 106  contains Line 107  contains
107      ! gnmax : borne superieure de gn (0.0233)      ! gnmax : borne superieure de gn (0.0233)
108      ! gninf : vrai si gn est en dessous de sa borne inferieure      ! gninf : vrai si gn est en dessous de sa borne inferieure
109      ! gnsup : vrai si gn est en dessus de sa borne superieure      ! gnsup : vrai si gn est en dessus de sa borne superieure
     ! gm : drole d'objet bien utile  
110      ! ri : nombre de Richardson      ! ri : nombre de Richardson
111      ! sn : coefficient de stabilite pour n      ! sn : coefficient de stabilite pour n
112      ! snq2 : premier terme du developement limite de sn en q2      ! snq2 : premier terme du developement limite de sn en q2
# Line 118  contains Line 118  contains
118      REAL gnmax      REAL gnmax
119      LOGICAL gninf      LOGICAL gninf
120      LOGICAL gnsup      LOGICAL gnsup
     REAL gm  
     ! REAL ri(klon, klev+1)  
121      REAL sn(klon, klev+1)      REAL sn(klon, klev+1)
122      REAL snq2(klon, klev+1)      REAL snq2(klon, klev+1)
123      REAL sm(klon, klev+1)      REAL sm(klon, klev+1)
# Line 196  contains Line 194  contains
194    
195      ! Initialisation de q2      ! Initialisation de q2
196    
197      call yamada(ngrid, g, rconst, plev, temp, zlev, zlay, u, v, teta, &      call yamada(ngrid, g, zlev, zlay, u, v, teta, q2diag, km, kn)
          q2diag, km, kn, ustar, l_mix)  
198      if (first.and.1.eq.1) then      if (first.and.1.eq.1) then
199         first=.false.         first=.false.
200         q2=q2diag         q2=q2diag
# Line 326  contains Line 323  contains
323    
324            gn=-long(igrid, ilev)**2 / q2(igrid, ilev) &            gn=-long(igrid, ilev)**2 / q2(igrid, ilev) &
325                 * n2(igrid, ilev)                 * n2(igrid, ilev)
           gm=long(igrid, ilev)**2 / q2(igrid, ilev) &  
                * m2(igrid, ilev)  
   
326            gninf=.false.            gninf=.false.
327            gnsup=.false.            gnsup=.false.
328            long(igrid, ilev)=long(igrid, ilev)            long(igrid, ilev)=long(igrid, ilev)
Legend:
Removed from v.105  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21