import arcpy
from arcpy import env
from arcpy.sa import *

#Parametros
estacion=arcpy.GetParameterAsText(0)
ID_interpol1=arcpy.GetParameterAsText(1)
mdt=Ra(arcpy.GetParameterAsText(2))
salida=arcpy.GetParameterAsText(3)


#Comienza a calcular presipitacion
p=Idw(estacion,ID_interpol1,30,2,12)

#Comienza caculo de evaporacion
A=-1.9595
B=3011.5

out_ras=((("mdt")*A)+B)



#Calculo de RESERVA DE AGUA

RS=p/E

#Comienza la reclasificacioN

Reclasificacion=Reclassify(RS,"Value",RemapValue.sa([0.81,0.80,1],[0.75,0.65,2]
                                                    [0.65,0.50,3],[0.50,0.45,4],
                                                    [0.45,0.25,5],[0.25,0.15,6],
                                                    [0.15,0.0,7]))

Reclasificacion.save("salida")

arcpy.AddMessage ("El proceso a finalizado")