from rsf.proj import *
import math

for case in (0,1):
    spike = 'spike%d' % case
    Flow(spike,None,
         '''
         spike n1=884 n2=201 d1=0.0025 d2=10 o1=0 o2=-1000
         label1=time unit1=s label2=midpoint unit2=m
         nsp=3 k1=%s k2=101,101,101 | bandpass fhi=60
         ''' % ('201,401,601','301,501,701')[case])
    rdmo = 'rdmo%d' % case
    Flow(rdmo,spike,
         '''
         dmo h=1000 velhalf=200 n=200 type=2 adj=n |
         spray axis=2 n=1 d=1 o=1000 |
         nmostretch v0=%g inv=%d half=n |
         window |
         dmo h=1000 velhalf=200 n=200 type=2 adj=y
         ''' % ((1./math.sqrt(1./(1000*1000)-1./(1333*1333)),
                 1./math.sqrt(1./(833*833)-1./(1000*1000)))[case],
                1-case))
    Plot(rdmo,
         '''
         window f2=50 n2=101 n1=800 |
         grey pclip=100 
         title="offset=1000 m gamma=%g"
         ''' % (1.333,0.833)[case])

Result('vlccps','rdmo0 rdmo1','SideBySideAniso')


End()