from rsf.proj import *
import fdmod,gfield
par = {
'nx':1200, 'ox':0, 'dx':1.0, 'lx':'x','ux':'m',
'nz':800, 'oz':0, 'dz':1.0, 'lz':'z','uz':'m',
'ru':30,'rv':6,'aa':1,'ff':45
}
fdmod.param(par)
gfield.execute('gg',1,par)
Result('gg',fdmod.cgrey('color=E',par))
Result('ggz','gg','stack norm=y axis=1 | graph')
Result('ggx','gg','stack norm=y axis=2 | graph')
Flow( 'vo','gg','math output="1500+x1"')
Result('vo',fdmod.cgrey('bias=1500 allpos=y',par))
Flow('nu','gg','scale rscale=0.1')
Flow('vv','vo nu',
'''
math v=${SOURCES[0]} n=${SOURCES[1]}
output="v/sqrt(1+n)"
''' % par)
Result('vv',fdmod.cgrey('bias=1500 allpos=y',par))
Flow('ivo','vo','math output=1/input^2 | stack')
Flow('ivv','vv','math output=1/input^2 | stack')
Result('ivall','ivo ivv',
'''
cat axis=2 space=n ${SOURCES[1]} |
graph
''')
End() |