Changeset 5185 for LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp_stats.F90

Timestamp:

Sep 11, 2024, 4:27:07 PM (2 months ago)

Author:

abarral

Message:

Replace REPROBUS CPP KEY by logical using handmade wonky wrapper

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp_stats.F90 (modified) (3 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp_stats.F90

@@ -217 +217 @@
           DO j=1,Nlevels
              sc_ratio = beta_tot(pr,i,j)/beta_mol(pr,j)
-             if ((sc_ratio .le. s_att) .and. (flag_sat .eq. 0)) flag_sat = j
+             if ((sc_ratio .le. s_att) .AND. (flag_sat .eq. 0)) flag_sat = j
              if (Ze_tot(pr,i,j) .lt. -30.) then  !radar can't detect cloud
                 if ( (sc_ratio .gt. s_cld) .or. (flag_sat .eq. j) ) then  !lidar sense cloud
@@ -267 +267 @@
     
     DO ij=2,Nbins+1
-       hist1D(ij-1) = count(var .ge. bins(ij-1) .and. var .lt. bins(ij))
+       hist1D(ij-1) = count(var .ge. bins(ij-1) .AND. var .lt. bins(ij))
        if (count(var .eq. R_GROUND) .ge. 1) hist1D(ij-1)=R_UNDEF
     enddo
@@ -300 +300 @@
     DO ij=2,nbin1+1
        DO ik=2,nbin2+1
-          jointHist(ij-1,ik-1)=count(var1 .ge. bin1(ij-1) .and. var1 .lt. bin1(ij) .and. &
-               var2 .ge. bin2(ik-1) .and. var2 .lt. bin2(ik))        
+          jointHist(ij-1,ik-1)=count(var1 .ge. bin1(ij-1) .AND. var1 .lt. bin1(ij) .AND. &
+               var2 .ge. bin2(ik-1) .AND. var2 .lt. bin2(ik))
        enddo
     enddo