Changeset 4143 for LMDZ6/trunk/libf/phylmdiso/cv_driver.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/cv_driver.F90 (modified) (14 diffs)

LMDZ6/trunk/libf/phylmdiso/cv_driver.F90

@@ -25 +25 @@
   USE dimphy
 #ifdef ISO
-  USE infotrac_phy, ONLY: ntraciso,niso,index_trac,ntraceurs_zone
+  USE infotrac_phy, ONLY: ntraciso=>ntiso,niso,itZonIso,nzone
   USE isotopes_mod, ONLY: iso_eau,iso_HDO,ridicule,bidouille_anti_divergence
 #ifdef ISOVERIF
@@ -511 +511 @@
     CALL cv_param(nd)
 #ifdef ISO
-       write(*,*) 'cv_driver 454: isos pas prévus ici'
+       write(*,*) 'cv_driver 454: isos pas prevus ici'
        stop
 #endif
@@ -687 +687 @@
 !c--debug
 #ifdef ISOVERIF
-       write(*,*) 'cv_driver 621: après cv3_undilute1'
+       write(*,*) 'cv_driver 621: apres cv3_undilute1'
        do k = 1, klev
         do i = 1, klon
@@ -752 +752 @@
         !write(*,*) 'xt1(iso_eau,1,1),q1(1,1)=',xt1(iso_eau,1,1),q1(1,1)
         !write(*,*) 'xt1(iso_eau,14,1),q1(14,1)=',xt1(iso_eau,14,1),q1(14,1) 
-        !write(*,*) 'iso_eau,use_iso=',iso_eau,use_iso
        do k = 1, klev
         do i = 1, nloc
@@ -783 +782 @@
 #ifdef ISO
 #ifdef ISOVERIF
-       write(*,*) 'cv_driver 720: après cv3_compress'           
+       write(*,*) 'cv_driver 720: apres cv3_compress'           
        do k = 1, klev
         do i = 1, ncum
@@ -883 +882 @@
                 ,cape,ep,hp,icb,inb,clw,nk,t,h,lv &
                 ,epmax_diag)
-        ! on écrase ep et recalcule hp
+        ! on écrase ep et recalcule hp
     END IF
 
@@ -910 +909 @@
 #ifdef ISO
 #ifdef ISOVERIF
-       write(*,*) 'cv_driver 837: après cv3_mixing'
+       write(*,*) 'cv_driver 837: apres cv3_mixing'
        do k = 1, klev
        do j = 1, klev
@@ -925 +924 @@
            call iso_verif_traceur_justmass(xtelij(1,i,j,k), &
      &           'cv_driver 847')
-           ! on ne vérfier pas le deltaD ici car peut dépasser le seuil
-           ! raisonable pour températures très froides.
+           ! on ne verifie pas le deltaD ici car peut depasser le seuil
+           ! raisonable pour temperatures tres froides.
 #endif               
         enddo
@@ -940 +939 @@
            call iso_verif_traceur(xt(1,i,k),'cv_driver 856')
            if (option_tmin.eq.1) then
-             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) &
      &           ,'cv_driver 909').eq.1) then
@@ -946 +945 @@
                write(*,*) 'xtclw=',xtclw(:,i,k)
                stop
-             endif !if (iso_verif_positif_nostop(xtclw(index_trac(
+             endif !if (iso_verif_positif_nostop(xtclw(itZonIso(
            endif !if ((option_traceurs.eq.17).or.
 #endif  
@@ -1000 +999 @@
        write(*,*) 'klev=',klev
 #ifdef ISOVERIF
-       write(*,*) 'cv_driver 930: après cv3_unsat'
+       write(*,*) 'cv_driver 930: apres cv3_unsat'
        do k = 1, klev
         do i = 1, ncum
@@ -1048 +1047 @@
             do i = 1, ncum
                do iiso=1,niso
-                  ixt_ddft=index_trac(izone_ddft,iiso)
-                  ixt_poubelle=index_trac(izone_poubelle,iiso)
+                  ixt_ddft=itZonIso(izone_ddft,iiso)
+                  ixt_poubelle=itZonIso(izone_poubelle,iiso)
                   xtp(ixt_ddft,i,k)=xtp(ixt_ddft,i,k) &
      &                    +xtp(ixt_poubelle,i,k)
@@ -1063 +1062 @@
           do k = 1, klev
             do i = 1, ncum
-               do izone=1,ntraceurs_zone
+               do izone=1,nzone
                  if (izone.eq.izone_ddft) then
                    do iiso=1,niso
-                     ixt_ddft=index_trac(izone,iiso)
-                     ixt_revap=index_trac(izone_revap,iiso)
+                     ixt_ddft=itZonIso(izone,iiso)
+                     ixt_revap=itZonIso(izone_revap,iiso)
                      xtp(ixt_ddft,i,k)=xtp(iiso,i,k)-xtp(ixt_revap,i,k)
                    enddo !do iiso=1,niso
                  elseif (izone.eq.izone_ddft) then
-                    ! rien à faire
+                    ! rien a faire
                  else !if (izone.eq.izone_ddft) then
                    do iiso=1,niso
-                     ixt=index_trac(izone,iiso)
+                     ixt=itZonIso(izone,iiso)
                      xtp(ixt,i,k)=0.0
                    enddo !do iiso=1,niso
                  endif !if (izone.eq.izone_ddft) then
-               enddo !do izone=1,ntraceurs_zone
+               enddo !do izone=1,nzone
 #ifdef ISOVERIF
                call iso_verif_traceur(xtp(1,i,k),'cv_driver 1059')
@@ -1247 +1246 @@
 ! si icvflag_Tpa=0, alors la fraction de glace dans l'ascendance adiabatique est
   ! fonction de la temperature de l'environnement et la temperature de l'ascendance est 
-  ! calculee en deux itérations, une en supposant qu'il n'y a pas de glace et l'autre 
-  ! en ajoutant la glace (ancien schéma d'Arnaud Jam).
+  ! calculee en deux iterations, une en supposant qu'il n'y a pas de glace et l'autre 
+  ! en ajoutant la glace (ancien schema d'Arnaud Jam).
 ! si icvflag_Tpa=1, alors la fraction de glace dans l'ascendance adiabatique est
   ! fonction de la temperature de l'environnement et la temperature de l'ascendance est