Python script:

import arcpy 
from arcpy.sa import * 
arcpy.env.workspace = "C:/Users/ThinkPad/PycharmProjects/Advanced_Gis_project/ski_slope_python.gdb" 
arcpy.env.overwriteOutput = True 

  
# Defining input variables 

  
Localitati = "Localitati_CJ.shp" 
Drumuri = "drumuri_CJ.shp" 
Ape = "Ape_CJ.shp" 
Partii = "Ski_slopes_CJ.shp" 
Zapada = Raster("snow_kriging") 
elev = Raster("dem_cj_stereo") 
dem = Raster("dem_cj_stereo") 
Output_cell_size = "25" 
Parallel_Processing_Factor = "0" 

  
# Extent = dem.extent 

  
Output_direction_raster = "" 
Output_back_direction_raster = "" 

  
# Process: Slope in percent from DEM 
Panta = Slope(dem, "PERCENT_RISE") 
Panta.save("Panta") 

  
# Process: Aspect from DEM 
Orientare = Aspect(dem) 
Orientare.save("Orientare") 

  
# Process: Euclidean distance Localitati, Drumuri, Ape, Partii 
Snap = arcpy.env.snapRaster 
arcpy.env.snapRaster = dem 
Parallel = arcpy.env.parallelProcessingFactor 
arcpy.env.parallelProcessingFactor = Parallel_Processing_Factor 
Extent = arcpy.env.extent 
arcpy.env.extent = dem 
Mask = arcpy.env.mask 
arcpy.env.mask = dem 
dist_loc = EucDistance(Localitati, "", Output_cell_size, Output_direction_raster) 
dist_loc.save("dist_loc") 
dist_drum = EucDistance(Drumuri, "", Output_cell_size, Output_direction_raster) 
dist_drum.save("dist_drum") 
dist_ape = EucDistance(Ape, "", Output_cell_size, Output_direction_raster) 
dist_ape.save("dist_ape") 
dist_partii = EucDistance(Partii, "", Output_cell_size, Output_direction_raster) 
dist_partii.save("dist_partii") 
arcpy.env.snapRaster = Snap 
arcpy.env.parallelProcessingFactor = Parallel 
arcpy.env.extent = Extent 
arcpy.env.mask = Mask 

  
# Process: Reclassify Panta, Orientare, dem, dist_loc, dist_drum, dist_ape, dist_partii 
remap_panta = RemapRange([[0, 6, 1], [6, 15, 6], [15, 25, 7], [25, 40, 9], [40, 45, 8], 
 [45, 55, 5], [55, 65, 4], [65, 75, 3], [75, 314, 2]]) 
reclass_panta = Reclassify(Panta, "VALUE", remap_panta) 
reclass_panta.save("reclass_panta")