sys.path.append("/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/00_Commons/03_Scripts/")
sys.path.append("/Users/gliu/Downloads/02_Research/01_Projects/03_SeaLevel/03_Scripts/cluster_ssh/")
from amv import proc,viz
import slutil
import yo_box as ybx
import tbx
import itertools
import matplotlib as mpl
#%%
# Set Paths
datpath = "/Users/gliu/Downloads/02_Research/01_Projects/03_SeaLevel/01_Data/01_Proc/"
outfigpath = "/Users/gliu/Downloads/02_Research/01_Projects/03_SeaLevel/02_Figures/20210610/"
proc.makedir(outfigpath)
# Experiment Names
#start = '1993-01'
#end = '2013-01'
#start = '1850-01'
#end = '2100-12'
nclusters = 6
rem_gmsl = False
maxiter = 5 # Number of iterations for elimiting points
minpts = 30 # Minimum points per cluster
# Other Toggles
debug = True
savesteps = True # Save Intermediate Variables
filteragain = False # Smooth variable again after coarsening
import matplotlib.pyplot as plt
import xarray as xr
import calendar as cal
import scm
import time
import cmocean
import copy
#%% Data Paths
datpath = "/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/01_hfdamping/01_Data/"
lipath = "/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/02_stochmod/01_Data/model_input/"
llpath = "/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/02_stochmod/01_Data/model_input/"
figpath = '/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/02_stochmod/02_Figures/20220415/'
proc.makedir(figpath)
mconfigs = ["PIC-SLAB", "PIC-FULL"]
dofs = [898 - 1 - 2 - 2, 1898 - 1 - 2 - 2]
lags = [1, 2, 3]
bboxplot = [-80, 0, 5, 60]
mons3 = [viz.return_mon_label(m, nletters=3) for m in np.arange(1, 13)]
#%% Functions
def load_dampraw(mconfig, datpath):
inpaths = datpath + "CESM-" + mconfig + "-Damping/"
damping = np.load(
inpaths +
import sys
projpath = "/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/02_stochmod/"
datpath = projpath + '01_Data/model_output/'
rawpath = projpath + '01_Data/model_input/'
outpathdat = datpath + '/proc/'
sys.path.append(
"/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/02_stochmod/03_Scripts/stochmod/model/"
)
sys.path.append(
"/Users/gliu/Downloads/02_Research/01_Projects/01_AMV/00_Commons/03_Scripts/"
)
outpath = projpath + '02_Figures/20210527/'
proc.makedir(outpath)
#%% Experiment parameters
# Options to determine the experiment ID
fscale = 1 # Number to scale NAO and other forcings by
nyrs = 1000 # Number of years to integrate over
applyfac = 2
mconfig = "SLAB_PIC"
# Do a stormtrackloop
runid = "303"
funiform = 1.5
# Experiment ID
expid = "%s_%iyr_funiform%i_run%s_fscale%03d_applyfac%i" % (
