merge_datasets.py
df1f2ca30d96master

merge_datasets.py
View Options

			import numpy as np
			import pandas as pd
			import xarray as xr
			import os
			import util

			target_path = '/work/hyenergy/output/meteo/results'

			var = 'SIS'
			model = 'ELM50_1000_CPU'
			queryname = 'dhm25_200_CH'
			datasets = ['prediction.nc', 'model_unc.nc', 'data_unc.nc']

			for dataset in datasets:
			# define lambda (path to all monthly datasets)
			# structure: dataset = lambda month : ("PATH_%02d_file" %month)
			pred = lambda month : os.path.join( '/scratch/walch/datasets', ('2001-%02d_all_%s_4D' %(month, var)),
			model, 'query', queryname, dataset)

			filename = os.path.join( target_path, '%s_%s_%s_%s' %(var, queryname, model, dataset))

			util.merge_models(pred, filename)

			print('Successfully saved merged dataset as %s' %filename)