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

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

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

revision 208 by guez, Wed Dec 7 16:44:53 2016 UTC revision 233 by guez, Tue Nov 7 10:52:46 2017 UTC
# Line 5  module coefkz_m Line 5  module coefkz_m
5  contains  contains
6    
7    SUBROUTINE coefkz(nsrf, paprs, pplay, ksta, ksta_ter, ts, rugos, u, v, t, &    SUBROUTINE coefkz(nsrf, paprs, pplay, ksta, ksta_ter, ts, rugos, u, v, t, &
8         q, qsurf, coefm, coefh)         q, qsurf, coefm, coefh, ycdragm, ycdragh)
9    
10      ! Authors: F. Hourdin, M. Forichon, Z. X. Li (LMD/CNRS)      ! Authors: F. Hourdin, M. Forichon, Z. X. Li (LMD/CNRS)
11      ! Date: September 22nd, 1993      ! Date: September 22nd, 1993
# Line 15  contains Line 15  contains
15      use clcdrag_m, only: clcdrag      use clcdrag_m, only: clcdrag
16      USE conf_phys_m, ONLY: iflag_pbl      USE conf_phys_m, ONLY: iflag_pbl
17      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
18      USE fcttre, ONLY: dqsatl, dqsats, foede, foeew, qsatl, qsats      USE fcttre, ONLY: foede, foeew
19      USE indicesol, ONLY: is_oce      USE indicesol, ONLY: is_oce
20      USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlstt, rlvtt, rtt      USE suphec_m, ONLY: rcpd, rd, retv, rg, rkappa, rlstt, rlvtt, rtt
21      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2      USE yoethf_m, ONLY: r2es, r5ies, r5les, rvtmp2
# Line 29  contains Line 29  contains
29      ! pression au milieu de chaque couche (en Pa)      ! pression au milieu de chaque couche (en Pa)
30    
31      REAL, intent(in):: ksta, ksta_ter      REAL, intent(in):: ksta, ksta_ter
32      REAL, intent(in):: ts(:) ! (klon) temperature du sol (en Kelvin)      REAL, intent(in):: ts(:) ! (knon) temperature du sol (en Kelvin)
33      REAL, intent(in):: rugos(:) ! (klon) longeur de rugosite (en m)      REAL, intent(in):: rugos(:) ! (klon) longeur de rugosite (en m)
34      REAL, intent(in):: u(:, :), v(:, :) ! (klon, klev) wind      REAL, intent(in):: u(:, :), v(:, :) ! (klon, klev) wind
35      REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K)      REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K)
36      real, intent(in):: q(:, :) ! (klon, klev) vapeur d'eau (kg/kg)      real, intent(in):: q(:, :) ! (klon, klev) vapeur d'eau (kg/kg)
37      real, intent(in):: qsurf(:) ! (klon)      real, intent(in):: qsurf(:) ! (knon)
38      REAL, intent(out):: coefm(:, :) ! (knon, klev) coefficient, vitesse      REAL, intent(out):: coefm(:, 2:) ! (knon, 2:klev) coefficient, vitesse
39    
40      real, intent(out):: coefh(:, :) ! (knon, klev)      real, intent(out):: coefh(:, 2:) ! (knon, 2:klev)
41      ! coefficient, chaleur et humidité      ! coefficient, chaleur et humidité
42    
43        real, intent(out):: ycdragm(:), ycdragh(:) ! (knon)
44    
45      ! Local:      ! Local:
46    
47      INTEGER knon ! nombre de points a traiter      INTEGER knon ! nombre de points a traiter
48        INTEGER itop(size(coefm, 1)) ! (knon) numero de couche du sommet
49      INTEGER itop(size(coefm, 1))                                   ! de la couche limite
     ! (knon) numero de couche du sommet de la couche limite  
50    
51      ! Quelques constantes et options:      ! Quelques constantes et options:
52    
# Line 133  contains Line 134  contains
134         z1(i) = zgeop(i, 1)         z1(i) = zgeop(i, 1)
135      ENDDO      ENDDO
136    
137      CALL clcdrag(klon, knon, nsrf, .false., u1, v1, t1, q1, z1, ts, qsurf, &      CALL clcdrag(nsrf, u1, v1, t1, q1, z1, ts, qsurf, rugos, ycdragm, ycdragh)
          rugos, coefm(:, 1), coefh(:, 1))  
138    
139      ! Calculer les coefficients turbulents dans l'atmosphere      ! Calculer les coefficients turbulents dans l'atmosphere
140    
Legend:
Removed from v.208  
changed lines
  Added in v.233
  ViewVC Help
Powered by ViewVC 1.1.21