import grass.script as gs

print("--- STARTING GSOC PROTOTYPE ---")
# Set the area we are working on
gs.run_command("g.region", raster="elev_lid792_1m")

# Run the calculation through Python
gs.run_command("r.slope.aspect", 
               elevation="elev_lid792_1m", 
               slope="gsoc_slope_final", 
               overwrite=True)

print("--- SUCCESS: gsoc_slope_final created! ---")
