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

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/dyn3dpar/vlsplt_p.F (modified) (7 diffs)

trunk/LMDZ.COMMON/libf/dyn3dpar/vlsplt_p.F

@@ -214 +214 @@
 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 +529 @@
            ! 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 +550 @@
          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 +602 @@
 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 +920 @@
            !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 +1051 @@
 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 +1207 @@
            !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