Stochastic Harmonic Autoregressive Parametric (SHArP) Weather Generator
|Authors:||Kimberly Smith Courtenay Strong Firas Rassoul-Agha|
|Owners:||Kimberly Smith Courtenay Strong iUTAH Data Manager|
|Resource type:||Script Resource|
|Storage:||The size of this resource is 45.9 KB|
|Created:||Jun 25, 2018 at 8:21 p.m.|
|Last updated:||Feb 20, 2019 at 2:33 a.m. by Kimberly Smith|
|Citation:||See how to cite this resource|
This collection of MATLAB scripts runs the temperature component of the stochastic harmonic autoregressive parametric (SHArP) weather generator, which simulates trended, nonstationary precipitation occurrence, precipitation amount, and temperature values. The precipitation component is not included here because it is based off of previous weather generators and doesn't use any new methods. The temperature component follows a concept introduced in previous studies, but SHArP simulates temperature values directly as opposed to simulating the temperature residuals. SHArP weather generator validation and illustrations at a single site have been published in the Journal of Applied Meteorology and Climatology (Smith et al. 2017; "A new method for generating stochastic simulations of daily air temperature for use in weather generators"; https://doi.org/10.1175/JAMC-D-16-0122.1), and the multisite generalization of SHArP is in press in the same journal (Smith et al. 2018; "Multisite generalization of the SHArP weather generator"; https://doi.org/10.1175/JAMC-D-17-0236.1).
How to use:
- All input data need to be combined into a single .mat structure file and contain the number of locations to be simulated, the number of days to simulate, a time matrix (in format: year month day), the day of year corresponding to each day, the precipitation occurrence 1s and 0s in logical format corresponding to "wet" or "dry" days (determined from observations or downscaled climate model output), the maximum temperature training data for each site (i.e. observations or downscaled climate model output), and the minimum temperature training data for each site.
- This .mat structure file (called "sett") is the input to fit_SHArP_1, which fits the parameters needed to run the model. fit_SHArP_1 internally uses fit_SHArP_var_1a to determine the coefficients for each of the 26 variables per site in the model; this latter script is not run separately. The output from fit_SHArP_1 is called the driver_t structure and needs to contain maximum and minimum temperature data (from observations or downscaled climate model output), the precipitation occurrence 1s and 0s in logical format corresponding to "wet" or "dry" days (determined from observations or downscaled climate model output), a time matrix (in format: year month day), a variable structure that contains the 26 coefficients for each variable at each site, b_k, c_k, and the results from the EOF analysis on the precipitation occurrence patterns.
- The driver_t structure is the input to simulate_SHArP_2, which gives maximum temperature and minimum temperature as output. This script can also be modified to output the temperature mean and noise but note that this will increase computation time significantly.
|The content of this resource serves as the data for:||https://doi.org/10.1175/JAMC-D-16-0122.1|
|The content of this resource serves as the data for:||https://doi.org/10.1175/JAMC-D-17-0236.1|
This resource was created using funding from the following sources:
|Agency Name||Award Title||Award Number|
|National Science Foundation||Random Polymer Measures||DMS-1407574|
|National Science Foundation||iUTAH-innovative Urban Transitions and Aridregion Hydro-sustainability||EPS-1208732|
|National Science Foundation||Collaborative Research: CI-WATER, Cyberinfrastructure to Advance High Performance Water Resource Modeling||EPS-1135483|
|National Science Foundation||Collaborative Research: CI-WATER, Cyberinfrastructure to Advance High Performance Water Resource Modeling||EPS-1135482|
How to Cite
This resource is shared under the Creative Commons Attribution CC BY.http://creativecommons.org/licenses/by/4.0/
Please wait for the process to complete.