Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr.F90

Timestamp:

Jul 19, 2024, 5:41:58 PM (4 months ago)

Author:

abarral

Message:

(lint) Fix obsolete boolean operators

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr.F90 (modified) (21 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cvltr.F90

@@ -49 +49 @@
   REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: sij        ! fraction dair de lenv
   REAL,DIMENSION(klon,klev),INTENT(IN)      :: wght_cvfd  ! weights of the layers feeding convection
-  REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij       ! contenu en eau condensée spécifique/conc deau condensée massique
+  REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: elij       ! contenu en eau condens�e sp�cifique/conc deau condens�e massique
   REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: epmlmMm    ! eau condensee precipitee dans mel masse dair sat
   REAL,DIMENSION(klon,klev),INTENT(IN)      :: eplaMm    ! eau condensee precipitee dans aa masse dair sat
 
-  REAL,DIMENSION(klon,klev),INTENT(IN)      :: clw        ! contenu en eau condensée dans lasc adiab
+  REAL,DIMENSION(klon,klev),INTENT(IN)      :: clw        ! contenu en eau condens�e dans lasc adiab
   REAL,DIMENSION(klon),INTENT(IN)           :: sigd
   INTEGER,DIMENSION(klon),INTENT(IN)        :: icb,inb
@@ -80 +80 @@
 
 ! RomP ! les variables sont nettoyees des valeurs aberrantes
-  REAL,DIMENSION(klon,klev)         :: Pa, Pm  ! pluie AA et mélanges, var temporaire
+  REAL,DIMENSION(klon,klev)         :: Pa, Pm  ! pluie AA et m�langes, var temporaire
   REAL,DIMENSION(klon,klev)         :: pmflxs,pmflxr ! pmflxrIN,pmflxsIN sans valeur aberante
   REAL,DIMENSION(klon,klev)         :: mp            ! flux de masse
@@ -155 +155 @@
 ! On prend la moyenne des precip entre le niveau i+1 et i
 ! I=3/4* (P(1+1)+P(i))/2 / (sigd*r*rho_l)
-!  1000kg/m3= densité de l'eau
+!  1000kg/m3= densit� de l'eau
 ! 0.75e-3 = 3/4 /1000
-! Par la suite, I est tout le temps multiplié par sig_d pour avoir l'impaction sur la surface de la maille
-! on le néglige ici pour simplifier le code
+! Par la suite, I est tout le temps multipli� par sig_d pour avoir l'impaction sur la surface de la maille
+! on le n�glige ici pour simplifier le code
   do j=1,klev-1
    do i=1,klon
@@ -201 +201 @@
   END DO
 
-! suppression des valeurs très faibles (~1e-320)
+! suppression des valeurs tr�s faibles (~1e-320)
 ! multiplication de levaporation pour lavoir par unite de temps
 ! et par unite de surface de la maille
@@ -207 +207 @@
   DO j=1,klev
    DO i=1,klon
-    if(ev(i,j).lt.1.e-16) then
+    if(ev(i,j)<1.e-16) then
      evap(i,j)=0.
     else
@@ -217 +217 @@
   DO j=1,klev
    DO i=1,klon
-   if(j.lt.klev) then
-    if(epIN(i,j).lt.1.e-32) then
+   if(j<klev) then
+    if(epIN(i,j)<1.e-32) then
      ep(i,j)=0.
     else
@@ -226 +226 @@
     ep(i,j)=epmax
    endif
-    if(mpIN(i,j).lt.1.e-32) then
+    if(mpIN(i,j)<1.e-32) then
      mp(i,j)=0.
     else
      mp(i,j)=mpIN(i,j)
     endif
-    if(pmflxsIN(i,j).lt.1.e-32) then
+    if(pmflxsIN(i,j)<1.e-32) then
      pmflxs(i,j)=0.
     else
      pmflxs(i,j)=pmflxsIN(i,j)
     endif
-    if(pmflxrIN(i,j).lt.1.e-32) then
+    if(pmflxrIN(i,j)<1.e-32) then
      pmflxr(i,j)=0.
     else
      pmflxr(i,j)=pmflxrIN(i,j)
     endif
-    if(wdtrainA(i,j).lt.1.e-32) then
+    if(wdtrainA(i,j)<1.e-32) then
      Pa(i,j)=0.
     else
      Pa(i,j)=wdtrainA(i,j)
     endif
-    if(wdtrainM(i,j).lt.1.e-32) then
+    if(wdtrainM(i,j)<1.e-32) then
      Pm(i,j)=0.
     else
@@ -257 +257 @@
   DO j = klev-1,1,-1
    DO i = 1,klon
-     NO_precip(i,j) = (pmflxr(i,j+1)+pmflxs(i,j+1)).lt.1.e-10&
-                    .and.Pa(i,j).lt.1.e-10.and.Pm(i,j).lt.1.e-10
+     NO_precip(i,j) = (pmflxr(i,j+1)+pmflxs(i,j+1))<1.e-10&
+                    .and.Pa(i,j)<1.e-10.and.Pm(i,j)<1.e-10
    END DO
   END DO
@@ -285 +285 @@
      DO j=k-1,1,-1
         DO i=1,klon
-           if(mp(i,j+1).gt.1.e-10) then
+           if(mp(i,j+1)>1.e-10) then
               zmd(i,j,k)=zmd(i,j+1,k)*min(1.,mp(i,j)/mp(i,j+1)) !det ~ mk(j)=mk(j+1)*mp(i,j)/mp(i,j+1)
            ENDif
@@ -302 +302 @@
      DO j=1,klev-1
         DO i=1,klon
-        if(mp(i,j+1).gt.1.e-10) then
+        if(mp(i,j+1)>1.e-10) then
           qTrdi(i,j+1,it)=qTrdi(i,j+1,it)+(zmd(i,j+1,k)/mp(i,j+1))*tr(i,k,it)
         else
@@ -397 +397 @@
   DO j=1,klev
      DO i=1,klon
-        if(j.ge.icb(i).and.j.le.inb(i)) then
-          if(clw(i,j).gt.1.e-16) then
+        if(j>=icb(i).and.j<=inb(i)) then
+          if(clw(i,j)>1.e-16) then
            qPa(i,j,it)=ccntrAA_coef*tr(i,1,it)/clw(i,j)
           else
@@ -413 +413 @@
       DO k=1,j-1
         DO i=1,klon
-        if(k.ge.icb(i).and.k.le.inb(i).and.& 
-           j.le.inb(i)) then
-            if(elij(i,k,j).gt.1.e-16) then
+        if(k>=icb(i).and.k<=inb(i).and.&
+           j<=inb(i)) then
+            if(elij(i,k,j)>1.e-16) then
              qMeltmp(i,j,it)=((1-ep(i,k))*ccntrAA_coef*tr(i,1,it)&
                          *(1.-sij(i,k,j))  +ccntrENV_coef&
@@ -431 +431 @@
      DO j=1,klev
         DO i=1,klon
-          if(Mint(i,j).gt.1.e-16) then
+          if(Mint(i,j)>1.e-16) then
            qMel(i,j,it)=qpmMint(i,j,it)/Mint(i,j)
           else
@@ -442 +442 @@
    DO j=klev-1,1,-1
     DO i=1,klon
-     if(mp(i,j+1).gt.mp(i,j).and.mp(i,j+1).gt.1.e-10) then  ! detrainement
+     if(mp(i,j+1)>mp(i,j).and.mp(i,j+1)>1.e-10) then  ! detrainement
        kappa(i,j)=((pmflxr(i,j+1)+pmflxs(i,j+1)+Pa(i,j)+Pm(i,j))*&
                 (-mp(i,j+1)-imp(i,j)/RG*dxpres(i,j))&
                 + (imp(i,j)/RG*dxpres(i,j))*(evap(i,j)/RG*dxpres(i,j)))
              
-     elseif(mp(i,j).gt.mp(i,j+1).and.mp(i,j).gt.1.e-10) then! entrainement
-       if(j.eq.1) then
+     elseif(mp(i,j)>mp(i,j+1).and.mp(i,j)>1.e-10) then! entrainement
+       if(j==1) then
         kappa(i,j)=((pmflxr(i,j+1)+pmflxs(i,j+1)+Pa(i,j)+Pm(i,j))*&
                  (-mp(i,2)-imp(i,j)/RG*dxpres(i,j))&
@@ -465 +465 @@
   DO j=klev-1,1,-1
    DO i=1,klon
-    if (abs(kappa(i,j)).lt.1.e-25) then    !si denominateur nul (il peut y avoir des mp!=0)
+    if (abs(kappa(i,j))<1.e-25) then    !si denominateur nul (il peut y avoir des mp!=0)
      kappa(i,j)=1.
-     if(j.eq.1) then
-       qDi(i,j,it)=qDi(i,j+1,it) !orig tr(i,j,it)   ! mp(1)=0 donc tout vient de la couche supérieure
-     elseif(mp(i,j+1).gt.mp(i,j).and.mp(i,j+1).gt.1.e-10) then
+     if(j==1) then
+       qDi(i,j,it)=qDi(i,j+1,it) !orig tr(i,j,it)   ! mp(1)=0 donc tout vient de la couche sup�rieure
+     elseif(mp(i,j+1)>mp(i,j).and.mp(i,j+1)>1.e-10) then
        qDi(i,j,it)=qDi(i,j+1,it)
-     elseif(mp(i,j).gt.mp(i,j+1).and.mp(i,j).gt.1.e-10) then! entrainement
+     elseif(mp(i,j)>mp(i,j+1).and.mp(i,j)>1.e-10) then! entrainement
        qDi(i,j,it)=(-mp(i,j+1)*(qDi(i,j+1,it)-tr(i,j,it))-mp(i,j)*tr(i,j,it))/(-mp(i,j))
      else  ! si  mp (i)=0 et mp(j+1)=0
@@ -490 +490 @@
 !!jyg  (20130119) correction pour le sommet du nuage
 !!     if(j.ge.inb(i)) then       !au-dessus du nuage, sommet inclu
-     if(j.gt.inb(i)) then       !au-dessus du nuage
+     if(j>inb(i)) then       !au-dessus du nuage
        qDi(i,j,it)=tr(i,j,it)   ! pas de descente => environnement = descente insaturee
        qPr(i,j,it)=0.
 
 !  vvv premiere couche du modele ou mp(1)=0  ! det tout le temps  vvv
-     elseif(j.eq.1) then 
-      if(mp(i,2).gt.1.e-10) then !mp(2) non nul -> detrainement (car mp(1) = 0) !ent pas possible
+     elseif(j==1) then
+      if(mp(i,2)>1.e-10) then !mp(2) non nul -> detrainement (car mp(1) = 0) !ent pas possible
        if(NO_precip(i,j)) then !pas de precip en (i)
         qDi(i,j,it)=qDi(i,j+1,it)
@@ -530 +530 @@
      else 
 !------------------------------------------------------------- detrainement 
-      if(mp(i,j+1).gt.mp(i,j).and.mp(i,j+1).gt.1.e-10) then !mp(i,j).gt.1.e-10) then
+      if(mp(i,j+1)>mp(i,j).and.mp(i,j+1)>1.e-10) then !mp(i,j).gt.1.e-10) then
          if(NO_precip(i,j)) then 
           qDi(i,j,it)=qDi(i,j+1,it)
@@ -548 +548 @@
          endif !precip
 !------------------------------------------------------------- entrainement 
-      elseif(mp(i,j).gt.mp(i,j+1).and.mp(i,j).gt.1.e-10) then 
+      elseif(mp(i,j)>mp(i,j+1).and.mp(i,j)>1.e-10) then
          if(NO_precip(i,j)) then
           qDi(i,j,it)=(-mp(i,j+1)*(qDi(i,j+1,it)-tr(i,j,it))-mp(i,j)*tr(i,j,it))/(-mp(i,j))
@@ -614 +614 @@
                -zmfdam(i,k,it)*ccntrAA_coef
 ! lessivage courants insatures
-   if(k.le.inb(i).and.k.gt.1) then   ! tendances dans le nuage
+   if(k<=inb(i).and.k>1) then   ! tendances dans le nuage
 !------------------------------------------------------------- detrainement
-      if(mp(i,k+1).gt.mp(i,k).and.mp(i,k+1).gt.1.e-10) then
+      if(mp(i,k+1)>mp(i,k).and.mp(i,k+1)>1.e-10) then
         uscavtrac= (-mp(i,k)+mp(i,k+1))*(qDi(i,k,it)-tr(i,k,it))&
                    + mp(i,k)*(tr(i,k-1,it)-tr(i,k,it))
@@ -625 +625 @@
 !                    'mp',mp(i,k)
 !------------------------------------------------------------- entrainement 
-      elseif(mp(i,k).gt.mp(i,k+1).and.mp(i,k).gt.1.e-10) then
+      elseif(mp(i,k)>mp(i,k+1).and.mp(i,k)>1.e-10) then
          uscavtrac= mp(i,k)*(tr(i,k-1,it)-tr(i,k,it))
 !
@@ -641 +641 @@
       endif  ! mp/det/ent
 !------------------------------------------------------------- premiere couche 
-   elseif(k.eq.1) then  !                                      mp(1)=0.
-      if(mp(i,2).gt.1.e-10) then  !detrainement
+   elseif(k==1) then  !                                      mp(1)=0.
+      if(mp(i,2)>1.e-10) then  !detrainement
          uscavtrac= (-0.+mp(i,2))*(qDi(i,k,it)-tr(i,k,it)) !&
 !                   + mp(i,2)*(0.-tr(i,k,it))