from rsf.proj import *

Fetch('hector.HH','ground')

Flow('dat','hector.HH','dd form=native | put o2=0.017')

Result('radialdat','dat',
       '''
       grey label2="Offset (km)" label1="time (s)" wanttitle=n
       crowd1=0.83 clip=300
       ''')

for nw in (2,4,8):
    rad = 'rad%d' %nw
    Flow(rad,'dat',
         'transp | radial vmin=0 vmax=0.65 nv=250 tp=0.2 nw=%d | transp' %
         nw)
    Plot(rad,
         '''
         grey label2="Radial velocity (km/s)" label1="time (s)"
         crowd1=0.83 clip=300 title="B-%d interpolation"
         ''' % nw)
    err = 'err%d' % nw
    Flow(err,[rad,'dat'],
         '''
         transp | radial tp=0.2 nw=%d inv=y | transp |
         add scale=1,-1 ${SOURCES[1]}
         ''' % nw)
    Plot(err,
         '''
         grey label2="Offset (km)" label1="time (s)" wanttitle=n
         crowd1=0.83 clip=300 title="B-%d error"
         ''' % nw)

Result('radial','rad2 rad4 err2 err4','TwoRows')

End()