Changeset 4050 for LMDZ6/trunk/libf/phylmdiso/cv30_routines.F90

Timestamp:

Dec 23, 2021, 6:54:17 PM (2 years ago)

Author:

dcugnet

Message:

Second commit for new tracers.

• include most of the keys in the tracers descriptor vector "tracers(:)".
• fix in phylmdiso/cv3_routines: fq_* variables were used where their fxt_* counterparts were expected.
• multiple IF(nqdesc(iq)>0) and IF(nqfils(iq)>0) tests suppressed, because they are not needed: "do ... enddo" loops with 0 upper bound are not executed.
• remove French accents from comments (encoding problem) in phylmdiso/cv3_routines and phylmdiso/cv30_routines.
• modifications in "isotopes_verif_mod", where the call to function "iso_verif_tag17_q_deltad_chn" in "iso_verif_tag17_q_deltad_chn" was not detected at linking stage, although defined in the same module (?).

File:

• LMDZ6/trunk/libf/phylmdiso/cv30_routines.F90 (modified) (38 diffs)

LMDZ6/trunk/libf/phylmdiso/cv30_routines.F90

@@ -449 +449 @@
 
 #ifdef ISOVERIF
-        write(*,*) 'cv30_routine undilute 1 413: entrée'
+        write(*,*) 'cv30_routine undilute 1 413: entree'
 #endif
 
@@ -602 +602 @@
           zfice(i) = MIN(MAX(zfice(i),0.0),1.0)          
        enddo
-       ! calcul de la composition du condensat glacé et liquide
+       ! calcul de la composition du condensat glace et liquide
 
        do i=1,len
@@ -647 +647 @@
 
 #ifdef ISOVERIF
-            write(*,*) 'cv30_routine undilute 1 598: après condiso'
+            write(*,*) 'cv30_routine undilute 1 598: apres condiso'
           
           if (iso_eau.gt.0) then
@@ -1012 +1012 @@
             else
               q(i,k)=0.0
-              clw(i,k)=0.0 ! mise en commentaire le 5 avril pour vérif
+              clw(i,k)=0.0 ! mise en commentaire le 5 avril pour verif
 !            convergence
             endif  !f (negation(essai_convergence)) then
@@ -1908 +1908 @@
       real xtrti(ntraciso,nloc)
       real xtres(ntraciso)
-      ! on ajoute la dimension nloc à xtrti pour vérifs dans les tags: 5 fev
+      ! on ajoute la dimension nloc a xtrti pour verifs dans les tags: 5 fev
       ! 2010
       real zfice(nloc),zxtliq(ntraciso,nloc),zxtice(ntraciso,nloc)
@@ -1923 +1923 @@
 #ifdef ISO
 #ifdef ISOVERIF
