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

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

@@ -57 +57 @@
             zzm=masse(ig,k)/ptimestep
             cfl=max(cfl,fm(ig,k)/zzm)
-            if (entr(ig,k).gt.zzm) then
+            if (entr(ig,k)>zzm) then
                print*,'entr*dt>m,1',k,entr(ig,k)*ptimestep,masse(ig,k)
                abort_message = 'entr dt > m, 1st'
@@ -68 +68 @@
 
 
-      if (prt_level.ge.1) print*,'Q2 THERMCEL_DQ 0'
+      if (prt_level>=1) print*,'Q2 THERMCEL_DQ 0'
 
 !   calcul du detrainement
@@ -76 +76 @@
 !           print*,'Q2 DQ ',detr(ig,k),fm(ig,k),entr(ig,k)
 !test
-            if (detr(ig,k).lt.0.) then
+            if (detr(ig,k)<0.) then
                entr(ig,k)=entr(ig,k)-detr(ig,k)
                detr(ig,k)=0.
@@ -82 +82 @@
 !     s         'f+1=',fm(ig,k+1),'e=',entr(ig,k),'d=',detr(ig,k)
             endif
-            if (fm(ig,k+1).lt.0.) then
+            if (fm(ig,k+1)<0.) then
 !               print*,'fm2<0!!!'
             endif
-            if (entr(ig,k).lt.0.) then
+            if (entr(ig,k)<0.) then
 !               print*,'entr2<0!!!'
             endif
@@ -98 +98 @@
       do k=2,nlay
          do ig=1,ngrid
-            if ((fm(ig,k+1)+detr(ig,k))*ptimestep.gt.  &
+            if ((fm(ig,k+1)+detr(ig,k))*ptimestep>  &
      &         1.e-5*masse(ig,k)) then
          qa(ig,k)=(fm(ig,k)*qa(ig,k-1)+entr(ig,k)*q(ig,k))  &
@@ -105 +105 @@
                qa(ig,k)=q(ig,k)
             endif
-            if (qa(ig,k).lt.0.) then
+            if (qa(ig,k)<0.) then
 !               print*,'qa<0!!!'
             endif
-            if (q(ig,k).lt.0.) then
+            if (q(ig,k)<0.) then
 !               print*,'q<0!!!'
             endif
@@ -200 +200 @@
             zzm=masse(ig,k)/ptimestep
             cfl=max(cfl,fm(ig,k)/zzm)
-            if (entr(ig,k).gt.zzm) then
+            if (entr(ig,k)>zzm) then
                print*,'entr*dt>m,2',k,entr(ig,k)*ptimestep,masse(ig,k)
                abort_message = 'entr dt > m, 2nd'
@@ -223 +223 @@
 
 do iter=1,niter
-      if (prt_level.ge.1) print*,'Q2 THERMCEL_DQ 0'
+      if (prt_level>=1) print*,'Q2 THERMCEL_DQ 0'
 
 !   calcul du detrainement
@@ -231 +231 @@
 !           print*,'Q2 DQ ',detr(ig,k),fm(ig,k),entr(ig,k)
 !test
-            if (detr(ig,k).lt.0.) then
+            if (detr(ig,k)<0.) then
                entr(ig,k)=entr(ig,k)-detr(ig,k)
                detr(ig,k)=0.
@@ -237 +237 @@
 !     s         'f+1=',fm(ig,k+1),'e=',entr(ig,k),'d=',detr(ig,k)
             endif
-            if (fm(ig,k+1).lt.0.) then
+            if (fm(ig,k+1)<0.) then
 !               print*,'fm2<0!!!'
             endif
-            if (entr(ig,k).lt.0.) then
+            if (entr(ig,k)<0.) then
 !               print*,'entr2<0!!!'
             endif
@@ -253 +253 @@
       do k=2,nlay
          do ig=1,ngrid
-            if ((fm(ig,k+1)+detr(ig,k))*ztimestep.gt.  &
+            if ((fm(ig,k+1)+detr(ig,k))*ztimestep>  &
      &         1.e-5*masse(ig,k)) then
          qa(ig,k)=(fm(ig,k)*qa(ig,k-1)+entr(ig,k)*q(ig,k))  &
@@ -260 +260 @@
                qa(ig,k)=q(ig,k)
             endif
-            if (qa(ig,k).lt.0.) then
+            if (qa(ig,k)<0.) then
 !               print*,'qa<0!!!'
             endif
-            if (q(ig,k).lt.0.) then
+            if (q(ig,k)<0.) then
 !               print*,'q<0!!!'
             endif
@@ -292 +292 @@
 #endif
 
-            if (wqd(ig,k).lt.0.) then
+            if (wqd(ig,k)<0.) then
 !               print*,'wqd<0!!!'
             endif