Changeset 3891 for LMDZ6/branches/LMDZ-tracers/libf/dyn3dmem/vlspltqs_loc.F

Timestamp:

May 11, 2021, 2:10:34 PM (4 years ago)

Author:

dcugnet

Message:

• Bugs corrections:
  • sequential gcm fixed
  • parallel gcm compilation fixed ; to be tested
• Some generic operations moved from infotrac to readTracFile
• Fixed algebrical reduction routine, used in the isotopes parameters file.
• Additional component "comp" in the tracers descriptor derived type "tra",

specifying the model component name(s) (cf. tracers sections) it belongs.

• isotopes class selection tool fixed.

File:

• LMDZ6/branches/LMDZ-tracers/libf/dyn3dmem/vlspltqs_loc.F (modified) (6 diffs)

LMDZ6/branches/LMDZ-tracers/libf/dyn3dmem/vlspltqs_loc.F

@@ -12 +12 @@
 c   --------------------------------------------------------------------
       USE parallel_lmdz
-      USE infotrac, ONLY : nqtot, tracers, tra,        ! CRisi                 &
-     &                     qperemin,masseqmin,ratiomin ! MVals et CRisi
+      USE infotrac, ONLY : nqtot, tracers, tra,         ! CRisi                &
+     &                     qprntmin, massqmin, ratiomin ! MVals et CRisi
       IMPLICIT NONE
 c
@@ -349 +349 @@
           DO ij=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),masseqmin)
-           if (q(ij,l,iq).gt.qperemin) then ! modif 13 nov 2020
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),massqmin)
+           if (q(ij,l,iq).gt.qprntmin) then ! modif 13 nov 2020
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else
@@ -374 +374 @@
          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),masseqmin)
+            new_m=max(masse(ij,l,iq)+u_m(ij-1,l)-u_m(ij,l),massqmin)
             q(ij,l,iq)=(q(ij,l,iq)*masse(ij,l,iq)+
      &      u_mq(ij-1,l)-u_mq(ij,l))
@@ -428 +428 @@
 c   --------------------------------------------------------------------
       USE parallel_lmdz
-      USE infotrac, ONLY : nqtot, tracers, tra,        ! CRisi                 &
-     &                     qperemin,masseqmin,ratiomin ! MVals et CRisi
+      USE infotrac, ONLY : nqtot, tracers, tra,         ! CRisi                &
+     &                     qprntmin, massqmin, ratiomin ! MVals et CRisi
       USE comconst_mod, ONLY: pi
       IMPLICIT NONE
@@ -759 +759 @@
           DO ij=ijbm,ijem 
            !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),masseqmin)
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),massqmin)
           enddo !DO ij=ijbm,ijem 
 
@@ -766 +766 @@
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
            !write(lunout,*) 'ij,l,q(ij,l,iq)=',ij,l,q(ij,l,iq)
-           if (q(ij,l,iq).gt.qperemin) then ! modif 13 nov 2020
+           if (q(ij,l,iq).gt.qprntmin) then ! modif 13 nov 2020
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else