#True for detrending data, False for raw data detr = True corr = False lterm = True rENSO = False drawbox = True ##360 is the central meridian (greenwich) #EDIT THIS FOR BOUNDS lonbounds = [0., 360.] latbounds = [-65, 65.] #detrend annual fields instead of monthly? prevents thf from blowing up for some reason... if detr: sst = detrend(sst) # ps = detrend(ps) # u = detrend(u) # v = detrend(v) field = detrend(field) lats = field.getLatitude()[:] lons = field.getLongitude()[:] #lats = sst.getLatitude()[:] #lons = sst.getLongitude()[:] #nt = sst.shape[0] #lons[0] = 0 #nlat = len(lats) #nlon = len(lons)
ps_mask = np.ma.mask_or(sst_mask[:ps.shape[0], :, :], ps_mask) ps = np.ma.array(ps, mask=ps_mask) #True for detrending data, False for raw data detr = True corr = False lterm = False #EDIT THIS FOR BOUNDS lonbounds = [290., 360.] latbounds = [-30, 70.] #detrend annual fields instead of monthly? prevents thf from blowing up for some reason... if detr: sst = detrend(sst) u = detrend(u) v = detrend(v) field = detrend(field) lats = field.getLatitude()[:] lons = field.getLongitude()[:] #lats = sst.getLatitude()[:] #lons = sst.getLongitude()[:] #nt = sst.shape[0] #lons[0] = 0 #nlat = len(lats) #nlon = len(lons) t = sst.getTime().asRelativeTime("months since 1980")