Geoclaw 1D / AMR

[paulinchin]

I'd like to run GEOCLAW (CLAWPACK 5.4.1) in 1D with Cartesian grid. As ndim=1 doesn't work, I set it as suggested here: https://www.clawpack.org/sphinx-versioning/amrclaw1d.html. However in this case it doesn't see the topography.

My setrun.py:

    clawdata.lower[0] = 0.
    clawdata.upper[0] = 5000000.

    clawdata.lower[1] = 0.
    clawdata.upper[1] = 0.

    # Number of grid cells: Coarsest grid
    clawdata.num_cells[0] = 5000
    clawdata.num_cells[1] = 1

    amrdata.refinement_ratios_x = [1]
    amrdata.refinement_ratios_y = [1]
    amrdata.refinement_ratios_t = [1]

My topography file:

5200  mx
1  my
0  xlower
0  ylower
1000  cellsize
-9999  nodataval
-4000
-4000
-4000
-4000
-4000
... 
5200 values

I get an error:
*** Error, grid cell does not overlap any topo grid xim = 0.00000000000000E+00 xip = 0.10000000000000E+04 yjm = 0.00000000000000E+00 yjp = 0.00000000000000E+00

But it doesn't give error about not overlapping area. I also tried to fully comment second dimension, but in that case it gives error about domain overlapping.

I found 1D version here: https://github.com/clawpack/geoclaw_1d, but not sure if it is ready to use.

Also, I found a problem setting any other than topotype=2 topography and qinit dtopo (as in example for "bowl-radial" simulation) in Cartesian grid simulation. GEOCLAW doesn't provide any error, but doesn't see dtopo in other formats.