def AgricultureWater(time):
pop = (Expo(Fit.PopulationFit(time),pPara)/Expo(Fit.PopulationFit(years[-1]),pPara))
pcg = (Expo(Fit.PCGDPFit(time),gPara)/Expo(Fit.PCGDPFit(years[-1]),gPara))
ira = (Expo(Fit.IrrigationAreaFit(time),iPara)/Expo(Fit.IrrigationAreaFit(years[-1]),iPara))
return Data.WaterUseAgriculture[years[-1]]*(pop/3+pcg/3+ira/3)
x = np.linspace(years[0],years[-1],100)
plt.figure()
plt.xlabel("Population/ 10k people")
plt.ylabel("Agreculture Water Usage/ 100m m^3")
plt.scatter([Data.Population[t] for t in years],[Data.WaterUseAgriculture[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.PopulationFit(x),Expo(Fit.PopulationFit(x),pPara),'g-',label="Fit")
plt.legend()
plt.show()
plt.close()
plt.figure()
plt.xlabel("PCGDP/ CNY")
plt.ylabel("Agreculture Water Usage/ 100m m^3")
plt.scatter([Data.PCGDP[t] for t in years],[Data.WaterUseAgriculture[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.PCGDPFit(x),Expo(Fit.PCGDPFit(x),gPara),'g-',label="Fit")
plt.legend(loc=0)
plt.show()
plt.close()
plt.figure()
plt.xlabel("Irrigation Area/ kha" )
plt.ylabel("Agreculture Water Usage/ 100m m^3")
plt.scatter([Data.IrrigationArea[t] for t in years],[Data.WaterUseAgriculture[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.IrrigationAreaFit(x),Expo(Fit.IrrigationAreaFit(x),iPara),'g-',label="Fit")
plt.legend(loc=0)
plt.show()
plt.close()
