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
