# -*- coding: utf-8 -*-
"""
Created on Sat Apr 25 14:36:29 2020

@author: eirik
"""
import numpy as np
import matplotlib.pyplot as plt

#Introduksjonstekst til den som bruker programmet
print("Prognosene viser at temperaturen i verden vil stige med rundt 2 grader frem mot 2050. Dette skyldes i stor grad utslipp av CO2, \
 og en av konsekvensene av økningen vil være at havnivået stiger og tar over hjemmene til opp mot 110 millioner mennesker. \
 Men tiltak kan gjøres. I denne simuleringen er du verdensleder og skal ta et valg hvert femte år for hvilket tiltak du vil gjennomføre. \
 Noen tiltak har større innvirkning enn andre, så valgene dine har stor konsekvens for hvordan verden ser ut i fremtiden. Lykke til! \n")

#Legger inn årlige besparelser for ulike tiltak
Tiltaksliste = ["Helse og utdannelse i fattige land", 85.42/30, "Plantebasert diett", 65.01/30, "Redusere matsvinn", 87.45/30, \
                "Forbedre aircondition", 57.75/30, "Forbedre kjøleskap", 43.53/30, "Regnskogsplanting", 54.45/30, "Vindkraft på land", 47.21/30]
Spartliste = [0]

#Prognose for maksveverdier 2050
TempSluttAr = 2
Flyktninger = 110
CO2prAr = 50


#La bruker velge et tiltak hvert femte år. Legg CO2-besparelsen for tiltaket til Spartliste
#og slett tiltaket fra Tiltakslite.
Spart = 0
for i in range(30, 0, -5):
    print("Mulige tiltak år", 2050-i)
    for n in range(0, len(Tiltaksliste),2):
        print(int((n+2)/2), Tiltaksliste[n])
    Tiltak = int(input("Skriv inn tallet på tiltaket du vil gjøre: "))
    Spart = Spart + (Tiltaksliste[Tiltak*2-1]*i)
    Spartliste.append(Tiltaksliste[Tiltak*2-1]*5)
    del Tiltaksliste[Tiltak*2-2]
    del Tiltaksliste[Tiltak*2-2]

#Regne ut kummulative tall til graf
Kum1 = np.cumsum(Spartliste)
Kum2 = np.cumsum(Kum1)
xliste = [2020, 2025, 2030, 2035, 2040, 2045, 2050]

#Graf over CO2-utslipp
xlist = xliste
MaxCO2list = [0, CO2prAr*5, CO2prAr*10, CO2prAr*15, CO2prAr*20, CO2prAr*25, CO2prAr*30]
fig = plt.figure()
ax1 = fig.add_subplot(111)
ax1.grid(True)
plt.title("CO2-utslipp")
plt.xlabel("Årstall")
plt.ylabel("Totalt utslipp siden 2020 i gigatonn")

ax1.plot(xlist,MaxCO2list, label="Forventet graf")
ax1.plot(xlist,MaxCO2list - Kum2, label="Graf etter tiltak")
plt.legend()
plt.show()

#Graf over temperaturstigning
xlist = xliste
MaxTmplist = [0, TempSluttAr/6, TempSluttAr*2/6, 3*TempSluttAr/6, 4*TempSluttAr/6, 5*TempSluttAr/6, TempSluttAr]
fig = plt.figure()
ax1 = fig.add_subplot(111)
ax1.grid(True)
plt.title("Økning i temperatur")
plt.xlabel("Årstall")
plt.ylabel("Grader Celsius")

ax1.plot(xlist,MaxTmplist, label="Forventet graf")
f=[0]
for i in range (1,len(MaxCO2list), 1):
    f.append((MaxCO2list[i] - Kum2[i])/MaxCO2list[i]*MaxTmplist[i])
ax1.plot(xlist, f, label="Graf etter tiltak")
plt.legend()
plt.show()

#Graf over økning i klimaflyktninger
xlist = xliste
MaxFlyktlist = [0, Flyktninger/6, 2*Flyktninger/6, 3*Flyktninger/6, 4*Flyktninger/6, 5*Flyktninger/6, Flyktninger]
fig = plt.figure()
ax1 = fig.add_subplot(111)
ax1.grid(True)
plt.title("Økning i klimaflyktninger")
plt.xlabel("Årstall")
plt.ylabel("Millioner mennesker")

ax1.plot(xlist, MaxFlyktlist, label="Forventet graf")
g=[0]
for i in range (1,len(MaxCO2list), 1):
    g.append((MaxCO2list[i] - Kum2[i])/MaxCO2list[i]*MaxFlyktlist[i])
ax1.plot(xlist, g, label="Graf etter tiltak")
plt.legend()
plt.show()

print("Gratulerer! Du klarte å redusere forventet utslipp med ", "%.2f" % Spart, "gigatonn CO2, temperaturen steg bare","%.2f" % f[6], \
      "grader og ", "%.2f" % (Flyktninger - g[6]), "millioner mennesker kan takke deg for at de ikke er blitt klimaflyktninger!")

if (Spart > 256.7):
    print("Veldig bra! Du valgte alltid de tiltakene som gav størst effekt!")
else:
    print("Godt forsøk! Du valgte godt, men det fantes mer effektive tiltak.")