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()