Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/so2_to_h2so4.F90

Timestamp:

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

Author:

abarral

Message:

(lint) Fix obsolete boolean operators

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/so2_to_h2so4.F90 (modified) (5 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/so2_to_h2so4.F90

@@ -37 +37 @@
      !only in the stratosphere
         IF (is_strato(ilon,ilev)) THEN
-           IF (SO2_lifetime(ilon,ilev).GT.0.0 .AND. SO2_lifetime(ilon,ilev).LT.1.E10) THEN
+           IF (SO2_lifetime(ilon,ilev)>0.0 .AND. SO2_lifetime(ilon,ilev)<1.E10) THEN
               IF (flag_OH_reduced) THEN
                  !--convert SO2 to H2SO4 (slimane)
@@ -83 +83 @@
                       & *pplay(ilon,ilev)/t_seri(ilon,ilev)/RD/1.E6/mSO2mol
                  
-                 IF (rrak1 .GE. 0.0) THEN
+                 IF (rrak1 >= 0.0) THEN
                     rreduce =( (rkho2o3+rkoho3)*rrak0 + rkso2oh*rrak1 ) / &
                          (  (rkho2o3+rkoho3)*rrak0 )
@@ -100 +100 @@
               ! Check lifetime rreduce < timestep*3 (such as SO2 loss > 0.28*SO2) with exp(-1/3)=0.72
               IF(flag_min_rreduce) THEN
-                 IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
+                 IF (rreduce < (3.*pdtphys)) rreduce = 3.*pdtphys
               ENDIF
               budg_3D_so2_to_h2so4(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/rreduce))
@@ -122 +122 @@
               ! test: quick sequential integration for SO2 to H2SO4 and reverse
               
-              IF (H2SO4_lifetime(ilon,ilev).GT.0.0 .AND. H2SO4_lifetime(ilon,ilev).LT.1.E10) THEN
+              IF (H2SO4_lifetime(ilon,ilev)>0.0 .AND. H2SO4_lifetime(ilon,ilev)<1.E10) THEN
                  
                  rreduce = H2SO4_lifetime(ilon,ilev)
@@ -130 +130 @@
                  ! Check lifetime rreduce < timestep*3 (such as H2SO4 loss > 0.28*H2SO4) with exp(-1/3)=0.72
                  IF(flag_min_rreduce) THEN
-                    IF (rreduce .LT. (3.*pdtphys)) rreduce = 3.*pdtphys
+                    IF (rreduce < (3.*pdtphys)) rreduce = 3.*pdtphys
                  ENDIF
                  dummyso4toso2 = (mSO2mol/mH2SO4mol)*tr_seri(ilon,ilev,id_H2SO4_strat)*(1.0-exp(-pdtphys/rreduce))