import numpy as N

ort='/storage/emulated/0/Download/'
#ort='./'
txt='tmax.txt'
htm='smartphone.html'

ja=N.genfromtxt(ort+txt,usecols=(0),dtype='i')
mo=N.genfromtxt(ort+txt,usecols=(1),dtype='i')
ta=N.genfromtxt(ort+txt,usecols=(2),dtype='i')
tg=N.genfromtxt(ort+txt,usecols=(3),dtype='f')

jo=N.arange(1893,2017,1);nj=len(jo)

to=N.zeros(nj,float)
tx=N.zeros(nj,float)
tn=N.zeros(nj,float)

for j in range(nj):

    print jo[j]

    w=N.ma.masked_where(ja!=jo[j],tg)
    w=w.compressed()
    
    to[j]=N.mean(w)
    
    tx[j]=len(w[w>30.])
    tn[j]=len(w[w<0.])
    
###########

xo=N.vstack([N.arange(nj),N.ones(nj)]).T

m,c=N.linalg.lstsq(xo,to)[0]

print m,c
    
#################

f=open(ort+htm,'w')

f.write('<!DOCTYPE html>\n')
f.write('<html>\n')
f.write('<head>\n')
f.write('<script type="text/javascript"\n')
f.write('src="'+ort+'dygraph-combined.js"></script>\n')
f.write('<link rel="stylesheet" src="'+ort+'dygraph.css"/>\n')
f.write('</head>\n')
f.write('<body>\n')

f.write('<div id="graphdiv2"\n')
f.write('style="width:500px;height:600px;background-color:lightgrey"></div>\n')
f.write('<script type="text/javascript">\n')

f.write('g2 = new Dygraph(document.getElementById("graphdiv2"),\n')
	
f.write('	"YEAR,HETA,EIST\\n"+\n')

for j in range(nj):

    if(j<nj-1):

       f.write('"%i,%5.1f,%5.2f\\n"+\n'%(jo[j],tx[j],tn[j]))#to[j],c+m*j))

    else:
    
       f.write('"%i,%5.1f,%5.2f\\n",\n'%(jo[j],tx[j],tn[j]))#to[j],c+m*j))

f.write('{\n')
f.write('title: "POTSDAM",\n')
f.write('ylabel: "Days",\n')
f.write('xlabel: "Years",\n')
f.write('rollPeriod: 7,\n')
f.write('showRoller: true,\n')    
f.write('costumBars: true,\n')
f.write('colors: ["red","blue"],\n')
f.write('legend: "always",\n')
 
f.write('}\n')
f.write(');\n')
f.write('</script>\n')
f.write('</body>\n')
f.write('</html>\n')
f.close()

print htm
print 'geschrieben in'
print ort
print 'ENDE'