from rsf.proj import *
from galilee import Galilee
from helderiv import Helderiv

def plot(title):
    return '''
    grey pclip=94 title="%s" yreverse=n transp=n
    label2="north-south (km)" label1="west-east (km)"
    ''' % title

Galilee('mesh')

Flow('fill','mesh','lapfill grad=y niter=20')
Plot('ruf2','fill','igrad | ' + plot('Filled and d/dx'))

Helderiv('hel',na=128,eps=0.00001)
Plot('held','fill hel',
     'helicon filt=${SOURCES[1]} | ' + plot('Filled and helix derivative'))

Result('helgal','ruf2 held','SideBySideAniso')

Flow('sfill','mesh','shapefill verb=y rect1=2 rect2=2 niter=20')
Plot('sfill','igrad | ' + plot('Filled and d/dx (Shaping)'))


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