Changeset 972 for LMDZ4/trunk/libf/phylmd/add_phys_tend.F90

Timestamp:

Jun 19, 2008, 12:24:22 PM (16 years ago)

Author:

lmdzadmin

Message:

Version thermique FH/CRio
Ajout tests cas physiques non pris en comptes et ajout/enleve prints
Nouvelle routine thermcell_flux2.F90
IM

File:

• LMDZ4/trunk/libf/phylmd/add_phys_tend.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/add_phys_tend.F90

@@ -18 +18 @@
 use phys_state_var_mod
 IMPLICIT none
+#include "iniprint.h"
 
 ! Arguments :
@@ -32 +33 @@
 INTEGER jadrs(klon*klev), jbad
 INTEGER jqadrs(klon*klev), jqbad
+INTEGER kadrs(klon*klev)
+INTEGER kqadrs(klon*klev)
 
 integer debug_level
-
+logical, save :: first=.true.
+INTEGER, SAVE :: itap
 !======================================================================
 ! Initialisations
 
 debug_level=10
-
+     if (first) then
+        itap=0
+        first=.false.
+     endif
+! Incrementer le compteur de la physique
+     itap   = itap + 1
 !======================================================================
 ! Ajout des tendances sur le vent et l'eau liquide
@@ -62 +71 @@
             jbad = jbad + 1
             jadrs(jbad) = i
+            kadrs(jbad) = k
             ENDIF
             IF ( zq<0. .or. zq>0.1 .or. abs(zdq(i,k))>1.e-2 ) then
             jqbad = jqbad + 1
             jqadrs(jqbad) = i
+            kqadrs(jqbad) = k
             ENDIF
             t_seri(i,k)=zt
@@ -97 +108 @@
          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)
+           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
+!              print*,' cas q_seri<1.e-15 i k q_seri zq zdq :',i,k,q_seri(i,k),zq,zdq(i,k)
+               q_seri(i,k)=1.e-15
+               zdq(i,k)=(1.e-15-q_seri(i,k))
+              endif
+           endif
 !           zq=q_seri(i,k)+zdq(i,k)
 !           if (zq.lt.1.e-15) then
@@ -103 +121 @@
 !           endif
          ENDDO
-         call print_debug_phys(i,debug_level,text)
       ENDDO
 ENDIF
 !
 
+!IM ajout memes tests pour reverifier les jbad, jqbad beg
+      jbad=0
+      jqbad=0
+      DO k = 1, klev
+         DO i = 1, klon
+            IF ( t_seri(i,k)>370. .or. t_seri(i,k)<130. .or. abs(zdt(i,k))>50. ) then
+            jbad = jbad + 1
+            jadrs(jbad) = i
+!           if(prt_level.ge.10) THEN 
+!             print*,'cas2 i k t_seri zdt',i,k,t_seri(i,k),zdt(i,k)
+!           endif
+            ENDIF
+            IF ( q_seri(i,k)<0. .or. q_seri(i,k)>0.1 .or. abs(zdq(i,k))>1.e-2 ) then
+            jqbad = jqbad + 1
+            jqadrs(jqbad) = i
+            kqadrs(jqbad) = k
+!           if(prt_level.ge.10) THEN 
+!             print*,'cas2 i k q_seri zdq',i,k,q_seri(i,k),zdq(i,k)
+!           endif
+            ENDIF
+         ENDDO
+      ENDDO
+IF (jbad .GT. 0) THEN
+      DO j = 1, jbad
+         i=jadrs(j)
+         k=kadrs(j)
+         print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
+       &        zdt(i,k),t_seri(i,k)-zdt(i,k)
+!        if(prt_level.ge.10) THEN 
+         if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 
+          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
+          call print_debug_phys(i,debug_level,text)
+         endif
+      ENDDO
+ENDIF
+!
+IF (jqbad .GT. 0) THEN
+      DO j = 1, jqbad
+         i=jqadrs(j)
+         k=kqadrs(j)
+         print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
+       &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
+!        if(prt_level.ge.10) THEN 
+         if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 
+          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
+          call print_debug_phys(i,debug_level,text)
+         endif
+      ENDDO
+ENDIF
+
+      CALL hgardfou(t_seri,ftsol,text)
       RETURN
       END