/[lmdze]/trunk/libf/phylmd/Conflx/flxmain.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/Conflx/flxmain.f90

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

revision 72 by guez, Mon Jul 8 18:12:18 2013 UTC revision 73 by guez, Fri Nov 15 17:48:30 2013 UTC
# Line 74  contains Line 74  contains
74         firstcal = .FALSE.         firstcal = .FALSE.
75      ENDIF      ENDIF
76    
77      DO i = 1, klon      ldcum = .FALSE.
78         ldcum(i) = .FALSE.      dt_con = 0.
79      ENDDO      dq_con = 0.
     DO k = 1, klev  
        DO i = 1, klon  
           dt_con(i, k) = 0.0  
           dq_con(i, k) = 0.0  
        ENDDO  
     ENDDO  
   
     ! initialiser les variables et faire l'interpolation verticale  
80    
81        ! Initialiser les variables et faire l'interpolation verticale :
82      CALL flxini(ten, qen, qsen, pgeo, paph, zgeoh, ztenh, zqenh, zqsenh, &      CALL flxini(ten, qen, qsen, pgeo, paph, zgeoh, ztenh, zqenh, zqsenh, &
83           ptu, pqu, ptd, pqd, mfd, zmfds, zmfdq, zdmfdp, mfu, mfus, mfuq, &           ptu, pqu, ptd, pqd, mfd, zmfds, zmfdq, zdmfdp, mfu, mfus, mfuq, &
84           zdmfup, zdpmel, plu, plude, ilab, pen_u, pde_u, pen_d, pde_d)           zdmfup, zdpmel, plu, plude, ilab, pen_u, pde_u, pen_d, pde_d)
85    
86      ! determiner les valeurs au niveau de base de la tour convective      ! Déterminer les valeurs au niveau de base de la tour convective :
87        CALL flxbase(ztenh, zqenh, zgeoh, paph, ptu, pqu, plu, ldcum, kcbot, ilab)
     CALL flxbase(ztenh, zqenh, zgeoh, paph, &  
          ptu, pqu, plu, ldcum, kcbot, ilab)  
   
     ! calculer la convergence totale de l'humidite et celle en provenance  
     ! de la couche limite, plus precisement, la convergence integree entre  
     ! le sol et la base de la convection. Cette derniere convergence est  
     ! comparee avec l'evaporation obtenue dans la couche limite pour  
     ! determiner le type de la convection  
88    
89      k=1      ! Calculer la convergence totale de l'humidité et celle en
90      DO i = 1, klon      ! provenance de la couche limite, plus précisément, la convergence
91         zdqcv(i) = pqte(i, k)*(paph(i, k+1)-paph(i, k))      ! intégrée entre le sol et la base de la convection. Cette
92         zdhpbl(i) = 0.0      ! dernière convergence est comparée avec l'&vaporation obtenue
93         zdqpbl(i) = 0.0      ! dans la couche limite pour déterminer le type de la convection.
94      ENDDO  
95        zdqcv = pqte(:, 1) * (paph(:, 2) - paph(:, 1))
96        zdhpbl = 0.
97        zdqpbl = 0.
98    
99      DO k=2, klev      DO k=2, klev
100         DO i = 1, klon         DO i = 1, klon
# Line 194  contains Line 182  contains
182            IF(ktype(i) == 2) zentr(i)=ENTRSCV            IF(ktype(i) == 2) zentr(i)=ENTRSCV
183         ENDDO         ENDDO
184    
185         IF (lmfdd) THEN ! si l'on considere le panache descendant         downdraft: IF (lmfdd) THEN
186              ! si l'on considere le panache descendant
187            ! calculer la precipitation issue du panache ascendant pour            ! calculer la precipitation issue du panache ascendant pour
188            ! determiner l'existence du panache descendant dans la convection            ! determiner l'existence du panache descendant dans la convection
189            DO i = 1, klon            DO i = 1, klon
# Line 266  contains Line 255  contains
255            DO i = 1, klon            DO i = 1, klon
256               IF (lddraf(i)) mfub(i)=mfub1(i)               IF (lddraf(i)) mfub(i)=mfub1(i)
257            ENDDO            ENDDO
258         ENDIF ! fin de test sur lmfdd         ENDIF downdraft
259    
260         ! calculer de nouveau le panache ascendant         ! calculer de nouveau le panache ascendant
261    
Legend:
Removed from v.72  
changed lines
  Added in v.73
  ViewVC Help
Powered by ViewVC 1.1.21