from rsfproj import *

Flow('mod3',None,
     '''sigmoid n1=200 d2=.008 n2=200 |
     smooth rect1=3 diff1=1 |
     smooth rect1=3''')
Flow('dki','mod3',
     'preconstkirch inv=y zero=y vel=1.5 nh=1 h0=0 dh=0')

Flow('mod1','mod3','window j2=4')
Flow('mod2','mod3','window j2=2')

for size in range(1,4):
    dki = 'dki%d' % size
    mod = 'mod%d' % size
    Flow(dki,mod,
         'preconstkirch inv=y zero=y vel=1.5 nh=1 h0=0 dh=0')

Flow('kir','dki3',
     'preconstkirch inv=n zero=y vel=1.5')

for case in ['mod3','dki3','kir','mod1','dki1','mod2','dki2']:
    Plot(case,'grey title="%s" labelsz=10' % case)

Result('kfgood','mod3 dki3','SideBySideAniso')
Result('skmig', 'mod3 kir ','SideBySideAniso')

Flow('smod','mod3','spectra all=y | scale axis=1')
Flow('skir','kir','spectra all=y | scale axis=1')
Flow('skir2','skir','pow pow1=1 | scale axis=1')
Result('kirspec','smod skir2 skir',
       '''
       cat axis=2 ${SOURCES[1:3]} |
       dots yreverse=1 labels="m=model:f * recon:recon"
       label1=Frequency unit1=Hz dots=0
       ''')

Result('skmod','mod1 dki1 mod2 dki2','TwoRows')

###########################################################################
End()