#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)
예제 #2
0
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")