source: LMDZ5/trunk/libf/phylmd/thermcell_height.F90 @ 2306

Last change on this file since 2306 was 1998, checked in by fhourdin, 11 years ago

Reactivation du calcul d'un zmax continu pour les thermiques

(thermcell_height.F90, thermcell_plume.F90, thermcell_dry.F90)
ouvre la voie à la réactivation d'une fermeture humide des thermiques
iflag_thermals_closure=2
(conf_phys_m.F90, thermcell.h, thermcell_main.F90)

Modification liée à la conservation de l'eau

(add_phys_tend.F90, add_pbl_tend.F90, physiq.F90)

Modifications liées au déclenchement stochastique

  1. possibilité de revenir à la Ale déterministe pour le criter ALE>|CIN| iflag_trig_bl=2, 1 par défaut)
  2. possibilité d'activer une fermeture statistique où ALP est divisé par la probabilité de déclenchement iflag_clos_bl=1 (0 par défaut, ancienne option 1 passée en =2)

Modification de l'entrainemement dans la version "stratocumulus" du

modèle du thermique (quand iflag_thermals_ed=8).
(modifie thermcell_plume.F90)

Catherine, Jean-Yves et Frédéric

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.7 KB
Line 
1      SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lmix,  &
2     &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax,lev_out)                           
3
4!-----------------------------------------------------------------------------
5!thermcell_height: calcul des caracteristiques du thermique: zmax,wmax,zmix
6!-----------------------------------------------------------------------------
7      IMPLICIT NONE
8#include "iniprint.h"
9#include "thermcell.h"
10
11      INTEGER ig,l
12      INTEGER ngrid,nlay
13      INTEGER lalim(ngrid),lmin(ngrid)
14      INTEGER lmix(ngrid)
15      REAL linter(ngrid)
16      integer lev_out                           ! niveau pour les print
17
18      REAL zw2(ngrid,nlay+1)
19      REAL zlev(ngrid,nlay+1)
20
21      REAL wmax(ngrid)
22      INTEGER lmax(ngrid)
23      REAL zmax(ngrid)
24      REAL zmax0(ngrid)
25      REAL zmix(ngrid)
26      REAL num(ngrid)
27      REAL denom(ngrid)
28
29      REAL zlevinter(ngrid)
30
31!calcul de la hauteur max du thermique
32      do ig=1,ngrid
33         lmax(ig)=lalim(ig)
34      enddo
35      do ig=1,ngrid
36         do l=nlay,lalim(ig)+1,-1
37            if (zw2(ig,l).le.1.e-10) then
38               lmax(ig)=l-1
39            endif
40         enddo
41      enddo
42
43! On traite le cas particulier qu'il faudrait éviter ou le thermique
44! atteind le haut du modele ...
45      do ig=1,ngrid
46      if ( zw2(ig,nlay) > 1.e-10 ) then
47          print*,'WARNING !!!!! W2 thermiques non nul derniere couche '
48          lmax(ig)=nlay
49      endif
50      enddo
51
52! pas de thermique si couche 1 stable
53      do ig=1,ngrid
54         if (lmin(ig).gt.1) then
55             lmax(ig)=1
56             lmin(ig)=1
57             lalim(ig)=1
58         endif
59      enddo
60!   
61! Determination de zw2 max
62      do ig=1,ngrid
63         wmax(ig)=0.
64      enddo
65
66      do l=1,nlay
67         do ig=1,ngrid
68            if (l.le.lmax(ig)) then
69                if (zw2(ig,l).lt.0.)then
70                  print*,'pb2 zw2<0'
71                endif
72                zw2(ig,l)=sqrt(zw2(ig,l))
73                wmax(ig)=max(wmax(ig),zw2(ig,l))
74            else
75                 zw2(ig,l)=0.
76            endif
77          enddo
78      enddo
79
80!   Longueur caracteristique correspondant a la hauteur des thermiques.
81      do  ig=1,ngrid
82         zmax(ig)=0.
83         zlevinter(ig)=zlev(ig,1)
84      enddo
85
86!     if (iflag_thermals_ed.ge.1) then
87      if (1==0) then
88!CR:date de quand le calcul du zmax continu etait buggue
89         num(:)=0.
90         denom(:)=0.
91         do ig=1,ngrid
92          do l=1,nlay
93             num(ig)=num(ig)+zw2(ig,l)*zlev(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
94             denom(ig)=denom(ig)+zw2(ig,l)*(zlev(ig,l+1)-zlev(ig,l))
95          enddo
96       enddo
97       do ig=1,ngrid
98       if (denom(ig).gt.1.e-10) then
99          zmax(ig)=2.*num(ig)/denom(ig)
100          zmax0(ig)=zmax(ig)
101       endif
102       enddo
103 
104      else
105!CR:Calcul de zmax continu via le linter     
106      do  ig=1,ngrid
107! calcul de zlevinter
108          zlevinter(ig)=(zlev(ig,lmax(ig)+1)-zlev(ig,lmax(ig)))*  &
109     &    linter(ig)+zlev(ig,lmax(ig))-lmax(ig)*(zlev(ig,lmax(ig)+1)  &
110     &    -zlev(ig,lmax(ig)))
111!pour le cas ou on prend tjs lmin=1
112!       zmax(ig)=max(zmax(ig),zlevinter(ig)-zlev(ig,lmin(ig)))
113       zmax(ig)=max(zmax(ig),zlevinter(ig)-zlev(ig,1))
114       zmax0(ig)=zmax(ig)
115      enddo
116
117
118      endif
119!endif iflag_thermals_ed
120!
121! def de  zmix continu (profil parabolique des vitesses)
122      do ig=1,ngrid
123           if (lmix(ig).gt.1) then
124! test
125              if (((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
126     &        *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
127     &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
128     &        *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))).gt.1e-10)  &
129     &        then
130!             
131            zmix(ig)=((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
132     &        *((zlev(ig,lmix(ig)))**2-(zlev(ig,lmix(ig)+1))**2)  &
133     &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
134     &        *((zlev(ig,lmix(ig)-1))**2-(zlev(ig,lmix(ig)))**2))  &
135     &        /(2.*((zw2(ig,lmix(ig)-1)-zw2(ig,lmix(ig)))  &
136     &        *((zlev(ig,lmix(ig)))-(zlev(ig,lmix(ig)+1)))  &
137     &        -(zw2(ig,lmix(ig))-zw2(ig,lmix(ig)+1))  &
138     &        *((zlev(ig,lmix(ig)-1))-(zlev(ig,lmix(ig))))))
139              else
140              zmix(ig)=zlev(ig,lmix(ig))
141              print*,'pb zmix'
142              endif
143          else
144              zmix(ig)=0.
145          endif
146!test
147         if ((zmax(ig)-zmix(ig)).le.0.) then
148            zmix(ig)=0.9*zmax(ig)
149!            print*,'pb zmix>zmax'
150         endif
151      enddo
152!
153! calcul du nouveau lmix correspondant
154      do ig=1,ngrid
155         do l=1,nlay
156            if (zmix(ig).ge.zlev(ig,l).and.  &
157     &          zmix(ig).lt.zlev(ig,l+1)) then
158              lmix(ig)=l
159             endif
160          enddo
161      enddo
162!
163      return
164      end
Note: See TracBrowser for help on using the repository browser.