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':
# 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')
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...'
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...'
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)