-      write(*,*) 'cv30_routines 1820: entrée dans cv3_mixing'
+      write(*,*) 'cv30_routines 1820: entree dans cv3_mixing'
       if (iso_eau.gt.0) then
       call iso_verif_egalite_vect2D( &
@@ -1965 +1965 @@
              xtelij(ixt,i,k,j)=0.0
             enddo !do ixt =1,niso
-            ! on initialise mieux que ça qent et elij, même si au final les
-            ! valeurs en nd=nl+1 ne sont pas utilisées
+            ! on initialise mieux que ca qent et elij, meme si au final les
+            ! valeurs en nd=nl+1 ne sont pas utilisees
             qent(i,k,j)=rr(i,j)
             elij(i,k,j)=0.0    
@@ -2121 +2121 @@
 !     :           'tcond(il),rs(il,j)=',
 !     :            il,i,j,xtent(:,il,i,j),tcond(il),rs(il,j)
-        ! colorier la vapeur résiduelle selon température de
-        ! condensation, et le condensat en un tag spécifique
+        ! colorier la vapeur residuelle selon temperature de
+        ! condensation, et le condensat en un tag spEcifique
           if ((elij(il,i,j).gt.0.0).and.(qent(il,i,j).gt.0.0)) then  
             if (option_traceurs.eq.17) then        
@@ -2241 +2241 @@
 #ifdef ISOTRAC          
         if (option_tmin.ge.1) then
-        ! colorier la vapeur résiduelle selon température de
-        ! condensation, et le condensat en un tag spécifique
+        ! colorier la vapeur residuelle selon temperature de
+        ! condensation, et le condensat en un tag specifique
 !        write(*,*) 'cv3 tmp 2095 il,i,j,xtent(:,il,i,j)=',
 !     :            il,i,j,xtent(:,il,i,j)
@@ -2475 +2475 @@
 #ifdef ISOTRAC          
         if (option_tmin.ge.1) then
-        ! colorier la vapeur résiduelle selon température de
-        ! condensation, et le condensat en un tag spécifique
+        ! colorier la vapeur residuelle selon temperature de
+        ! condensation, et le condensat en un tag specifique
 !        write(*,*) 'cv3 tmp 2314 il,i,j,xtent(:,il,i,j)=',
 !     :            il,i,j,xtent(:,il,i,j)
@@ -2579 +2579 @@
 #ifdef ISO
 #ifdef ISOTRAC
-        ! seulement à la fin on taggue le condensat
+        ! seulement a la fin on taggue le condensat
         if (option_cond.ge.1) then
          do im = 1, nd
          do jm = 1, nd
          do il = 1, ncum    
-           ! colorier le condensat en un tag spécifique
+           ! colorier le condensat en un tag specifique
            do ixt=niso+1,ntraciso
              if (index_zone(ixt).eq.izone_cond) then
@@ -2603 +2603 @@
          do im = 1, nd
          do il = 1, ncum    
-           ! colorier le condensat en un tag spécifique
+           ! colorier le condensat en un tag specifique
            do ixt=niso+1,ntraciso
              if (index_zone(ixt).eq.izone_cond) then
@@ -2739 +2739 @@
   ! ------------------------------------------------------
 !#ifdef ISOVERIF
-!        write(*,*) 'cv30_routines 2382: entrée dans cv3_unsat'
+!        write(*,*) 'cv30_routines 2382: entree dans cv3_unsat'
 !#endif
 
@@ -2747 +2747 @@
   mp(:, :) = 0.
 #ifdef ISO
-  ! initialisation plus complète de water et rp
+  ! initialisation plus complete de water et rp
   water(:,:)=0.0
   xtwater(:,:,:)=0.0
@@ -2936 +2936 @@
         call iso_verif_traceur(xtwdtrain(1,il),'cv30_routine 2540')
         if (option_cond.ge.1) then
-           ! on vérifie que tout le détrainement est taggé condensat
+           ! on verifie que tout le detrainement est tagge condensat
            if (iso_verif_positif_nostop( &
      &          xtwdtrain(index_trac(izone_cond,iso_eau),il) &
@@ -3032 +3032 @@
 
 #ifdef ISO
-      ! ajout cam: éviter les evaporations ou eaux négatives
-!      water(il,i)=max(0.0,water(il,i)) ! ceci est toujours vérifié
+      ! ajout cam: eviter les evaporations ou eaux negatives
+!      water(il,i)=max(0.0,water(il,i)) ! ceci est toujours verifie
 #ifdef ISOVERIF
           call iso_verif_positif(water(il,i),'cv30_unsat 2376')
@@ -3189 +3189 @@
 #ifdef ISO
 #ifdef ISOVERIF
-! verif des inputs à appel stewart
+! verif des inputs a appel stewart
 !        write(*,*) 'cv30_routines 2842 tmp: appel de appel_stewart'
       do il=1,ncum 
@@ -3208 +3208 @@
        enddo
 #endif
-        ! appel de appel_stewart_vectorisé
+        ! appel de appel_stewart_vectorise
         call appel_stewart_vectall(lwork,ncum, &
      &                   ph,t,evap,xtwdtrain, &
@@ -3268 +3268 @@
 #endif
         
-! équivalent isotopique de rp(il,i)=amin1(rp(il,i),rs(il,i))
+! equivalent isotopique de rp(il,i)=amin1(rp(il,i),rs(il,i))
        do il=1,ncum
         if (i.lt.inb(il) .and. lwork(il)) then
@@ -3463 +3463 @@
       real xtbx(ntraciso), xtawat(ntraciso)
       ! cam debug
-      ! pour l'homogénéisation sous le nuage:
+      ! pour l'homogeneisation sous le nuage:
       real frsum(nloc), bxtsum(ntraciso,nloc), fxtsum(ntraciso,nloc)
-      ! correction dans calcul tendance liée à Am:
+      ! correction dans calcul tendance liee a Am:
       real dq_tmp,k_tmp,dx_tmp,R_tmp,dqreste_tmp,dxreste_tmp,kad_tmp
       logical correction_excess_aberrant
       parameter (correction_excess_aberrant=.false.)
-        ! correction qui permettait d'éviter deltas et dexcess aberrants. Mais
+        ! correction qui permettait d'eviter deltas et dexcess aberrants. Mais
         ! pb: ne conserve pas la masse d'isotopes!
 #ifdef DIAGISO
-        ! diagnostiques juste: tendance des différents processus
+        ! diagnostiques juste: tendance des differents processus
       real fxt_detrainement(ntraciso,nloc,nd)
       real fxt_fluxmasse(ntraciso,nloc,nd)
@@ -3517 +3517 @@
 #ifdef ISO
        ! cam debug
-!       write(*,*) 'cv30_routines 3082: entrée dans cv3_yield'
+!       write(*,*) 'cv30_routines 3082: entree dans cv3_yield'
        ! en cam debug
        do ixt = 1, ntraciso
@@ -3749 +3749 @@
         do ixt = 1, ntraciso
 !        fxt_fluxmasse(ixt,il,1)=fxt_fluxmasse(ixt,il,1) &
-!     &      +0.01*grav*am(il)*(xt(ixt,il,2)-xt(ixt,il,1))*work(il) ! déplacé
-!     plus haut car il existe différents cas
+!     &      +0.01*grav*am(il)*(xt(ixt,il,2)-xt(ixt,il,1))*work(il) ! deplace
+!     plus haut car il existe differents cas
         fxt_ddft(ixt,il,1)=fxt_ddft(ixt,il,1) &
      &      +0.01*grav*mp(il,2)*(xtp(ixt,il,2)-xt(ixt,il,1))*work(il)
@@ -3759 +3759 @@
 
 
-        ! pour l'ajout de la tendance liée au flux de masse Am, il faut être
+        ! pour l'ajout de la tendance liee au flux de masse Am, il faut etre
         ! prudent.
         ! On a dq1=k*(q2-q1) avec k=dt*0.01*grav*am(il)*work(il)
-        ! Pour les isotopes, la formule utilisée depuis 2006 et qui avait toujours marché est:
+        ! Pour les isotopes, la formule utilisee depuis 2006 et qui avait toujours marche est:
         ! dx1=k*(x2-x1)
-        ! Mais on plante dans un cas pathologique en décembre 2017 lors du test
-        ! d'un cas d'Anne Cozic: les isotopes deviennent négatifs.
+        ! Mais on plante dans un cas pathologique en decembre 2017 lors du test
+        ! d'un cas d'Anne Cozic: les isotopes deviennent negatifs.
         ! C'est un cas pas physique: on perd 99% de la masse de vapeur d'eau! 
         ! q2=1.01e-3 et q1=1.25e-3 kg/kg
-        ! et dq=-1.24e-3: comment est-ce possible qu'un flux venant d'un air à
-        ! q2= 1.01e-3 assèche q1 jusqu'à 0.01e-3kg/kg!
-        ! Pour les isotopes, ça donne des x1+dx négatifs.
-        ! Ce n'est pas physique mais il faut quand même s'adapter.
-        ! Pour cela, on considère que d'abord on fait rentrer le flux de masse
+        ! et dq=-1.24e-3: comment est-ce possible qu'un flux venant d'un air a
+        ! q2= 1.01e-3 asseche q1 jusqu'a 0.01e-3kg/kg!
+        ! Pour les isotopes, ca donne des x1+dx negatifs.
+        ! Ce n'est pas physique mais il faut quand meme s'adapter.
+        ! Pour cela, on considere que d'abord on fait rentrer le flux de masse
         ! descendant, et ensuite seulement on fait sortir le flux de masse
         ! sortant.
@@ -3778 +3778 @@
         ! isotopique de la vapeur d'eau q1. 
         ! A la fin, on a R=(x1+dx)/(q1+dq)=(x1+k*x2)/(q1+k*q2)
-        ! On vérifie que quand k est petit, on tend vers la formulation
+        ! On verifie que quand k est petit, on tend vers la formulation
         ! habituelle.
-        ! Comme on est habitués à la formulation habituelle, qu'elle a fait ses
-        ! preuves, on la garde sauf dans le cas où dq/q<-0.9 où on utilise la
+        ! Comme on est habitues a la formulation habituelle, qu'elle a fait ses
+        ! preuves, on la garde sauf dans le cas ou dq/q<-0.9 ou on utilise la
         ! nouvelle formulation.
         ! rappel: dq_tmp=0.01*grav*am(il)*(rr(il,2)-rr(il,1))*work(il)*delt
-        ! Même avec cette nouvelle foirmulation, on a encore des isotopes
-        ! négatifs, cette fois à cause des ddfts
-        ! On considère donc les tendances et série et non en parallèle quand on
+        ! Meme avec cette nouvelle foirmulation, on a encore des isotopes
+        ! negatifs, cette fois a cause des ddfts
+        ! On considere donc les tendances et serie et non en parallele quand on
         ! calcule R_tmp.
         dq_tmp=0.01*grav*am(il)*(rr(il,2)-rr(il,1))*work(il)*delt ! utile ci-dessous
         if ((dq_tmp/rr(il,1).lt.-0.9).and.correction_excess_aberrant) then
-                ! nouvelle formulation où on fait d'abord entrer k*q2 et ensuite
+                ! nouvelle formulation ou on fait d'abord entrer k*q2 et ensuite
                 ! seulement on fait sortir k*q1 sans changement de composition
                 ! isotopique
@@ -3828 +3828 @@
            enddo ! do ixt = 1, ntraciso
         else !if (dq_tmp/rr(il,1).lt.-0.9) then
-                ! formulation habituelle qui avait toujours marché de 2006 à
-                ! décembre 2017.
+                ! formulation habituelle qui avait toujours marche de 2006 a
+                ! decembre 2017.
            do ixt = 1, ntraciso     
                 fxt(ixt,il,1)=fxt(ixt,il,1) &
@@ -4232 +4232 @@
         ! ad.
 #endif
-       ! ici, on sépare 2 cas, pour éviter le cas pathologique décrit plus haut
-       ! pour la tendance liée à Am en i=1, qui peut conduire à des isotopes
-       ! négatifs dans les cas où les flux de masse soustrait plus de 90% de la
-       ! vapeur de la couche. Voir plus haut le détail des équations.
-       ! La différence ici est qu'on considère les flux de masse amp1 et ad en
-       ! même temps.
+       ! ici, on separe 2 cas, pour eviter le cas pathologique decrit plus haut
+       ! pour la tendance liee a Am en i=1, qui peut conduire a des isotopes
+       ! negatifs dans les cas ou les flux de masse soustrait plus de 90% de la
+       ! vapeur de la couche. Voir plus haut le detail des equations.
+       ! La difference ici est qu'on considere les flux de masse amp1 et ad en
+       ! meme temps.
        dq_tmp= 0.01*grav*dpinv*(amp1(il)*(rr(il,i+1)-rr(il,i)) &
     &            -ad(il)*(rr(il,i)-rr(il,i-1)))*delt
-       ! c'est équivalent à dqi= kamp1*qip1+kad*qim1-(kamp1+kad)*qi
+       ! c'est equivalent a dqi= kamp1*qip1+kad*qim1-(kamp1+kad)*qi
        if ((dq_tmp/rr(il,i).lt.-0.9).and.correction_excess_aberrant) then 
         ! nouvelle formulation
@@ -4430 +4430 @@
         ! on change le traitement de cette ligne le 8 mai 2009:
         ! avant, on avait: xtawat=xtelij(il,k,i)-(1.-xtep(il,i))*xtclw(il,i)
-        ! c'est à dire que Rawat=Relij+(1-ep)*clw/awat*(Relij-Rclw)
-        ! si Relij!=Rclw, alors un fractionnement isotopique non physique était
+        ! c'est a dire que Rawat=Relij+(1-ep)*clw/awat*(Relij-Rclw)
+        ! si Relij!=Rclw, alors un fractionnement isotopique non physique etait
         ! introduit.
-        ! En fait, awat représente le surplus de condensat dans le mélange par
-        ! rapport à celui restant dans la colonne adiabatique
-        ! ce surplus à la même compo que le elij, sans fractionnement.
-        ! d'où le nouveau traitement ci-dessous.
+        ! En fait, awat represente le surplus de condensat dans le melange par
+        ! rapport a celui restant dans la colonne adiabatique
+        ! ce surplus a la meme compo que le elij, sans fractionnement.
+        ! d'ou le nouveau traitement ci-dessous.
       if (elij(il,k,i).gt.0.0) then
         do ixt = 1, ntraciso
           xtawat(ixt)=awat*(xtelij(ixt,il,k,i)/elij(il,k,i))
-!          xtawat(ixt)=amax1(xtawat(ixt),0.0) ! pas nécessaire
+!          xtawat(ixt)=amax1(xtawat(ixt),0.0) ! pas necessaire
         enddo
       else !if (elij(il,k,i).gt.0.0) then
           ! normalement, si elij(il,k,i)<=0, alors awat=0
-          ! on le vérifie. Si c'est vrai -> xtawat=0 aussi
+          ! on le verifie. Si c'est vrai -> xtawat=0 aussi
 #ifdef ISOVERIF
         call iso_verif_egalite(awat,0.0,'cv30_yield 3779')
@@ -4942 +4942 @@
      &       'cv30_yield 5029,O18, evap')
           if ((il.eq.1636).and.(i.eq.9)) then
-            write(*,*) 'cv30_yield 5057: ici, on vérifie deltaD_nobx'
+            write(*,*) 'cv30_yield 5057: ici, on verifie deltaD_nobx'
             write(*,*) 'il,i=',il,i
             write(*,*) 'fr(il,i),bx,fr(il,i)-bx=',fr(il,i),bx,fr(il,i)-bx
@@ -4973 +4973 @@
         else ! taggage des ddfts:
         ! la formule pour fq_ddft suppose que le ddft est en RP. Ce n'est pas le
-        ! cas pour le water tagging puisqu'il y a conversion des molécules
-        ! blances entrainées en molécule rouges.
+        ! cas pour le water tagging puisqu'il y a conversion des molecules
+        ! blances entrainees en molecule rouges.
         ! Il faut donc prendre en compte ce taux de conversion quand
         ! entrainement d'env vers ddft
@@ -4983 +4983 @@
 !     :           -conversion(iiso)   
 
-        ! Pb: quand on discretise, dqp/dt n'est pas vérifée numériquement.
-        ! on se retrouve donc avec des d Ye/dt différents de 0 même si ye=0 ( on
-        ! note X les molécules poubelles et Y les molécules ddfts).
+        ! Pb: quand on discretise, dqp/dt n'est pas verifee numeriquement.
+        ! on se retrouve donc avec des d Ye/dt differents de 0 meme si ye=0 ( on
+        ! note X les molecules poubelles et Y les molecules ddfts).
 
         ! Solution alternative: Dans le cas entrainant, Ye ne varie que par
         ! ascendance compensatoire des ddfts et par perte de Ye vers le ddft. On
-        ! calcule donc ce terme directement avec schéma amont: 
-
-        ! ajout déjà de l'évap
+        ! calcule donc ce terme directement avec schema amont: 
+
+        ! ajout deja de l'evap
         do ixt = 1+niso,ntraciso
              fxt(ixt,il,i)=fxt(ixt,il,i) &
@@ -5069 +5069 @@
 #endif
                 else !if (abs(dXe).gt.ridicule) then
-                    ! dans ce cas, fxtXe doit être faible
+                    ! dans ce cas, fxtXe doit etre faible
                     
 #ifdef ISOVERIF
@@ -5085 +5085 @@
                       fxt(ixt,il,i)=fxt(ixt,il,i)+fxtXe(iiso)
                     else !if (izone.eq.izone_poubelle) then
-                        ! pas de tendance pour ce tag là
+                        ! pas de tendance pour ce tag la
                     endif !if (izone.eq.izone_poubelle) then
                    endif !if ((izone.ne.izone_revap).and.
@@ -5099 +5099 @@
                
             else !if (mp(il,i).gt.mp(il,i+1)) then
-                ! cas détrainant: pas de problèmes
+                ! cas detrainant: pas de problemes
                 do ixt=1+niso,ntraciso
                 fxt(ixt,il,i)=fxt(ixt,il,i) &
@@ -5389 +5389 @@
   DO il = 1, ncum
 
-! attention, on corrige un problème C Risi
+! attention, on corrige un probleme C Risi
       IF (cvflag_grav) then
 
@@ -5722 +5722 @@
 !             write(*,*) 'cv30_routine 3990: fin des il pour i=',i
           enddo !do i=1,nl
-!          write(*,*) 'cv30_routine 3990: fin des vérifs sur homogen'
+!          write(*,*) 'cv30_routine 3990: fin des verifs sur homogen'
 #endif
 
@@ -6027 +6027 @@
 
   ! fraction deau condensee dans les melanges convertie en precip : epm
-  ! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
+  ! et eau condensee precipitee dans masse d'air sature : l_m*dM_m/dzdz.dzdz
   DO j = 1, nam1
     DO k = 1, j - 1
@@ -6226 +6226 @@
 
 #ifdef ISOVERIF
-        write(*,*) 'cv30_routines 4293: entrée dans cv3_uncompress'
+        write(*,*) 'cv30_routines 4293: entree dans cv3_uncompress'
 #endif
   DO i = 1, ncum
@@ -6346 +6346 @@
 
         ! On fait varier epmax en fn de la cape
-        ! Il faut donc recalculer ep, et hp qui a déjà été calculé et
-        ! qui en dépend
-        ! Toutes les autres variables fn de ep sont calculées plus bas.
+        ! Il faut donc recalculer ep, et hp qui a deja ete calcule et
+        ! qui en depend
+        ! Toutes les autres variables fn de ep sont calculees plus bas.
 
 #include "cvthermo.h"