Index: /trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F
===================================================================
--- /trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F	(revision 2306)
+++ /trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F	(revision 2307)
@@ -139,5 +139,6 @@
       END
       RECURSIVE SUBROUTINE vlx(q,pente_max,masse,u_m,iq)
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
 
 c     Auteurs:   P.Le Van, F.Hourdin, F.Forget 
@@ -459,6 +460,6 @@
            !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
@@ -481,5 +482,5 @@
          DO ij=iip2+1,ip1jm
             !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),1e-16)
+            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),masseqmin)
             q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+
      &      u_mq(ij-1,l)-u_mq(ij,l))
@@ -518,5 +519,6 @@
       END
       RECURSIVE SUBROUTINE vly(q,pente_max,masse,masse_adv_v,iq)
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
       USE comconst_mod, ONLY: pi
 c
@@ -788,6 +790,6 @@
            !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
@@ -886,5 +888,6 @@
       END
       RECURSIVE SUBROUTINE vlz(q,pente_max,masse,w,iq)
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi 
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
 c
 c     Auteurs:   P.Le Van, F.Hourdin, F.Forget 
@@ -1015,6 +1018,6 @@
            !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masseq(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
Index: /trunk/LMDZ.COMMON/libf/dyn3d_common/infotrac.F90
===================================================================
--- /trunk/LMDZ.COMMON/libf/dyn3d_common/infotrac.F90	(revision 2306)
+++ /trunk/LMDZ.COMMON/libf/dyn3d_common/infotrac.F90	(revision 2307)
@@ -32,4 +32,7 @@
   INTEGER, ALLOCATABLE, DIMENSION(:,:), SAVE    :: iqfils
   INTEGER, ALLOCATABLE, DIMENSION(:), SAVE    :: iqpere
+  REAL :: qperemin,masseqmin ! MVals: thresholds for zq(pere) and masseq in the transport of the isotopic Ratio (vlsplt_p.F) 
+  PARAMETER (qperemin=1.e-16)
+  PARAMETER (masseqmin=1.e-16)
 
 ! conv_flg(it)=0 : convection desactivated for tracer number it 
Index: /trunk/LMDZ.COMMON/libf/dyn3dpar/vlsplt_p.F
===================================================================
--- /trunk/LMDZ.COMMON/libf/dyn3dpar/vlsplt_p.F	(revision 2306)
+++ /trunk/LMDZ.COMMON/libf/dyn3dpar/vlsplt_p.F	(revision 2307)
@@ -214,5 +214,6 @@
 c   --------------------------------------------------------------------
       USE Parallel_lmdz
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
       IMPLICIT NONE
 c
@@ -528,6 +529,6 @@
            ! les calcule donc que de ijb à ije
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
@@ -549,5 +550,5 @@
          DO ij=ijb+1,ije
             !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),1e-16)
+            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),masseqmin)
             q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+
      &      u_mq(ij-1,l)-u_mq(ij,l))
@@ -601,5 +602,6 @@
 c   --------------------------------------------------------------------
       USE parallel_lmdz
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
       USE comconst_mod, ONLY: pi
 
@@ -918,6 +920,6 @@
            !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
@@ -1049,5 +1051,6 @@
 c   --------------------------------------------------------------------
       USE Parallel_lmdz
-      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils ! CRisi
+      USE infotrac, ONLY : nqtot,nqfils,nqdesc,iqfils, ! CRisi
+     &                     qperemin,masseqmin ! MVals
       IMPLICIT NONE
 c
@@ -1204,6 +1207,6 @@
            !Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),1e-16)
-           if (q(ij,l,iq).gt.1e-16) then
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),masseqmin)
+           if (q(ij,l,iq).gt.qperemin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else 
Index: /trunk/LMDZ.MARS/README
===================================================================
--- /trunk/LMDZ.MARS/README	(revision 2306)
+++ /trunk/LMDZ.MARS/README	(revision 2307)
@@ -2986,2 +2986,5 @@
 Some code tidying: use getin_p() instead of getin() and use "call abort_physic"
 instead of "stop" or "call abort"
+
+== 30/04/2020 == MV
+follow-up of the commit regarding the dynamical transport of isotopes: new variables for the thresholds for zq(pere) and masseq in the transport of the isotopic Ratio