Randall J. LeVeque
Randall J. LeVeque
The [tutorials](https://github.com/geohackweek/tutorial_contents) from a recent [Geohack week](https://geohackweek.github.io/ghw2017/) contain lots of some hints on reading in raster files via `rasterio`, applying map projections, etc. For those using conda environments, this [environment.yml](https://github.com/geohackweek/vizhack/blob/master/environment.yml)...
This is not documented except in http://www.clawpack.org/release_5_3_1.html where it says: Improved the checkpoint and restart capabilities, adding the option to set checkpt_style < 0 in setrun.py if you want the...
Often a script `make_fgmax.py` is used to generate the text file that specifies an fgmax grid. This can be executed automatically from setrun.py (add to the main program at the...
When creating topo `B = aux(1,:,;)` in `setaux.f90` on a fine grid along the coast, there may be some cells where the average topo is just slightly negative, `-dry_tolerance <...
Arrival times are computed `src/2d/shallow/fgmax_frompatch.f90` in a way that assumes the sea level and initial shoreline did not change due to subsidence or uplift during an earthquake, which is often...
Currently geoclaw provides output like what's below, which is confusing. This is also different from what is provided by amrclaw. Clean this up since most users don't care about `Bufnst`...
Should have example testing bc2amr in this case, as mentioned in #122.
See discussion in #134 and #137
Both compute the great circle distance, but the old routine took `x1,y1,x1,y2` and computed distance from `(x1,y1)` to `(x2,y2)` while the new one takes `x,y` and computes distance from `(x[0],x[1])`...
The read below gives `Z` with the wrong shape because the file is really topo_type 2. We should check that shapes end up being consistent... ``` from clawpack.geoclaw import topotools...