from rsf.proj import *

Fetch('dynamics.H','toldi')

Result('toldi','dynamics.H',
       '''
       window n1=375 n2=1 f2=9 |
       dd form=native type=float |
       grey wanttitle=n allpos=y
       ''')

Flow('cmp','dynamics.H',
     '''
     window f3=5 n3=1 |
     dd form=native type=float |
     add add=-128 |
     nmostretch v0=1.5 half=n inv=y
     ''')

near = [0,5,9]
far = map(lambda x: 27-x,near)

for f2 in near+far:
    win = 'win%d' % f2
    Flow(win,'cmp','window f2=%d n2=1' % f2)

for case in range(2):
    plot = 'case%d' % case
    Plot(plot,
         'win%d win%d win%d' % map(tuple,[near,far])[case],
         '''
         cat axis=2 ${SOURCES[1:3]} |
         wiggle wanttitle=n plotcol=%d zplot=.4
         wantaxis2=n xll=0 xur=13 yll=.5 yur=10
         ''' % (6-case))

Result('reciptrace','case0 case1','Overlay')

for s in range(2):
    slice = 'slice%d' % s
    Plot(slice,'dynamics.H',
         '''
         window min1=%g n1=1 |
         dd form=native type=float |
         grey wanttitle=n allpos=y
         ''' % (1,2.5-0.004)[s])

Result('recipslice','slice0 slice1','OverUnderAniso')

End()