def test_time_slice(directory, sourcefiles, sd, ed):
ctt.loadtestfiles(directory, sourcefiles)
filepattern = "*.nc"
ens = cd.mkensemble(filepattern)
ens = cd.time_slice(ens, start_date=sd, end_date=ed)
data = ctt.sha(ens)
return data
def test_ensembleenvelope(directory, sourcefiles, var, sd, ed, cwd):
""" Tests ensemble_envelope_timeseries
Parameters
----------
directory : string
the directory to link the files to
sourcefiles : list
a list of the locations of the files to link
var : string
the variable plotted
sd : int
start date for time slice
ed : int
end date for time slice
cwd : string
directory of the test modules
"""
ctt.loadtestfiles(directory, sourcefiles)
filepattern = '*.nc'
ens = cd.mkensemble(filepattern)
ens = cd.time_slice(ens, sd, ed)
ens = cd.remap(ens, 'r360x180', delete=True)
ens = cd.my_operator(ens, my_cdo_str='cdo -fldmean -yearmean {infile} {outfile}', output_prefix='processed_', delete=True)
cd.ens_stats(ens, var)
cd.ensemble_envelope_timeseries(ens, 'ENS-MEAN_processed_remap_ts_Amon_CanESM2_historical_R-MEAN_195001-200012.nc', 'ENS-STD_processed_remap_ts_Amon_CanESM2_historical_STD_195001-200012.nc', var, pdf='ensenv')
plotcompare(directory, cwd, 'ensenv')
os.system('rm ensenv.pdf')
def test_time_slice(directory, sourcefiles, sd, ed):
ctt.loadtestfiles(directory, sourcefiles)
filepattern = '*.nc'
ens = cd.mkensemble(filepattern)
ens = cd.time_slice(ens, start_date=sd, end_date=ed)
data = ctt.sha(ens)
return data
def process_local_cmip5_data(var):
"""time-merge
"""
# 2. Join the time-slices
# First build a cmipdata ensemble object
filepattern = var + '*.nc'
ens = cd.mkensemble(filepattern)
# Join the time-slices and limit to years between 1880 and 2012
ens = cd.cat_experiments(ens, var, 'historical', 'rcp45')
ens = cd.time_slice(ens, '1881-01-01', '2012-12-31')
def test_ensembleenvelope(directory, sourcefiles, var, sd, ed, cwd):
""" Tests ensemble_envelope_timeseries
Parameters
----------
directory : string
the directory to link the files to
sourcefiles : list
a list of the locations of the files to link
var : string
the variable plotted
sd : int
start date for time slice
ed : int
end date for time slice
cwd : string
directory of the test modules
"""
ctt.loadtestfiles(directory, sourcefiles)
filepattern = '*.nc'
ens = cd.mkensemble(filepattern)
ens = cd.time_slice(ens, sd, ed)
ens = cd.remap(ens, 'r360x180', delete=True)
ens = cd.my_operator(
ens,
my_cdo_str='cdo -fldmean -yearmean {infile} {outfile}',
output_prefix='processed_',
delete=True)
cd.ens_stats(ens, var)
cd.ensemble_envelope_timeseries(
ens,
'ENS-MEAN_processed_remap_ts_Amon_CanESM2_historical_R-MEAN_195001-200012.nc',
'ENS-STD_processed_remap_ts_Amon_CanESM2_historical_STD_195001-200012.nc',
var,
pdf='ensenv')
plotcompare(directory, cwd, 'ensenv')
os.system('rm ensenv.pdf')
