#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")
