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