-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtemis-uvi.py
55 lines (37 loc) · 1.11 KB
/
temis-uvi.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
# Analysis of interannual variability of incident solar UV index radiation from TEMIS data.
# Original data source and documentation:
# https://www.temis.nl/uvradiation/product/uvi-uvd.html
# Instructions from:
# https://realpython.com/python-csv/
#import csv
print('começa aqui...')
import numpy as np
uvi_path = '/home/santiago/Projetos/temis-uvi/TEMIS_UVI_Acarau.csv'
uvi_data = np.genfromtxt(uvi_path, delimiter=';', names=True, dtype=None, encoding=None)
print(uvi_data.dtype.names)
yy = uvi_data['yyyy']
mm = uvi_data['mm']
dd = uvi_data['dd']
dt = uvi_data['date']
ief = uvi_data['uvief']
uef = uvi_data['uvdef']
dvf = uvi_data['uvdvf']
ddf = uvi_data['uvddf']
ozo = uvi_data['ozone']
clim_ozo = []
for m in range(1,13):
print(m)
print(np.average(ozo[np.where(mm == m)]))
np.append(clim_ozo, np.average(ozo[np.where(mm == m)]))
# clim_ozo = np.average(ozo[np.where(mm == m)])
#print(clim_ozo)
#ijan = np.where(mm == 1)
#print(ijan)
#print(ozo[ijan])
#print(np.average(ozo[ijan]))
#print(range(1,12))
#range(12)
#for i in range(1, 12):
# print(i)
print(clim_ozo)
print('chegamos ao fim!')