module mis_mod {  
  use tcai1+cgstep_mod+solver_mod
contains
# fill in missing data on 1-axis by minimizing power out of a given filter.
  subroutine mis1 ( niter, mm, aa) {
    integer,             intent (in)      :: niter   # number of iterations
    real, dimension (:), pointer          :: aa      # roughening filter
    real, dimension (:), intent (in out)  :: mm      # in - data with zeroes
                                                     # out - interpolated
    real, dimension (size(mm) + size(aa)) :: zero    # filter output

    call tcai1_init( aa, size( mm));         zero = 0.
		#%<  zero = tcail_lop mm  : x0=mm :  known=(mm!=0.) : niter=niter
    call cgstep_close ()
    }
}