ind = 1
# Note 21499 cdd Na
for x in range(0, len(lats)): #len(lats[0:253])):
for y in range(0, len(lons)): #len(lons[0:191])):
o = feature.objects.get(pk=ind)
clim = climatecal(year[0], year[1], lats[x], lons[y], str_date)
pr = clim.cal_val('pr')
tmax = clim.cal_val('tasmax')
tmin = clim.cal_val('tasmin')
value = {}
# res_su
value[eindex[0]] = powerpack.su(year[0], year[1], len(pr), str_date,
tmax, tmin, pr)
# res_fd
value[eindex[1]] = powerpack.fd(year[0], year[1], len(pr), str_date,
tmax, tmin, pr)
# res_id
value[eindex[2]] = powerpack.id(year[0], year[1], len(pr), str_date,
tmax, tmin, pr)
# res_tr
value[eindex[3]] = powerpack.tr(year[0], year[1], len(pr), str_date,
tmax, tmin, pr)
# res_gsl
value[eindex[4]] = powerpack.gsl(year[0], year[1], len(pr), str_date,
tmax, tmin, pr)
# res_wsdi
# value[eindex[5]] = powerpack.wsdi(year[0],year[1],len(pr),str_date,tmax,tmin,pr)
# res_csdi
# value[eindex[6]] = powerpack.csdi(year[0],year[1],len(pr),str_date,tmax,tmin,pr)
# res_dtr
value[eindex[5]] = powerpack.dtr(year[0], year[1], len(pr), str_date,
