# import sys import os import numpy as nmp from netCDF4 import Dataset import barakuda_tool as bt import barakuda_orca as bo import barakuda_plot as bp venv_needed = {'ORCA','RUN','DIAG_D','NN_T','NN_S'} vdic = bt.check_env_var(sys.argv[0], venv_needed) CONFRUN = vdic['ORCA']+'-'+vdic['RUN'] cname_temp = vdic['NN_T'] cname_sali = vdic['NN_S'] #if len(sys.argv) != 4 and len(sys.argv) != 6 : # print 'Usage: '+sys.argv[0]+' <YYYY1> <YYYY2> <Nb. Levels> (<name temp.> <name salin.>)' # sys.exit(0) #cy1 = sys.argv[1] ; cy2 = sys.argv[2] ; jy1=int(cy1); jy2=int(cy2)
from netCDF4 import Dataset import barakuda_tool as bt import barakuda_orca as bo import barakuda_plot as bp lfig0 = True lfig1 = True lfig2 = True venv_needed = {'ORCA','EXP','DIAG_D','COMP2D','i_do_sect','MM_FILE','ANNUAL_3D', 'NN_SST','NN_T','NN_SSS','NN_S','F_T_OBS_3D_12','F_S_OBS_3D_12', 'F_SST_OBS_12','NN_SST_OBS','NN_T_OBS','NN_S_OBS','NM_TS_OBS'} vdic = bt.check_env_var(sys.argv[0], venv_needed) CONFEXP = vdic['ORCA']+'-'+vdic['EXP'] CC = vdic['COMP2D'] i_do_sect = int(vdic['i_do_sect']) l_3df_are_annual = False ; # we expect 3D fields (T & S) to be in monthly 'grid_T' mclim file if len(vdic['ANNUAL_3D']) != 0: l_3df_are_annual = True ; # 3D fields (T & S) are in annual 'grid_T' aclim file # Bounds and increment for comparison maps: if CC == 'OBS': cname_obs = vdic['NM_TS_OBS'] tmin=-5. ; tmax=5. ; dtemp = 0.5
import barakuda_plot as bp Socean = 363. ; # Surface of the ocean in 10^6 km^2 Lt_1y = 365.*24.*3600.*1E-6 ; # Length of 1 year in 10^6 seconds (=31.536) csn = sys.argv[0] cv_evb = 'evap_ao_cea' ; # debug evap in ec-earth... DEFAULT_LEGEND_LOC = 'lower left' venv_needed = {'ORCA','EXP','NN_SST','NN_SSS','NN_SSH','NN_T','NN_S','NN_MLD','LMOCLAT', 'TRANSPORT_SECTION_FILE','FIG_FORM','BM_FILE'} vdic = bt.check_env_var(csn, venv_needed) CONFEXP = vdic['ORCA']+'-'+vdic['EXP'] ff = vdic['FIG_FORM'] ; # format for figures (usually "png" or "svg") narg = len(sys.argv) if narg != 2: print 'Usage: {} <diag>'.format(csn) sys.exit(0) cdiag = sys.argv[1] print '\n '+csn+': diag => "'+cdiag+'"' if cdiag == 'mean_tos':
import barakuda_tool as bt import barakuda_ncio as bn import barakuda_orca as bo import barakuda_plot as bp cv_evb = 'evap_ao_cea' # debug evap in ec-earth... DEFAULT_LEGEND_LOC = 'lower left' venv_needed = { 'ORCA', 'RUN', 'NN_SST', 'NN_SSS', 'NN_SSH', 'NN_T', 'NN_S', 'NN_MLD', 'LMOCLAT', 'TRANSPORT_SECTION_FILE', 'FIG_FORM' } vdic = bt.check_env_var(sys.argv[0], venv_needed) CONFRUN = vdic['ORCA'] + '-' + vdic['RUN'] ff = vdic['FIG_FORM'] # format for figures (usually "png" or "svg") narg = len(sys.argv) if narg != 2: print 'Usage: {} <diag>'.format(sys.argv[0]) sys.exit(0) cdiag = sys.argv[1] print '\n plot_time_series.py: diag => "' + cdiag + '"' if cdiag == 'mean_tos':
from netCDF4 import Dataset import barakuda_tool as bt import barakuda_ncio as bn import barakuda_orca as bo import barakuda_plot as bp csn = sys.argv[0] cv_evb = 'evap_ao_cea' ; # debug evap in ec-earth... DEFAULT_LEGEND_LOC = 'lower left' venv_needed = {'ORCA','RUN','NN_SST','NN_SSS','NN_SSH','NN_T','NN_S','NN_MLD','LMOCLAT','TRANSPORT_SECTION_FILE','FIG_FORM'} vdic = bt.check_env_var(csn, venv_needed) CONFRUN = vdic['ORCA']+'-'+vdic['RUN'] ff = vdic['FIG_FORM'] ; # format for figures (usually "png" or "svg") def __test_nb_years__(nb_months, cd): if nb_months%12 != 0: print 'ERROR: '+csn+' for diag='+cd print ' => number of "supposed" monthly time records not a multiple of 12!' print ' => Nb. rec. = '+str(nb_months) sys.exit(0) return nb_months/12 narg = len(sys.argv) if narg != 2: