clawpack / geoclaw

Version of Clawpack for geophysical waves and flows
http://www.clawpack.org/geoclaw
BSD 3-Clause "New" or "Revised" License
76 stars 87 forks source link

Running a forward model - reading in bathymetry file failure #552

Open saeedmohanna opened 1 year ago

saeedmohanna commented 1 year ago

Hello,

I am trying to run a forward model for a tsunami and I get the following error when running make .plots:

Screen Shot 2023-07-26 at 12 35 57 PM

I'm not sure what could be causing this as the file is not that large (44.9 MB)

Please advise.

Thanks.

Best, Saeed.

mandli commented 1 year ago

Can you provide more of the output? There's nothing that should cause that to occur off-hand.

saeedmohanna commented 1 year ago

Sure, here's the rest:

Screenshot from 2023-07-26 16-29-07

mandli commented 1 year ago

My best guess is that there's something wrong with your topography file. I do not suppose you could post it somewhere so one of us could download the file and see if we can figure out what is wrong?

saeedmohanna commented 1 year ago

Thanks for the prompt response. Attached is the topo file (type 3). The original version was in type 1 (too large to be allowed to be uploaded by github) - I read it in and converted it to type 3. type3topo.zip

mandli commented 1 year ago

I can read the topography file fine so that is not it. Attached is a zip file with a working Makefileand setrun.py that runs to completion. Can you send your setrun.py maybe? scripts.zip

saeedmohanna commented 1 year ago

Here's my setrun.py

setrun.zip

mandli commented 1 year ago

Can you also send the dtopo file?

saeedmohanna commented 1 year ago

I was able to fix it! The dtopo file type I specified was type 1 which I thought was correct as my dtopo file had 3 columns (x,y,z). However, I forgot that geoclaw's type 1 dtopo file actually had 4 columns which I think are level, x, y , and z. When I used the correct format of dtopo_type=1, the code ran successfully! Thank you for your help

mandli commented 1 year ago

The dtopo type 1 is t, x, y, dz I believe, level is not part of it.