Masking fix in cosp_diag_warmrain? about cospv2.0 HOT 4 CLOSED

The wrinkle to this is cloud-free scenes and dark scenes are both set to R_UNDEF in the modis-simulator. We want to include to former, but not the latter.

from cospv2.0.

@dustinswales I apporogize for this bug and late reply. The warm rain diagnostics did intend to use only daytime cloudy scenes, but actually the current code did not explicitly exclude dark-scene columns from the analysis.
To fix this bug, I will change the following two places as suggested by @brhillman and @dustinswales at the previous PR from 2019 (#36).

https://github.com/CFMIP/COSPv2.0/blob/master/src/cosp_stats.F90#L338

!! initialize
 do i = 1, Npoints
    if ( lwp(i) .eq. R_UNDEF ) then  ! for non-sunlit columns
       cfodd_ntotal(i,:,:,:) = R_UNDEF
       wr_occfreq_ntotal(i,:) = R_UNDEF
       icod(i,:,:) = R_UNDEF
    else
       cfodd_ntotal(i,:,:,:)  = 0._wp
       wr_occfreq_ntotal(i,:) = 0._wp
       icod(i,:,:) = 0._wp
    endif
 enddo

https://github.com/CFMIP/COSPv2.0/blob/master/src/cosp_stats.F90#L345

    if ( ( lwp(i)   .le.  CWP_THRESHOLD  .and. lwp(i)  .ne.  R_UNDEF  ) .or.  &

I think this change will work appropriately for the issue. I push the bug-fix later via pull request.

from cospv2.0.

@dustinswales I think that we can now close this thread, as we have corrected this bug in the recent PR (#73). Thank you for your support.

from cospv2.0.

@takmichibata
Excellent!
Thanks for your work on this!

from cospv2.0.