#Defining area of interest in order to save computational time
g.region -s n=4814000 s=4810000 w=769000 e=775000

#Infinite slope stability model, factor of safety, only one lithological class, constant soil depth
r.slope.stability -v prefix=slx_inf_fos cellsize=10 model=i elevation=dem soilclass=lithology_simple seepage=1 observed=landslides numlayers=2 depthvals=2,100 geotech=1,1,15000,2000,30,40,1,2,27000,100000,89,0 overlap=100

#Infinite slope stability model, factor of safety, five lithological classes, constant soil depth
r.slope.stability -v prefix=sly_inf_fos cellsize=10 model=i elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthvals=2,100,2,100,2,100,2,100,2,100 geotech=1,1,15000,4000,30,40,1,2,27000,100000,89,0,2,1,15000,500,15,35,2,2,27000,100000,89,0,3,1,15000,2500,40,40,3,2,27000,100000,89,0,4,1,15000,5000,30,40,4,2,27000,100000,89,0,5,1,15000,0,20,40,5,2,27000,100000,89,0

#Infinite slope stability model, factor of safety, five lithological classes, varying soil depth
r.slope.stability -v prefix=sl_inf_fos cellsize=10 model=i elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthmaps=depth1,depth2 geotech=1,1,15000,4000,30,40,1,2,27000,100000,89,0,2,1,15000,500,15,35,2,2,27000,100000,89,0,3,1,15000,2500,40,40,3,2,27000,100000,89,0,4,1,15000,5000,30,40,4,2,27000,100000,89,0,5,1,15000,0,20,40,5,2,27000,100000,89,0

#Sliding surface model, soil class mode, truncated ellipsoids, factor of safety, five lithological classes
r.slope.stability -t -v prefix=sl_ell_fos cellsize=10 model=c elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthmaps=depth1,depth2 geotech=1,1,15000,4000,30,40,1,2,27000,100000,89,0,2,1,15000,500,15,35,2,2,27000,100000,89,0,3,1,15000,2500,40,40,3,2,27000,100000,89,0,4,1,15000,5000,30,40,4,2,27000,100000,89,0,5,1,15000,0,20,40,5,2,27000,100000,89,0 elldens=100 ellips=100,20,50,15,4,2,0.5,0.0,3,1,-9999,-9999,-9999,-9999 overlap=100

#Sliding surface model, soil class mode, truncated ellipsoids, Newmark seismic slope stability analysis, five lithological classes
r.slope.stability -t -v prefix=sl_ell_newmark cellsize=10 model=c elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthmaps=depth1,depth2 geotech=1,1,15000,4000,30,40,1,2,27000,100000,89,0,2,1,15000,500,15,35,2,2,27000,100000,89,0,3,1,15000,2500,40,40,3,2,27000,100000,89,0,4,1,15000,5000,30,40,4,2,27000,100000,89,0,5,1,15000,0,20,40,5,2,27000,100000,89,0 elldens=100 ellips=100,20,50,15,4,2,0.5,0.0,3,1,-9999,-9999,-9999,-9999 overlap=100 newmark=distance,10,7.0,0,1 newmarkref=newmarkref

#Sliding surface model, soil class mode, truncated ellipsoids, slope failure probability, five lithological classes, additional output
r.slope.stability -a -p -t -v prefix=sl_ell_pf cellsize=10 model=c elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthmaps=depth1,depth2 geotech=1,1,15000,4000,30,40,1000,1000,7000,5,20,40,1,2,27000,97500,88,0,1,95000,100000,1,87,89,2,1,15000,500,15,35,500,0,1000,5,10,20,2,2,27000,97500,88,0,1,95000,100000,1,87,89,3,1,15000,2500,40,40,750,1000,4000,5,35,45,3,2,27000,97500,88,0,1,95000,100000,1,87,89,4,1,15000,5000,30,40,1000,2000,8000,5,25,35,4,2,27000,97500,88,0,1,95000,100000,1,87,89,5,1,15000,50,20,40,25,0,100,5,15,25,5,2,27000,97500,88,0,1,95000,100000,1,87,89 sampling=2,3,3,0,1,1,1 elldens=100 ellips=100,20,50,15,4,2,0.5,0.0,3,1,-9999,-9999,-9999,-9999 overlap=100

#Sliding surface model, soil class mode, truncated ellipsoids, slope failure probability, division into tiles, five lithological classes
r.slope.stability -m -p -t -v prefix=sl_ell_pf_tiles cellsize=10 model=c elevation=dem soilclass=lithology seepage=1 observed=landslides numlayers=2,2,2,2,2 depthmaps=depth1,depth2 geotech=1,1,15000,4000,30,40,1000,1000,7000,5,20,40,1,2,27000,97500,88,0,1,95000,100000,1,87,89,2,1,15000,500,15,35,500,0,1000,5,10,20,2,2,27000,97500,88,0,1,95000,100000,1,87,89,3,1,15000,2500,40,40,750,1000,4000,5,35,45,3,2,27000,97500,88,0,1,95000,100000,1,87,89,4,1,15000,5000,30,40,1000,2000,8000,5,25,35,4,2,27000,97500,88,0,1,95000,100000,1,87,89,5,1,15000,50,20,40,25,0,100,5,15,25,5,2,27000,97500,88,0,1,95000,100000,1,87,89 sampling=2,3,3,0,1,1,1 elldens=100 ellips=100,20,50,15,4,2,0.5,0.0,3,1,-9999,-9999,-9999,-9999 ncores=16 tilesx=5 tilesy=5 overlap=100

#Restoring original region
g.region -s rast=dem