Changeset 2307 for trunk/LMDZ.COMMON/libf/dyn3d

Timestamp:

Apr 30, 2020, 5:14:24 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
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
MV

File:

• trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F (modified) (7 diffs)

trunk/LMDZ.COMMON/libf/dyn3d/vlsplt.F

@@ -139 +139 @@
       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 +460 @@
            !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 +482 @@
          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 +519 @@
       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 +790 @@
            !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 +888 @@
       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 +1018 @@
            !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