Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_dry.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/lmdz_thermcell_dry.F90 (modified) (9 diffs)

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

@@ -15 +15 @@
 ! Le calcul pourrait etre sans doute simplifier.
 ! La temperature potentielle virtuelle dans la panache ascendant est
-! la temperature potentielle virtuelle pondérée par alim_star.
+! la temperature potentielle virtuelle pondérée par alim_star.
 !--------------------------------------------------------------------------
        USE lmdz_thermcell_ini, ONLY: prt_level, RG
@@ -73 +73 @@
        do l=1,nlay-2
          do ig=1,ngrid
-            if (l.eq.lmin(ig).and.lalim(ig).gt.1) then
+            if (l==lmin(ig).and.lalim(ig)>1) then
 
 !------------------------------------------------------------------------
@@ -90 +90 @@
 ! 1. le flux de masse en haut  f_star(ig,l+1)
 ! 2. la temperature potentielle virtuelle dans la couche ztva(ig,l)
-! 3. la vitesse au carré en haut zw2(ig,l+1)
+! 3. la vitesse au carr� en haut zw2(ig,l+1)
 !------------------------------------------------------------------------
 
-            else if (zw2(ig,l).ge.1e-10) then
+            else if (zw2(ig,l)>=1e-10) then
 
                ztva(ig,l)=(f_star(ig,l)*ztva(ig,l-1)+alim_star(ig,l)  &
@@ -104 +104 @@
 !------------------------------------------------------------------------
 
-            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1).lt.1.e-10) then
+            if (zw2(ig,l+1)>0. .and. zw2(ig,l+1)<1.e-10) then
 !               stop'On tombe sur le cas particulier de thermcell_dry'
 !               print*,'On tombe sur le cas particulier de thermcell_dry'
@@ -112 +112 @@
             endif
 
-            if (zw2(ig,l+1).lt.0.) then
+            if (zw2(ig,l+1)<0.) then
                linter(ig)=(l*(zw2(ig,l+1)-zw2(ig,l))  &
      &           -zw2(ig,l))/(zw2(ig,l+1)-zw2(ig,l))
@@ -119 +119 @@
 !            endif
 !CR:zmax continu 06/05/12: calcul de linter quand le thermique est stoppe par le detrainement
-            elseif (f_star(ig,l+1).lt.0.) then
+            elseif (f_star(ig,l+1)<0.) then
                linter(ig)=(l*(f_star(ig,l+1)-f_star(ig,l))  &
      &           -f_star(ig,l))/(f_star(ig,l+1)-f_star(ig,l))
@@ -128 +128 @@
                wa_moy(ig,l+1)=sqrt(zw2(ig,l+1))
 
-            if (wa_moy(ig,l+1).gt.wmaxa(ig)) then
+            if (wa_moy(ig,l+1)>wmaxa(ig)) then
 !   lmix est le niveau de la couche ou w (wa_moy) est maximum
                lmix(ig)=l+1
@@ -135 +135 @@
          enddo
       enddo
-       if (prt_level.ge.1) print*,'fin calcul zw2'
+       if (prt_level>=1) print*,'fin calcul zw2'
 !
 ! Determination de zw2 max
@@ -144 +144 @@
       do l=1,nlay
          do ig=1,ngrid
-            if (l.le.lmax(ig)) then
+            if (l<=lmax(ig)) then
                 zw2(ig,l)=sqrt(zw2(ig,l))
                 wmax(ig)=max(wmax(ig),zw2(ig,l))