source: LMDZ5/trunk/libf/phylmd/add_phys_tend.F90 @ 1998

Last change on this file since 1998 was 1998, checked in by fhourdin, 10 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: 7.4 KB
Line 
1!
2! $Id: add_phys_tend.F90 1998 2014-03-19 14:14:30Z fhourdin $
3!
4SUBROUTINE add_phys_tend (zdu,zdv,zdt,zdq,zdql,paprs,text)
5!======================================================================
6! Ajoute les tendances des variables physiques aux variables
7! d'etat de la dynamique t_seri, q_seri ...
8! On en profite pour faire des tests sur les tendances en question.
9!======================================================================
10
11
12!======================================================================
13! Declarations
14!======================================================================
15
16use dimphy
17use phys_local_var_mod
18use phys_state_var_mod
19IMPLICIT none
20#include "iniprint.h"
21#include "YOMCST.h"
22
23! Arguments :
24!------------
25REAL zdu(klon,klev),zdv(klon,klev)
26REAL zdt(klon,klev),zdq(klon,klev),zdql(klon,klev)
27REAL paprs(klon,klev+1)
28CHARACTER*(*) text
29
30! Local :
31!--------
32REAL zt,zq
33REAL zq_int, zqp_int, zq_new
34
35REAL zqp(klev)
36
37INTEGER i, k,j
38INTEGER jadrs(klon*klev), jbad
39INTEGER jqadrs(klon*klev), jqbad
40INTEGER kadrs(klon*klev)
41INTEGER kqadrs(klon*klev)
42
43LOGICAL done(klon)
44
45integer debug_level
46logical, save :: first=.true.
47!$OMP THREADPRIVATE(first)
48INTEGER, SAVE :: itap
49!$OMP THREADPRIVATE(itap)
50!======================================================================
51! Initialisations
52
53debug_level=10
54     if (first) then
55        itap=0
56        first=.false.
57     endif
58! Incrementer le compteur de la physique
59     itap   = itap + 1
60!======================================================================
61! Ajout des tendances sur le vent et l'eau liquide
62!======================================================================
63
64     u_seri(:,:)=u_seri(:,:)+zdu(:,:)
65     v_seri(:,:)=v_seri(:,:)+zdv(:,:)
66     ql_seri(:,:)=ql_seri(:,:)+zdql(:,:)
67
68!======================================================================
69! On ajoute les tendances de la temperature et de la vapeur d'eau
70! en verifiant que ca ne part pas dans les choux
71!======================================================================
72
73      jbad=0
74      jqbad=0
75      DO k = 1, klev
76         DO i = 1, klon
77            zt=t_seri(i,k)+zdt(i,k)
78            zq=q_seri(i,k)+zdq(i,k)
79            IF ( zt>370. .or. zt<130. .or. abs(zdt(i,k))>50. ) then
80            jbad = jbad + 1
81            jadrs(jbad) = i
82            kadrs(jbad) = k
83            ENDIF
84            IF ( zq<0. .or. zq>0.1 .or. abs(zdq(i,k))>1.e-2 ) then
85            jqbad = jqbad + 1
86            jqadrs(jqbad) = i
87            kqadrs(jqbad) = k
88            ENDIF
89            t_seri(i,k)=zt
90            q_seri(i,k)=zq
91         ENDDO
92      ENDDO
93
94!=====================================================================================
95! Impression et stop en cas de probleme important
96!=====================================================================================
97
98IF (jbad .GT. 0) THEN
99      DO j = 1, jbad
100         i=jadrs(j)
101         if(prt_level.ge.debug_level) THEN
102          print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
103          print*,'l    T     dT       Q     dQ    '
104          DO k = 1, klev
105             write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
106          ENDDO
107          call print_debug_phys(i,debug_level,text)
108         endif
109      ENDDO
110ENDIF
111!
112!=====================================================================================
113! Impression, warning et correction en cas de probleme moins important
114!=====================================================================================
115IF (jqbad .GT. 0) THEN
116      done(:) = .false.                         !jyg
117      DO j = 1, jqbad
118        i=jqadrs(j)
119        IF (.NOT.done(i)) THEN                  !jyg
120          if(prt_level.ge.debug_level) THEN
121           print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
122           print*,'l    T     dT       Q     dQ    '
123           DO k = 1, klev
124              write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
125           ENDDO
126          endif
127!jyg<20140228 Corrections pour conservation de l'eau
128          DO k = 1, klev
129            zqp(k) = max(q_seri(i,k),1.e-15)
130          ENDDO
131          zq_int  = 0.
132          zqp_int = 0.
133          DO k = 1, klev
134            zq_int  = zq_int  + q_seri(i,k)*(paprs(i,k)-paprs(i,k+1))/Rg
135            zqp_int = zqp_int + zqp(k)     *(paprs(i,k)-paprs(i,k+1))/Rg
136          ENDDO
137          if(prt_level.ge.debug_level) THEN
138           print*,' cas q_seri<1.e-15 i k zq_int zqp_int zq_int/zqp_int :', &
139                                i, kqadrs(j), zq_int, zqp_int, zq_int/zqp_int
140          endif
141          DO k = 1, klev
142            zq_new = zqp(k)*zq_int/zqp_int
143            zdq(i,k) = zdq(i,k) + zq_new - q_seri(i,k)
144            q_seri(i,k) = zq_new
145          ENDDO
146          done(i) = .true.
147        ENDIF !(.NOT.done(i))
148
149!         DO k = 1, klev
150!           zq=q_seri(i,k)+zdq(i,k)
151!           if (zq.lt.1.e-15) then
152!              if (q_seri(i,k).lt.1.e-15) then
153!               if(prt_level.ge.debug_level) THEN
154!                print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
155!               endif
156!               q_seri(i,k)=1.e-15
157!               zdq(i,k)=(1.e-15-q_seri(i,k))
158!              endif
159!           endif
160!!           zq=q_seri(i,k)+zdq(i,k)
161!!           if (zq.lt.1.e-15) then
162!!              zdq(i,k)=(1.e-15-q_seri(i,k))
163!!           endif
164!         ENDDO
165!jyg>
166      ENDDO
167ENDIF
168!
169
170!IM ajout memes tests pour reverifier les jbad, jqbad beg
171      jbad=0
172      jqbad=0
173      DO k = 1, klev
174         DO i = 1, klon
175            IF ( t_seri(i,k)>370. .or. t_seri(i,k)<130. .or. abs(zdt(i,k))>50. ) then
176            jbad = jbad + 1
177            jadrs(jbad) = i
178!            if(prt_level.ge.debug_level) THEN
179!             print*,'cas2 i k t_seri zdt',i,k,t_seri(i,k),zdt(i,k)
180!            endif
181            ENDIF
182            IF ( q_seri(i,k)<0. .or. q_seri(i,k)>0.1 .or. abs(zdq(i,k))>1.e-2 ) then
183            jqbad = jqbad + 1
184            jqadrs(jqbad) = i
185            kqadrs(jqbad) = k
186!            if(prt_level.ge.debug_level) THEN
187!             print*,'cas2 i k q_seri zdq',i,k,q_seri(i,k),zdq(i,k)
188!            endif
189            ENDIF
190         ENDDO
191      ENDDO
192IF (jbad .GT. 0) THEN
193      DO j = 1, jbad
194         i=jadrs(j)
195         k=kadrs(j)
196         if(prt_level.ge.debug_level) THEN
197          print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
198       &        zdt(i,k),t_seri(i,k)-zdt(i,k)
199!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
200          print*,'l    T     dT       Q     dQ    '
201          DO k = 1, klev
202             write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
203          ENDDO
204          call print_debug_phys(i,debug_level,text)
205         endif
206      ENDDO
207ENDIF
208!
209IF (jqbad .GT. 0) THEN
210      DO j = 1, jqbad
211         i=jqadrs(j)
212         k=kqadrs(j)
213         if(prt_level.ge.debug_level) THEN
214          print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
215       &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
216!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
217          print*,'l    T     dT       Q     dQ    '
218          DO k = 1, klev
219            write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
220          ENDDO
221          call print_debug_phys(i,debug_level,text)
222         endif
223      ENDDO
224ENDIF
225
226      CALL hgardfou(t_seri,ftsol,text)
227      RETURN
228      END
Note: See TracBrowser for help on using the repository browser.