Python get_sections_names_from_file Examples

Example #1

    idfig = 'amoc'
    cyu = r'Sv'
    ym = 3.5
    yp = 24.5

elif cdiag == 'mean_mldr10_1':
    cvar = vdic['NN_MLD']
    idfig = 'mld'
    clnm = 'Mean mixed-layer depth, '
    cyu = r'm'
    ym = yp = 0.

elif cdiag == 'transport_sections':
    idfig = 'transport'
    print '  Using TRANSPORT_SECTION_FILE = ' + vdic['TRANSPORT_SECTION_FILE']
    list_sections = bo.get_sections_names_from_file(
        vdic['TRANSPORT_SECTION_FILE'])
    print 'List of sections to treat: ', list_sections

elif cdiag == 'seaice':
    idfig = 'ice'
    cyu = r'10$^6$km$^2$'

else:
    print 'ERROR: plot_time_series.py => diagnostic ' + cdiag + ' unknown!'
    sys.exit(0)

##########################################
# Basic temp., sali. and SSH time series #
##########################################

if idfig == 'simple':

Example #2

File: compare_time_series.py Project: brodeau/barakuda

# Transport through sections
############################

if itrsp == 1:

    print '\nUsing TRANSPORT_SECTION_FILE = '+vdic['TRANSPORT_SECTION_FILE']
    list_sections = bo.get_sections_names_from_file(vdic['TRANSPORT_SECTION_FILE'])
    print 'List of sections to treat: ', list_sections
    nbsect = len(list_sections)

    vstuff = [ 'volume', 'heat' , 'salt' ]
    vunit  = [ 'Sv'    , 'PW'   , 'kt/s' ]


    jrun = 0
    for confrun in clist_confruns:

        jsect=0
        for csect in list_sections:
            print '\n Treating transports through '+csect

            cf_in = cd_diag+'/'+confrun+'/transport_sect_'+csect+'.nc' ; bt.chck4f(cf_in, script_name='compare_time_series.py')

Example #3

File: plot_trsp_sigma.py Project: OscarEReyes/barakuda

DENSITY_SECTION_FILE = os.getenv('DENSITY_SECTION_FILE')
if DENSITY_SECTION_FILE == None: print 'The DENSITY_SECTION_FILE environement variable is no set'; sys.exit(0)


rsigdense0 = bphy.rsigma_dense

l_force_lim = True

cfig_type = 'png'

path_fig = DIAG_D+'/'

CONFRUN = CONF+'-'+RUN

print '  Using DENSITY_SECTION_FILE = '+DENSITY_SECTION_FILE
list_sections = bo.get_sections_names_from_file(DENSITY_SECTION_FILE)
print 'List of sections to treat: ', list_sections
nbsec = len(list_sections)



cf_in = DIAG_D+'/transport_by_sigma_class.nc' ; bt.chck4f(cf_in)



id_in = Dataset(cf_in)

vtime  = id_in.variables['time'][:]       ; nbm  = len(vtime)
vsigma = id_in.variables['sigma_bins'][:] ; nbin = len(vsigma)

print '      => '+str(nbin)+' sigma-density bins and '+str(nbm)+' time snapshots...'

Example #4

File: plot_trsp_sigma.py Project: jvegasbsc/barakuda

vdic = bt.check_env_var(sys.argv[0], venv_needed)

CONFRUN = vdic['ORCA'] + '-' + vdic['RUN']

rsigdense0 = bphy.rsigma_dense

l_force_lim = True

cfig_type = 'png'

path_fig = vdic['DIAG_D'] + '/'

cf_dens_sect = vdic['DENSITY_SECTION_FILE']
print '  Using cf_dens_sect = ' + cf_dens_sect
list_sections = bo.get_sections_names_from_file(cf_dens_sect)
print 'List of sections to treat: ', list_sections
nbsec = len(list_sections)

cf_in = vdic['DIAG_D'] + '/transport_by_sigma_class.nc'
bt.chck4f(cf_in)

id_in = Dataset(cf_in)

vtime = id_in.variables['time'][:]
nbm = len(vtime)
vsigma = id_in.variables['sigma_bins'][:]
nbin = len(vsigma)

print '      => ' + str(nbin) + ' sigma-density bins and ' + str(
    nbm) + ' time snapshots...'

Example #5

File: plot_trsp_sigma.py Project: brodeau/barakuda

vdic = bt.check_env_var(sys.argv[0], venv_needed)

CONFRUN = vdic["ORCA"] + "-" + vdic["RUN"]

rsigdense0 = bphy.rsigma_dense

l_force_lim = True

cfig_type = "png"

path_fig = vdic["DIAG_D"] + "/"

cf_dens_sect = vdic["DENSITY_SECTION_FILE"]
print "  Using cf_dens_sect = " + cf_dens_sect
list_sections = bo.get_sections_names_from_file(cf_dens_sect)
print "List of sections to treat: ", list_sections
nbsec = len(list_sections)


cf_in = vdic["DIAG_D"] + "/transport_by_sigma_class.nc"
bt.chck4f(cf_in)


id_in = Dataset(cf_in)

vtime = id_in.variables["time"][:]
nbm = len(vtime)
vsigma = id_in.variables["sigma_bins"][:]
nbin = len(vsigma)