from rsf.proj import *

Flow('qdome',None,'qdome n1=200 d1=0.008 | window n3=1 f3=45')
#Flow('qdome',None,'qdome n1=200 d1=0.008')
Flow('slope','qdome','dip rect1=5 rect2=5')

Plot('qdome','grey title=Input')
Plot('slope','grey scalebar=y title=Slope')

Result('reflectors','qdome slope','SideBySideAniso')

Flow('mask1',None,
     'spike n1=200 n2=100 nsp=9 k2=40,52,65,10,75,88,15,20,24')
Flow('mask4',None,
     'spike n1=200 n2=100 nsp=4 k2=5,25,40,70 k1=1,1,1,1 l1=103,190,44,187')

for case in ('1','4'):
    mask = 'mask'+case
    well = 'well'+case
    fill = 'fill'+case

    Flow(well,['qdome',mask],'add mode=p ${SOURCES[1]}')
    Plot(well,'grey allpos=y title="Well Logs" ')

    Flow(fill,[well,mask,'slope'],
         '''
         planemis2 dip=${SOURCES[2]} mask=${SOURCES[1]} verb=y |
         bandpass fhi=30
         ''')

    Plot(fill,'grey title=Interpolated')
    Result('combo'+case,['qdome',well,fill],'SideBySideAniso')

End()