Ignore:
Timestamp:
Jul 19, 2024, 5:41:58 PM (16 months ago)
Author:
abarral
Message:

(lint) Fix obsolete boolean operators

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/inlandsis.F

    r3900 r5082  
    409409        IF (BloMod) THEN
    410410 
    411         if (klonv.eq.1) then
    412           if(isnoSV(1).ge.2                   .and.
     411        if (klonv==1) then
     412          if(isnoSV(1)>=2                   .and.
    413413     .         TsisSV(1,max(1,isnoSV(1)))<273.  .and.
    414414     .         ro__SV(1,max(1,isnoSV(1)))<500.  .and.
    … …  
    518518C +     Gallee et al., 2001    eq 5, p5
    519519 
    520           if (usth_param .eq. "gal") then
     520          if (usth_param == "gal") then
    521521            Salt_us   =   (log(2.868) - log(1 + SaltMo)) * rCd10n/0.085
    522522            Salt_us   = Salt_us * Fac_Mo
    … …  
    525525          endif
    526526 
    527           if (usth_param .eq. "lis") then !Liston et al. 2007
     527          if (usth_param == "lis") then !Liston et al. 2007
    528528            if(ro__SV(ikl,isn)>300.) then
    529529              Salt_us   = 0.005*exp(0.013*ro__SV(ikl,isn))
    … …  
    607607        hSalSV(ikl) = 8.436e-2  * us__SV(ikl)**SblPom
    608608 
    609         if (qsalt_param .eq. "pom") then
     609        if (qsalt_param == "pom") then
    610610          qSalSV(ikl) = (us__SV(ikl)**2 - usthSV(ikl)**2) *signus
    611611     .               / (hSalSV(ikl) * gravit * us__SV(ikl) * 3.25)
    612612        endif
    613613 
    614         if (qsalt_param .eq. "bin") then
     614        if (qsalt_param == "bin") then
    615615          qSalSV(ikl) = (us__SV(ikl) * us__SV(ikl)
    616616     .                -usthSV(ikl) * usthSV(ikl))*signus
    … …  
    668668C + ...         [ A compromise between
    669669C + ...           Kotlyakov (1961) and Lenaerts (2012, JGR, Part1) ]
    670                 if (tt_tmp.ge.-10) then
     670                if (tt_tmp>=-10) then
    671671                  BufPro   =  max( rosMin,
    672672     .            104. *sqrt( max( vv_tmp-6.0,0.0))) ! Kotlyakov (1961)
    … …  
    894894
    895895         
    896       if (discret_xf.AND.klonv.eq.1) then
    897 
    898        if(isnoSV(1).ge.1.or.NLaysv(1).ge.1) then
     896      if (discret_xf.AND.klonv==1) then
     897
     898       if(isnoSV(1)>=1.or.NLaysv(1)>=1) then
    899899C +          **********
    900900         call SISVAT_zSn
    … …  
    10811081
    10821082
    1083       if (iflag_temp_inlandsis .eq. 0) then
     1083      if (iflag_temp_inlandsis == 0) then
    10841084
    10851085       call SISVAT_TSo
    … …  
    11981198! ^^^^^^^^^^^^^^^^^^^^
    11991199
    1200           IF (iflag_tsurf_inlandsis .EQ. 0) THEN   
     1200          IF (iflag_tsurf_inlandsis == 0) THEN
    12011201
    12021202            Tsrfsv(ikl) =TsisSV(ikl,isnoSV(ikl))
    12031203
    1204           ELSE IF (iflag_tsurf_inlandsis .GT. 0) THEN
     1204          ELSE IF (iflag_tsurf_inlandsis > 0) THEN
    12051205! Etienne: extrapolation from the two uppermost levels:
    12061206
    … …  
    12081208           zm1=-dzsnSV(ikl,isnoSV(ikl))/2.
    12091209           zm2=-(dzsnSV(ikl,isnoSV(ikl)) + dzsnSV(ikl,isnoSV(ikl)-1)/2.)
    1210          else if (isnoSV(ikl) .EQ. 1) then
     1210         else if (isnoSV(ikl) == 1) then
    12111211           zm1=-dzsnSV(ikl,isnoSV(ikl))/2.
    12121212           zm2=-(dzsnSV(ikl,isnoSV(ikl))+dz_dSV(0)/2.)
    … …  
    12361236      IF (SnoMod)                                                 THEN
    12371237
    1238       if (discret_xf .AND. klonv.eq.1) then
    1239       if(isnoSV(1).ge.1) then
     1238      if (discret_xf .AND. klonv==1) then
     1239      if(isnoSV(1)>=1) then
    12401240C +          **********
    12411241      call SISVAT_GSn
    … …  
    13051305! Calculation of snow roughness length
    13061306!=====================================
    1307           IF (iflag_z0m_snow .EQ. 0) THEN
     1307          IF (iflag_z0m_snow == 0) THEN
    13081308
    13091309          Z0m_Sn=prescribed_z0m_snow
    13101310
    1311           ELSE IF (iflag_z0m_snow .EQ. 1) THEN
     1311          ELSE IF (iflag_z0m_snow == 1) THEN
    13121312
    13131313          Z0m_Sn=Z0enSV(ikl)
    13141314
    1315           ELSE IF (iflag_z0m_snow .EQ. 2) THEN                             
     1315          ELSE IF (iflag_z0m_snow == 2) THEN
    13161316
    13171317C +--Rough   Snow Surface Roughness Length (Variable Sastrugi Height)
    … …  
    13361336          coefd = log(z03)-coefc*ta3
    13371337
    1338           if (TaT_SV(ikl) .lt. ta1) then
     1338          if (TaT_SV(ikl) < ta1) then
    13391339            Z0_obs = z01
    1340           else if (TaT_SV(ikl).ge.ta1 .and. TaT_SV(ikl).lt.ta2) then
     1340          else if (TaT_SV(ikl)>=ta1 .and. TaT_SV(ikl)<ta2) then
    13411341            Z0_obs = exp(coefa*TaT_SV(ikl) + coefb)
    1342           else if (TaT_SV(ikl).ge.ta2 .and. TaT_SV(ikl).lt.ta3) then
     1342          else if (TaT_SV(ikl)>=ta2 .and. TaT_SV(ikl)<ta3) then
    13431343            ! if st > 0, melting induce smooth surface
    13441344            Z0_obs = exp(coefc*TaT_SV(ikl) + coefd)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_gsn.F

    r3792 r5082  
    805805        DO ikl=1,knonv
    806806        DO isn=iiceSV(ikl),isnoSV(ikl)
    807           IF  (G1snSV(ikl,isn).ge.0.)                               THEN
    808             IF(G1snSV(ikl,isn).lt.vsphe4.and.istoSV(ikl,isn).eq.0)  THEN
     807          IF  (G1snSV(ikl,isn)>=0.)                               THEN
     808            IF(G1snSV(ikl,isn)<vsphe4.and.istoSV(ikl,isn)==0)  THEN
    809809                   istoSV(ikl,isn)=istdSV(1)
    810             ELSEIF(G1_dSV-G1snSV(ikl,isn)         .lt.vsphe4.and.
    811      .             etaSno(ikl,isn)/dzsnSV(ikl,isn).gt.vtelv1)       THEN
    812               IF  (istoSV(ikl,isn).eq.0)
     810            ELSEIF(G1_dSV-G1snSV(ikl,isn)         <vsphe4.and.
     811     .             etaSno(ikl,isn)/dzsnSV(ikl,isn)>vtelv1)       THEN
     812              IF  (istoSV(ikl,isn)==0)
    813813     .             istoSV(ikl,isn)=   istdSV(2)
    814               IF  (istoSV(ikl,isn).eq.istdSV(1))
     814              IF  (istoSV(ikl,isn)==istdSV(1))
    815815     .             istoSV(ikl,isn)=   istdSV(3)
    816             ELSEIF(TsisSV(ikl,isn).lt.TfSnow)                       THEN
    817               IF  (istoSV(ikl,isn).eq.istdSV(2))
     816            ELSEIF(TsisSV(ikl,isn)<TfSnow)                       THEN
     817              IF  (istoSV(ikl,isn)==istdSV(2))
    818818     .             istoSV(ikl,isn)=   istdSV(4)
    819               IF  (istoSV(ikl,isn).eq.istdSV(3))
     819              IF  (istoSV(ikl,isn)==istdSV(3))
    820820     .             istoSV(ikl,isn)=   istdSV(5)
    821821            END IF

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_qsn.F

    r3900 r5082  
    438438      END DO
    439439
    440       IF  (isnitr.GT.0)                                       GO TO 1000
     440      IF  (isnitr>0)                                       GO TO 1000
    441441 
    442442 

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_sno_albedo.F

    r3980 r5082  
    218218
    219219
    220           IF (iflag_albcalc .GE. 1) THEN  ! Albedo calculation according to Kokhanovsky and Zege 2004
     220          IF (iflag_albcalc >= 1) THEN  ! Albedo calculation according to Kokhanovsky and Zege 2004
    221221
    222222          dalbed = 0.0
    … …  
    752752          alb30 = exp(-(alpha**0.5)* 3.0 / 7.0 * (1.0 + 20 * cos30))
    753753
    754           IF ((wavelengths(i).GE.lambdamin).AND.
    755      .       (wavelengths(i).LT.lambdamax)) THEN
     754          IF ((wavelengths(i)>=lambdamin).AND.
     755     .       (wavelengths(i)<lambdamax)) THEN
    756756          albint=albint+R*Pas  ! rectangle integration -> can be improved with trapezintegration
    757757          norm=norm+Pas

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_ts2.F

    r3792 r5082  
    207207
    208208      DO ig=1,knonv   
    209         IF (isnoSV(ig).GT.0) THEN
    210           IF (isnoSV(ig).GT.1) THEN
     209        IF (isnoSV(ig)>0) THEN
     210          IF (isnoSV(ig)>1) THEN
    211211           mug(ig)=1./(1.+dzsnSV(ig,isnoSV(ig)-1)/dzsnSV(ig,isnoSV(ig))) !mu
    212212          ELSE
    … …  
    217217        ENDIF
    218218
    219         IF (mug(ig)  .LE. 0.05) THEN
     219        IF (mug(ig)  <= 0.05) THEN
    220220         write(*,*)'Attention mu low', mug(ig)
    221221        ENDIF
    222         IF (mug(ig)  .GE. 0.98) THEN
     222        IF (mug(ig)  >= 0.98) THEN
    223223         write(*,*)'Attention mu high', mug(ig)
    224224        ENDIF
    … …  
    331331
    332332      DO ig = 1,knonv
    333         IF (ps__SV(ig).LT.1.) THEN
     333        IF (ps__SV(ig)<1.) THEN
    334334!          write(*,*)'ig',ig,'ps',ps__SV(ig)
    335335          ps__SV(ig)=max(ps__SV(ig),1.e-8)
    336336        ENDIF
    337         IF (p1l_SV(ig).LT.1.) THEN
     337        IF (p1l_SV(ig)<1.) THEN
    338338!          write(*,*)'ig',ig,'p1l',p1l_SV(ig)
    339339          p1l_SV(ig)=max(p1l_SV(ig),1.e-8)
    340340        ENDIF
    341         IF (TaT_SV(ig).LT.180.) THEN
     341        IF (TaT_SV(ig)<180.) THEN
    342342!          write(*,*)'ig',ig,'TaT',TaT_SV(ig)
    343343          TaT_SV(ig)=max(TaT_SV(ig),180.)
    344344        ENDIF
    345         IF (QaT_SV(ig).LT.1.e-8) THEN
     345        IF (QaT_SV(ig)<1.e-8) THEN
    346346!          write(*,*)'ig',ig,'QaT',QaT_SV(ig)
    347347          QaT_SV(ig)=max(QaT_SV(ig),1.e-8)
    348348        ENDIF
    349         IF (Tsf_SV(ig).LT.100.) THEN
     349        IF (Tsf_SV(ig)<100.) THEN
    350350!          write(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
    351351          Tsf_SV(ig)=max(Tsf_SV(ig),180.)
    352352        ENDIF
    353         IF (Tsf_SV(ig).GT.500.) THEN
     353        IF (Tsf_SV(ig)>500.) THEN
    354354!          write(*,*)'ig',ig,'Tsf',Tsf_SV(ig)
    355355          Tsf_SV(ig)=min(Tsf_SV(ig),400.)
    … …  
    359359!           Tsrf(ig)=max(Tsrf(ig),TaT_SV(ig)-20.)
    360360!         ENDIF
    361          IF (cdH_SV(ig).LT.1.e-10) THEN
     361         IF (cdH_SV(ig)<1.e-10) THEN
    362362!           IF (ig.le.3)   write(*,*)'ig',ig,'cdH',cdH_SV(ig)
    363363           cdH_SV(ig)=.5
    … …  
    380380     &          /LhvH2O / zx_pkh(ig)
    381381        ELSE
    382           IF (Tsf_SV(ig).LT.t_coup) THEN
     382          IF (Tsf_SV(ig)<t_coup) THEN
    383383             zx_qs = qsats(Tsf_SV(ig)) / ps__SV(ig)
    384384             zx_dq_s_dh = dqsats(Tsf_SV(ig),zx_qs)/LhvH2O               &
    … …  
    408408      DO ig = 1,knonv
    409409        zx_sl(ig) = LhvH2O
    410         IF (Tsf_SV(ig) .LT. RTT) zx_sl(ig) = LhsH2O
     410        IF (Tsf_SV(ig) < RTT) zx_sl(ig) = LhsH2O
    411411        zx_k1(ig) = zx_coef(ig)
    412412      ENDDO

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_weq.F

    r3792 r5082  
    5858 
    5959           ikl   = 1
    60       IF          (isnoSV(ikl).gt.iiceSV(ikl))                      THEN
     60      IF          (isnoSV(ikl)>iiceSV(ikl))                      THEN
    6161 
    6262          SnoWEQ = 0.
    … …  
    7171C +  =====================
    7272 
    73       IF        (iiceSV(1).gt.0)                                    THEN
     73      IF        (iiceSV(1)>0)                                    THEN
    7474 
    7575          IceWEQ = 0.

  • LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_zsn.F

    r3900 r5082  
    459459      DO isn=nsno,2,-1
    460460      DO ikl=1,knonv
    461         IF (Agrege(ikl).gt.0..AND.i_thin(ikl).lt.isnoSV(ikl))       THEN
     461        IF (Agrege(ikl)>0..AND.i_thin(ikl)<isnoSV(ikl))       THEN
    462462          staggr          =  min(1,max(0,isn-i_thin(ikl)    -1))
    463463     .                    *  min(1,max(0,    isnoSV(ikl)-isn+2))
Note: See TracChangeset for help on using the changeset viewer.