Changeset 4050 for LMDZ6/trunk/libf/phylmdiso/cv3_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/cv3_routines.F90 (modified) (34 diffs)

LMDZ6/trunk/libf/phylmdiso/cv3_routines.F90

@@ -403 +403 @@
         enddo !do i=1,len
 #endif          
-! initialiser quelques variables oubliées
+! initialiser quelques variables oubliees
        do i=1,len
           plcllo(i)=0.0
@@ -900 +900 @@
           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
@@ -959 +959 @@
 
 #ifdef ISOVERIF
-            write(*,*) 'cv3_routine undilute 1 598: après condiso'
+            write(*,*) 'cv3_routine undilute 1 598: apres condiso'
           
           if (iso_eau.gt.0) then
@@ -1435 +1435 @@
 
 !JAM--------------------------------------------------------------------
-! Calcul de la quantité d'eau sous forme de glace
+! Calcul de la quantite d'eau sous forme de glace
 ! --------------------------------------------------------------------
   INTEGER nl, len
@@ -2856 +2856 @@
       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)
@@ -2873 +2873 @@
 #ifdef ISO
 #ifdef ISOVERIF
-!       write(*,*) 'cv3_routines 1820: entrée dans cv3_mixing'
+!       write(*,*) 'cv3_routines 1820: entree dans cv3_mixing'
        do i=minorig+1,nl
         do il=1,ncum
@@ -3083 +3083 @@
 !     :           '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        
@@ -3194 +3194 @@
 #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)
@@ -3431 +3431 @@
 #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)
@@ -3543 +3543 @@
 #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
@@ -3567 +3567 @@
          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
@@ -3991 +3991 @@
         call iso_verif_traceur(xtwdtrain(1,il),'cv3_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) &
@@ -4156 +4156 @@
 !!---end jyg---
 
-! --------retour à la formulation originale d'Emanuel.
+! --------retour a la formulation originale d'Emanuel.
         IF (cvflag_ice) THEN
 
@@ -4170 +4170 @@
 
 !JAM  Attention: evap=sigt*E
-!    Modification: evap devient l'évaporation en milieu de couche
-!    car nécessaire dans cv3_yield
-!    Du coup, il faut modifier pas mal d'équations...
+!    Modification: evap devient l'evaporation en milieu de couche
+!    car necessaire dans cv3_yield
+!    Du coup, il faut modifier pas mal d'equations...
 !    et l'expression de afac qui devient afac1
 !    revap=sqrt((prec(i+1)+prec(i))/2)
@@ -4191 +4191 @@
 !JYG    Dans sa formulation originale, Emanuel calcule l'evaporation par:
 ! c             evap(il,i)=sigt*afac*revap
-! ce qui n'est pas correct. Dans cv_routines, la formulation a été modifiee.
+! ce qui n'est pas correct. Dans cv_routines, la formulation a ete modifiee.
 ! Ici,l'evaporation evap est simplement calculee par l'equation de
 ! conservation.
@@ -4525 +4525 @@
 #ifdef ISO
 #ifdef ISOVERIF
-! verif des inputs à appel stewart
+! verif des inputs a appel stewart
       do il=1,ncum 
        if (i.le.inb(il) .and. lwork(il)) then
@@ -4543 +4543 @@
        enddo
 #endif
-        ! appel de appel_stewart_vectorisé
+        ! appel de appel_stewart_vectorise
         call appel_stewart_vectall_np(lwork,ncum, &
      &                   ph,t,evap,xtwdtrain, &
@@ -4611 +4611 @@
 #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
@@ -4651 +4651 @@
 #endif
           rpprec(il,i)=rs(il,i)      
-         ! sous cas rajouté le 11dec 2011. Normalement, pas utile 
+         ! sous cas rajoute le 11dec 2011. Normalement, pas utile 
          else if (rp(il,i).eq.0.0) then                 
             do ixt=1,ntraciso
@@ -4864 +4864 @@
       real xtbx(ntraciso), xtawat(ntraciso,nloc)
       ! cam debug
-      ! pour l'homogénéisation sous le nuage:
+      ! pour l'homogeneisation sous le nuage:
       real bxtsum(ntraciso,nloc), fxtsum(ntraciso,nloc)
 #ifdef DIAGISO
-        ! diagnostiques juste: tendance des différents processus
+        ! diagnostiques juste: tendance des differents processus
       real fxt_detrainement(niso,nloc,nd)
       real fxt_fluxmasse(niso,nloc,nd)
@@ -4917 +4917 @@
 #ifdef ISO
        ! cam debug
-!       write(*,*) 'cv3_routines 3082: entrée dans cv3_yield'
+!       write(*,*) 'cv3_routines 3082: entree dans cv3_yield'
        ! en cam debug
        do ixt = 1, ntraciso
@@ -4994 +4994 @@
   END DO
 #ifdef ISO
-! on initialise mieux fr et fxt par securité
+! on initialise mieux fr et fxt par securite
   fr(:,:)=0.0
   fxt(:,:,:)=0.0
@@ -5845 +5845 @@
         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
@@ -5855 +5855 @@
 !     :           -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) &
@@ -5941 +5941 @@
 #endif
                 else !if (abs(dXe).gt.ridicule) then
-                    ! dans ce cas, fxtXe doit être faible
+                    ! dans ce cas, fxtXe doit etre faible
                     
 #ifdef ISOVERIF
@@ -5957 +5957 @@
                       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.
@@ -5971 +5971 @@
                
             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) &
