pro radar_mask_post_mmcr common mask_variables, dif_sig, count_pro, levels_t, mask mask2=mask ; assume the SNR comes in. set mask=1 where snr > -20 dbz and mask=-1 where snr < -50 for l=0,1 do begin for j=0,count_pro-1 do begin for k=0,levels_t-1 do begin if mask[j,k] eq 1 then begin count=0 for jj=j-1,j+1 do begin for kk=k-1,k+1 do begin if jj ge 0 and jj lt count_pro-1 and kk ge 0 and kk le levels_t-1 then begin if mask[jj,kk] eq 1 then count=count+1 endif endfor endfor if count le 5 then begin mask2[j,k]=0 endif endif endfor endfor endfor mask=mask2 end