Changeset 4033 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Nov 30, 2021, 2:00:13 PM (3 years ago)

Author:

crisi

Message:

nettoyage de sorties ecran

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -2421 +2421 @@
        do itr=1,nqtottr
          iq=itr_indice(itr)
-#ifdef ISOVERIF
-         write(*,*) 'physiq 1973: itr,iq=',itr,iq
-         write(*,*) 'qx(1,1,iq)=',qx(1,1,iq)
-#endif
+!#ifdef ISOVERIF
+!         write(*,*) 'physiq 1973: itr,iq=',itr,iq
+!         write(*,*) 'qx(1,1,iq)=',qx(1,1,iq)
+!#endif
          DO  k = 1, klev
              DO  i = 1, klon
@@ -3523 +3523 @@
      &           'physiq 1456, avant concvl')
             endif   
-#endif      
+#endif     
       enddo !do k=1,nlev
-   enddo  !do i=1,klon
-        if (iso_eau.gt.0) then
-        i=1
-        k=1
-        write(*,*) 'physic 2376: xt_seri(iso_eau,i,k),q_seri(i,k)=',xt_seri(iso_eau,i,k),q_seri(i,k)
-        write(*,*) 'xt_seri(:,i,k)=',xt_seri(:,i,k)
-        write(*,*) 'physic 2376: xt_x(iso_eau,i,k),q_x(i,k)=',xt_x(iso_eau,i,k),q_x(i,k)
-        write(*,*) 'xt_x(:,i,k)=',xt_x(:,i,k)
-        endif
-#endif 
+   enddo  !do i=1,klon   
+#endif
 !ISOVERIF         
           if ((bidouille_anti_divergence).and. &
@@ -3610 +3602 @@
 #ifdef ISOVERIF
 !          write(*,*) 'q_detrainement(1,:)=',q_detrainement(1,:)
-          call iso_verif_noNaN_vect2D(d_xt_con, &
+      call iso_verif_noNaN_vect2D(d_xt_con, &
      &           'physiq 3203a apres conv',ntraciso,klon,klev) 
-          call iso_verif_noNaN_vect2D(xt_seri, &
-     &           'physiq 3203b apres conv',ntraciso,klon,klev)          
+      call iso_verif_noNaN_vect2D(xt_seri, &
+     &           'physiq 3203b apres conv',ntraciso,klon,klev) 
+      if (iso_HDO.gt.0) then
+      call iso_verif_aberrant_enc_vect2D( &
+     &           xt_seri,q_seri, &
+     &           'physiq 3619a',ntraciso,klon,klev)
+      endif
+      if ((iso_HDO.gt.0).and.(iso_O18.gt.0)) then
+      call iso_verif_O18_aberrant_enc_vect2D( &
+     &           xt_seri,q_seri, &
+     &           'physiq 3619b',ntraciso,klon,klev)
+      endif         
 #endif
 #ifdef ISOVERIF
@@ -3906 +3908 @@
 #ifdef ISO
 #ifdef ISOVERIF
-     write(*,*) 'physiq 3425'
+     write(*,*) 'physiq 3425: apres convection'
       if (iso_HDO.gt.0) then
       call iso_verif_aberrant_enc_vect2D( &
@@ -3915 +3917 @@
       call iso_verif_O18_aberrant_enc_vect2D( &
      &           xt_seri,q_seri, &
-     &           'physiq 3691a',ntraciso,klon,klev)
+     &           'physiq 3691b',ntraciso,klon,klev)
       endif
 #ifdef ISOTRAC