up [pdf]
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')

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

sfdd
sfwindow
sfmask
sfheaderwindow
sfmath
sfbin
sflapfill
sfigrad
sfgrey
sfwilson
sfhelicon

data/galilee/galilee.h