Melt causing elevations to drop to zero

jtunnicl / ThawScape

0 stars 2 forks source link

Melt causing elevations to drop to zero #13

Open chrisdjscott opened 5 years ago

chrisdjscott commented 5 years ago

Change ordering of melt operation to be from low to high elevations
Compute melt potential / incoming raster first, then apply it during the avalanche routine
Make sure there are no pit/flats in DEM before computing the low to high ordering
Could try changing melt parameter

chrisdjscott commented 5 years ago

Running for 1 day with all components on ends with:

erosion

Turning off melt and running for 1 day gives:

erosion

So it seems something to do with melt is causing the elevations to drop to zero very quickly.

Changing the ordering to low to high looks worse:

erosion

chrisdjscott commented 5 years ago

I added an option debug_melt that outputs some rasters.

During the first step they look like:

debug_2010_78_13_shade_raster asc

debug_2010_78_13_I_P asc

debug_2010_78_13_incoming asc

chrisdjscott commented 5 years ago

If I set the melt parameter to 25000 (was 250) then it behaves better over the 1 day run.

erosion

This is for the version where we first compute incoming in a separate loop and then apply it ordered from low to high pixels. I will now try moving the application of the melt into the avalanche routine.