from rsf.proj import *

Fetch('input.HH','pwd')

Flow('input','input.HH',
     'dd form=native | put label2=Offset unit2=km label1=Time unit1=s')
Plot('input','grey wanttitle=n')

Flow('mask','input',
     '''
     math output="((x1-1.4)/0.3)^2+((x2-0.7)/0.4)^2" |
     mask min=1 | dd type=float
     ''')
Flow('hole','input mask','add mode=p ${SOURCES[1]}')
Plot('hole','grey wanttitle=n')

Flow('dip','hole mask','dip mask=${SOURCES[1]} rect1=50 rect2=10 order=3')

Flow('fill','hole mask dip',
     'planemis2 dip=${SOURCES[2]} mask=${SOURCES[1]} verb=y order=3')
Plot('fill','grey wanttitle=n')

Result('combo','input hole fill','SideBySideAniso')

End()