Index: /LMDZ5/trunk/libf/dyn3d/advtrac.F90
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/advtrac.F90	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/advtrac.F90	(revision 2286)
@@ -238,5 +238,5 @@
            ! CRisi: on fait passer tout q pour avoir acces aux fils
            
-           write(*,*) 'advtrac 239: iq,q(1721,19,:)=',iq,q(1721,19,:)     
+           !write(*,*) 'advtrac 239: iq,q(1721,19,:)=',iq,q(1721,19,:)     
            call vlsplt(q,2.,massem,wg,pbarug,pbarvg,dtvr,iq)
            
@@ -247,5 +247,5 @@
         else if(iadv(iq).eq.14) then
            !
-           write(*,*) 'advtrac 248: iq,q(1721,19,:)=',iq,q(1721,19,:)
+           !write(*,*) 'advtrac 248: iq,q(1721,19,:)=',iq,q(1721,19,:)
            CALL vlspltqs( q, 2., massem, wg , &
                 pbarug,pbarvg,dtvr,p,pk,teta,iq)
Index: /LMDZ5/trunk/libf/dyn3d/caladvtrac.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/caladvtrac.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/caladvtrac.F	(revision 2286)
@@ -85,7 +85,7 @@
           ENDDO
 
-          write(*,*) 'caladvtrac 87'
+          !write(*,*) 'caladvtrac 87'
           CALL qminimum( q, nqtot, finmasse )
-          write(*,*) 'caladvtrac 89'
+          !write(*,*) 'caladvtrac 89'
 
           CALL SCOPY   ( ip1jmp1*llm, masse, 1, finmasse,       1 )
Index: /LMDZ5/trunk/libf/dyn3d/leapfrog.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/leapfrog.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/leapfrog.F	(revision 2286)
@@ -344,5 +344,5 @@
      *        p, masse, dq,  teta,
      .        flxw, pk)
-          write(*,*) 'caladvtrac 346'
+          !write(*,*) 'caladvtrac 346'
 
          
Index: /LMDZ5/trunk/libf/dyn3d/qminimum.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/qminimum.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/qminimum.F	(revision 2286)
@@ -55,5 +55,5 @@
 
       zx_defau_diag(:,:,:)=0.0
-      q_follow(:,:,:)=q(:,:,:)  
+      q_follow(:,:,1:2)=q(:,:,1:2)  
       DO 1000 k = 1, llm
         DO 1040 i = 1, ip1jmp1
@@ -110,5 +110,5 @@
       ENDIF
 
-      write(*,*) 'qminimum 128'
+      !write(*,*) 'qminimum 128'
       if (ok_isotopes) then
       ! CRisi: traiter de même les traceurs d'eau
@@ -158,5 +158,5 @@
       
         ! 3) transfert d'eau de la vapeur au liquide
-        write(*,*) 'qminimum 164'
+        !write(*,*) 'qminimum 164'
         do k=1,llm
         DO i = 1,ip1jmp1
@@ -186,5 +186,5 @@
 
       endif !if (ok_isotopes) then
-      write(*,*) 'qminimum 188'
+      !write(*,*) 'qminimum 188'
       
 c
Index: /LMDZ5/trunk/libf/dyn3d/vlsplt.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/vlsplt.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/vlsplt.F	(revision 2286)
@@ -454,5 +454,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlsplt 326: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vlsplt 326: iq,nqfils(iq)=',iq,nqfils(iq) 
       
       if (nqdesc(iq).gt.0) then  
@@ -579,5 +579,5 @@
       DATA temps0,temps1,temps2,temps3,temps4,temps5/0.,0.,0.,0.,0.,0./
 
-      write(*,*) 'vly 578: entree, iq=',iq
+      !write(*,*) 'vly 578: entree, iq=',iq
 
       IF(first) THEN
@@ -758,5 +758,5 @@
       ENDDO
 
