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

Timestamp:

May 11, 2021, 2:10:34 PM (3 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/vlsplt_loc.F (modified) (9 diffs)

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

@@ -14 +14 @@
 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
@@ -346 +346 @@
            ! 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),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
@@ -369 +369 @@
          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))
@@ -425 +425 @@
 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 
 
@@ -765 +765 @@
          DO ij=ijb,ije
            !MVals: veiller a ce qu'on n'ait pas de denominateur nul
-           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
@@ -901 +901 @@
       USE parallel_lmdz
       USE vlz_mod
-      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
@@ -949 +949 @@
 !$OMP THREADPRIVATE(first)
 
-      !REAL masseq(ijb_u:ije_u,llm,nqtot),Ratio(ijb_u:ije_u,llm,nqtot) ! CRisi
+      !REAL massq(ijb_u:ije_u,llm,nqtot),Ratio(ijb_u:ije_u,llm,nqtot) ! CRisi
       ! Ces varibles doivent être déclarées en pointer et en save dans
       ! vlz_loc si on veut qu'elles soient vues par tous les threads.  
       INTEGER ichld,iq2 ! CRisi
-
+      TYPE(tra), POINTER :: tr
+      tr => tracers(iq)
 
       IF (first) THEN
@@ -1174 +1175 @@
           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
+           masse(ij,l,iq2)=max(masse(ij,l,iq)*q(ij,l,iq),massqmin)
+           if (q(ij,l,iq).gt.qprntmin) then
              Ratio(ij,l,iq2)=q(ij,l,iq2)/q(ij,l,iq)
            else