Changeset 1998 for LMDZ5/trunk/libf/phylmd/add_phys_tend.F90

Timestamp:

Mar 19, 2014, 3:14:30 PM (10 years ago)

Author:

fhourdin

Message:

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

File:

• LMDZ5/trunk/libf/phylmd/add_phys_tend.F90 (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/add_phys_tend.F90

@@ -2 +2 @@
 ! $Id$
 !
-SUBROUTINE add_phys_tend (zdu,zdv,zdt,zdq,zdql,text)
+SUBROUTINE add_phys_tend (zdu,zdv,zdt,zdq,zdql,paprs,text)
 !======================================================================
 ! Ajoute les tendances des variables physiques aux variables 
@@ -19 +19 @@
 IMPLICIT none
 #include "iniprint.h"
+#include "YOMCST.h"
 
 ! Arguments :
@@ -24 +25 @@
 REAL zdu(klon,klev),zdv(klon,klev)
 REAL zdt(klon,klev),zdq(klon,klev),zdql(klon,klev)
+REAL paprs(klon,klev+1)
 CHARACTER*(*) text
 
@@ -29 +31 @@
 !--------
 REAL zt,zq
+REAL zq_int, zqp_int, zq_new
+
+REAL zqp(klev)
 
 INTEGER i, k,j
@@ -35 +40 @@
 INTEGER kadrs(klon*klev)
 INTEGER kqadrs(klon*klev)
+
+LOGICAL done(klon)
 
 integer debug_level
@@ -107 +114 @@
 !=====================================================================================
 IF (jqbad .GT. 0) THEN
+      done(:) = .false.                         !jyg
       DO j = 1, jqbad
-         i=jqadrs(j)
-         if(prt_level.ge.debug_level) THEN
-          print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
-          print*,'l    T     dT       Q     dQ    '
-         endif
-         DO k = 1, klev
-           zq=q_seri(i,k)+zdq(i,k)
-           if (zq.lt.1.e-15) then
-              if (q_seri(i,k).lt.1.e-15) then
-               if(prt_level.ge.debug_level) THEN
-                print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
-               endif
-               q_seri(i,k)=1.e-15
-               zdq(i,k)=(1.e-15-q_seri(i,k))
-              endif
-           endif
+        i=jqadrs(j)
+        IF (.NOT.done(i)) THEN                  !jyg
+          if(prt_level.ge.debug_level) THEN
+           print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
+           print*,'l    T     dT       Q     dQ    '
+           DO k = 1, klev
+              write(*,'(i3,2f14.4,2e14.2)') k,t_seri(i,k),zdt(i,k),q_seri(i,k),zdq(i,k)
+           ENDDO
+          endif
+!jyg<20140228 Corrections pour conservation de l'eau
+          DO k = 1, klev
+            zqp(k) = max(q_seri(i,k),1.e-15)
+          ENDDO
+          zq_int  = 0.
+          zqp_int = 0.
+          DO k = 1, klev
+            zq_int  = zq_int  + q_seri(i,k)*(paprs(i,k)-paprs(i,k+1))/Rg
+            zqp_int = zqp_int + zqp(k)     *(paprs(i,k)-paprs(i,k+1))/Rg
+          ENDDO
+          if(prt_level.ge.debug_level) THEN
+           print*,' cas q_seri<1.e-15 i k zq_int zqp_int zq_int/zqp_int :', &
+                                i, kqadrs(j), zq_int, zqp_int, zq_int/zqp_int
+          endif
+          DO k = 1, klev
+            zq_new = zqp(k)*zq_int/zqp_int
+            zdq(i,k) = zdq(i,k) + zq_new - q_seri(i,k)
+            q_seri(i,k) = zq_new
+          ENDDO
+          done(i) = .true.
+        ENDIF !(.NOT.done(i))
+
+!         DO k = 1, klev
 !           zq=q_seri(i,k)+zdq(i,k)
 !           if (zq.lt.1.e-15) then
-!              zdq(i,k)=(1.e-15-q_seri(i,k))
+!              if (q_seri(i,k).lt.1.e-15) then
+!               if(prt_level.ge.debug_level) THEN
+!                print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
+!               endif
+!               q_seri(i,k)=1.e-15
+!               zdq(i,k)=(1.e-15-q_seri(i,k))
+!              endif
 !           endif
-         ENDDO
+!!           zq=q_seri(i,k)+zdq(i,k)
+!!           if (zq.lt.1.e-15) then
+!!              zdq(i,k)=(1.e-15-q_seri(i,k))
+!!           endif
+!         ENDDO
+!jyg>
       ENDDO
 ENDIF