geoclaw icon indicating copy to clipboard operation
geoclaw copied to clipboard

Running a forward model - reading in bathymetry file failure

Open saeedmohanna opened this issue 1 year ago • 9 comments

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.

saeedmohanna avatar Jul 26 '23 19:07 saeedmohanna

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

mandli avatar Jul 26 '23 23:07 mandli

Sure, here's the rest:

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

saeedmohanna avatar Jul 26 '23 23:07 saeedmohanna

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?

mandli avatar Jul 27 '23 13:07 mandli

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

saeedmohanna avatar Jul 27 '23 17:07 saeedmohanna

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

mandli avatar Jul 27 '23 20:07 mandli

Here's my setrun.py

setrun.zip

saeedmohanna avatar Jul 27 '23 22:07 saeedmohanna

Can you also send the dtopo file?

mandli avatar Jul 28 '23 02:07 mandli

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

saeedmohanna avatar Jul 28 '23 19:07 saeedmohanna

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

mandli avatar Jul 28 '23 21:07 mandli