-      write(*,*) 'vly 756'
+      !write(*,*) 'vly 756'
       DO l=1,llm
        DO ij=1,ip1jm
@@ -776,5 +776,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
    
       if (nqfils(iq).gt.0) then  
@@ -874,5 +874,5 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vly 853: sortie'
+      !write(*,*) 'vly 853: sortie'
 
       RETURN
@@ -933,5 +933,5 @@
 c    sens de W
 
-      write(*,*) 'vlz 923: entree'
+      !write(*,*) 'vlz 923: entree'
 
 #ifdef BIDON
@@ -964,5 +964,5 @@
       ENDDO
 
-      write(*,*) 'vlz 954'
+      !write(*,*) 'vlz 954'
       DO ij=1,ip1jmp1
          dzq(ij,1)=0.
@@ -981,5 +981,5 @@
 c calcul de  - d( q   * w )/ d(sigma)    qu'on ajoute a  dq pour calculer dq
 
-       write(*,*) 'vlz 969'
+       !write(*,*) 'vlz 969'
        DO l = 1,llm-1
          do  ij = 1,ip1jmp1
@@ -1002,5 +1002,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
+      !write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
       if (nqfils(iq).gt.0) then  
        do ifils=1,nqdesc(iq)
@@ -1041,5 +1041,5 @@
         enddo !do ifils=1,nqdesc(iq)
       endif !if (nqfils(iq).gt.0) then
-      write(*,*) 'vlsplt 1032'
+      !write(*,*) 'vlsplt 1032'
 
       RETURN
Index: /LMDZ5/trunk/libf/dyn3d/vlspltqs.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3d/vlspltqs.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3d/vlspltqs.F	(revision 2286)
@@ -175,5 +175,5 @@
       enddo !do ifils=1,nqdesc(iq)  
       endif ! if (nqfils(iq).gt.0) then
-      write(*,*) 'vlspltqs 183: fin de la routine'
+      !write(*,*) 'vlspltqs 183: fin de la routine'
 
       RETURN
@@ -486,5 +486,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlspltqs 326: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vlspltqs 326: iq,nqfils(iq)=',iq,nqfils(iq) 
       
       if (nqfils(iq).gt.0) then  
@@ -797,5 +797,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq) 
    
       if (nqfils(iq).gt.0) then  
@@ -812,5 +812,5 @@
         do ifils=1,nqfils(iq)
          iq2=iqfils(ifils,iq)
-         write(*,*) 'vlyqs 783: appel rec de vly, iq2=',iq2
+         !write(*,*) 'vlyqs 783: appel rec de vly, iq2=',iq2
          call vly(Ratio,pente_max,masseq,qbyv,iq2)
         enddo !do ifils=1,nqfils(iq)
@@ -868,5 +868,5 @@
       ENDDO
 
-      write(*,*) 'vly 866'
+      !write(*,*) 'vly 866'
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -881,5 +881,5 @@
         enddo !do ifils=1,nqdesc(iq)
       endif !if (nqfils(iq).gt.0) then
-      write(*,*) 'vly 879'
+      !write(*,*) 'vly 879'
 
       RETURN
Index: /LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/advtrac_loc.F	(revision 2286)
@@ -155,10 +155,10 @@
 c$OMP BARRIER
                  
-          write(*,*) 'advtrac 157: appel de vlspltgen_loc'
+          !write(*,*) 'advtrac 157: appel de vlspltgen_loc'
           call vlspltgen_loc( q,iadv, 2., massem, wg ,
      *                        pbarug,pbarvg,dtvr,p,
      *                        pk,teta )
 
-          write(*,*) 'advtrac 162: apres appel vlspltgen_loc'
+          !write(*,*) 'advtrac 162: apres appel vlspltgen_loc'
       if (ok_iso_verif) then
            call check_isotopes(q,ijb_u,ije_u,'advtrac 162')
