# 1 "" # 1 "" # 1 "" # 1 "" module fixbad ! Given a PEF, find bad data and restore it. use mis2_mod use helicon use quantile_mod implicit none contains subroutine fixbad1 (niter, aa, yy) integer, intent (in) :: niter type( filter), intent (in) :: aa real, dimension (:) :: yy real, dimension (size (yy)) :: rr logical, dimension (size (yy)) :: known integer :: stat call helicon_init( aa) stat = helicon_lop (.false., .false., yy, rr) rr = abs (rr) known = ( yy > 0.) .and. ( rr < 4. * quantile( size(rr)/2, rr)) call mis2 (niter, yy, aa, known, 2) end subroutine end module