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

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

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

trunk/libf/phylmd/coefkzmin.f90 revision 47 by guez, Fri Jul 1 15:00:48 2011 UTC trunk/Sources/phylmd/coefkzmin.f revision 233 by guez, Tue Nov 7 10:52:46 2017 UTC
# Line 19  contains Line 19  contains
19      REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1      REAL, intent(in):: yu(klon, klev), yv(klon, klev) ! wind, in m s-1
20      REAL, intent(in):: yt(klon, klev) ! temperature, in K      REAL, intent(in):: yt(klon, klev) ! temperature, in K
21      REAL, intent(in):: yq(klon, klev)      REAL, intent(in):: yq(klon, klev)
22        REAL, intent(in):: ycoefm(:) ! (ngrid) drag coefficient
23    
24      REAL, intent(in):: ycoefm(klon) ! drag coefficient      REAL, intent(inout):: km(:, 2:) ! (knon, 2:klev)
   
     REAL, intent(inout):: km(klon, klev)  
25      ! coefficient de diffusion turbulente de quantité de mouvement (au      ! coefficient de diffusion turbulente de quantité de mouvement (au
26      ! bas de chaque couche) (en sortie : la valeur à la fin du pas de      ! bas de chaque couche) (en sortie : la valeur à la fin du pas de
27      ! temps), m2 s-1      ! temps), m2 s-1
28    
29      REAL, intent(inout):: kn(klon, klev)      REAL, intent(inout):: kn(:, 2:) ! (knon, 2:klev)
30      ! coefficient de diffusion turbulente des scalaires (au bas de      ! coefficient de diffusion turbulente des scalaires (au bas de
31      ! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1      ! chaque couche) (en sortie : la valeur à la fin du pas de temps), m2 s-1
32    
# Line 74  contains Line 73  contains
73      enddo      enddo
74    
75      forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))      forall (k = 2: klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))
76      ustar = SQRT(ycoefm(:ngrid) * (yu(:ngrid, 1)**2 + yv(:ngrid, 1)**2))      ustar = SQRT(ycoefm * (yu(:ngrid, 1)**2 + yv(:ngrid, 1)**2))
77    
78      ! Fin de la partie qui doit être incluse à terme dans clmain      ! Fin de la partie qui doit être incluse à terme dans clmain
79    
Legend:
Removed from v.47  
changed lines
  Added in v.233
  ViewVC Help
Powered by ViewVC 1.1.21