clawpack/pyclaw — Issue #107: Stegaton Application Missing mthlim
Author: mandli
Created: 2011-12-15 22:56:42+00:00
State: open

This example is not working in $PYCLAW/apps/elasticity/1d/stegoton giving a storage complaint about not having the correct number of arguments to the fortran module:

```
=============================================
apps/./elasticity/1d/stegoton
=============================================
Traceback (most recent call last):
File "./stegoton.py", line 185, in
output = run_app_from_main(stegoton)
File "/Users/mandli/src/clawpack_dev/pyclaw/src/pyclaw/util.py", line 51, in run_app_from_main
output=application(**app_kwargs)
File "./stegoton.py", line 177, in stegoton
status = claw.run()
File "/Users/mandli/src/clawpack_dev/pyclaw/src/pyclaw/controller.py", line 268, in run
status = self.solver.evolve_to_time(self.solution,t)
File "/Users/mandli/src/clawpack_dev/pyclaw/src/pyclaw/solver.py", line 671, in evolve_to_time
self.step(solution)
File "/Users/mandli/src/clawpack_dev/pyclaw/src/pyclaw/clawpack.py", line 148, in step
self.step_hyperbolic(solution)
File "/Users/mandli/src/clawpack_dev/pyclaw/src/pyclaw/clawpack.py", line 323, in step_hyperbolic
self.qbc,cfl = classic.step1(mbc,mx,self.qbc,self.auxbc,dx,dt,self.method,self.mthlim)
TypeError: Required argument 'mthlim' (pos 9) not found
```
It looks like mthlim is never set in the app script in $PYCLAW/apps/elasticity/1d/stegoton/stegaton.py so this may be an easy fix for someone who knows what the appropriate limiters are. If this is the problem we may want to add a check to make sure that mthlim is set in the solver `is_valid` function.