ref_subregion_mean, ref_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std([ref_dataset], subregions)
model_subregion_mean, model_subregion_std, subregion_array = utils.calc_subregion_area_mean_and_std(model_datasets, subregions)
""" Step 7: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
workdir = workdir+'/'
print 'Writing a netcdf file: ',workdir+config['output_netcdf_filename']
if not os.path.exists(workdir):
os.system("mkdir "+workdir)
if config['use_subregions']:
dsp.write_netcdf_multiple_datasets_with_subregions(ref_dataset, ref_name, model_datasets, model_names,
path=workdir+config['output_netcdf_filename'],
subregions=subregions, subregion_array = subregion_array,
ref_subregion_mean=ref_subregion_mean, ref_subregion_std=ref_subregion_std,
model_subregion_mean=model_subregion_mean, model_subregion_std=model_subregion_std)
else:
dsp.write_netcdf_multiple_datasets_with_subregions(ref_dataset, ref_name, model_datasets, model_names,
path=workdir+config['output_netcdf_filename'])
""" Step 8: Calculate metrics and draw plots """
nmetrics = config['number_of_metrics_and_plots']
if config['use_subregions']:
Map_plot_subregion(subregions, ref_dataset, workdir)
if nmetrics > 0:
print 'Calculating metrics and generating plots'
for imetric in np.arange(nmetrics)+1:
metrics_name = config['metrics'+'%1d' %imetric]
utils.calc_subregion_area_mean_and_std(model_datasets, subregions))
""" Step 7: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
workdir = workdir + '/'
print 'Writing a netcdf file: ', workdir + config['output_netcdf_filename']
if not os.path.exists(workdir):
os.system("mkdir -p " + workdir)
if config['use_subregions']:
dsp.write_netcdf_multiple_datasets_with_subregions(
obs_dataset,
obs_name,
model_datasets,
model_names,
path=workdir + config['output_netcdf_filename'],
subregions=subregions,
subregion_array=subregion_array,
ref_subregion_mean=obs_subregion_mean,
ref_subregion_std=obs_subregion_std,
model_subregion_mean=model_subregion_mean,
model_subregion_std=model_subregion_std)
else:
dsp.write_netcdf_multiple_datasets_with_subregions(
obs_dataset,
obs_name,
model_datasets,
model_names,
path=workdir + config['output_netcdf_filename'])
""" Step 8: Calculate metrics and draw plots """
nmetrics = config['number_of_metrics_and_plots']
if config['use_subregions']:
utils.calc_subregion_area_mean_and_std(target_datasets, subregions))
""" Step 6: Write a netCDF file """
workdir = config['workdir']
if workdir[-1] != '/':
workdir = workdir + '/'
print('Writing a netcdf file: ', workdir + config['output_netcdf_filename'])
if not os.path.exists(workdir):
os.system("mkdir -p " + workdir)
if config['use_subregions']:
dsp.write_netcdf_multiple_datasets_with_subregions(
reference_dataset,
reference_name,
target_datasets,
target_names,
path=workdir + config['output_netcdf_filename'],
subregions=subregions,
subregion_array=subregion_array,
ref_subregion_mean=reference_subregion_mean,
ref_subregion_std=reference_subregion_std,
model_subregion_mean=target_subregion_mean,
model_subregion_std=target_subregion_std)
else:
dsp.write_netcdf_multiple_datasets_with_subregions(
reference_dataset,
reference_name,
target_datasets,
target_names,
path=workdir + config['output_netcdf_filename'])
""" Step 7: Calculate metrics and draw plots """
nmetrics = config['number_of_metrics_and_plots']
if config['use_subregions']: