Index: LMDZ6/trunk/libf/phylmd/wake.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/wake.F90	(revision 4294)
+++ LMDZ6/trunk/libf/phylmd/wake.F90	(revision 4295)
@@ -291,5 +291,4 @@
   REAL, DIMENSION(klon)                                 :: awdens_in, wdens_in   ! pour les prints
 
-   write(81) klon,klev,znatsurf,p,pi,ph,omgb,dtime,tenv0,qe0,dtdwn,dqdwn,amdwn,amup,dta,dqa,wgen,sigd_con,cin,deltatw,deltaqw,sigmaw,awdens,wdens
   ! -------------------------------------------------------------------------
   ! Initialisations
@@ -319,5 +318,5 @@
   ! alpk = 0.5
   ! alpk = 0.05
-print *,'XXXX dtime input ', dtime
+!print *,'XXXX dtime input ', dtime
  igout = klon/2+1/klon
 
@@ -451,5 +450,5 @@
     d_sig_bnd2(i) = sigmaw_targ - sigmaw(i)
     d_sigmaw2(i) = sigmaw_targ - sigmaw(i)
-  print *,'XXXX1 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX1 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
     sigmaw(i) = sigmaw_targ
 !>jyg
@@ -836,5 +835,5 @@
       d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i)
       d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
-  print *,'XXXX2 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX2 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
       sigmaw(i) = sigmaw_targ
 !>jyg
@@ -950,5 +949,5 @@
         d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i)
         d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
-  print *,'XXXX3 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX3 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
         sigmaw(i) = sigmaw_targ
 !>jyg
@@ -1496,5 +1495,5 @@
         sigmaw(i) = sigmaw(i) + d_sigmaw(i)
         d_sigmaw2(i) = d_sigmaw2(i) + d_sigmaw(i)
-  print *,'XXXX4 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX4 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
       END IF
     END DO
@@ -1518,5 +1517,5 @@
           d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i)
           d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
-  print *,'XXXX5 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX5 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
           sigmaw(i) = sigmaw_targ
         END IF
@@ -1822,5 +1821,5 @@
           d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i)
           d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
-  print *,'XXXX6 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX6 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
           sigmaw(i) = sigmaw_targ
 !>jyg
@@ -2041,5 +2040,5 @@
       d_sig_bnd2(i) = d_sig_bnd2(i) + sigmaw_targ - sigmaw(i)
       d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
-  print *,'XXXX7 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX7 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
       sigmaw(i) = sigmaw_targ
 !>jyg
@@ -2141,5 +2140,5 @@
 !!      d_sigmaw2(i) = d_sigmaw2(i) + sigmaw_targ - sigmaw(i)
       d_sigmaw2(i) = sigmaw_targ - sigmaw_in(i)      ! _in = correction jyg 20220124
-  print *,'XXXX8 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
+!  print *,'XXXX8 d_sigmaw2(i), sigmaw(i) ', d_sigmaw2(i), sigmaw(i)
       sigmaw(i) = sigmaw_targ
       IF (iflag_wk_pop_dyn >= 1) THEN
@@ -2209,5 +2208,5 @@
     d_sig_bnd2(i) = d_sig_bnd2(i)/dtime
     d_sigmaw2(i) = d_sigmaw2(i)/dtime
-  print *,'XXXX9 d_sigmaw2(i), sigmaw(i), dtime ', d_sigmaw2(i), sigmaw(i), dtime
+!  print *,'XXXX9 d_sigmaw2(i), sigmaw(i), dtime ', d_sigmaw2(i), sigmaw(i), dtime
 !
     d_dens_gen2(i) = d_dens_gen2(i)/dtime