Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
saeedmohanna opened this issue Jul 26, 2023 · 9 comments
Open

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

saeedmohanna opened this issue Jul 26, 2023 · 9 comments

Comments

@saeedmohanna
Copy link

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
Copy link
Member

mandli commented Jul 26, 2023

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

@saeedmohanna
Copy link
Author

Sure, here's the rest:

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

@mandli
Copy link
Member

mandli commented Jul 27, 2023

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
Copy link
Author

saeedmohanna commented Jul 27, 2023

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
Copy link
Member

mandli commented Jul 27, 2023

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
Copy link
Author

Here's my setrun.py

setrun.zip

@mandli
Copy link
Member

mandli commented Jul 28, 2023

Can you also send the dtopo file?

@saeedmohanna
Copy link
Author

saeedmohanna commented Jul 28, 2023

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
Copy link
Member

mandli commented Jul 28, 2023

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants