@make_skill_metrics
pro make_mask_metrics, phase_true, binary_mask, pc, missed, false, skill


;--- make input binary (0/1) masks
idx = where(finite(phase_true) and finite(binary_mask) and phase_true ge 0 and binary_mask ge 0, cc)
if (cc gt 0) then begin
   forecast = binary_mask[idx]
   observed = phase_true[idx]
endif else begin
   print, "ERROR:  No valid data in make_mask_metrics, stopping"
endelse

idx = where(phase_true gt 0,cc) & if (cc gt 0) then observed[idx] = 1

make_skill_metrics, observed, forecast, pc, missed, false, skill

end