In fact there are two possibilities to play a 1D simulations:

  • 1) you have a DTM and you want to run multiple 1D simulations in different points of the domain: in this case you NEED the DTM and in the 2nd block of the file "options_point.txt"

inside the folder /param you should set the coordinates of interest:

/** 2 block - COORDINATES of the points for which the simulations is run

if state_px_coord==1 coordinate are (E, N) --- if state_px_coord==0 coordinate are (row,col) (max 9999 points)

#1 E(row) #2 N(col) #3 Z[m] #4 landuse #5 soiltype #6 slopes[rad] #7 aspect[rad] #8 sky[-] #9 dz/dE(net) #10 dz/dS(net) write -99 if you want to read the values from the distributed files

if the files are not present, set E=1 and N=1 and set the other values to the desired ones (default values: Z=0m, landuse=1, soiltype=1, slopes=0, aspect=0, sky=1)

#1 #2 #3 #4 #5 #6 #7 #8 #9 #10 */

2: double matrix local_results {1,10}

558175 7626807 -99 -99 -99 -99 -99 -99 -99 -99

In this case, putting -99, you are "asking" GEOtop to go and automatically set the parameters according to the files (slope, DEM, aspect...) inside the folder /morpho.

  • 2) you don't have a DTM and you want to run the simulation in a "virtual" point with specified topographical characteristics. Then you DON'T NEED the DTM and you have to set the parameters differently as:

/** 2 block - COORDINATES of the points for which the simulations is run

if state_px_coord==1 coordinate are (E, N) --- if state_px_coord==0 coordinate are (row,col) (max 9999 points)

#1 E(row) #2 N(col) #3 Z[m] #4 landuse #5 soiltype #6 slopes[rad] #7 aspect[rad] #8 sky[-] #9 dz/dE(net) #10 dz/dS(net) write -99 if you want to read the values from the distributed files

if the files are not present, set E=1 and N=1 and set the other values to the desired ones (default values: Z=0m, landuse=1, soiltype=1, slopes=0, aspect=0, sky=1)

#1 #2 #3 #4 #5 #6 #7 #8 #9 #10 */

2: double matrix local_results {1,10}

1 1 85.917 1 1 2.319 8.576 0.971 0.0156 -0.0617