Changeset 8527 for trunk/anuga_core/validation_tests/Tests/Simple/deep_wave
- Timestamp:
- Aug 21, 2012, 7:43:27 AM (13 years ago)
- Location:
- trunk/anuga_core/validation_tests/Tests/Simple/deep_wave
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/validation_tests/Tests/Simple/deep_wave
-
Property
svn:ignore
set to
.deps
-
Property
svn:ignore
set to
-
trunk/anuga_core/validation_tests/Tests/Simple/deep_wave/produce_results.py
r8468 r8527 1 2 from anuga.utilities.argparsing import run_validation_script3 1 4 2 5 run_validation_script('run_wave.py') 6 run_validation_script('plotme.py') 3 from fabricate import * 7 4 5 # Get the values of alg and cfl and any other 6 # paarameters set for these validation tests 7 from validation_tests.parameters import * 8 9 10 # Setup the python scripts which produce the output for this 11 # validation test 12 def build(): 13 run('python', 'run_wave.py', '-alg', alg, '-cfl', cfl) 14 run('python', 'plotme.py', '-alg', alg, '-cfl', cfl ) 15 16 def clean(): 17 autoclean() 18 19 main() -
trunk/anuga_core/validation_tests/Tests/Simple/deep_wave/run_wave.py
r8468 r8527 90 90 Bw = anuga.Time_boundary(domain=domain, # Time dependent boundary 91 91 ## Sine wave 92 f =lambda t: [(-amplitude*sin((1./wave_length)*t*2*pi)), 0.0, 0.0])92 function=lambda t: [(-amplitude*sin((1./wave_length)*t*2*pi)), 0.0, 0.0]) 93 93 ## Sawtooth? 94 # f =lambda t: [(-8.0*(sin((1./180.)*t*2*pi))+(1./2.)*sin((2./180.)*t*2*pi)+(1./3.)*sin((3./180.)*t*2*pi)), 0.0, 0.0])94 # function=lambda t: [(-8.0*(sin((1./180.)*t*2*pi))+(1./2.)*sin((2./180.)*t*2*pi)+(1./3.)*sin((3./180.)*t*2*pi)), 0.0, 0.0]) 95 95 ## Sharp rise, linear fall 96 # f =lambda t: [(5.0*(-((t-0.)/300.)*(t<300.)-cos((t-300.)*2.*pi*(1./240.))*(t>=300. and t<420.)+(1.-(t-420.)/300.)*(t>=420. and t <720.))), 0.0, 0.0])97 # f =lambda t: [amplitude*(1.-2.*(pi*(1./720.)*(t-720.))**2)/exp((pi*(1./720.)*(t-720.))**2) , 0.0, 0.0])98 # f =lambda t: [(-8.0*sin((1./720.)*t*2*pi))*((t<720.)-0.5*(t<360.)), 0.0, 0.0])96 # function=lambda t: [(5.0*(-((t-0.)/300.)*(t<300.)-cos((t-300.)*2.*pi*(1./240.))*(t>=300. and t<420.)+(1.-(t-420.)/300.)*(t>=420. and t <720.))), 0.0, 0.0]) 97 # function=lambda t: [amplitude*(1.-2.*(pi*(1./720.)*(t-720.))**2)/exp((pi*(1./720.)*(t-720.))**2) , 0.0, 0.0]) 98 # function=lambda t: [(-8.0*sin((1./720.)*t*2*pi))*((t<720.)-0.5*(t<360.)), 0.0, 0.0]) 99 99 def waveform(t): 100 100 return -amplitude*sin((1./wave_length)*t*2*pi)
Note: See TracChangeset
for help on using the changeset viewer.