Python tseries1の例

コード例 #1

ファイル: tseries_sea_ice_fw_flux.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

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)

コード例 #2

ファイル: tseries1_batch_gyre.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

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

コード例 #3

#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',

コード例 #4

ファイル: tseries_westshelf.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

#             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,

コード例 #5

ファイル: B-ISMF_figs.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

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)

コード例 #6

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,

コード例 #7

ファイル: tseries_batch_shelf.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

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)

コード例 #8

#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,

コード例 #9

ファイル: tseries_batch_cross_trough.py プロジェクト: cbegeman/e3sm-cryo-analysis-scripts

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)