for i in range(1979,2019): filename = str(i) + '.nc' print(filename) nc = Dataset(filename) temp = nc.variables['t2m'][:] for m in range(12): print(m) ver = temp[m,282,1040] - 273.15 print(ver) f = open('veracruz.txt', 'a') f.write(str(i)+ ' ' + str(m)+ ' '+ str(ver) + '\n') clim = np.mean(t_ver[0:360]) for i in range(1979,2018): print(i) nc = Dataset(str(i)+'.nc') var = nc.variables['t2m'][:] tver = (np.mean(var[:,275,1093])) - 273.15 anom = tver - clim print(anom)