Index: /LMDZ5/trunk/libf/dyn3dmem/caladvtrac_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/caladvtrac_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/caladvtrac_loc.F	(revision 2286)
@@ -56,5 +56,5 @@
 !$OMP THREADPRIVATE(Request_vanleer)
 
-      write(*,*) 'caladvtrac 58: entree'     
+      !write(*,*) 'caladvtrac 58: entree'     
       ijbu=ij_begin
       ijeu=ij_end
@@ -109,5 +109,5 @@
 
       IF ( iadvtr.EQ.iapp_tracvl ) THEN
-      write(*,*) 'caladvtrac 133'
+      !write(*,*) 'caladvtrac 133'
 c$OMP MASTER
         call suspend_timer(timer_caldyn)
@@ -184,9 +184,9 @@
          CALL WriteField_u('wg1',wg_adv)
 #endif        
-      write(*,*) 'caladvtrac 185'  
+      !write(*,*) 'caladvtrac 185'  
       CALL advtrac_loc( pbarug_adv,pbarvg_adv,wg_adv, 
      *             p_adv,  massem_adv,q_adv, teta_adv,
      .             pk_adv)     
-      write(*,*) 'caladvtrac 189'
+      !write(*,*) 'caladvtrac 189'
 
 
Index: /LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F	(revision 2286)
@@ -164,5 +164,5 @@
 c$OMP END MASTER
 c$OMP BARRIER
-        write(*,*) 'integrd 170'
+        !write(*,*) 'integrd 170'
       IF (.NOT. Checksum_all) THEN
         call WriteField_v('int_vcov',vcov)
@@ -194,5 +194,5 @@
        
 c
-        write(*,*) 'integrd 200'
+        !write(*,*) 'integrd 200'
 C$OMP MASTER
       if (pole_nord) THEN
@@ -285,5 +285,5 @@
 c
 c
-        write(*,*) 'integrd 291'
+        !write(*,*) 'integrd 291'
       IF (pole_nord) THEN
        
@@ -352,7 +352,7 @@
         endif !if (ok_iso_verif) then
 
-        write(*,*) 'integrd 341'
+        !write(*,*) 'integrd 341'
         CALL qminimum_loc( q, nq, deltap )
-        write(*,*) 'integrd 343'
+        !write(*,*) 'integrd 343'
 
         if (ok_iso_verif) then
Index: /LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F	(revision 2286)
@@ -308,5 +308,5 @@
 
    1  CONTINUE ! Matsuno Forward step begins here
-      write(*,*) 'leapfrog 298: itau=',itau
+      !write(*,*) 'leapfrog 298: itau=',itau
       jD_cur = jD_ref + day_ini - day_ref +                             &
      &          itau/day_step
@@ -713,10 +713,10 @@
       IF( forward. OR . leapf )  THEN
 ! Ehouarn: NB: fields sent to advtrac are those at the beginning of the time step
-        write(*,*) 'leapfrog 679: avant CALL caladvtrac_loc'
+        !write(*,*) 'leapfrog 679: avant CALL caladvtrac_loc'
          CALL caladvtrac_loc(q,pbaru,pbarv,
      *        p, masse, dq,  teta,
      .        flxw,pk, iapptrac)
 
-         write(*,*) 'leapfrog 719'
+         !write(*,*) 'leapfrog 719'
          if (ok_iso_verif) then
            call check_isotopes(q,ijb_u,ije_u,
@@ -755,5 +755,5 @@
 !       CALL FTRACE_REGION_BEGIN("integrd")
 
-       write(*,*) 'leapfrog 720' 
+       !write(*,*) 'leapfrog 720' 
         if (ok_iso_verif) then
            call check_isotopes(q,ijb_u,ije_u,'leapfrog 756')
@@ -765,5 +765,5 @@
 !     $              finvmaold                                    )
 
