import netCDF4 import datetime import numpy as np import matplotlib.pyplot as plt import scipy.signal from matplotlib import cm from math import pi import weddell_mod as wed #run_incr = ['ISMF','ISMF-noEAIS','ISMF-noDIB'] run_incr = ['ISMF', 'ISMF-noEAIS', 'ISMF-3dGM', 'ISMF-noDIB'] placename = 'frisEAcoast' var_incr = ['sea_ice_fw_flux'] year_range = [10, 180] wed.tseries1( run_incr, var_incr, year_range=year_range, placename=placename, input_filename= 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_sea_ice_fw_flux_frisEAcoast_t010-180', print_to_file=False, create_figure=True, operation='area_sum', #reference_run = 'ISMF', #tav = 12, #apply_filter = True, cutoff = 1/4, year_overlay=True, overwrite=False)
# placename = 'wed_pyc_Ryan', # apply_filter = True, cutoff = 1/4, # print_to_file=True,create_figure=True, # input_filename = 'ISMF_zpyc_wed_pyc_Ryan_70-101_zlim-4500--3500') #wed.tseries1(run_incr,['mean'],year_range=year_range, # placename = 'wed_pyc_Ryan', # apply_filter = True, cutoff = 1/4, # print_to_file=True,create_figure=True, # input_filename = 'ISMF_zpyc_wed_pyc_Ryan_70-101_zlim-2000--500') #wed.tseries1(run_incr,var_incr,year_range=year_range, # placename = 'gyre_interior', # print_to_file=True,create_figure=True, # apply_filter = True, cutoff = 1/4, # ztop_pyc = [True], zbottom_pyc = [False], # year_overlay=False,overwrite=True) #wed.tseries1(run_incr,var_incr,year_range=year_range, # placename = placename, # print_to_file=True,create_figure=True, # varlim = rholim, # #apply_filter = True, cutoff = 1/4, # ztop_pyc = [False], zbottom_pyc = [True], # year_overlay=False,overwrite=True) wed.tseries1(run_incr,var_incr,year_range=year_range, placename = placename, print_to_file=False,create_figure=True, varlim = rholim, input_filename = filename, ztop_pyc = [True,False], zbottom_pyc = [False,True], year_overlay=False,overwrite=True)
#wed.tseries1(run_incr,var_incr,year_range=year_range, # input_filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_M31W_z20.000000mab_t010-180', # placename = placename,zeval=20.,zab=True, # print_to_file=False,create_figure=True, # show_obs = True, obs = M31W_T_obs, # year_overlay=False,overwrite=False) # var_incr = ['S'] wed.tseries1( run_incr, var_incr, year_range=year_range, input_filename= 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_M31W_z20.000000mab_t010-180', placename=placename, zeval=20., zab=True, print_to_file=False, create_figure=True, show_obs=True, obs=M31W_S_obs, year_overlay=False, overwrite=False) var_incr = ['rho'] wed.tseries1( run_incr, var_incr, year_range=year_range, input_filename= 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_M31W_z20.000000mab_t010-180',
# ztop_pyc = [True], zbottom_pyc = [False], # year_overlay=False,overwrite=True) #wed.tseries1(run_incr,var_incr,year_range=year_range, # placename = placename, # print_to_file=True,create_figure=True, # varlim = rholim, # #apply_filter = True, cutoff = 1/4, # ztop_pyc = [False], zbottom_pyc = [True], # year_overlay=False,overwrite=True) wed.tseries1( run_incr, var_incr, year_range=year_range, placename=placename, #reference_run='ISMF', #print_to_file=True,create_figure=False, print_to_file=False, create_figure=True, input_filename=filename, #apply_filter = True, cutoff = 1/4, #ztop_pyc = [True], zbottom_pyc = [False], ztop_pyc=[False], zbottom_pyc=[True], year_overlay=False, overwrite=True) #zeval = [-100,-400], #zrange = [-100,-500], #lat=lat, lon=lon, #option = 'coord',placename = 'S4E',#'M31W', #zrange=[0,20],zab=True,
filename = '_transect_u_' + locname + '_10-100' #wed.fluxgate(locname,yrrange=year_range,morange = [1,12+1],run_incr=run_incr, # mode = 'barotropic-baroclinic', #mode = 'pos-neg', # plot_map = False, plot_transect = False, # overwrite=False,savepath = savepath) #wed.tseries1(run_incr, # ['F_barotropic','F_baroclinic_pos'], # year_range=year_range, # placename=locname, # input_filename = run_incr[0] + filename, # savepath = savepath) #placename = 'M31W' #wed.tseries1(run_incr,['T','S'],year_range=year_range, # placename = placename,zeval=20.,zab=True, # print_to_file=True,create_figure=True, # savepath=savepath,year_overlay=False,overwrite=False) placename = 'wed_pyc_Ryan_shelf' wed.tseries1( run_incr, ['rho', 'rho'], year_range=year_range, placename=placename, #input_filename = 'B-ISMF_rho_wed_pyc_Ryan_shelf_abovepyc_belowpyc_t010-011', print_to_file=True, create_figure=True, ztop_pyc=[True, False], zbottom_pyc=[False, True], savepath=savepath, year_overlay=False, overwrite=False)
import weddell_mod as wed from plot_config import wedwang_c_btr_flux_obs run_incr = ['ISMF', 'ISMF-noEAIS', 'ISMF-3dGM', 'ISMF-noDIB'] placename = 'wedwang_c' #var_incr = ['F_baroclinic_pos'] var_incr = ['F_barotropic'] savepath_flux = '/global/homes/c/cbegeman/weddell_output/fluxgate/' filename = 'ISMF_transect_u_wedwang_c_10-180' year_range = [10, 180] wed.tseries1(run_incr, var_incr, year_range=year_range, placename=placename, print_to_file=False, create_figure=True, input_filename=filename, savepath=savepath_flux, show_obs=True, obs=wedwang_c_btr_flux_obs, year_overlay=False, overwrite=True) #ratio_barotropic = [True], #apply_filter = True, cutoff = 1/4, #zeval = [-100,-400], #zrange = [-100,-500], #lat=lat, lon=lon, #option = 'coord',placename = 'S4E',#'M31W', #zrange=[0,20],zab=True,
import matplotlib.pyplot as plt import scipy.signal from matplotlib import cm from math import pi import weddell_mod as wed run_incr = ['ISMF', 'ISMF-noEAIS', 'ISMF-3dGM', 'ISMF-noDIB'] #placename = 'wed_pyc_Ryan_shallow' #placename = 'gyre_interior' placename = 'wed_pyc_Ryan_shelf' #var_incr = ['taux','tauy'] #var_incr = ['T','S','rho','u','v'] #var_incr = ['T','S'] #var_incr = ['unormal'] var_incr = ['rho'] #filename = 'ISMF_ISMF-noEAIS_rho_wed_pyc_Ryan_shallow_abovepyc_t070-101' year_range = [10, 180] wed.tseries1( run_incr, var_incr, year_range=year_range, placename=placename, print_to_file=False, create_figure=True, #input_filename = filename, #apply_filter = True, cutoff = 1/4, ztop_pyc=[True for i in run_incr], year_overlay=False, overwrite=False)
#var_incr = ['unormal'] var_incr = ['rho'] #var_incr = ['rho','rho'] rholim = [1027.2, 1027.9] lat = -70 lon = 340 #filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_gyre_interior_abovepyc_t010-180' filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_' + placename + '_abovepyc_t010-180' year_range = [10, 180] wed.tseries1(run_incr, var_incr, year_range=year_range, placename=placename, print_to_file=False, create_figure=True, input_filename=filename, ztop_pyc=[True], diff_pyc=[True], show_obs=True, obs=gyre_interior_drho_obs, overwrite=False) #wed.tseries1(run_incr,['mean'],year_range=year_range, # placename = 'wed_pyc_Ryan', # apply_filter = True, cutoff = 1/4, # print_to_file=True,create_figure=True, # input_filename = 'ISMF_zpyc_wed_pyc_Ryan_70-101_zlim-4500--3500') #wed.tseries1(run_incr,['mean'],year_range=year_range, # placename = 'wed_pyc_Ryan', # apply_filter = True, cutoff = 1/4, # print_to_file=True,create_figure=True,
import numpy as np import matplotlib.pyplot as plt import scipy.signal from matplotlib import cm from math import pi import weddell_mod as wed run_incr = ['ISMF','ISMF-noEAIS','ISMF-3dGM','ISMF-noDIB'] #placename = 'S4E' #placename = 'M31W' placename = 'N31W' #var_incr = ['T','S','rho']#,'u','v'] #var_incr = ['T']#,'u','v'] var_incr = ['rho']#,'u','v'] #filename = 'ISMF_ISMF-noEAIS_rho_wed_pyc_Ryan_shallow_abovepyc_t070-101' year_range = [10,180] #year_range = [10,11] wed.tseries1(run_incr,var_incr,year_range=year_range, #input_filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_S4E_z20.000000mab_t010-180', input_filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_M31W_z20.000000mab_t010-180', input_filename2 = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_N31W_z500.000000m_t010-180', #input_filename = 'ISMF_ISMF-noEAIS_ISMF-3dGM_ISMF-noDIB_rho_M31W_W_z20.000000mab_t010-180', var2 = 'rho_z500', placename = placename,zeval=20.,zab=True, #placename = placename,zeval=500.,zab=False, print_to_file=False,create_figure=True, #print_to_file=True,create_figure=False, year_overlay=False,overwrite=False)