@@ -0,0 +1,1558 @@
+%22%22%22A Dakotathon uncertainty quantification experiment with Hydrotrend.%0A%0AThis experiment requires a WMT executor with PyMT installed. It also%0Arequires Dakotathon and Hydrotrend installed as CSDMS components.%0A%0A%22%22%22%0Aimport os%0Aimport numpy as np%0Afrom pymt.components import PolynomialChaos, Hydrotrend%0Afrom dakotathon.utils import configure_parameters%0A%0A%0Amodel, dakota = Hydrotrend(), PolynomialChaos()%0A%0Aexperiment = %7B%0A    'component': type(model).__name__,%0A    'run_duration': 365,               # days%0A    'auxiliary_files': 'HYDRO0.HYPS',  # the default Waipaoa hypsometry%0A    'quadrature_order': 4,%0A    'samples': 10000,%0A    'seed': 17,%0A    'probability_levels': %5B0.05, 0.10, 0.33, 0.50, 0.67, 0.90, 0.95%5D,%0A    'variance_based_decomp': True,%0A    'descriptors': %5B'starting_mean_annual_temperature',%0A                    'total_annual_precipitation'%5D,%0A    'variable_type': 'uniform_uncertain',%0A    'lower_bounds': %5B12.8, 1.4%5D,%0A    'upper_bounds': %5B15.8, 1.8%5D,%0A    'response_descriptors': 'channel_exit_water_sediment~suspended__mass_flow_rate',%0A    'response_statistics': 'median',%0A    %7D%0Adakota_parameters, model_parameters = configure_parameters(experiment)%0A%0Adakota_parameters%5B'run_directory'%5D = model.setup(os.getcwd(), **model_parameters)%0A%0Acfg_file = 'HYDRO.IN'  # get from pymt eventually%0Adakota_tmpl_file = cfg_file + '.dtmpl'%0Aos.rename(cfg_file, dakota_tmpl_file)%0Adakota_parameters%5B'template_file'%5D = dakota_tmpl_file%0A%0Adakota.setup(dakota_parameters%5B'run_directory'%5D, **dakota_parameters)%0A%0Adakota.initialize('dakota.yaml')%0Adakota.update()%0Adakota.finalize()%0A

Add script for UQ experiment