-       write(*,*) 'leapfrog 724'       
+       !write(*,*) 'leapfrog 724'       
         if (ok_iso_verif) then
            call check_isotopes(q,ijb_u,ije_u,'leapfrog 762')
@@ -1148,5 +1148,5 @@
            call check_isotopes(q,ijb_u,ije_u,'leapfrog 1132')
         endif !if (ok_iso_verif) then
-        write(*,*) 'leapfrog 1134: iflag_phys=',iflag_phys
+        !write(*,*) 'leapfrog 1134: iflag_phys=',iflag_phys
 
       IF(iflag_phys.EQ.2) THEN ! "Newtonian" case
@@ -1221,5 +1221,5 @@
 c   dissipation horizontale et verticale  des petites echelles:
 c   ----------------------------------------------------------
-      write(*,*) 'leapfrog 1163: apdiss=',apdiss
+      !write(*,*) 'leapfrog 1163: apdiss=',apdiss
       IF(apdiss) THEN
       
Index: /LMDZ5/trunk/libf/dyn3dmem/qminimum_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/qminimum_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/qminimum_loc.F	(revision 2286)
@@ -50,5 +50,5 @@
 c
 
-        write(*,*) 'qminimum 52: entree'
+        !write(*,*) 'qminimum 52: entree'
         if (ok_iso_verif) then
            call check_isotopes(q,ij_begin,ij_end,'qminimum 52')   
Index: /LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/vlsplt_loc.F	(revision 2286)
@@ -56,6 +56,6 @@
       INTEGER ijb,ije,ijb_x,ije_x
       
-      write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
-     &   iq,ijb_x
+      !write(*,*) 'vlsplt 58: entree dans vlx_loc, iq,ijb_x=',
+!     &   iq,ijb_x
 c   calcul de la pente a droite et a gauche de la maille
 
@@ -154,5 +154,5 @@
       ENDIF ! (pente_max.lt.-1.e-5)
 
-      write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlx 156: iq,ijb_x=',iq,ijb_x
 
 c   bouclage de la pente en iip1:
@@ -273,5 +273,5 @@
                ENDDO
                niju=iju
-               PRINT*,'vlx 278, niju,nl',niju,nl(l)
+               !PRINT*,'vlx 278, niju,nl',niju,nl(l)
 
 c  traitement des mailles
@@ -331,5 +331,5 @@
 ! Il faut faire ça avant d'avoir mis à jour q et masse
 
-      write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq) 
+      !write(*,*) 'vlsplt 326: iq,ijb_x,nqfils(iq)=',iq,ijb_x,nqfils(iq) 
 
       if (nqfils(iq).gt.0) then  
@@ -354,5 +354,5 @@
 ! end CRisi
 
-      write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 360: iq,ijb_x=',iq,ijb_x
 
 c   calcul des tENDances
@@ -373,5 +373,5 @@
       ENDDO
 c$OMP END DO NOWAIT
-      write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 380: iq,ijb_x=',iq,ijb_x
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -394,5 +394,5 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlsplt 399: iq,ijb_x=',iq,ijb_x
 c     CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1)
 c     CALL SCOPY((jjm-1)*llm,masse(iip1+iip1,1),iip1,masse(iip2,1),iip1)
@@ -726,5 +726,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
+      !write(*,*) 'vly 689: iq,nqfils(iq)=',iq,nqfils(iq) 
 
       ijb=ij_begin-2*iip1
@@ -926,5 +926,5 @@
 c    sens de W
 
-      write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
+      !write(*,*) 'vlsplt 926: entree dans vlz_loc, iq=',iq
 #ifdef BIDON
       IF(testcpu) THEN
@@ -982,5 +982,5 @@
 c calcul de  - d( q   * w )/ d(sigma)    qu'on ajoute a  dq pour calculer dq
 
-       write(*,*) 'vlz 982,ijb,ije=',ijb,ije
+       !write(*,*) 'vlz 982,ijb,ije=',ijb,ije
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
        DO l = 1,llm-1
