Pellegrino

import math
import random

def calcv(latitud,t):
      return (23.5*math.sin(math.pi*(t-80)/180)+90-latitud)/90

def decide(v):
    if 0<v<1:
        return v**2
    if v>=1:
        return 1
    if v<=0:
        return 0

randomlist=[]
def calcenergy(f):
    t=random.random()
    randomlist.append(t)
    return area*f*sundigit*t
print randomlist
area=500
sundigit=10
t=int(input("How many latitudes do you want to calculate?"))

x=0
listalatitud=[]
while x<t:
      latitud=int(input("latitud"))
      if latitud>0 and latitud<90:
          listalatitud.append(latitud)
          x=x+1
      else:
          print("Next time say a latitude between 0 and 90")

lista=[]
total=0
for latitud in listalatitud:
    tmp = []
    for t in xrange(365):
        v=calcv(latitud,t)
        f=decide(v)
        l=calcenergy(f)
        tmp.append(l)
    lista.append(tmp)

latitudlista=["45","25"]

print "Area", "soltal", "latitud", "dag", "solighetsfaktor","f(t,latitud)","W(t)"

Month=["Januari","Februari"]
for x in listalatitud:
    for t in range(0,2):
        print Month[t]
        for i in range (1,31):
            print 500,10,x,i,randomlist[i],decide(calcv(x,i)),calcenergy(decide(calcv(x,i)))