@@ -6176 +6176 @@
         ! 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,il)=awat(il)*(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 !do ixt = 1, ntraciso
       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(il),0.0,'cv3_yield 3779')
@@ -6850 +6850 @@
         fq_detrainement(il, i) = fq_detrainement(il, i)/alpha_qpos(il)
         do ixt=1,ntraciso
-          fq_ddft(ixt,il, i) = fq_ddft(ixt,il, i)/alpha_qpos(il)
-          fq_evapprecip(ixt,il, i) = fq_evapprecip(ixt,il, i)/alpha_qpos(il)
-          fq_fluxmasse(ixt,il, i) = fq_fluxmasse(ixt,il, i)/alpha_qpos(il)
-          fq_detrainement(ixt,il, i) = fq_detrainement(ixt,il, i)/alpha_qpos(il)
+          fxt_ddft(ixt,il, i) = fxt_ddft(ixt,il, i)/alpha_qpos(il)
+          fxt_evapprecip(ixt,il, i) = fxt_evapprecip(ixt,il, i)/alpha_qpos(il)
+          fxt_fluxmasse(ixt,il, i) = fxt_fluxmasse(ixt,il, i)/alpha_qpos(il)
+          fxt_detrainement(ixt,il, i) = fxt_detrainement(ixt,il, i)/alpha_qpos(il)
         enddo ! do ixt=1,ntraciso
 #endif       
@@ -7179 +7179 @@
     ENDDO       ! k
 
-! 14/01/15 AJ delta n'a rien à faire là...                                                 
+! 14/01/15 AJ delta n'a rien a faire la...                                                 
     DO il = 1, ncum                                                  ! cld
 !!      IF (wa(il,i)>0.0 .AND. iflag(il)<=1) &                         ! cld
@@ -7195 +7195 @@
 
 ! IM cf. FH 
-! 14/01/15 AJ ne correspond pas à ce qui a été codé par JYG et SB           
+! 14/01/15 AJ ne correspond pas a ce qui a ete code par JYG et SB           
                                                          
       IF (iflag_clw==0) THEN                                         ! cld
@@ -7290 +7290 @@
 
 ! 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, nl
     DO k = 1, nl
@@ -7576 +7576 @@
 
         ! 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"
@@ -7613 +7613 @@
 
         ! il faut calculer la cape: on fait un calcule simple car tant qu'on ne
-        ! connait pas ep, on ne connait pas les mélanges, ddfts etc... qui sont
+        ! connait pas ep, on ne connait pas les melanges, ddfts etc... qui sont
         ! necessaires au calcul de la cape dans la nouvelle physique