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

Last change on this file since 1998 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: 7.4 KB
RevLine 
[904]1!
[1163]2! $Id: add_phys_tend.F90 1998 2014-03-19 14:14:30Z fhourdin $
[904]3!
[1998]4SUBROUTINE add_phys_tend (zdu,zdv,zdt,zdq,zdql,paprs,text)
[904]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
[941]18use phys_state_var_mod
[904]19IMPLICIT none
[972]20#include "iniprint.h"
[1998]21#include "YOMCST.h"
[904]22
23! Arguments :
24!------------
25REAL zdu(klon,klev),zdv(klon,klev)
26REAL zdt(klon,klev),zdq(klon,klev),zdql(klon,klev)
[1998]27REAL paprs(klon,klev+1)
[904]28CHARACTER*(*) text
29
30! Local :
31!--------
32REAL zt,zq
[1998]33REAL zq_int, zqp_int, zq_new
[904]34
[1998]35REAL zqp(klev)
36
[904]37INTEGER i, k,j
[924]38INTEGER jadrs(klon*klev), jbad
39INTEGER jqadrs(klon*klev), jqbad
[972]40INTEGER kadrs(klon*klev)
41INTEGER kqadrs(klon*klev)
[904]42
[1998]43LOGICAL done(klon)
44
[904]45integer debug_level
[972]46logical, save :: first=.true.
[987]47!$OMP THREADPRIVATE(first)
[972]48INTEGER, SAVE :: itap
[987]49!$OMP THREADPRIVATE(itap)
[904]50!======================================================================
51! Initialisations
52
53debug_level=10
[972]54     if (first) then
55        itap=0
56        first=.false.
57     endif
58! Incrementer le compteur de la physique
59     itap   = itap + 1
[904]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
[972]82            kadrs(jbad) = k
[904]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
[972]87            kqadrs(jqbad) = k
[904]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)
[1047]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
[904]105             write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
[1047]106          ENDDO
107          call print_debug_phys(i,debug_level,text)
108         endif
[904]109      ENDDO
110ENDIF
111!
112!=====================================================================================
113! Impression, warning et correction en cas de probleme moins important
114!=====================================================================================
115IF (jqbad .GT. 0) THEN
[1998]116      done(:) = .false.                         !jyg
[904]117      DO j = 1, jqbad
[1998]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
[904]150!           zq=q_seri(i,k)+zdq(i,k)
151!           if (zq.lt.1.e-15) then
[1998]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
[904]159!           endif
[1998]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>
[904]166      ENDDO
167ENDIF
168!
169
[972]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
[1163]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
[972]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
[1163]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
[972]189            ENDIF
190         ENDDO
191      ENDDO
192IF (jbad .GT. 0) THEN
193      DO j = 1, jbad
194         i=jadrs(j)
195         k=kadrs(j)
[1047]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, &
[972]198       &        zdt(i,k),t_seri(i,k)-zdt(i,k)
[1047]199!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
[972]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
[1163]207ENDIF
[972]208!
209IF (jqbad .GT. 0) THEN
210      DO j = 1, jqbad
211         i=jqadrs(j)
212         k=kqadrs(j)
[1047]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,&
[972]215       &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
[1047]216!!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN
[972]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)
[904]227      RETURN
228      END
Note: See TracBrowser for help on using the repository browser.