Index: LMDZ5/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2606)
+++ LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2611)
@@ -216,4 +216,6 @@
 
     USE paramLMDZ_phy_mod
+
+    USE cmp_seri_mod
 
     !IM stations CFMIP
@@ -801,7 +803,16 @@
     ! eva: evaporation de l'eau liquide nuageuse
     ! vdf: couche limite (Vertical DiFfusion)
-
+    !
     ! tendance nulles
     REAL, dimension(klon,klev):: du0, dv0, dt0, dq0, dql0, dqi0
+    !
+    ! Flag pour pouvoir ne pas ajouter les tendances. 
+    ! Par defaut, les tendances doivente etre ajoutees et
+    ! flag_inhib_tend = 0
+    ! flag_inhib_tend > 0 : tendances non ajoutees, avec un nombre 
+    ! croissant de print quand la valeur du flag augmente
+    !!! attention, ce flag doit etre change avec prudence !!!
+    INTEGER :: flag_inhib_tend = 0 !  0 is the default value
+!!    INTEGER :: flag_inhib_tend = 2
 
     !
@@ -1772,4 +1783,6 @@
        ENDDO
     ENDDO
+    ! Initialize variables used for diagnostic purpose
+    if (flag_inhib_tend .ne. 0) call init_cmp_seri
     !IM
     IF (ip_ebil_phy.ge.1) THEN 
@@ -2217,9 +2230,9 @@
           CALL add_pbl_tend &
                (d_u_vdf,d_v_vdf,d_t_vdf+d_t_diss,d_q_vdf,dql0,dqi0,paprs,&
-               'vdf',abortphy)
+               'vdf',abortphy,flag_inhib_tend)
        ELSE
           CALL add_phys_tend &
                (d_u_vdf,d_v_vdf,d_t_vdf+d_t_diss,d_q_vdf,dql0,dqi0,paprs,&
-               'vdf',abortphy)
+               'vdf',abortphy,flag_inhib_tend)
        ENDIF
        !--------------------------------------------------------------------
@@ -2599,5 +2612,5 @@
 
     CALL add_phys_tend(d_u_con, d_v_con, d_t_con, d_q_con, dql0, dqi0, paprs, &
-         'convection',abortphy)
+         'convection',abortphy,flag_inhib_tend)
 
     !-------------------------------------------------------------------------
@@ -2772,5 +2785,5 @@
        d_q_wake(:,:)=dq_wake(:,:)*dtime
        CALL add_phys_tend(du0,dv0,d_t_wake,d_q_wake,dql0,dqi0,paprs,'wake', &
-            abortphy)
+            abortphy,flag_inhib_tend)
        !------------------------------------------------------------------------
 
@@ -2919,5 +2932,5 @@
           !
           CALL add_phys_tend(d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs,  &
-                             dql0,dqi0,paprs,'thermals', abortphy)
+                             dql0,dqi0,paprs,'thermals', abortphy,flag_inhib_tend)
           !
           !>jyg
@@ -2983,5 +2996,5 @@
           ! ajout des tendances de l'ajustement sec ou des thermiques
           CALL add_phys_tend(du0,dv0,d_t_ajsb,d_q_ajsb,dql0,dqi0,paprs, &
-               'ajsb',abortphy)
+               'ajsb',abortphy,flag_inhib_tend)
           d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_ajsb(:,:)
           d_q_ajs(:,:)=d_q_ajs(:,:)+d_q_ajsb(:,:)
@@ -3042,5 +3055,5 @@
 
     CALL add_phys_tend(du0,dv0,d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc,paprs, &
-         'lsc',abortphy)
+         'lsc',abortphy,flag_inhib_tend)
     !---------------------------------------------------------------------------
     DO k = 1, klev
@@ -3804,6 +3817,6 @@
     ENDDO
 
-    CALL add_phys_tend(du0,dv0,d_t_swr,dq0,dql0,dqi0,paprs,'SW',abortphy)
-    CALL add_phys_tend(du0,dv0,d_t_lwr,dq0,dql0,dqi0,paprs,'LW',abortphy)
+    CALL add_phys_tend(du0,dv0,d_t_swr,dq0,dql0,dqi0,paprs,'SW',abortphy,flag_inhib_tend)
+    CALL add_phys_tend(du0,dv0,d_t_lwr,dq0,dql0,dqi0,paprs,'LW',abortphy,flag_inhib_tend)
 
     !
@@ -3889,5 +3902,5 @@
        ! ajout des tendances de la trainee de l'orographie
        CALL add_phys_tend(d_u_oro,d_v_oro,d_t_oro,dq0,dql0,dqi0,paprs,'oro', &
-            abortphy)
+            abortphy,flag_inhib_tend)
        !----------------------------------------------------------------------
        !
@@ -3935,5 +3948,5 @@
        ! ajout des tendances de la portance de l'orographie
        CALL add_phys_tend(d_u_lif, d_v_lif, d_t_lif, dq0, dql0, dqi0, paprs, &
-            'lif', abortphy)
+            'lif', abortphy,flag_inhib_tend)
     ENDIF ! fin de test sur ok_orolf
 
@@ -3958,5 +3971,5 @@
        d_t_hin(:, :)=0.
        CALL add_phys_tend(du_gwd_hines, dv_gwd_hines, d_t_hin, dq0, dql0, &
-            dqi0, paprs, 'hin', abortphy)
+            dqi0, paprs, 'hin', abortphy,flag_inhib_tend)
     ENDIF
 
@@ -3975,5 +3988,5 @@
 
        CALL add_phys_tend(du_gwd_front, dv_gwd_front, dt0, dq0, dql0, dqi0, &
-            paprs, 'front_gwd_rando', abortphy)
+            paprs, 'front_gwd_rando', abortphy,flag_inhib_tend)
     ENDIF
 
@@ -3983,5 +3996,5 @@
             du_gwd_rando, dv_gwd_rando, east_gwstress, west_gwstress)
        CALL add_phys_tend(du_gwd_rando, dv_gwd_rando, dt0, dq0, dql0, dqi0, &
-            paprs, 'flott_gwd_rando', abortphy)
+            paprs, 'flott_gwd_rando', abortphy,flag_inhib_tend)
        zustr_gwd_rando=0.
        zvstr_gwd_rando=0.
@@ -4046,5 +4059,5 @@
        ! ajout de la tendance d'humidite due au methane
        CALL add_phys_tend(du0, dv0, dt0, d_q_ch4*dtime, dql0, dqi0, paprs, &
-            'q_ch4', abortphy)
+            'q_ch4', abortphy,flag_inhib_tend)
     END IF
     !
