from rsfproj 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()