rjleveque on dev
more changes_to_master for 5.8.0 update setrun descriptions fo a… (compare)
batchmodule is an add-on you can get from my github account. I use it to run parameter sweeps or testing easily with Clawpack setups.
And on that note, I am having trouble using a wall in the simple advection example. When I simply change the bc to
extrap, it works fine, but when I change it to
wallI get the following traceback:
File "/home/tmabbot/anaconda/lib/python2.7/site-packages/clawpack/pyclaw/solver.py", line 474, in _bc_lower array[self.reflect_index[idim],i,...] = -array[self.reflect_index[idim],2*self.num_ghost-1-i,...] # Negate normal velocity IndexError: index 1 is out of bounds for axis 0 with size 1
Still trying to get the feel for everything Clawpack is doing, so any guidance would be appreciated.
So if I'm reading
solver.py correctly, the wall is a reflecting boundary, fixing the gradient to be zero along the normal vector to the boundary?
If so, my question then is why this function throws an error. I understand why mathematically, but am curious what the mechanics of the program that cause this error.
I've been trying to use the BoussClaw add-on model to GeoClaw (https://www.sciencedirect.com/science/article/pii/S0378383917300546#s0010). However, it doesn't seem to gel very well with the version of ClawPack that I have installed - it gives a lot of errors for undefined references to modules in AMR.
Could this be an issue with the version I'm using (5.4.1), where BoussClaw might be developed for an older version? Does anyone have experience with this?
The error messages I am getting are all to do with the amr module, where I get undefined reference errors for amr_module_MOD _listofgrids, _amr_module_MOD _liststart and amr_module_MOD _bndlist and a few others sprinkled throughout the code.
Did those functions get replaced? And if so, by what?