all sorted! it seems to have fixed onto Mt Hood when I used Ailsa’s trick with the name selection
Thank you!
1 Like
Hi! I’ve been trying to run TITAN2D with my own volcano DEM, and am hitting a wall. I’m using a GeoTIFF file which loads and plots correctly, and am trying to model a pile collapse. I think I’ve specified all units correctly, but at ‘Run Titan here’ the code line completes instantly - suggesting there’s no actual flow being modelled. Any insights appreciated!
Hi Ailsa,
Could you please upload a screenshot of your input.py file? That will help me trying to figure out what is wrong. Also, please use the notebook located in shared/Victor Courses/Spring 2025/Week 5/ instead of the one located in the Titan2d folder from ‘victor setup’. Sam was supposed to update it because there were some bugs but not sure he did…
Hi guys! I was checking the assignment for today on Gradescope. However, I can´t find any instructions, and there isn´t any folder for week 4 (lava flows). I ran the four models that were explained in class, but what are the next instructions? Thanks!
Hi Omar,
The pdf file with the instructions is located in Gradescope, named Week4_TakeHome_Lava2D. I attached it here too. All you have to do is upload the outputs of each model into Gradescope as a word or pdf file. Hope this helps.
2 Likes
Sure! Attached is a copy of my input file. I have copied the notebook located in shared/Victor Courses/Spring 2025/Week 5/ to my own folder /ailsa/Online Course 2025/Week 5/, and run this - but it produces the same error. I’ve also attached a screenshot of the cell with the problem. Under Run Titan here!, when I execute the line it immediately returns the message:
CompletedProcess(args='./titan -nt 4 input.py', returncode=126) .
Googling ‘return code 126’ suggests this may have to do with permissions? E.g., “Error code 126 seems to imply you do not have permission to execute the script or a command inside it.”
Output when running line after Run Titan here!:
Screenshot of input.py file:
Hi Ailsa,
Must be an issue with the DEM. Try first to put it in a path that does not have any space. Then remove the argument ‘stderr=…’ to be able to see what the Titan shell while running. It should show you what the error is.
Hmm, could be a DEM-specific issue, but I have run the code with the Merapi DEM in > shared > DEMs and have the same error. Have renamed folders so file path has no whitespaces, but has not resolved the issue. Removing the argument ‘stderr=…’ reveals it is a permission issue (screenshot attached).
Interestingly, when I use the notebook in the TITAN2D folder, I get a different error message saying TITAN2D not found … even though it is installed and up to date according to ‘victor setup’ and ‘update’ in the terminal?
Ok, please make sure your notebook is in the same folder as your input.py file and your titan executable (files created with you run ‘victor setup’ for Titan). If the issue still persists after that, let us know.
Had partial success! I had issues running the notebook in my folder, so deleted and reinstalled ‘TITAN2D’ using victor setup and ran the notebook in there. Also reprojected DEM to UTM coordinates (first DEM actually in WGS84). I modelled a single pile source with initial velocity and direction to simulate a blast. Here’s the output, which I think is valid (the pile is called “pile_colima.asc”, but I think that’s just a naming relic from a previous Colima run?). Chose a low bed friction angle, hence the large runout (I think)?
Error generated by earlier cell, Plot output here, where I should be able to plot the pile at a given timestep. This generates a very long error message, which concludes: ValueError: 'vertices' must be 2D with shape (N, 2), but your input has shape (0,). I wondered if it was to do with the line, ax = plt.axes(projection=ccrs.epsg(32628)) (because my volcano’s UTM zone is 50S or 32750), but changing that did not fix. Not sure of what vertices is - insight appreciated!
Great job! Footprint of the PDCs looks good. Regarding the issue with getting a snapshot at a particular timestep, I will turn to Sam for this because I thought we fixed it but there is still maybe something wrong with this cell. I see in your input file that you put both a diter and dtime. Please set one to ‘None’ (generally the diter because we prefer having outputs every number of seconds) and try again.
1 Like