@@ -998,5 +998,5 @@
        ENDDO
 c$OMP END DO NOWAIT    
-       write(*,*) 'vlz 1001'    
+       !write(*,*) 'vlz 1001'    
 
 c$OMP MASTER
@@ -1010,5 +1010,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
+      !write(*,*) 'vlsplt 942: iq,nqfils(iq)=',iq,nqfils(iq)
       if (nqfils(iq).gt.0) then  
        do ifils=1,nqdesc(iq)
Index: /LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/vlspltgen_loc.F	(revision 2286)
@@ -202,5 +202,5 @@
 !      DO iq=1,nqtot
       DO iq=1,nqperes ! CRisi: on ne boucle que sur les pères= ceux qui sont transportés directement par l'air
-       write(*,*) 'vlspltgen 192: iq,iadv=',iq,iadv(iq)
+       !write(*,*) 'vlspltgen 192: iq,iadv=',iq,iadv(iq)
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -298,5 +298,5 @@
 
       do iq=1,nqperes
-        write(*,*) 'vlspltgen 279: iq=',iq
+        !write(*,*) 'vlspltgen 279: iq=',iq
 
         if(iadv(iq) == 0) then
@@ -351,5 +351,5 @@
 
       do iq=1,nqperes
-       write(*,*) 'vlspltgen 321: iq=',iq
+       !write(*,*) 'vlspltgen 321: iq=',iq
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -383,5 +383,5 @@
 
       do iq=1,nqperes
-      write(*,*) 'vlspltgen 349: iq=',iq
+      !write(*,*) 'vlspltgen 349: iq=',iq
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -448,5 +448,5 @@
 c$OMP BARRIER
       do iq=1,nqperes
-      write(*,*) 'vlspltgen 409: iq=',iq
+      !write(*,*) 'vlspltgen 409: iq=',iq
 
         if(iadv(iq) == 0) then
@@ -470,8 +470,8 @@
       
       enddo
-      write(*,*) 'vlspltgen_loc 476'
-
-c$OMP BARRIER
-      write(*,*) 'vlspltgen_loc 477'
+      !write(*,*) 'vlspltgen_loc 476'
+
+c$OMP BARRIER
+      !write(*,*) 'vlspltgen_loc 477'
 c$OMP MASTER
       call VTb(VTHallo)
@@ -489,5 +489,5 @@
 
 
-      write(*,*) 'vlspltgen_loc 494'
+      !write(*,*) 'vlspltgen_loc 494'
       if (ok_iso_verif) then
            call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 461')
@@ -495,5 +495,5 @@
 
       do iq=1,nqperes
-      write(*,*) 'vlspltgen 449: iq=',iq
+      !write(*,*) 'vlspltgen 449: iq=',iq
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -526,5 +526,5 @@
 
       do iq=1,nqperes
-      write(*,*) 'vlspltgen 477: iq=',iq
+      !write(*,*) 'vlspltgen 477: iq=',iq
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -553,5 +553,5 @@
        enddo !do iq=1,nqperes
 
-      write(*,*) 'vlspltgen 550: apres derniere serie de call vlx'
+      !write(*,*) 'vlspltgen 550: apres derniere serie de call vlx'
       if (ok_iso_verif) then
            call check_isotopes(zq,ij_begin,ij_end,'vlspltgen_loc 521')
@@ -560,10 +560,10 @@
       ijb=ij_begin
       ije=ij_end
-      write(*,*) 'vlspltgen_loc 557'
+      !write(*,*) 'vlspltgen_loc 557'
 c$OMP BARRIER      
 
-      write(*,*) 'vlspltgen_loc 559'  
+      !write(*,*) 'vlspltgen_loc 559'  
       DO iq=1,nqtot
