- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp_stats.F90
r3358 r5082 203 203 do j=1,Nlevels 204 204 sc_ratio = beta_tot(pr,i,j)/beta_mol(pr,j) 205 if ((sc_ratio .le. s_att) .and. (flag_sat .eq.0)) flag_sat = j206 if (Ze_tot(pr,i,j) .lt.-30.) then !radar can't detect cloud207 if ( (sc_ratio .gt. s_cld) .or. (flag_sat .eq.j) ) then !lidar sense cloud205 if ((sc_ratio <= s_att) .and. (flag_sat == 0)) flag_sat = j 206 if (Ze_tot(pr,i,j) < -30.) then !radar can't detect cloud 207 if ( (sc_ratio > s_cld) .or. (flag_sat == j) ) then !lidar sense cloud 208 208 lidar_only_freq_cloud(pr,j)=lidar_only_freq_cloud(pr,j)+1. !top->surf 209 209 flag_cld=1 … … 213 213 endif 214 214 enddo !levels 215 if (flag_cld .eq.1) tcc(pr)=tcc(pr)+1._wp215 if (flag_cld == 1) tcc(pr)=tcc(pr)+1._wp 216 216 enddo !columns 217 217 enddo !points … … 245 245 246 246 do ij=2,Nbins+1 247 hist1D(ij-1) = count(var .ge. bins(ij-1) .and. var .lt.bins(ij))248 if (count(var .eq. R_GROUND) .ge.1) hist1D(ij-1)=R_UNDEF247 hist1D(ij-1) = count(var >= bins(ij-1) .and. var < bins(ij)) 248 if (count(var == R_GROUND) >= 1) hist1D(ij-1)=R_UNDEF 249 249 enddo 250 250 … … 278 278 do ij=2,nbin1+1 279 279 do ik=2,nbin2+1 280 jointHist(ij-1,ik-1)=count(var1 .ge. bin1(ij-1) .and. var1 .lt.bin1(ij) .and. &281 var2 .ge. bin2(ik-1) .and. var2 .lt. bin2(ik))280 jointHist(ij-1,ik-1)=count(var1 >= bin1(ij-1) .and. var1 < bin1(ij) .and. & 281 var2 >= bin2(ik-1) .and. var2 < bin2(ik)) 282 282 enddo 283 283 enddo
Note: See TracChangeset
for help on using the changeset viewer.