Changeset 4143 for LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90

Timestamp:

May 9, 2022, 12:35:40 PM (2 years ago)

Author:

dcugnet

Message:

• Some variables are renamed or replaced by direct equivalents:
  • iso_indnum -> tracers(:)%iso_iName
  • niso_possibles -> niso
  • iqiso -> iqIsoPha ; index_trac -> itZonIso
  • ok_iso_verif -> isoCheck
  • ntraceurs_zone -> nzone ; ntraciso -> ntiso
  • qperemin -> min_qparent ; masseqmin -> min_qmass ; ratiomin -> min_ratio
• Some renamed variables are only aliased with the older name (using USE <module>, ONLY: <oldName> => <newName>) in routines where they are repeated many times.
• Few hard-coded indexes are now computed (examples: ilic, iso, ivap, irneb, iq_vap, iq_liq, iso_H2O, iso_HDO, iso_HTO, iso_O17, iso_O18).
• The IF(isoCheck) test is now embedded in the check_isotopes_seq and check_isotopes_loc routines (lighter calling).

File:

• LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90 (modified) (20 diffs)

LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90

@@ -314 +314 @@
      &   )
 #ifdef ISO
-    use infotrac_phy, ONLY: ntraciso
+    use infotrac_phy, ONLY: ntraciso=>ntiso
 #ifdef ISOVERIF
     use isotopes_verif_mod, ONLY: iso_verif_positif,iso_verif_noNaN,iso_verif_egalite
@@ -685 +685 @@
      &                   )
 #ifdef ISO
-USE infotrac_phy, ONLY: ntraciso
+USE infotrac_phy, ONLY: ntraciso=>ntiso
 USE isotopes_mod, ONLY: pxtmelt,pxtice,pxtmin,pxtmax,cond_temp_env, &
         iso_eau,iso_HDO,ridicule
@@ -1274 +1274 @@
   USE print_control_mod, ONLY: lunout
 #ifdef ISO
-    use infotrac_phy, ONLY: ntraciso
+    use infotrac_phy, ONLY: ntraciso=>ntiso
     use isotopes_mod, ONLY: essai_convergence, iso_eau,iso_HDO
 #ifdef ISOVERIF
@@ -1474 +1474 @@
   USE print_control_mod, ONLY: prt_level
 #ifdef ISO
-use infotrac_phy, ONLY: ntraciso
+use infotrac_phy, ONLY: ntraciso=>ntiso
 USE isotopes_mod, ONLY: pxtmelt,pxtice,pxtmin,pxtmax,cond_temp_env, &
         iso_eau,iso_HDO
@@ -2777 +2777 @@
 
 #ifdef ISO
-use infotrac_phy, ONLY: ntraciso,niso,index_trac
+use infotrac_phy, ONLY: ntraciso=>ntiso,niso,itZonIso
 USE isotopes_mod, ONLY: pxtmelt,pxtice,pxtmin,pxtmax, iso_eau,iso_HDO, &
         ridicule
@@ -3580 +3580 @@
         call iso_verif_traceur(xtclw(1,il,im), &
      &          'condiso_liq_ice_vectiso_trac 358')