-       write(*,*) 'vlspltgen_loc 561, iq=',iq  
+       !write(*,*) 'vlspltgen_loc 561, iq=',iq  
 #ifdef DEBUG_IO    
        CALL WriteField_u('zq',zq(:,:,iq))
@@ -579,5 +579,5 @@
         ENDDO
 c$OMP END DO NOWAIT    
-      write(*,*) 'vlspltgen_loc 575'      
+      !write(*,*) 'vlspltgen_loc 575'      
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
@@ -588,5 +588,5 @@
         ENDDO
 c$OMP END DO NOWAIT  
-      write(*,*) 'vlspltgen_loc 583'  
+      !write(*,*) 'vlspltgen_loc 583'  
       ENDDO !DO iq=1,nqtot
         
@@ -603,5 +603,5 @@
 cc$OMP BARRIER
 
-      write(*,*) 'vlspltgen 597: sortie'  
+      !write(*,*) 'vlspltgen 597: sortie'  
       RETURN
       END
Index: /LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F
===================================================================
--- /LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F	(revision 2285)
+++ /LMDZ5/trunk/libf/dyn3dmem/vlspltqs_loc.F	(revision 2286)
@@ -47,6 +47,6 @@
       INTEGER ijb,ije,ijb_x,ije_x
       
-      write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
-     &   iq,ijb_x
+      !write(*,*) 'vlspltqs 58: entree vlxqs_loc, iq,ijb_x=',
+!     &   iq,ijb_x
 
 c   calcul de la pente a droite et a gauche de la maille
@@ -280,5 +280,5 @@
                ENDDO
                niju=iju
-               PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
+               !PRINT*,'vlxqs 280: niju,nl',niju,nl(l)
 
 c  traitement des mailles
@@ -336,6 +336,6 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=',
-     &     iq,ijb_x,nqfils(iq)  
+      !write(*,*) 'vlspltqs 336: iq,ijb_x,nqfils(iq)=',
+!     &     iq,ijb_x,nqfils(iq)  
 
       if (nqfils(iq).gt.0) then  
@@ -353,5 +353,5 @@
         do ifils=1,nqfils(iq)
          iq2=iqfils(ifils,iq)
-         write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
+         !write(*,*) 'vlxqs 349: on appelle vlx pour iq2=',iq2
          call vlx_loc(Ratio,pente_max,masse,u_mq,ijb_x,ije_x,iq2)
         enddo !do ifils=1,nqfils(iq)
@@ -359,5 +359,5 @@
 ! end CRisi
 
-      write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x   
+      !write(*,*) 'vlspltqs 360: iq,ijb_x=',iq,ijb_x   
 
 c   calcul des tendances
@@ -379,5 +379,5 @@
 c$OMP END DO NOWAIT
 
-      write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlspltqs 380: iq,ijb_x=',iq,ijb_x
 
 ! retablir les fils en rapport de melange par rapport a l'air:
@@ -398,5 +398,5 @@
       endif !if (nqfils(iq).gt.0) then
 
-      write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
+      !write(*,*) 'vlspltqs 399: iq,ijb_x=',iq,ijb_x
 
 c     CALL SCOPY((jjm-1)*llm,q(iip1+iip1,1),iip1,q(iip2,1),iip1)
@@ -477,6 +477,6 @@
       l=3
       if ((ij.ge.ijb).and.(ij.le.ije)) then
-        write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
-     &             ij,l,iq,ijb,q(ij,l,:)
+        !write(*,*) 'vlyqs 480: ij,l,iq,ijb,q(ij,l,:)=',
+!     &             ij,l,iq,ijb,q(ij,l,:)
       endif  
 
@@ -725,5 +725,5 @@
 ! CRisi: appel récursif de l'advection sur les fils.
 ! Il faut faire ça avant d'avoir mis à jour q et masse
-      write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq)
+      !write(*,*) 'vlyqs 689: iq,nqfils(iq)=',iq,nqfils(iq)
       
       ijb=ij_begin-2*iip1