from rsf.proj import *
## If you change the order of lines here, do not forget to modify ../rsfplot.tex
# Graph, x linear example
Flow('line',None,'math n1=10 output="x1"')
Result('line','graph title=Line')

# Graph, x squared example
Flow('xsq',None,'math n1=10 output="x1*x1"')
Flow('sqline','xsq line','cmplx ${SOURCES[1]}')
Result('sqline','''
       graph symbol="+" title="Line, x\\^2" symbolsz=12''')

# Multiple lines example
Flow('mline',None,'math n1=10 n2=3 output="(x2+1)*x1"')
Result('mlines','mline','graph title="Multiple lines, n3 > 1"')
Result('mplines','mline','''
       graph title="With symbol=" symbol=+*o symbolsz=8''')

# Dots example
Result('dplines','mline','''
       dots title="Lines with bars" gaineach=n radius=0.15''')

# Box example
Plot('label1',None,'''
     box y0=5.6 x0=8 xt=-1 yt=1 label="Line \F10 a"''')
Plot('mline','mline','graph title="Lines with a label (overlay)"')
Result('mblines','mline label1','Overlay')

# Surface and contours
Flow('sinsq',None,'''
     math n1=201 n2=201 o1=-1.0 o2=-1.0 d1=0.01 d2=0.01
          output="sin(3.14*x1)*sin(3.14*x2)" label1=x label2=y''')
Result('sqsurf','sinsq','grey color=j title=Surface scalebar=y')
Result('csqsurf','sinsq','contour title=Isolines allpos=n')

# Surface and contours, overlay
Plot('psinsq','sinsq','''grey color=j title=Overlay scalebar=y
                              min1=-1 max1=1 min2=-1 max2=1''')
Plot('pcsinsq','sinsq','contour title=Overlay allpos=n scalebar=y')
Result('csurf','psinsq pcsinsq','Overlay')

# Wiggle example
Flow('signal',None,'''spike nsp1 n1=301 d1=0.004 o1=0.0 k1=101 l1=101 title=
                            label1=Time unit1=s | ricker1 frequency=15 |
                      spray axis=2 d=0.01 n=51 o=0 label=Distance unit=km''')
Flow('vel','signal','math output=2')
Flow('gath',['signal','vel'],'inmo velocity=${SOURCES[1]}')
Plot('wgath','gath','''wiggle poly=y min2=0 max2=0.5 plotcol=7
                       yreverse=y wherexlabel=top gridcol=7''')
Plot('rgath','gath','grey transp=n min2=0 max2=0.5')
Result('sgath','rgath wgath','Overlay')

# Cubeplot example
Flow('ganom',None,'''math n1=41 n2=41 n3=41 o1=-1 o2=-1 o3=-1
     d1=0.05 d2=0.05 d3=0.05 output="exp(-2*(x1^2+x2^2+x3^2))"''')
Result('fganom','ganom','''byte gainpanel=all bar=gbar.rsf
       allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf
       title="Slices, flat" frame1=20 frame2=20 frame3=20''')

# Cubeplot, non-flat example
Result('ganom3','ganom','''byte gainpanel=all bar=gbar.rsf
       allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf
       title="Cube" frame1=20 frame2=20 frame3=20 flat=n''')

# SideBySide example
Plot('pfganom','ganom','''byte gainpanel=all bar=gbar.rsf
     allpos=y | cubeplot color=j scalebar=y bar=gbar.rsf
     title="Surface" frame1=20 frame2=20 frame3=20''')
Plot('pcganom','ganom','''contour3 color=j scalebar=y
     title=Isolines nc=30 frame1=20 frame2=20 frame3=20''')
Result('sganom','pfganom pcganom','SideBySideAniso')

# PLplot example
Result('sqplsin','sinsq','''
       plsurf mesh=n color=e title="PLplot surface"''')

End()