-        if (iso_verif_positif_nostop(xtclw(index_trac( &
+        if (iso_verif_positif_nostop(xtclw(itZonIso( &
      &           izone_cond,iso_eau),i,k)-xtclw(iso_eau,i,k) &
      &           ,'cv3_routines 909').eq.1) then
@@ -3588 +3588 @@
      &             niso,ntraciso,index_zone,izone_cond      
                stop
-         endif !if (iso_verif_positif_nostop(xtclw(index_trac(
+         endif !if (iso_verif_positif_nostop(xtclw(itZonIso(
 #endif             
          enddo !do il = 1, ncum   
@@ -3615 +3615 @@
   USE print_control_mod, ONLY: prt_level, lunout
 #ifdef ISO
-    use infotrac_phy, ONLY: ntraciso
+    use infotrac_phy, ONLY: ntraciso=>ntiso
     use isotopes_mod, ONLY: essai_convergence, iso_eau,iso_HDO, &
         ridicule
@@ -3628 +3628 @@
 #ifdef ISOTRAC
     use isotrac_mod, only: option_cond,izone_cond
-    use infotrac_phy, ONLY: index_trac
+    use infotrac_phy, ONLY: itZonIso
 #ifdef ISOVERIF
     use isotopes_verif_mod, ONLY: iso_verif_traceur_justmass, &
@@ -3993 +3993 @@
           ! on verifie que tout le detrainement est tagge condensat
           if (iso_verif_positif_nostop( &
-     &          xtwdtrain(index_trac(izone_cond,iso_eau),il) &
+     &          xtwdtrain(itZonIso(izone_cond,iso_eau),il) &
      &          -xtwdtrain(iso_eau,il), &
      &          'cv3_routines 2795').eq.1) then
@@ -4535 +4535 @@
 !        if (option_tmin.ge.1) then
 !           call iso_verif_positif(xtwater(
-!     :           index_trac(izone_cond,iso_eau),il,i+1)
+!     :           itZonIso(izone_cond,iso_eau),il,i+1)
 !     :           -xtwater(iso_eau,il,i+1),
 !     :          'cv3_routines 3083')
@@ -4602 +4602 @@
 !        if (option_tmin.ge.1) then
 !         call iso_verif_positif(xtwater(
-!     :           index_trac(izone_cond,iso_eau),il,i)
+!     :           itZonIso(izone_cond,iso_eau),il,i)
 !     :           -xtwater(iso_eau,il,i),
 !     :          'cv3_routines 3143')
@@ -4741 +4741 @@
 
 #ifdef ISO
-    use infotrac_phy, ONLY: ntraciso,niso, &
-&       ntraceurs_zone,index_trac
+    use infotrac_phy, ONLY: ntraciso=>ntiso, niso, nzone, itZonIso
     use isotopes_mod, ONLY: essai_convergence,ridicule,iso_eau,iso_HDO,iso_O18
 #ifdef ISOVERIF
@@ -5889 +5888 @@
           do iiso = 1, niso
              
-             ixt_ddft=index_trac(izone_ddft,iiso)  
+             ixt_ddft=itZonIso(izone_ddft,iiso)  
              if (mp(il,i).gt.mp(il,i+1)) then
                 fxtYe(iiso)=0.01*grav*dpinv*mp(il,i) &
@@ -5902 +5901 @@
      &              -mp(il,i)*(xtp(iiso,il,i)-xt(iiso,il,i-1))) 
         
-             ixt_poubelle=index_trac(izone_poubelle,iiso)
+             ixt_poubelle=itZonIso(izone_poubelle,iiso)
              fxt(ixt_ddft,il,i)=fxt(ixt_ddft,il,i)+fxtYe(iiso)
              fxt(ixt_poubelle,il,i)=fxt(ixt_poubelle,il,i) &
@@ -5919 +5918 @@
      &              -mp(il,i)*(xtp(iiso,il,i)-xt(iiso,il,i-1)))
 
-                ixt_ddft=index_trac(izone_ddft,iiso) 
+                ixt_ddft=itZonIso(izone_ddft,iiso) 
                 fxtYe(iiso)=0.01*grav*dpinv*mp(il,i) &
      &           *(xt(ixt_ddft,il,i-1)-xt(ixt_ddft,il,i))
                 fxt(ixt_ddft,il,i)=fxt(ixt_ddft,il,i)+fxtYe(iiso)  
 
-               ixt_revap=index_trac(izone_revap,iiso)  
+               ixt_revap=itZonIso(izone_revap,iiso)  
                fxt_revap(iiso)=0.01*grav*dpinv*(mp(il,i+1)* &
      &                  (xtp(ixt_revap,il,i+1)-xt(ixt_revap,il,i)) &
@@ -5935 +5934 @@
      &                   -xt(ixt_ddft,il,i)-xt(ixt_revap,il,i)
                 if (Xe(iiso).gt.ridicule) then
-                  do izone=1,ntraceurs_zone
+                  do izone=1,nzone
                    if ((izone.ne.izone_revap).and. &
      &                   (izone.ne.izone_ddft)) then
-                    ixt=index_trac(izone,iiso) 
+                    ixt=itZonIso(izone,iiso) 
                     fxt(ixt,il,i)=fxt(ixt,il,i) &
      &                   +xt(ixt,il,i)/Xe(iiso)*fxtXe(iiso)
                    endif !if ((izone.ne.izone_revap).and.
-                  enddo !do izone=1,ntraceurs_zone   
+                  enddo !do izone=1,nzone   
 #ifdef ISOVERIF
 !                write(*,*) 'iiso=',iiso
@@ -5964 +5963 @@
                 endif
 #endif                   
-                do izone=1,ntraceurs_zone
+                do izone=1,nzone
                    if ((izone.ne.izone_revap).and. &
      &                   (izone.ne.izone_ddft)) then                    
-                    ixt=index_trac(izone,iiso) 
+                    ixt=itZonIso(izone,iiso) 
                     if (izone.eq.izone_poubelle) then
                       fxt(ixt,il,i)=fxt(ixt,il,i)+fxtXe(iiso)
@@ -5974 +5973 @@
                     endif !if (izone.eq.izone_poubelle) then
                    endif !if ((izone.ne.izone_revap).and.
-                enddo !do izone=1,ntraceurs_zone
+                enddo !do izone=1,nzone
 #ifdef ISOVERIF
                   call iso_verif_traceur_justmass(fxt(1,il,i), &
@@ -7459 +7458 @@
      &     )    
 #ifdef ISO
-    use infotrac_phy, ONLY: ntraciso
+    use infotrac_phy, ONLY: ntraciso=>ntiso
 #ifdef ISOVERIF
     use isotopes_verif_mod, ONLY: Tmin_verif,iso_verif_aberrant, &