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

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

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

trunk/phylmd/stdlevvar.f revision 108 by guez, Tue Sep 16 14:00:41 2014 UTC trunk/Sources/phylmd/stdlevvar.f revision 188 by guez, Tue Mar 22 16:31:39 2016 UTC
# Line 11  contains Line 11  contains
11    
12      use coefcdrag_m, only: coefcdrag      use coefcdrag_m, only: coefcdrag
13      USE suphec_m, ONLY: rg, rkappa      USE suphec_m, ONLY: rg, rkappa
14        use screenp_m, only: screenp
15    
16      ! Objet : calcul de la température et de l'humidité relative à 2 m      ! Objet : calcul de la température et de l'humidité relative à 2 m
17      ! et du module du vent à 10 m à partir des relations de      ! et du module du vent à 10 m à partir des relations de
# Line 67  contains Line 68  contains
68      ! RKAR : constante de von Karman      ! RKAR : constante de von Karman
69      REAL, PARAMETER:: RKAR=0.40      REAL, PARAMETER:: RKAR=0.40
70      ! niter : nombre iterations calcul "corrector"      ! niter : nombre iterations calcul "corrector"
71      INTEGER, parameter:: niter=2, ncon=niter-1      INTEGER, parameter:: niter=2
72    
73      ! Variables locales      ! Variables locales
74      INTEGER i, n      INTEGER i, n
# Line 83  contains Line 84  contains
84      REAL, dimension(klon):: zdte, zdq      REAL, dimension(klon):: zdte, zdq
85      ! lmon : longueur de Monin-Obukhov selon Hess, Colman and McAvaney      ! lmon : longueur de Monin-Obukhov selon Hess, Colman and McAvaney
86      DOUBLE PRECISION, dimension(klon):: lmon      DOUBLE PRECISION, dimension(klon):: lmon
     DOUBLE PRECISION, parameter:: eps=1.0D-20  
87      REAL, dimension(klon):: delu, delte, delq      REAL, dimension(klon):: delu, delte, delq
88      REAL, dimension(klon):: u_zref, te_zref, q_zref      REAL, dimension(klon):: u_zref, te_zref, q_zref
89      REAL, dimension(klon):: temp, pref      REAL, dimension(klon):: temp, pref
90      LOGICAL okri      LOGICAL okri
91      REAL, dimension(klon):: u_zref_p, temp_p, q_zref_p      REAL, dimension(klon):: u_zref_p, temp_p, q_zref_p
92      !convertgence      !convertgence
     REAL, dimension(klon):: te_zref_con, q_zref_con  
93      REAL, dimension(klon):: u_zref_c, temp_c, q_zref_c      REAL, dimension(klon):: u_zref_c, temp_c, q_zref_c
94      REAL, dimension(klon):: ok_pred, ok_corr      REAL, dimension(klon):: ok_pred, ok_corr
95    
# Line 124  contains Line 123  contains
123    
124      ! First aproximation of variables at zref        ! First aproximation of variables at zref  
125      zref = 2.0      zref = 2.0
126      CALL screenp(klon, knon, nsrf, speed, tpot, q1, &      CALL screenp(klon, knon, speed, tpot, q1, &
127           ts1, qsurf, rugos, lmon, &           ts1, qsurf, rugos, lmon, &
128           ustar, testar, qstar, zref, &           ustar, testar, qstar, zref, &
129           delu, delte, delq)           delu, delte, delq)
# Line 157  contains Line 156  contains
156            ! return to normal temperature            ! return to normal temperature
157    
158            temp(i) = te_zref(i) * (psol(i)/pref(i))**(-RKAPPA)            temp(i) = te_zref(i) * (psol(i)/pref(i))**(-RKAPPA)
   
           IF(n == ncon) THEN  
              te_zref_con(i) = te_zref(i)  
              q_zref_con(i) = q_zref(i)  
           ENDIF  
159         ENDDO         ENDDO
160      ENDDO      ENDDO
161    
# Line 181  contains Line 175  contains
175      ! First aproximation of variables at zref        ! First aproximation of variables at zref  
176    
177      zref = 10.0      zref = 10.0
178      CALL screenp(klon, knon, nsrf, speed, tpot, q1, &      CALL screenp(klon, knon, speed, tpot, q1, &
179           ts1, qsurf, rugos, lmon, &           ts1, qsurf, rugos, lmon, &
180           ustar, testar, qstar, zref, &           ustar, testar, qstar, zref, &
181           delu, delte, delq)           delu, delte, delq)
Legend:
Removed from v.108  
changed lines
  Added in v.188
  ViewVC Help
Powered by ViewVC 1.1.21