├── .gitignore ├── EXAMPLE.md ├── LICENSE ├── README.md ├── docs └── SHAPEFILE_FIELD_EXPLANATIONS.md ├── examples ├── README.md ├── hru_param_explanations.xlsx ├── output_compare │ ├── gsflow_prms.out │ └── sagehen_statvar.dat ├── sagehen │ ├── BUILDING.md │ ├── Troubleshoot_notes.txt │ ├── dem │ │ └── ned10m_nad83.img │ ├── hru_params │ │ ├── prms_dimensions.csv │ │ └── prms_parameters.csv │ ├── impervious │ │ └── nlcd2011_imp.img │ ├── model │ │ ├── input │ │ │ └── prms │ │ │ │ └── sagehen_datafile_1sta.data │ │ └── windows │ │ │ ├── gsflow.bat │ │ │ └── sagehen.control │ ├── nhd │ │ ├── NHDFlowline.dbf │ │ ├── NHDFlowline.prj │ │ ├── NHDFlowline.sbn │ │ ├── NHDFlowline.sbx │ │ ├── NHDFlowline.shp │ │ ├── NHDFlowline.shx │ │ ├── WBDHU12.dbf │ │ ├── WBDHU12.prj │ │ ├── WBDHU12.sbn │ │ ├── WBDHU12.sbx │ │ ├── WBDHU12.shp │ │ └── WBDHU12.shx │ ├── prism │ │ ├── ppt │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_01_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_02_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_03_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_04_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_05_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_06_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_07_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_08_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_09_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_10_bil.bil │ │ │ ├── PRISM_ppt_30yr_normal_800mM2_11_bil.bil │ │ │ └── PRISM_ppt_30yr_normal_800mM2_12_bil.bil │ │ ├── tmax │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_01_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_02_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_03_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_04_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_05_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_06_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_07_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_08_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_09_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_10_bil.bil │ │ │ ├── PRISM_tmax_30yr_normal_800mM2_11_bil.bil │ │ │ └── PRISM_tmax_30yr_normal_800mM2_12_bil.bil │ │ └── tmin │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_01_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_02_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_03_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_04_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_05_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_06_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_07_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_08_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_09_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_10_bil.bil │ │ │ ├── PRISM_tmin_30yr_normal_800mM2_11_bil.bil │ │ │ └── PRISM_tmin_30yr_normal_800mM2_12_bil.bil │ ├── readme.txt │ ├── runscripts.bat │ ├── runscripts_python_removed.bat │ ├── sagehen_parameters.ini │ ├── save │ │ ├── calibration_parameters.param.save │ │ ├── cascade.param.save │ │ ├── default_values.param.save │ │ ├── gis_derived_parameters.param.save │ │ ├── sagehen.control.save │ │ └── sagehen_datafile_1sta.data.save │ ├── shapefiles │ │ ├── model_points.dbf │ │ ├── model_points.prj │ │ ├── model_points.sbn │ │ ├── model_points.sbx │ │ ├── model_points.shp │ │ ├── model_points.shx │ │ ├── watershed.dbf │ │ ├── watershed.prj │ │ ├── watershed.sbn │ │ ├── watershed.sbx │ │ ├── watershed.shp │ │ └── watershed.shx │ ├── soils │ │ ├── awc.img │ │ ├── clay.img │ │ ├── ksat.img │ │ └── sand.img │ └── veg │ │ ├── us_140evc.img │ │ ├── us_140evc.img.vat.dbf │ │ ├── us_140evt.img │ │ └── us_140evt.img.vat.dbf ├── sagehen_tests │ ├── README.md │ ├── active_fill │ │ ├── dem │ │ │ └── ned90m_nad83z11.img │ │ ├── shapefiles │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ │ └── test_parameters_fill.ini │ ├── clipped_elev │ │ ├── dem │ │ │ └── ned90m_nad83z11_clip.img │ │ ├── shapefiles │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ │ └── test_parameters_elev_clip.ini │ ├── lake │ │ ├── dem │ │ │ └── ned90m_nad83z11.img │ │ ├── lake_parameters.ini │ │ ├── runscripts.bat │ │ └── shapefiles │ │ │ ├── lakes.dbf │ │ │ ├── lakes.prj │ │ │ ├── lakes.sbn │ │ │ ├── lakes.sbx │ │ │ ├── lakes.shp │ │ │ ├── lakes.shx │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ ├── precipitation_zones │ │ ├── dem │ │ │ └── ned90m_nad83z11.img │ │ ├── hru_params │ │ │ ├── prms_dimensions.csv │ │ │ └── prms_parameters.csv │ │ ├── ppt_zones_parameters.ini │ │ ├── runscripts.bat │ │ └── shapefiles │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── ppt_zones.dbf │ │ │ ├── ppt_zones.prj │ │ │ ├── ppt_zones.shp │ │ │ ├── ppt_zones.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ ├── temperature_zones │ │ ├── dem │ │ │ └── ned90m_nad83z11.img │ │ ├── hru_params │ │ │ ├── prms_dimensions.csv │ │ │ └── prms_parameters.csv │ │ ├── runscripts.bat │ │ ├── shapefiles │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── temp_zones.dbf │ │ │ ├── temp_zones.prj │ │ │ ├── temp_zones.sbn │ │ │ ├── temp_zones.sbx │ │ │ ├── temp_zones.shp │ │ │ ├── temp_zones.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ │ ├── temp_1sta_parameters.ini │ │ └── temp_zones_parameters.ini │ ├── terminal_basin │ │ ├── dem │ │ │ └── ned90m_nad83z11_terminal_basin.img │ │ ├── hru_params │ │ │ ├── prms_dimensions.csv │ │ │ └── prms_parameters.csv │ │ ├── runscripts.bat │ │ ├── shapefiles │ │ │ ├── model_points.dbf │ │ │ ├── model_points.prj │ │ │ ├── model_points.sbn │ │ │ ├── model_points.sbx │ │ │ ├── model_points.shp │ │ │ ├── model_points.shx │ │ │ ├── watershed.dbf │ │ │ ├── watershed.prj │ │ │ ├── watershed.sbn │ │ │ ├── watershed.sbx │ │ │ ├── watershed.shp │ │ │ └── watershed.shx │ │ └── terminal_basin_parameters.ini │ └── terminal_lake │ │ ├── dem │ │ └── ned90m_nad83z11_terminal_lake.img │ │ ├── hru_params │ │ ├── prms_dimensions.csv │ │ └── prms_parameters.csv │ │ ├── runscripts.bat │ │ ├── shapefiles │ │ ├── lakes.dbf │ │ ├── lakes.prj │ │ ├── lakes.sbn │ │ ├── lakes.sbx │ │ ├── lakes.shp │ │ ├── lakes.shx │ │ ├── model_points.dbf │ │ ├── model_points.prj │ │ ├── model_points.sbn │ │ ├── model_points.sbx │ │ ├── model_points.shp │ │ ├── model_points.shx │ │ ├── watershed.dbf │ │ ├── watershed.prj │ │ ├── watershed.sbn │ │ ├── watershed.sbx │ │ ├── watershed.shp │ │ └── watershed.shx │ │ └── terminal_lake_parameters.ini └── template │ └── template_parameters.ini ├── gsflow_arcpy_tutorial.pdf ├── misc ├── MISC.md ├── dem_adj.lyr └── flow_dir_points.lyr ├── remaps ├── REMAPS.md ├── arc10p1 │ ├── aspect.rmp │ ├── convert_remap_arc10p2.py │ ├── covdensum.rmp │ ├── covdenwin.rmp │ ├── covtype.rmp │ ├── rtdepth.rmp │ ├── snow_intcp.rmp │ ├── srain_intcp.rmp │ ├── temp_adj_x10.rmp │ └── wrain_intcp.rmp ├── aspect.rmp ├── covdensum.rmp ├── covdenwin.rmp ├── covtype.rmp ├── covtype_to_snow_intcp.rmp ├── covtype_to_srain_intcp.rmp ├── covtype_to_wrain_intcp.rmp ├── rtdepth.rmp ├── snow_intcp.rmp ├── srain_intcp.rmp ├── temp_adj_x10.rmp └── wrain_intcp.rmp └── scripts ├── crt_fill_parameters.py ├── daymet_normals.py ├── dem_2_streams.py ├── dem_parameters.py ├── field_list.ini ├── fishnet_generator.py ├── hru_parameters.py ├── impervious_parameters.py ├── ppt_ratio_parameters.py ├── prism_4km_normals.py ├── prism_800m_normals.py ├── prms_template_fill.py ├── soil_parameters.py ├── soil_raster_prep.py ├── stream_parameters.py ├── support_functions.py ├── temp_adjust_parameters.py ├── template_parameters.ini ├── template_prms_dimensions.csv ├── template_prms_parameters.csv └── veg_parameters.py /.gitignore: -------------------------------------------------------------------------------- 1 | # Python specific 2 | *.pyc 3 | 4 | # Generic files to ignore 5 | *~ 6 | *.lock 7 | *.bak 8 | log 9 | 10 | # PyCharm 11 | .idea 12 | 13 | # CRT 14 | CRT*.exe 15 | 16 | # PRMS/GSFlow model files and folders 17 | *.control 18 | *.dat 19 | *.data 20 | *.out 21 | *.log 22 | 23 | # Model parameter files, folders, and shapefiles 24 | # hru_params 25 | hru_params.* 26 | hru_params_label.* 27 | cascade_work 28 | cascade_gw_work 29 | dem_rasters 30 | fill_work 31 | fishnet_temp 32 | flow_rasters 33 | hru_temp 34 | impervious_rasters 35 | logs 36 | ppt_rasters 37 | ppt_ratio 38 | soil_rasters 39 | stream_rasters 40 | temp_adjust 41 | tmax_rasters 42 | tmin_rasters 43 | veg_rasters 44 | *.param 45 | 46 | # Project specific files and file types 47 | *.mxd 48 | *.img.aux.xml 49 | *.shp.xml 50 | *.cpg 51 | 52 | # Script folders 53 | scripts/custom 54 | scripts/old 55 | 56 | # Ancillary data folders 57 | impervious 58 | nlcd 59 | prism 60 | soils 61 | veg 62 | vegetation 63 | 64 | # Exclude all executables 65 | *.exe 66 | -------------------------------------------------------------------------------- /EXAMPLE.md: -------------------------------------------------------------------------------- 1 | # Example 2 | Create the project folder if necessary (i.e. D:\Projects\sagehen) and create an hru_params sub-folder (i.e. D:\Projects\sagehen\hru_params). 3 | 4 | #### Clone the repository 5 | 6 | If you already have a local repository of the gsflow-arcpy-tools, pull the latest version from GitHub. If you don't already have a local repository, either clone the repository locally or download a zip file of the scripts from Github. For this example, the local copy will be cloned directly in the project folder (i.e. D:\Projects\sagehen\scripts and D:\Projects\sagehen\remaps) 7 | 8 | #### Study Area 9 | 10 | Create a "shapefiles" sub-folder (i.e. D:\Projects\sagehen\shapefiles). 11 | Place the study area shapefile in the shapefiles sub-folder (ie.e. D:\Projects\sagehen\shapefiles\watershed.shp). The study area is typically a single watershed (like a HUC8) or a collection of watersheds that will be modeled together. Check the projection of the study area shapefile and make sure it is in a projected coordinate system (like NAD83 UTM Zone 10N). 12 | 13 | ##### Input File 14 | 15 | Copy the template input file (scripts\template_parameters.ini) to the hru_params folder and rename (i.e. D:\Projects\sagehen\hru_params\example_params.ini) 16 | Within the input file, all of the file and folder paths need to be set to the project folder. A simple find and replace of "D:\Projects\gsflow-arcpy-example" to your project folder (i.e. "D:\Projects\sagehen") should work, but doublecheck all of the paths. Eventually, all of the paths will be set as relative paths to the main project folder, but for now they must be absolute paths. 17 | 18 | ##### Fishnet / Model Grid 19 | 20 | The following parameters must be explicitly set whether you are building a new fishnet/grid shapefile or reading in an existing one. If reading an existing grid, the paramters must match the grid properties exactly. 21 | - hru_cellsize - cellsize (units will depend on projection) 22 | - hru_ref_x: snap x coordinate (units will depend on projection, 0 is a good choice for a new grid) 23 | - hru_ref_x: snap y corodinate (units will depend on projection, 0 is a good choice for a new grid) 24 | - hru_projection: EPSG code (see spatialreference.org) 25 | - study_area_path: The full filepath of your study area shapefile 26 | 27 | ##### Running the Scripts 28 | 29 | Currently, the scripts must be run from the windows command prompt so that the input file can be passed as an argument directly to the script. Eventually, the scipts will be modified so they can be executed by double clicking on the script in Explorer and the input file will be set through a dialog box.
30 | First open the windows command prompt by pressing the "windows" key and "r" or clicking the "Start" button, "Accessories", and then "Command Prompt". 31 | 32 | Within the command prompt, change to the target drive if necessary: 33 | ``` 34 | > D: 35 | ``` 36 | 37 | Navigate to the project folder: 38 | ``` 39 | > cd D:\Projects\sagehen 40 | ``` 41 | Run the fisnet_generator.py script if building a new fishnet. Note, this will overwrite an existing fishnet shapefile if one already exists at the path specified in the input file. 42 | ``` 43 | > python scripts\fishnet_generator.py -i hru_params\example_params.ini 44 | ``` 45 | 46 | ##### Digitial Elevation Model (DEM) 47 | 48 | If you don't already have a DEM raster of the study area, download the DEM tiles using the download_ned.py script. You must set the study area shapefile path using the "--extent" argument and the output folder using the "--output" argument. The script will project the shapefile to a geographic coordinate system and download all NED 1x1 degree tiles that intersect the study area. 49 | ``` 50 | > python tools\download_ned.py --extent shapefiles\watershed.shp --output dem\tiles 51 | ``` 52 | After downloading the tiles, they need to be merged together into a single image and then projected to the fishnet (or watershed) spatial reference. This could be done in ArcGIS or using the GDAL utilities (http://www.gdal.org/gdal_utilities.html) at the command line. 53 | 54 | First, the gdal_merge tool can be used to merge the tiles into a single raster. 55 | ``` 56 | > gdal_merge.py -of HFA -co COMPRESSED=YES -o dem\ned_30m_merge.img dem\tiles\imgn39w107_1.img dem\tiles\imgn39w108_1.img dem\tiles\imgn40w107_1.img dem\tiles\imgn40w108_1.img 57 | ``` 58 | 59 | It is important that the elevation units match the linear unit of the coordinate system. If the fishnet is in a coordinate system with feet as the linear unit, scale the NED rasters by 0.3048 to convert meters to feet. This can be done using the gdal_calc utility 60 | ``` 61 | gdal_calc.py -A dem\ned_30m_merge.img --outfile=dem\ned_nad83_feet.img --calc="0.3048*A" --format HFA --co COMPRESSED=YES 62 | ``` 63 | 64 | Project the raster to the fishnet (or study area) spatial reference using the "gdalwarp" utility. 65 | ``` 66 | gdalwarp -r "bilinear" -tr 30 30 -s_srs "EPSG:4269" -t_srs "EPSG:26910" -overwrite -ot Float32 -srcnodata None -dstnodata -3.4028234663852886e+38 -of HFA -co COMPRESSED=YES dem\ned_30m_merge.img dem\ned_30m.img 67 | ``` 68 | The output cellsize may need to be adjusted using "-tr" argument for coordinate systems that are not in meters, otherwise 30 should work well for most applications. The output EPGS code (set with "-t_srs") needs to match the fishnet EPSG set in the input file. To subset the raster at this step, use the "-te xmin ymin xmax ymax" argument with coordinates in the projected coordinate system of the fishnet. 69 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | GSFlow ArcPy 2 | ============ 3 | 4 | Series of Python/ArcPy (ArcGIS) scripts for developing inputs for a GSFLOW model. 5 | 6 | ## Script Execution Order 7 | - fishnet_generator.py 8 | - hru_parameters.py 9 | - dem_parameters.py 10 | - veg_parameters.py 11 | - soil_raster_prep.py 12 | - soil_parameters.py 13 | - impervious_parameters.py 14 | - prism_4km_normals.py / prism_800m_normals.py 15 | - ppt_ratio_parameters.py 16 | - *Iterate to define the stream network* 17 | - dem_2_streams.py 18 | - crt_fill_parameters.py 19 | - stream_parameters.py 20 | - prms_template_fill.py 21 | 22 | ## Ancillary Data 23 | 24 | Almost all of the following data can be downloaded for a study area using the [USGS Geospatial Data Gateway](http://datagateway.nrcs.usda.gov/). They could probably also be downloaded using the [National Map Viewer](http://viewer.nationalmap.gov/viewer/), but this has not been tested. Specific download instructions are provided for each dataset below. 25 | 26 | #### Elevation 27 | 28 | Elevation data is set using the 10m (1/3 arc-second) or 30m (1 arc-second) National Elevation Dataset (NED) rasters. These can be easily downloaded in 1x1 degree tiles for the CONUS from the [USGS FTP](rockyftp.cr.usgs.gov) in the folder vdelivery/Datasets/Staged/Elevation. 29 | 30 | #### LANDFIRE 31 | 32 | Vegetation type and cover percent is set using [LANDFIRE](http://www.landfire.gov/). The data can be downloaded in [tiles](http://www.landfire.gov/viewer/) or for the entire [CONUS](http://www.landfire.gov/lf_mosaics.php) (but only version 1.3.0?). The default remap files were developed and tested for version 1.2.0 LANDFIRE 2010 version 1.2.0. 33 | 34 | #### Soils 35 | 36 | Available water capacity (AWC), percent sand, percent clay, and saturated hydraulic conductivity (Ksat) can be set using SSURGO or STATSGO. The easiest way to acquire these data are through the [USGS Geospatial Data Gateway](http://datagateway.nrcs.usda.gov/). Shapefiles of the soil properties can be extracted using the [NRCS Soil Data Viewer](http://www.nrcs.usda.gov/wps/portal/nrcs/detailfull/soils/home/?cid=nrcs142p2_053620). The shapefiles then need to be converted to raster. 37 | 38 | #### PRISM 39 | 40 | PRISM precipitation, minimum temperature, and maximum temperature 30 year normals for the CONUS can be downloaded from the [PRISM site](http://www.prism.oregonstate.edu/normals/). 41 | 42 | #### CRT 43 | 44 | User must have [Cascade Routing Tool](http://water.usgs.gov/ogw/CRT/) (CRT) version 1.3.1 45 | 46 | #### Remap files 47 | 48 | Example ASCII remap files are provided, although it may be necessary to modify these to include new or missing LANDFIRE vegetation types. In versions of ArcGIS before 10.2, you could have comments after the values (indicated by a /*) but this was removed in 10.2. Now, comments must be on a separate line and begin with the "#" symbol. The convert_remap_10p2.py script will convert the ArcGIS 10.1 style ASCII remap files to the ArcGIS 10.2 style. 49 | 50 | ## Requirements 51 | 52 | + Python 2.7 53 | + ArcPy 54 | 55 | ## References 56 | 57 | Gardner, M., Morton, C., Huntington, J., Niswonger, R., Henson, W. (2018). Input data processing tools for the integrated hydrologic model GSFLOW. Environmental *Modelling & Software*, 109. [https://doi.org/10.1016/j.envsoft.2018.07.020](https://doi.org/10.1016/j.envsoft.2018.07.020) 58 | -------------------------------------------------------------------------------- /examples/README.md: -------------------------------------------------------------------------------- 1 | # Example Models 2 | 3 | ## Sagehen 4 | 5 | ## Template 6 | -------------------------------------------------------------------------------- /examples/hru_param_explanations.xlsx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/hru_param_explanations.xlsx -------------------------------------------------------------------------------- /examples/output_compare/gsflow_prms.out: -------------------------------------------------------------------------------- 1 | 2 | Precipitation-Runoff Modeling System (PRMS) 3 | Version 5.0.0 06/01/2018 4 | 5 | Process Available Modules 6 | -------------------------------------------------------------------- 7 | Basin Definition: basin 8 | Cascading Flow: cascade 9 | Time Series Data: obs, water_use_read, dynamic_param_read 10 | Potet Solar Rad: soltab 11 | Temperature Dist: temp_1sta, temp_laps, temp_dist2, climate_hru 12 | Precip Dist: precip_1sta, precip_laps, precip_dist2, 13 | climate_hru 14 | Temp & Precip Dist: xyz_dist, ide_dist 15 | Solar Rad Dist: ccsolrad, ddsolrad, climate_hru 16 | Transpiration Dist: transp_tindex, climate_hru, transp_frost 17 | Potential ET: potet_hamon, potet_jh, potet_pan, climate_hru, 18 | potet_hs, potet_pt, potet_pm, potet_pm_sta 19 | Interception: intcp 20 | Snow Dynamics: snowcomp 21 | Surface Runoff: srunoff_smidx, srunoff_carea 22 | Soil Zone: soilzone 23 | Groundwater: gwflow 24 | Streamflow Routing: strmflow, strmflow_in_out, muskingum, 25 | muskingum_lake 26 | Stream Temperature: stream_temp 27 | Output Summary: basin_sum, subbasin, map_results, prms_summary, 28 | nhru_summary, nsub_summary, water_balance 29 | basin_summary, nsegment_summary 30 | Preprocessing: write_climate_hru, frost_date 31 | -------------------------------------------------------------------- 32 | 33 | 34 | Active modules listed in the order in which they are called 35 | 36 | Process Module Version Date 37 | ==================================================================== 38 | GSFLOW Computation Order gsflow_prms 2018-05-25 39 | Basin Definition basin 2018-04-06 40 | Cascading Flow cascade 2018-05-25 41 | Common States and Fluxes climateflow 2018-04-25 42 | Potential Solar Radiation soltab 2016-09-09 43 | PRMS Set Time Variables prms_time 2015-03-31 44 | Time Series Data obs 2018-02-23 45 | Temperature Distribution temp_1sta 2018-01-16 46 | Precipitation Distribution precip_1sta 2016-10-21 47 | Solar Radiation Distribution ddsolrad 2016-11-03 48 | Transpiration Distribution transp_tindex 2015-01-06 49 | Potential Evapotranspiration potet_jh 2016-05-10 50 | Canopy Interception intcp 2018-02-26 51 | Snow Dynamics snowcomp 2018-05-04 52 | Surface Runoff srunoff_smidx 2018-04-25 53 | Soil Zone Computations soilzone 2018-04-25 54 | Groundwater gwflow 2018-02-26 55 | Streamflow Routing strmflow 2014-12-02 56 | Summary basin_sum 2017-10-21 57 | Output Summary map_results 2018-04-25 58 | Output Summary subbasin 2018-04-25 59 | ==================================================================== 60 | 61 | Using Control File: sagehen.control 62 | 63 | Using Parameter File: ./inputs/calibration_parameters.param 64 | 65 | 66 | Start time: 1982/08/01 00:00:00 67 | End time: 1997/03/31 00:00:00 68 | 69 | 70 | Model domain area: 13114.18 Active basin area: 6855.32 71 | Fraction impervious: 0.0011 Fraction pervious: 0.9989 72 | 73 | 74 | 75 | 76 | Year Month Day Precip ET Storage S-Runoff M-Runoff 77 | (inches) (inches) (inches) (inches) (inches) 78 | ------------------------------------------------------------- 79 | initial 2.417 80 | 1982 5.670 1.232 5.531 0.469 1.188 81 | 1983 65.257 22.482 8.598 31.477 37.966 82 | 1984 50.765 21.767 4.239 22.405 25.890 83 | 1985 29.253 15.762 3.209 8.376 11.141 84 | 1986 63.422 22.498 6.975 30.017 26.756 85 | 1987 18.430 13.804 1.608 3.776 6.370 86 | 1988 19.422 13.229 1.609 2.344 4.096 87 | 1989 47.413 19.449 5.563 17.748 12.548 88 | 1990 29.125 16.063 4.029 7.297 6.422 89 | 1991 29.454 17.970 3.267 6.927 5.515 90 | 1992 21.883 13.612 1.889 4.393 3.806 91 | 1993 49.247 18.694 4.280 22.158 18.685 92 | 1994 20.541 13.459 2.200 3.341 4.138 93 | 1995 65.579 20.940 7.136 33.138 29.179 94 | 1996 52.437 20.360 4.606 23.576 22.527 95 | ************************************************************* 96 | Total for run 604.014 253.953 26.152 224.944 232.203 97 | ************************************************************* 98 | Execution elapsed time 0 minutes 12.44 seconds 99 | 100 | -------------------------------------------------------------------------------- /examples/sagehen/Troubleshoot_notes.txt: -------------------------------------------------------------------------------- 1 | If getting the following error then you need to add 'python' in your envrionmental variables. 2 | 3 | error: 'python' is not recognized as an internal or external command, operable program or batch file. 4 | 5 | To enter the path for python to your environmental variables go to: 6 | 7 | control panel > system > advanced > |Environmental Variables| > system variables -> Path 8 | 9 | The path needs to include: C:\Python26; (or equivalent). 10 | 11 | Then restart your command file prompt, and try typing 'python'. 12 | 13 | The other option is to remove the 'python' from each command line in the batch file. 14 | 15 | This has been done for the alternate batch file called 'runscripts_python_removed.bat. 16 | 17 | 18 | If there are lakes present in a model then then 'nlakes' and 'nlake_hrus' need to be 19 | set to 'calculated' in the file prms_dimensions.csv 20 | 21 | If this is not done correctly then the following error occurs: 22 | 23 | [int(row[1]) for row in s_cursor if int(row[1]) > 0])) 24 | ValueError: max() arg is an empty sequence -------------------------------------------------------------------------------- /examples/sagehen/dem/ned10m_nad83.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/dem/ned10m_nad83.img -------------------------------------------------------------------------------- /examples/sagehen/hru_params/prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,1 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,0 18 | nsnow,calibration_parameters,0 19 | nsol,calibration_parameters,0 20 | nssr,calibration_parameters,calculated 21 | nsub,calibration_parameters,calculated 22 | ntemp,calibration_parameters,config_file 23 | one,calibration_parameters,1 24 | nlake_hrus,calibration_parameters,0 25 | -------------------------------------------------------------------------------- /examples/sagehen/impervious/nlcd2011_imp.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/impervious/nlcd2011_imp.img -------------------------------------------------------------------------------- /examples/sagehen/model/windows/gsflow.bat: -------------------------------------------------------------------------------- 1 | @ECHO OFF 2 | 3 | ..\..\..\..\gsflow\gsflow.exe .\sagehen.control 4 | ECHO. 5 | ECHO Run complete. Please press enter when you want to continue. 6 | PAUSE>NUL -------------------------------------------------------------------------------- /examples/sagehen/model/windows/sagehen.control: -------------------------------------------------------------------------------- 1 | PRMS Control File for Sagehen 2 | #### 3 | mapOutVar_names 4 | 1 5 | 4 6 | ./output/recharge 7 | #### 8 | nmapOutVars 9 | 1 10 | 1 11 | 0 12 | #### 13 | mapOutON_OFF 14 | 1 15 | 1 16 | 1 17 | #### 18 | parameter_check_flag 19 | 1 20 | 1 21 | 0 22 | #### 23 | prms_warmup 24 | 1 25 | 1 26 | 3 27 | #### 28 | aniOutON_OFF 29 | 1 30 | 1 31 | 0 32 | #### 33 | print_debug 34 | 1 35 | 1 36 | -1 37 | #### 38 | aniOutVar_names 39 | 5 40 | 4 41 | upslope_interflow 42 | upslope_dunnianflow 43 | soil_moist_tot 44 | upslope_hortonian 45 | sroff 46 | #### 47 | ani_output_file 48 | 1 49 | 4 50 | animation.out 51 | #### 52 | csv_output_file 53 | 1 54 | 4 55 | gsflow.csv 56 | #### 57 | model_output_file 58 | 1 59 | 4 60 | ../output/gsflow_prms.out 61 | #### 62 | data_file 63 | 1 64 | 4 65 | ../input/prms/sagehen_datafile_1sta.data 66 | #### 67 | dispGraphsBuffSize 68 | 1 69 | 4 70 | 15 71 | #### 72 | dispVar_element 73 | 15 74 | 4 75 | 2 76 | 2 77 | 2 78 | 2 79 | 2 80 | 1 81 | 2 82 | 1 83 | 1 84 | 1 85 | 2 86 | 1 87 | 1 88 | 1 89 | 1 90 | #### 91 | dispVar_names 92 | 15 93 | 4 94 | runoff 95 | sub_cfs 96 | sub_gwflow 97 | sub_interflow 98 | sub_sroff 99 | basin_pweqv 100 | subinc_pkweqv 101 | basin_dunnian 102 | basin_prefflow 103 | basin_slowflow 104 | subinc_precip 105 | basin_rain 106 | basin_snow 107 | basin_tmax 108 | basin_tmin 109 | #### 110 | dispVar_plot 111 | 15 112 | 1 113 | 1 114 | 1 115 | 2 116 | 2 117 | 2 118 | 3 119 | 3 120 | 4 121 | 4 122 | 4 123 | 5 124 | 5 125 | 5 126 | 6 127 | 6 128 | #### 129 | end_time 130 | 6 131 | 1 132 | 1997 133 | 3 134 | 31 135 | 0 136 | 0 137 | 0 138 | #### 139 | et_module 140 | 1 141 | 4 142 | potet_jh 143 | #### 144 | executable_desc 145 | 1 146 | 4 147 | GSFLOW model 148 | #### 149 | executable_model 150 | 1 151 | 4 152 | prmsIV.exe 153 | #### 154 | gsflow_output_file 155 | 1 156 | 4 157 | ../output/gsflow.out 158 | #### 159 | init_vars_from_file 160 | 1 161 | 1 162 | 0 163 | #### 164 | initial_deltat 165 | 1 166 | 2 167 | 24.0 168 | #### 169 | mms_user_dir 170 | 1 171 | 4 172 | ./ 173 | #### 174 | mms_user_out_dir 175 | 1 176 | 4 177 | ./ 178 | #### 179 | model_mode 180 | 1 181 | 4 182 | PRMS 183 | #### 184 | modflow_name 185 | 1 186 | 4 187 | sagehen.nam 188 | #### 189 | naniOutVars 190 | 1 191 | 1 192 | 5 193 | #### 194 | ndispGraphs 195 | 1 196 | 1 197 | 6 198 | #### 199 | param_file 200 | 4 201 | 4 202 | ../input/prms/calibration_parameters.param 203 | ../input/prms/cascade.param 204 | ../input/prms/default_values.param 205 | ../input/prms/gis_derived_parameters.param 206 | #### 207 | precip_module 208 | 1 209 | 4 210 | precip_1sta 211 | #### 212 | param_print_file 213 | 1 214 | 4 215 | ../output/sagehen.parprt 216 | #### 217 | save_vars_to_file 218 | 1 219 | 1 220 | 1 221 | #### 222 | solrad_module 223 | 1 224 | 4 225 | ddsolrad 226 | #### 227 | srunoff_module 228 | 1 229 | 4 230 | srunoff_smidx 231 | #### 232 | start_time 233 | 6 234 | 1 235 | 1982 236 | 8 237 | 1 238 | 0 239 | 0 240 | 0 241 | #### 242 | nstatVars 243 | 1 244 | 1 245 | 12 246 | #### 247 | statVar_element 248 | 12 249 | 4 250 | 1 251 | 1 252 | 1 253 | 1 254 | 1 255 | 1 256 | 1 257 | 1 258 | 1 259 | 1 260 | 1 261 | 1 262 | #### 263 | statVar_names 264 | 12 265 | 4 266 | basin_cfs 267 | basin_gwflow_cfs 268 | basin_sroff_cfs 269 | basin_dunnian 270 | basin_ssflow_cfs 271 | basin_ppt 272 | basin_snow 273 | basin_pweqv 274 | basin_snowcov 275 | basin_actet 276 | basin_potsw 277 | basin_recharge 278 | #### 279 | stat_var_file 280 | 1 281 | 4 282 | ../output/sagehen_statvar.dat 283 | #### 284 | statsON_OFF 285 | 1 286 | 1 287 | 1 288 | #### 289 | cascade_flag 290 | 1 291 | 1 292 | 1 293 | #### 294 | cascadegw_flag 295 | 1 296 | 1 297 | 1 298 | #### 299 | temp_module 300 | 1 301 | 4 302 | temp_1sta 303 | #### 304 | var_save_file 305 | 1 306 | 4 307 | ../output/prms_ic.out 308 | #### 309 | stats_output_file 310 | 1 311 | 4 312 | ../output/stats1.out 313 | #### 314 | var_init_file 315 | 1 316 | 4 317 | vars 318 | #### 319 | rpt_days 320 | 1 321 | 1 322 | 1 323 | #### 324 | capillary_module 325 | 1 326 | 4 327 | soilzone_prms -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.dbf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/NHDFlowline.dbf -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.prj: -------------------------------------------------------------------------------- 1 | GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/NHDFlowline.sbn -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/NHDFlowline.sbx -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/NHDFlowline.shp -------------------------------------------------------------------------------- /examples/sagehen/nhd/NHDFlowline.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/NHDFlowline.shx -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.dbf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/WBDHU12.dbf -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.prj: -------------------------------------------------------------------------------- 1 | GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/WBDHU12.sbn -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/WBDHU12.sbx -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/WBDHU12.shp -------------------------------------------------------------------------------- /examples/sagehen/nhd/WBDHU12.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/nhd/WBDHU12.shx -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_01_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_01_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_02_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_02_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_03_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_03_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_04_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_04_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_05_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_05_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_06_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_06_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_07_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_07_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_08_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_08_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_09_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_09_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_10_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_10_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_11_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_11_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_12_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/ppt/PRISM_ppt_30yr_normal_800mM2_12_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_01_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_01_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_02_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_02_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_03_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_03_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_04_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_04_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_05_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_05_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_06_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_06_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_07_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_07_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_08_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_08_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_09_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_09_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_10_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_10_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_11_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_11_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_12_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmax/PRISM_tmax_30yr_normal_800mM2_12_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_01_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_01_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_02_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_02_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_03_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_03_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_04_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_04_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_05_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_05_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_06_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_06_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_07_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_07_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_08_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_08_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_09_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_09_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_10_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_10_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_11_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_11_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_12_bil.bil: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/prism/tmin/PRISM_tmin_30yr_normal_800mM2_12_bil.bil -------------------------------------------------------------------------------- /examples/sagehen/readme.txt: -------------------------------------------------------------------------------- 1 | Gsflow-Arcpy example model readme 2 | 3 | Included with the Sagehen example problem: 4 | 5 | DEM - National elevation dataset 10m in NAD83 projection 6 | Boundary shapefile - Sagehen watershed polygon shapefile - watershed.shp 7 | Outlet point shapefile - Sagehen outlet point shapefile - model_points.shp 8 | Flowlines - NHD stream lines clipped to Sagehen model boundary 9 | CRT executable (CRT1.4_beta.exe) 10 | Vegetation data - LANDFIRE existing vegetation cover and type - us_140evc.img and us_140evt.img 11 | Soil data - Rasters of available water capacity, clay percent, sand percent, and saturated hydraulic conductivity - awc.img, clay.img, sand.img, ksat.img 12 | Climate data - PRISM 30 year monthly normals for precipitation and max/min temps 13 | Impervious dataset - National land coverage database, impervious cover 2011 - nlcd2011_imp.img 14 | Configuration file - Configuration file specific to Sagehen example model with all folders and fileneames 15 | Batch file to run all Gsflow-Arcpy scripts - runscripts.bat 16 | 17 | Upon running the batch file, the Sagehen example model will run through Gsflow-Arcpy and output parameter files in the hru_params folder 18 | A control file is provided in the prms folder, and a data file is provided in the inputs folder, where the parameter files will need to be transferred to run GSFLOW 19 | Control file - sagehen.control 20 | Data file - sagehen_datafile_1sta.data 21 | Output data will be stored in the outputs folder, and users can compare their results to the original outputs found in the output_compare folder 22 | GSFLOW executable (gsflow.exe) 23 | 24 | -------------------------------------------------------------------------------- /examples/sagehen/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\scripts\fishnet_generator.py -i sagehen_parameters.ini --overwrite 2 | 3 | python ..\..\scripts\hru_parameters.py -i sagehen_parameters.ini 4 | python ..\..\scripts\dem_parameters.py -i sagehen_parameters.ini 5 | python ..\..\scripts\dem_2_streams.py -i sagehen_parameters.ini 6 | python ..\..\scripts\crt_fill_parameters.py -i sagehen_parameters.ini 7 | python ..\..\scripts\dem_2_streams.py -i sagehen_parameters.ini 8 | python ..\..\scripts\crt_fill_parameters.py -i sagehen_parameters.ini 9 | python ..\..\scripts\stream_parameters.py -i sagehen_parameters.ini 10 | python ..\..\scripts\veg_parameters.py -i sagehen_parameters.ini 11 | python ..\..\scripts\soil_raster_prep.py -i sagehen_parameters.ini 12 | python ..\..\scripts\soil_parameters.py -i sagehen_parameters.ini 13 | python ..\..\scripts\prism_800m_normals.py -i sagehen_parameters.ini 14 | python ..\..\scripts\ppt_ratio_parameters.py -i sagehen_parameters.ini 15 | python ..\..\scripts\impervious_parameters.py -i sagehen_parameters.ini 16 | python ..\..\scripts\prms_template_fill.py -i sagehen_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen/runscripts_python_removed.bat: -------------------------------------------------------------------------------- 1 | ..\..\scripts\fishnet_generator.py -i sagehen_parameters.ini 2 | 3 | ..\..\scripts\hru_parameters.py -i sagehen_parameters.ini 4 | ..\..\scripts\dem_parameters.py -i sagehen_parameters.ini 5 | ..\..\scripts\dem_2_streams.py -i sagehen_parameters.ini 6 | ..\..\scripts\crt_fill_parameters.py -i sagehen_parameters.ini 7 | ..\..\scripts\dem_2_streams.py -i sagehen_parameters.ini 8 | ..\..\scripts\crt_fill_parameters.py -i sagehen_parameters.ini 9 | ..\..\scripts\stream_parameters.py -i sagehen_parameters.ini 10 | ..\..\scripts\veg_parameters.py -i sagehen_parameters.ini 11 | ..\..\scripts\soil_raster_prep.py -i sagehen_parameters.ini 12 | ..\..\scripts\soil_parameters.py -i sagehen_parameters.ini 13 | ..\..\scripts\prism_800m_normals.py -i sagehen_parameters.ini 14 | ..\..\scripts\ppt_ratio_parameters.py -i sagehen_parameters.ini 15 | ..\..\scripts\impervious_parameters.py -i sagehen_parameters.ini 16 | ..\..\scripts\prms_template_fill.py -i sagehen_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen/sagehen_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned10m_nad83.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = BILINEAR 48 | dem_cellsize = 10 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = True 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_FLOWAC 54 | reset_dem_adj_flag = False 55 | # Round DEM_ADJ to the given number of decimals 56 | dem_adj_decimals = 2 57 | 58 | 59 | ## DEM 2 Streams 60 | # Cells with flow accumulations >= threshold will be designated as stream cells 61 | flow_acc_threshold = 100 62 | # All 1st order streams with a length below threshold will be removed 63 | flow_length_threshold = 3 64 | 65 | calc_flow_dir_points_flag = True 66 | 67 | 68 | ## CRT Parameters 69 | crt_exe_path = ..\..\crt\CRT_1.3.1.exe 70 | crt_hruflg = 0 71 | crt_flowflg = 3 72 | crt_dpit = 0.01 73 | crt_outitmax = 100000 74 | 75 | 76 | ## CRT Fill Parameters 77 | use_crt_fill_flag = True 78 | 79 | 80 | ## Vegetation Parameters 81 | # Vegetation Type 82 | # Assume NEAREST Resampling 83 | veg_type_orig_path = .\veg\us_140evt.img 84 | veg_type_field = VALUE 85 | veg_type_cellsize = 10 86 | 87 | # Vegetation Cover 88 | # Assume NEAREST Resampling 89 | veg_cover_orig_path = .\veg\us_140evc.img 90 | veg_cover_cellsize = 10 91 | 92 | 93 | ## Soils Parameters 94 | soil_orig_folder = .\soils 95 | ksat_name = ksat.img 96 | awc_name = awc.img 97 | clay_pct_name = clay.img 98 | sand_pct_name = sand.img 99 | soil_cellsize = 50 100 | # Soil rasters are percent (i.e. 25%) 101 | # If false interpret as decimals (i.e. 0.25) 102 | soil_pct_flag = True 103 | # Fill nodata in soil rasters using nibble technique 104 | fill_soil_nodata_flag = True 105 | # Initial soil moisture and recharge are calculated as a fraction of the max 106 | moist_init_ratio = 0.1 107 | rechr_init_ratio = 0.1 108 | # Geology based multiplier raster for ssr2gw rate 109 | ssr2gw_mult_flag = False 110 | ssr2gw_mult_name = ssr2g_mult.img 111 | # If no raster provided, a default value must be set 112 | ssr2gw_k_default = 0.001 113 | # Read and apply soil depth raster 114 | # Otherwise soil depth will only be derived from rooting depth 115 | soil_depth_flag = False 116 | # soil_depth_name = soil_depth.img 117 | 118 | 119 | ## Impervious Cover Parameters 120 | impervious_orig_path = .\impervious\nlcd2011_imp.img 121 | # Resampling method: BILINEAR, CUBIC, NEAREST 122 | impervious_projection_method = BILINEAR 123 | impervious_cellsize = 10 124 | # Impervious cover rasters are percent (i.e. 25%) 125 | # If false interpret as decimals (i.e. 0.25) 126 | impervious_pct_flag = True 127 | 128 | 129 | ## Remap Files 130 | remap_folder = ..\..\remaps 131 | aspect_remap = aspect.rmp 132 | # Output values are 10 * value 133 | temp_adj_remap = temp_adj_x10.rmp 134 | # Output values are floats 135 | cov_type_remap = covtype.rmp 136 | covden_sum_remap = covdensum.rmp 137 | covden_win_remap = covdenwin.rmp 138 | snow_intcp_remap = covtype_to_snow_intcp.rmp 139 | srain_intcp_remap = covtype_to_srain_intcp.rmp 140 | wrain_intcp_remap = covtype_to_wrain_intcp.rmp 141 | root_depth_remap = rtdepth.rmp 142 | 143 | 144 | ## PRISM Parameters 145 | prism_folder = .\prism 146 | # Resampling method: BILINEAR, CUBIC, NEAREST 147 | prism_projection_method = BILINEAR 148 | # Output projected cellsize, not PRISM input cellsize 149 | prism_cellsize = 90 150 | # Recalculate JH coefficient with PRISM temperature values 151 | calc_prism_jh_coef_flag = True 152 | 153 | # PPT Ratios 154 | set_ppt_zones_flag = False 155 | ppt_obs_units = mm 156 | 157 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 158 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 159 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 160 | ppt_hru_id = 0 161 | 162 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 163 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 164 | ppt_zone_path = .\shapefiles\ppt_zones.shp 165 | ppt_zone_id_field = PPT_ZONE 166 | ppt_hru_id_field = PPT_HRU_ID 167 | ppt_obs_field_format = PPT_{:02d} 168 | 169 | 170 | ## PRMS Parameter 171 | prms_parameter_folder = .\model\input\prms 172 | 173 | single_param_file_flag = False 174 | # single_param_file_name = prms_inputs.param 175 | param_column_flag = True 176 | 177 | # Default/template values 178 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 179 | prms_param_csv_path = .\hru_params\prms_parameters.csv 180 | 181 | # Parameters passed through to PRMS param file(s) 182 | # Number of air temperature measurement stations in the data file 183 | ntemp = 1 184 | # Elevation units (0=feet, 1=meter) 185 | elev_units = 0 186 | # Index of the temperature station used to compute basin temperature values (1's based) 187 | basin_tsta = 1 188 | # Index of the base temperature station used for lapse rate calculations (1's based) 189 | hru_tsta = 1 190 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 191 | hru_tlaps = 1 192 | # Elevation of each air temperature measurment station 193 | tsta_elev = 6340 194 | 195 | # If using temp_1sta module, manually define monthly lapse rates 196 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 197 | tmin_lapse = 4.5, 4.5, 3.5, 3.5, 3, 2.5, 2.5, 2, 3, 3.5, 4, 4.5 198 | 199 | 200 | ## Miscellaneous 201 | # Scratch Workspace ('in_memory' or 'scratch') 202 | # in_memory may not work for computers with limited RAM 203 | # scratch will write all intermediate files to a scratch folder on the disk 204 | scratch_name = in_memory 205 | # scratch_name = scratch 206 | 207 | # Scale floating point values before converting to Int and calculating Median 208 | int_factor = 1 209 | 210 | # Don't change this unless it conflicts with an existing field 211 | orig_fid_field = ORIG_FID 212 | -------------------------------------------------------------------------------- /examples/sagehen/save/sagehen.control.save: -------------------------------------------------------------------------------- 1 | PRMS Control File for Sagehen 2 | #### 3 | mapOutVar_names 4 | 1 5 | 4 6 | ./output/recharge 7 | #### 8 | nmapOutVars 9 | 1 10 | 1 11 | 0 12 | #### 13 | mapOutON_OFF 14 | 1 15 | 1 16 | 1 17 | #### 18 | parameter_check_flag 19 | 1 20 | 1 21 | 0 22 | #### 23 | prms_warmup 24 | 1 25 | 1 26 | 3 27 | #### 28 | aniOutON_OFF 29 | 1 30 | 1 31 | 0 32 | #### 33 | print_debug 34 | 1 35 | 1 36 | -1 37 | #### 38 | aniOutVar_names 39 | 5 40 | 4 41 | upslope_interflow 42 | upslope_dunnianflow 43 | soil_moist_tot 44 | upslope_hortonian 45 | sroff 46 | #### 47 | ani_output_file 48 | 1 49 | 4 50 | animation.out 51 | #### 52 | csv_output_file 53 | 1 54 | 4 55 | gsflow.csv 56 | #### 57 | model_output_file 58 | 1 59 | 4 60 | ../output/gsflow_prms.out 61 | #### 62 | data_file 63 | 1 64 | 4 65 | ../input/prms/sagehen_datafile_1sta.data 66 | #### 67 | dispGraphsBuffSize 68 | 1 69 | 4 70 | 15 71 | #### 72 | dispVar_element 73 | 15 74 | 4 75 | 2 76 | 2 77 | 2 78 | 2 79 | 2 80 | 1 81 | 2 82 | 1 83 | 1 84 | 1 85 | 2 86 | 1 87 | 1 88 | 1 89 | 1 90 | #### 91 | dispVar_names 92 | 15 93 | 4 94 | runoff 95 | sub_cfs 96 | sub_gwflow 97 | sub_interflow 98 | sub_sroff 99 | basin_pweqv 100 | subinc_pkweqv 101 | basin_dunnian 102 | basin_prefflow 103 | basin_slowflow 104 | subinc_precip 105 | basin_rain 106 | basin_snow 107 | basin_tmax 108 | basin_tmin 109 | #### 110 | dispVar_plot 111 | 15 112 | 1 113 | 1 114 | 1 115 | 2 116 | 2 117 | 2 118 | 3 119 | 3 120 | 4 121 | 4 122 | 4 123 | 5 124 | 5 125 | 5 126 | 6 127 | 6 128 | #### 129 | end_time 130 | 6 131 | 1 132 | 1997 133 | 3 134 | 31 135 | 0 136 | 0 137 | 0 138 | #### 139 | et_module 140 | 1 141 | 4 142 | potet_jh 143 | #### 144 | executable_desc 145 | 1 146 | 4 147 | GSFLOW model 148 | #### 149 | executable_model 150 | 1 151 | 4 152 | prmsIV.exe 153 | #### 154 | gsflow_output_file 155 | 1 156 | 4 157 | ../output/gsflow.out 158 | #### 159 | init_vars_from_file 160 | 1 161 | 1 162 | 0 163 | #### 164 | initial_deltat 165 | 1 166 | 2 167 | 24.0 168 | #### 169 | mms_user_dir 170 | 1 171 | 4 172 | ./ 173 | #### 174 | mms_user_out_dir 175 | 1 176 | 4 177 | ./ 178 | #### 179 | model_mode 180 | 1 181 | 4 182 | PRMS 183 | #### 184 | modflow_name 185 | 1 186 | 4 187 | sagehen.nam 188 | #### 189 | naniOutVars 190 | 1 191 | 1 192 | 5 193 | #### 194 | ndispGraphs 195 | 1 196 | 1 197 | 6 198 | #### 199 | param_file 200 | 4 201 | 4 202 | ../input/prms/calibration_parameters.param 203 | ../input/prms/cascade.param 204 | ../input/prms/default_values.param 205 | ../input/prms/gis_derived_parameters.param 206 | #### 207 | precip_module 208 | 1 209 | 4 210 | precip_1sta 211 | #### 212 | param_print_file 213 | 1 214 | 4 215 | ../output/sagehen.parprt 216 | #### 217 | save_vars_to_file 218 | 1 219 | 1 220 | 1 221 | #### 222 | solrad_module 223 | 1 224 | 4 225 | ddsolrad 226 | #### 227 | srunoff_module 228 | 1 229 | 4 230 | srunoff_smidx 231 | #### 232 | start_time 233 | 6 234 | 1 235 | 1982 236 | 8 237 | 1 238 | 0 239 | 0 240 | 0 241 | #### 242 | nstatVars 243 | 1 244 | 1 245 | 12 246 | #### 247 | statVar_element 248 | 12 249 | 4 250 | 1 251 | 1 252 | 1 253 | 1 254 | 1 255 | 1 256 | 1 257 | 1 258 | 1 259 | 1 260 | 1 261 | 1 262 | #### 263 | statVar_names 264 | 12 265 | 4 266 | basin_cfs 267 | basin_gwflow_cfs 268 | basin_sroff_cfs 269 | basin_dunnian 270 | basin_ssflow_cfs 271 | basin_ppt 272 | basin_snow 273 | basin_pweqv 274 | basin_snowcov 275 | basin_actet 276 | basin_potsw 277 | basin_recharge 278 | #### 279 | stat_var_file 280 | 1 281 | 4 282 | ../output/sagehen_statvar.dat 283 | #### 284 | statsON_OFF 285 | 1 286 | 1 287 | 1 288 | #### 289 | cascade_flag 290 | 1 291 | 1 292 | 1 293 | #### 294 | cascadegw_flag 295 | 1 296 | 1 297 | 1 298 | #### 299 | temp_module 300 | 1 301 | 4 302 | temp_1sta 303 | #### 304 | var_save_file 305 | 1 306 | 4 307 | ../output/prms_ic.out 308 | #### 309 | stats_output_file 310 | 1 311 | 4 312 | ../output/stats1.out 313 | #### 314 | var_init_file 315 | 1 316 | 4 317 | vars 318 | #### 319 | rpt_days 320 | 1 321 | 1 322 | 1 323 | #### 324 | capillary_module 325 | 1 326 | 4 327 | soilzone_prms -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | vA3TYPEC2 OUTLET  -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/watershed.dbf -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen/soils/awc.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/soils/awc.img -------------------------------------------------------------------------------- /examples/sagehen/soils/clay.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/soils/clay.img -------------------------------------------------------------------------------- /examples/sagehen/soils/ksat.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/soils/ksat.img -------------------------------------------------------------------------------- /examples/sagehen/soils/sand.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/soils/sand.img -------------------------------------------------------------------------------- /examples/sagehen/veg/us_140evc.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/veg/us_140evc.img -------------------------------------------------------------------------------- /examples/sagehen/veg/us_140evc.img.vat.dbf: -------------------------------------------------------------------------------- 1 | uaValueN 2 | CountF 11 3.13300000000e+003 14 2.00000000000e+000 31 4.70000000000e+001 100 6.85000000000e+002 101 5.93000000000e+002 102 3.27300000000e+003 103 1.87850000000e+004 104 3.19540000000e+004 105 2.11240000000e+004 106 8.76800000000e+003 107 2.33000000000e+002 108 9.00000000000e+000 111 2.96000000000e+002 112 4.70000000000e+001 113 3.60000000000e+001 114 1.91800000000e+003 115 1.48100000000e+003 116 7.71000000000e+002 117 5.20000000000e+001 118 8.00000000000e+000 121 1.00000000000e+000 122 2.61000000000e+002 123 2.00000000000e+000 124 4.00000000000e+000 125 1.20000000000e+001 126 1.20000000000e+001 127 7.00000000000e+000 128 1.00000000000e+000 -------------------------------------------------------------------------------- /examples/sagehen/veg/us_140evt.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen/veg/us_140evt.img -------------------------------------------------------------------------------- /examples/sagehen/veg/us_140evt.img.vat.dbf: -------------------------------------------------------------------------------- 1 | u!aValueN 2 | CountF 3001 1.00000000000e+000 3002 5.00000000000e+001 3011 1.80000000000e+001 3017 6.00000000000e+000 3027 5.40000000000e+001 3028 1.56490000000e+004 3029 2.10000000000e+001 3031 3.58280000000e+004 3032 2.92330000000e+004 3033 3.90000000000e+001 3034 5.50000000000e+001 3044 2.60000000000e+001 3058 5.10000000000e+001 3071 9.00000000000e+000 3080 7.10000000000e+001 3098 3.52500000000e+003 3103 2.00000000000e+000 3126 2.65000000000e+002 3135 6.00000000000e+000 3138 2.91000000000e+002 3152 2.69000000000e+003 3154 1.06000000000e+003 3181 1.00000000000e+000 3219 4.00000000000e+000 3220 6.50000000000e+001 3221 6.30000000000e+002 3231 9.00000000000e+000 3255 6.42000000000e+002 3292 3.13300000000e+003 3294 4.70000000000e+001 3901 2.00000000000e+000 3923 3.00000000000e+001 3924 2.00000000000e+000 -------------------------------------------------------------------------------- /examples/sagehen_tests/README.md: -------------------------------------------------------------------------------- 1 | # Tests Models 2 | 3 | There are currently seven test models that have been developed from the Sagehen model to test different aspects of the scripts. 4 | 5 | ## active_fill 6 | 7 | Define the study area to include cells that are well outside the basin. 8 | These cells will be made active and filled so that they flow to the outlet point. 9 | 10 | ## clipped_elev 11 | 12 | Clip the input elevation grid so that all inactive cells have no elevation data. 13 | Check that the model still correctly identifies the outlet point(s). 14 | 15 | ## lake 16 | 17 | Add a single lake to the model. The elevations in the lake have been manually adjusted to the same value in the provided DEM. 18 | 19 | ## precipitation_zones 20 | 21 | ## temperatures_zones 22 | 23 | ## terminal_basin 24 | The DEM was modified to "dam" the eastern portion of the watershed and the outlet point was changed to a "SWALE" type to simulate a closed/terminal basin. 25 | This test was added because the scripts were not running when the model did not include an "OUTLET" point. 26 | 27 | ## terminal_lake 28 | This is the same DEM as the other terminal basin example but a lake was included around the "SWALE" point. 29 | This test was added to ensure that a "SWALE" point in a lake works the same as one outside. -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/dem/ned90m_nad83z11.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/dem/ned90m_nad83z11.img -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | va=TYPEC2ZONE_VALUEN 2 | OUTLET 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | vA HRU_TYPEN 2 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/active_fill/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/active_fill/test_parameters_fill.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned90m_nad83z11.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = NEAREST 48 | dem_cellsize = 90 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = False 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_MEAN 54 | reset_dem_adj_flag = False 55 | 56 | 57 | ## DEM 2 Streams 58 | # Cells with flow accumulations >= threshold will be designated as stream cells 59 | flow_acc_threshold = 100 60 | # All 1st order streams with a length below threshold will be removed 61 | flow_length_threshold = 3 62 | 63 | calc_flow_dir_points_flag = True 64 | 65 | 66 | ## CRT Parameters 67 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 68 | crt_hruflg = 0 69 | crt_flowflg = 3 70 | crt_dpit = 0.01 71 | crt_outitmax = 100000 72 | 73 | 74 | ## CRT Fill Parameters 75 | use_crt_fill_flag = True 76 | 77 | 78 | ## Vegetation Parameters 79 | # Vegetation Type 80 | # Assume NEAREST Resampling 81 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 82 | veg_type_field = VALUE 83 | veg_type_cellsize = 10 84 | 85 | # Vegetation Cover 86 | # Assume NEAREST Resampling 87 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 88 | veg_cover_cellsize = 10 89 | 90 | 91 | ## Soils Parameters 92 | soil_orig_folder = ..\..\sagehen\soils 93 | ksat_name = ksat.img 94 | awc_name = awc.img 95 | clay_pct_name = clay.img 96 | sand_pct_name = sand.img 97 | soil_cellsize = 50 98 | # Soil rasters are percent (i.e. 25%) 99 | # If false interpret as decimals (i.e. 0.25) 100 | soil_pct_flag = True 101 | # Fill nodata in soil rasters using nibble technique 102 | fill_soil_nodata_flag = True 103 | # Initial soil moisture and recharge are calculated as a fraction of the max 104 | moist_init_ratio = 0.1 105 | rechr_init_ratio = 0.1 106 | ssr2gw_k_default = 0.001 107 | # Setting true will clip the root depth to the soil depth 108 | clip_root_depth_flag = False 109 | # Soil depth will only be read if clip_root_depth_flag = True 110 | # Soil depth must be in inches (?) 111 | soil_depth_name = avg_depth 112 | 113 | 114 | ## Impervious Cover Parameters 115 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 116 | # Resampling method: BILINEAR, CUBIC, NEAREST 117 | impervious_projection_method = BILINEAR 118 | impervious_cellsize = 10 119 | # Impervious cover rasters are percent (i.e. 25%) 120 | # If false interpret as decimals (i.e. 0.25) 121 | impervious_pct_flag = True 122 | 123 | 124 | ## Remap Files 125 | remap_folder = ..\..\..\remaps 126 | aspect_remap = aspect.rmp 127 | # Output values are 10 * value 128 | temp_adj_remap = temp_adj_x10.rmp 129 | # Output values are floats 130 | cov_type_remap = covtype.rmp 131 | covden_sum_remap = covdensum.rmp 132 | covden_win_remap = covdenwin.rmp 133 | snow_intcp_remap = snow_intcp.rmp 134 | srain_intcp_remap = srain_intcp.rmp 135 | wrain_intcp_remap = wrain_intcp.rmp 136 | root_depth_remap = rtdepth.rmp 137 | 138 | 139 | ## PRISM Parameters 140 | prism_folder = ..\..\sagehen\prism 141 | # Resampling method: BILINEAR, CUBIC, NEAREST 142 | prism_projection_method = BILINEAR 143 | # Output projected cellsize, not PRISM input cellsize 144 | prism_cellsize = 90 145 | # Recalculate JH coefficient with PRISM temperature values 146 | calc_prism_jh_coef_flag = True 147 | 148 | # PPT Ratios 149 | set_ppt_zones_flag = False 150 | ppt_obs_units = mm 151 | 152 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 153 | # NOTE - These are made up values! 154 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 155 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 156 | ppt_hru_id = 0 157 | 158 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 159 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 160 | ppt_zone_path = .\shapefiles\ppt_zones.shp 161 | ppt_zone_id_field = PPT_ZONE 162 | ppt_hru_id_field = PPT_HRU_ID 163 | ppt_obs_field_format = PPT_{:02d} 164 | 165 | 166 | ## PRMS Parameter 167 | prms_parameter_folder = .\hru_params 168 | 169 | single_param_file_flag = False 170 | # single_param_file_name = prms_inputs.param 171 | param_column_flag = True 172 | 173 | # Default/template values 174 | prms_dimen_csv_path = ..\..\sagehen\hru_params\prms_dimensions.csv 175 | prms_param_csv_path = ..\..\sagehen\hru_params\prms_parameters.csv 176 | 177 | # Parameters passed through to PRMS param file(s) 178 | # Number of air temperature measurement stations in the data file 179 | ntemp = 1 180 | # Elevation units (0=feet, 1=meter) 181 | elev_units = 0 182 | # Index of the temperature station used to compute basin temperature values (1's based) 183 | basin_tsta = 1 184 | # Index of the base temperature station used for lapse rate calculations (1's based) 185 | hru_tsta = 1 186 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 187 | hru_tlaps = 1 188 | # Elevation of each air temperature measurment station 189 | tsta_elev = 6340 190 | 191 | # If using temp_1sta module, manually define monthly lapse rates 192 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 193 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 194 | 195 | 196 | ## Miscellaneous 197 | # Scratch Workspace ('in_memory' or 'scratch') 198 | # in_memory may not work for computers with limited RAM 199 | # scratch will write all intermediate files to a scratch folder on the disk 200 | scratch_name = in_memory 201 | # scratch_name = scratch 202 | 203 | # Scale floating point values before converting to Int and calculating Median 204 | int_factor = 1 205 | 206 | # Don't change this unless it conflicts with an existing field 207 | orig_fid_field = ORIG_FID 208 | -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/dem/ned90m_nad83z11_clip.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/dem/ned90m_nad83z11_clip.img -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | va=TYPEC2ZONE_VALUEN 2 | OUTLET 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | vA HRU_TYPEN 2 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/clipped_elev/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/clipped_elev/test_parameters_elev_clip.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned90m_nad83z11_clip.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = NEAREST 48 | dem_cellsize = 90 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = False 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_MEAN 54 | reset_dem_adj_flag = False 55 | 56 | 57 | ## DEM 2 Streams 58 | # Cells with flow accumulations >= threshold will be designated as stream cells 59 | flow_acc_threshold = 100 60 | # All 1st order streams with a length below threshold will be removed 61 | flow_length_threshold = 3 62 | 63 | calc_flow_dir_points_flag = True 64 | 65 | 66 | ## CRT Parameters 67 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 68 | crt_hruflg = 0 69 | crt_flowflg = 3 70 | crt_dpit = 0.01 71 | crt_outitmax = 100000 72 | 73 | 74 | ## CRT Fill Parameters 75 | use_crt_fill_flag = True 76 | 77 | 78 | ## Vegetation Parameters 79 | # Vegetation Type 80 | # Assume NEAREST Resampling 81 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 82 | veg_type_field = VALUE 83 | veg_type_cellsize = 10 84 | 85 | # Vegetation Cover 86 | # Assume NEAREST Resampling 87 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 88 | veg_cover_cellsize = 10 89 | 90 | 91 | ## Soils Parameters 92 | soil_orig_folder = ..\..\sagehen\soils 93 | ksat_name = ksat.img 94 | awc_name = awc.img 95 | clay_pct_name = clay.img 96 | sand_pct_name = sand.img 97 | soil_cellsize = 50 98 | # Soil rasters are percent (i.e. 25%) 99 | # If false interpret as decimals (i.e. 0.25) 100 | soil_pct_flag = True 101 | # Fill nodata in soil rasters using nibble technique 102 | fill_soil_nodata_flag = True 103 | # Initial soil moisture and recharge are calculated as a fraction of the max 104 | moist_init_ratio = 0.1 105 | rechr_init_ratio = 0.1 106 | ssr2gw_k_default = 0.001 107 | # Setting true will clip the root depth to the soil depth 108 | clip_root_depth_flag = False 109 | # Soil depth will only be read if clip_root_depth_flag = True 110 | # Soil depth must be in inches (?) 111 | soil_depth_name = avg_depth 112 | 113 | 114 | ## Impervious Cover Parameters 115 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 116 | # Resampling method: BILINEAR, CUBIC, NEAREST 117 | impervious_projection_method = BILINEAR 118 | impervious_cellsize = 10 119 | # Impervious cover rasters are percent (i.e. 25%) 120 | # If false interpret as decimals (i.e. 0.25) 121 | impervious_pct_flag = True 122 | 123 | 124 | ## Remap Files 125 | remap_folder = ..\..\..\remaps 126 | aspect_remap = aspect.rmp 127 | # Output values are 10 * value 128 | temp_adj_remap = temp_adj_x10.rmp 129 | # Output values are floats 130 | cov_type_remap = covtype.rmp 131 | covden_sum_remap = covdensum.rmp 132 | covden_win_remap = covdenwin.rmp 133 | snow_intcp_remap = snow_intcp.rmp 134 | srain_intcp_remap = srain_intcp.rmp 135 | wrain_intcp_remap = wrain_intcp.rmp 136 | root_depth_remap = rtdepth.rmp 137 | 138 | 139 | ## PRISM Parameters 140 | prism_folder = ..\..\sagehen\prism 141 | # Resampling method: BILINEAR, CUBIC, NEAREST 142 | prism_projection_method = BILINEAR 143 | # Output projected cellsize, not PRISM input cellsize 144 | prism_cellsize = 90 145 | # Recalculate JH coefficient with PRISM temperature values 146 | calc_prism_jh_coef_flag = True 147 | 148 | # PPT Ratios 149 | set_ppt_zones_flag = False 150 | ppt_obs_units = mm 151 | 152 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 153 | # NOTE - These are made up values! 154 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 155 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 156 | ppt_hru_id = 0 157 | 158 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 159 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 160 | ppt_zone_path = .\shapefiles\ppt_zones.shp 161 | ppt_zone_id_field = PPT_ZONE 162 | ppt_hru_id_field = PPT_HRU_ID 163 | ppt_obs_field_format = PPT_{:02d} 164 | 165 | 166 | ## PRMS Parameter 167 | prms_parameter_folder = .\hru_params 168 | 169 | single_param_file_flag = False 170 | # single_param_file_name = prms_inputs.param 171 | param_column_flag = True 172 | 173 | # Default/template values 174 | prms_dimen_csv_path = ..\..\sagehen\hru_params\prms_dimensions.csv 175 | prms_param_csv_path = ..\..\sagehen\hru_params\prms_parameters.csv 176 | 177 | # Parameters passed through to PRMS param file(s) 178 | # Number of air temperature measurement stations in the data file 179 | ntemp = 1 180 | # Elevation units (0=feet, 1=meter) 181 | elev_units = 0 182 | # Index of the temperature station used to compute basin temperature values (1's based) 183 | basin_tsta = 1 184 | # Index of the base temperature station used for lapse rate calculations (1's based) 185 | hru_tsta = 1 186 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 187 | hru_tlaps = 1 188 | # Elevation of each air temperature measurment station 189 | tsta_elev = 6340 190 | 191 | # If using temp_1sta module, manually define monthly lapse rates 192 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 193 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 194 | 195 | 196 | ## Miscellaneous 197 | # Scratch Workspace ('in_memory' or 'scratch') 198 | # in_memory may not work for computers with limited RAM 199 | # scratch will write all intermediate files to a scratch folder on the disk 200 | scratch_name = in_memory 201 | # scratch_name = scratch 202 | 203 | # Scale floating point values before converting to Int and calculating Median 204 | int_factor = 1 205 | 206 | # Don't change this unless it conflicts with an existing field 207 | orig_fid_field = ORIG_FID 208 | -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/dem/ned90m_nad83z11.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/dem/ned90m_nad83z11.img -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/lake_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = True 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned90m_nad83z11.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = NEAREST 48 | dem_cellsize = 90 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = False 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_MEAN 54 | reset_dem_adj_flag = False 55 | 56 | 57 | ## DEM 2 Streams 58 | # Cells with flow accumulations >= threshold will be designated as stream cells 59 | flow_acc_threshold = 100 60 | # All 1st order streams with a length below threshold will be removed 61 | flow_length_threshold = 3 62 | 63 | calc_flow_dir_points_flag = True 64 | 65 | 66 | ## CRT Parameters 67 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 68 | crt_hruflg = 0 69 | crt_flowflg = 3 70 | crt_dpit = 0.01 71 | crt_outitmax = 100000 72 | 73 | 74 | ## CRT Fill Parameters 75 | use_crt_fill_flag = True 76 | 77 | 78 | ## Vegetation Parameters 79 | # Vegetation Type 80 | # Assume NEAREST Resampling 81 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 82 | veg_type_field = VALUE 83 | veg_type_cellsize = 10 84 | 85 | # Vegetation Cover 86 | # Assume NEAREST Resampling 87 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 88 | veg_cover_cellsize = 10 89 | 90 | 91 | ## Soils Parameters 92 | soil_orig_folder = ..\..\sagehen\soils 93 | ksat_name = ksat.img 94 | awc_name = awc.img 95 | clay_pct_name = clay.img 96 | sand_pct_name = sand.img 97 | soil_cellsize = 50 98 | # Soil rasters are percent (i.e. 25%) 99 | # If false interpret as decimals (i.e. 0.25) 100 | soil_pct_flag = True 101 | # Fill nodata in soil rasters using nibble technique 102 | fill_soil_nodata_flag = True 103 | # Initial soil moisture and recharge are calculated as a fraction of the max 104 | moist_init_ratio = 0.1 105 | rechr_init_ratio = 0.1 106 | ssr2gw_k_default = 0.001 107 | # Setting true will clip the root depth to the soil depth 108 | clip_root_depth_flag = False 109 | # Soil depth will only be read if clip_root_depth_flag = True 110 | # Soil depth must be in inches (?) 111 | soil_depth_name = avg_depth 112 | 113 | 114 | ## Impervious Cover Parameters 115 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 116 | # Resampling method: BILINEAR, CUBIC, NEAREST 117 | impervious_projection_method = BILINEAR 118 | impervious_cellsize = 10 119 | # Impervious cover rasters are percent (i.e. 25%) 120 | # If false interpret as decimals (i.e. 0.25) 121 | impervious_pct_flag = True 122 | 123 | 124 | ## Remap Files 125 | remap_folder = ..\..\..\remaps 126 | aspect_remap = aspect.rmp 127 | # Output values are 10 * value 128 | temp_adj_remap = temp_adj_x10.rmp 129 | # Output values are floats 130 | cov_type_remap = covtype.rmp 131 | covden_sum_remap = covdensum.rmp 132 | covden_win_remap = covdenwin.rmp 133 | snow_intcp_remap = snow_intcp.rmp 134 | srain_intcp_remap = srain_intcp.rmp 135 | wrain_intcp_remap = wrain_intcp.rmp 136 | root_depth_remap = rtdepth.rmp 137 | 138 | 139 | ## PRISM Parameters 140 | prism_folder = ..\..\sagehen\prism 141 | # Resampling method: BILINEAR, CUBIC, NEAREST 142 | prism_projection_method = BILINEAR 143 | # Output projected cellsize, not PRISM input cellsize 144 | prism_cellsize = 90 145 | # Recalculate JH coefficient with PRISM temperature values 146 | calc_prism_jh_coef_flag = True 147 | 148 | # PPT Ratios 149 | set_ppt_zones_flag = False 150 | ppt_obs_units = mm 151 | 152 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 153 | # NOTE - These are made up values! 154 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 155 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 156 | ppt_hru_id = 0 157 | 158 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 159 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 160 | ppt_zone_path = .\shapefiles\ppt_zones.shp 161 | ppt_zone_id_field = PPT_ZONE 162 | ppt_hru_id_field = PPT_HRU_ID 163 | ppt_obs_field_format = PPT_{:02d} 164 | 165 | 166 | ## PRMS Parameter 167 | prms_parameter_folder = .\hru_params 168 | 169 | single_param_file_flag = False 170 | # single_param_file_name = prms_inputs.param 171 | param_column_flag = True 172 | 173 | # Default/template values 174 | prms_dimen_csv_path = ..\..\sagehen\hru_params\prms_dimensions.csv 175 | prms_param_csv_path = ..\..\sagehen\hru_params\prms_parameters.csv 176 | 177 | # Parameters passed through to PRMS param file(s) 178 | # Number of air temperature measurement stations in the data file 179 | ntemp = 1 180 | # Elevation units (0=feet, 1=meter) 181 | elev_units = 0 182 | # Index of the temperature station used to compute basin temperature values (1's based) 183 | basin_tsta = 1 184 | # Index of the base temperature station used for lapse rate calculations (1's based) 185 | hru_tsta = 1 186 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 187 | hru_tlaps = 1 188 | # Elevation of each air temperature measurment station 189 | tsta_elev = 6340 190 | 191 | # If using temp_1sta module, manually define monthly lapse rates 192 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 193 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 194 | 195 | 196 | ## Miscellaneous 197 | # Scratch Workspace ('in_memory' or 'scratch') 198 | # in_memory may not work for computers with limited RAM 199 | # scratch will write all intermediate files to a scratch folder on the disk 200 | scratch_name = in_memory 201 | # scratch_name = scratch 202 | 203 | # Scale floating point values before converting to Int and calculating Median 204 | int_factor = 1 205 | 206 | # Don't change this unless it conflicts with an existing field 207 | orig_fid_field = ORIG_FID 208 | -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\..\scripts\fishnet_generator.py -i lake_parameters.ini --overwrite 2 | 3 | python ..\..\..\scripts\hru_parameters.py -i lake_parameters.ini 4 | python ..\..\..\scripts\dem_parameters.py -i lake_parameters.ini 5 | python ..\..\..\scripts\dem_2_streams.py -i lake_parameters.ini 6 | python ..\..\..\scripts\crt_fill_parameters.py -i lake_parameters.ini 7 | python ..\..\..\scripts\dem_2_streams.py -i lake_parameters.ini 8 | python ..\..\..\scripts\crt_fill_parameters.py -i lake_parameters.ini 9 | python ..\..\..\scripts\stream_parameters.py -i lake_parameters.ini 10 | python ..\..\..\scripts\veg_parameters.py -i lake_parameters.ini 11 | python ..\..\..\scripts\soil_raster_prep.py -i lake_parameters.ini 12 | python ..\..\..\scripts\soil_parameters.py -i lake_parameters.ini 13 | python ..\..\..\scripts\prism_800m_normals.py -i lake_parameters.ini 14 | python ..\..\..\scripts\ppt_ratio_parameters.py -i lake_parameters.ini 15 | python ..\..\..\scripts\impervious_parameters.py -i lake_parameters.ini 16 | python ..\..\..\scripts\prms_template_fill.py -i lake_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.dbf: -------------------------------------------------------------------------------- 1 | vANAMEC Fake Lake  -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/lakes.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/lakes.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/lakes.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/lakes.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/lakes.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | va=TYPEC2ZONE_VALUEN 2 | OUTLET 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | vA HRU_TYPEN 2 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/lake/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/lake/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/dem/ned90m_nad83z11.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/dem/ned90m_nad83z11.img -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/hru_params/prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,1 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,0 18 | nsnow,calibration_parameters,0 19 | nsol,calibration_parameters,0 20 | nssr,calibration_parameters,calculated 21 | nsub,calibration_parameters,calculated 22 | ntemp,calibration_parameters,config_file 23 | one,calibration_parameters,1 24 | nlake_hrus,calibration_parameters,0 25 | -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/ppt_zones_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs - Sagehen test model - PPT zones 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | 38 | 39 | ## DEM Parameters 40 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 41 | dem_orig_path = .\dem\ned90m_nad83z11.img 42 | dem_units = meters 43 | # Resampling method: BILINEAR, CUBIC, NEAREST 44 | dem_projection_method = NEAREST 45 | dem_cellsize = 90 46 | # Calculate flow accumulation weighted elevation 47 | calc_flow_acc_dem_flag = False 48 | # flow_acc_dem_factor = 0.001 49 | # Field to initially set DEM_ADJ 50 | dem_adj_copy_field = DEM_MEAN 51 | reset_dem_adj_flag = False 52 | 53 | 54 | ## DEM 2 Streams 55 | # Cells with flow accumulations >= threshold will be designated as stream cells 56 | flow_acc_threshold = 100 57 | # All 1st order streams with a length below threshold will be removed 58 | flow_length_threshold = 3 59 | 60 | calc_flow_dir_points_flag = True 61 | 62 | 63 | ## CRT Parameters 64 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 65 | crt_hruflg = 0 66 | crt_flowflg = 3 67 | crt_dpit = 0.01 68 | crt_outitmax = 100000 69 | 70 | 71 | ## CRT Fill Parameters 72 | use_crt_fill_flag = True 73 | 74 | 75 | ## Vegetation Parameters 76 | # Vegetation Type 77 | # Assume NEAREST Resampling 78 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 79 | veg_type_field = VALUE 80 | veg_type_cellsize = 10 81 | 82 | # Vegetation Cover 83 | # Assume NEAREST Resampling 84 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 85 | veg_cover_cellsize = 10 86 | 87 | 88 | ## Soils Parameters 89 | soil_orig_folder = ..\..\sagehen\soils 90 | ksat_name = ksat.img 91 | awc_name = awc.img 92 | clay_pct_name = clay.img 93 | sand_pct_name = sand.img 94 | soil_cellsize = 50 95 | # Soil rasters are percent (i.e. 25%) 96 | # If false interpret as decimals (i.e. 0.25) 97 | soil_pct_flag = True 98 | # Fill nodata in soil rasters using nibble technique 99 | fill_soil_nodata_flag = True 100 | # Initial soil moisture and recharge are calculated as a fraction of the max 101 | moist_init_ratio = 0.1 102 | rechr_init_ratio = 0.1 103 | ssr2gw_k_default = 0.001 104 | # Setting true will clip the root depth to the soil depth 105 | clip_root_depth_flag = False 106 | # Soil depth will only be read if clip_root_depth_flag = True 107 | # Soil depth must be in inches (?) 108 | soil_depth_name = avg_depth 109 | 110 | 111 | ## Impervious Cover Parameters 112 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 113 | # Resampling method: BILINEAR, CUBIC, NEAREST 114 | impervious_projection_method = BILINEAR 115 | impervious_cellsize = 10 116 | # Impervious cover rasters are percent (i.e. 25%) 117 | # If false interpret as decimals (i.e. 0.25) 118 | impervious_pct_flag = True 119 | 120 | 121 | ## Remap Files 122 | remap_folder = ..\..\..\remaps 123 | aspect_remap = aspect.rmp 124 | # Output values are 10 * value 125 | temp_adj_remap = temp_adj_x10.rmp 126 | # Output values are floats 127 | cov_type_remap = covtype.rmp 128 | covden_sum_remap = covdensum.rmp 129 | covden_win_remap = covdenwin.rmp 130 | snow_intcp_remap = snow_intcp.rmp 131 | srain_intcp_remap = srain_intcp.rmp 132 | wrain_intcp_remap = wrain_intcp.rmp 133 | root_depth_remap = rtdepth.rmp 134 | 135 | 136 | ## PRISM Parameters 137 | prism_folder = ..\..\sagehen\prism 138 | # Resampling method: BILINEAR, CUBIC, NEAREST 139 | prism_projection_method = BILINEAR 140 | # Output projected cellsize, not PRISM input cellsize 141 | prism_cellsize = 90 142 | # Recalculate JH coefficient with PRISM temperature values 143 | calc_prism_jh_coef_flag = True 144 | 145 | 146 | # Precipitation Ratios 147 | set_ppt_zones_flag = True 148 | # ppt_calc_method = ZONES 149 | ppt_obs_units = mm 150 | 151 | # If set_ppt_zones_flag is false, mean monthly precipitation must be set manually 152 | # NOTE - These are made up values! 153 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 154 | # If a valid non-zero HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 155 | ppt_hru_id = 0 156 | 157 | # If set_ppt_zones_flag is True, mean monthly precipitation will be read from 158 | # the precipitation zones shapefile 159 | # Precipitation ratios will be adjusted to be 1 at each zone HRU_ID 160 | # (set in ppt_hru_id_field) 161 | ppt_zone_path = .\shapefiles\ppt_zones.shp 162 | ppt_zone_id_field = PPT_ZONE 163 | ppt_hru_id_field = PPT_HRU_ID 164 | ppt_obs_field_format = PPT_{:02d} 165 | 166 | 167 | ## PRMS Parameter 168 | prms_parameter_folder = .\hru_params 169 | 170 | single_param_file_flag = False 171 | # single_param_file_name = prms_inputs.param 172 | param_column_flag = True 173 | 174 | # Default/template values 175 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 176 | prms_param_csv_path = .\hru_params\prms_parameters.csv 177 | 178 | # Parameters passed through to PRMS param file(s) 179 | # Number of air temperature measurement stations in the data file 180 | ntemp = 1 181 | # Elevation units (0=feet, 1=meter) 182 | elev_units = 0 183 | # Index of the temperature station used to compute basin temperature values (1's based) 184 | basin_tsta = 1 185 | # Index of the base temperature station used for lapse rate calculations (1's based) 186 | hru_tsta = 1 187 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 188 | hru_tlaps = 1 189 | # Elevation of each air temperature measurment station 190 | tsta_elev = 6340 191 | 192 | # If using temp_1sta module, manually define monthly lapse rates 193 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 194 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 195 | 196 | 197 | ## Miscellaneous 198 | # Scratch Workspace ('in_memory' or 'scratch') 199 | # in_memory may not work for computers with limited RAM 200 | # scratch will write all intermediate files to a scratch folder on the disk 201 | scratch_name = in_memory 202 | # scratch_name = scratch 203 | 204 | # Scale floating point values before converting to Int and calculating Median 205 | int_factor = 1 206 | 207 | # Don't change this unless it conflicts with an existing field 208 | orig_fid_field = ORIG_FID 209 | -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\..\scripts\fishnet_generator.py -i ppt_zones_parameters.ini --overwrite 2 | 3 | python ..\..\..\scripts\hru_parameters.py -i ppt_zones_parameters.ini 4 | python ..\..\..\scripts\dem_parameters.py -i ppt_zones_parameters.ini 5 | python ..\..\..\scripts\dem_2_streams.py -i ppt_zones_parameters.ini 6 | python ..\..\..\scripts\crt_fill_parameters.py -i ppt_zones_parameters.ini 7 | python ..\..\..\scripts\dem_2_streams.py -i ppt_zones_parameters.ini 8 | python ..\..\..\scripts\crt_fill_parameters.py -i ppt_zones_parameters.ini 9 | python ..\..\..\scripts\stream_parameters.py -i ppt_zones_parameters.ini 10 | python ..\..\..\scripts\veg_parameters.py -i ppt_zones_parameters.ini 11 | python ..\..\..\scripts\soil_raster_prep.py -i ppt_zones_parameters.ini 12 | python ..\..\..\scripts\soil_parameters.py -i ppt_zones_parameters.ini 13 | python ..\..\..\scripts\prism_800m_normals.py -i ppt_zones_parameters.ini 14 | python ..\..\..\scripts\ppt_ratio_parameters.py -i ppt_zones_parameters.ini 15 | python ..\..\..\scripts\impervious_parameters.py -i ppt_zones_parameters.ini 16 | python ..\..\..\scripts\prms_template_fill.py -i ppt_zones_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | va=TYPEC2ZONE_VALUEN 2 | OUTLET 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.dbf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.dbf -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]] -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/ppt_zones.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | v 2 | A HRU_TYPEN 3 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/precipitation_zones/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/precipitation_zones/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/dem/ned90m_nad83z11.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/dem/ned90m_nad83z11.img -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/hru_params/prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,1 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,0 18 | nsnow,calibration_parameters,0 19 | nsol,calibration_parameters,0 20 | nssr,calibration_parameters,calculated 21 | nsub,calibration_parameters,calculated 22 | ntemp,calibration_parameters,config_file 23 | one,calibration_parameters,1 24 | nlake_hrus,calibration_parameters,0 25 | -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\..\scripts\fishnet_generator.py -i temp_zones_parameters.ini --overwrite 2 | 3 | python ..\..\..\scripts\hru_parameters.py -i temp_zones_parameters.ini 4 | python ..\..\..\scripts\dem_parameters.py -i temp_zones_parameters.ini 5 | python ..\..\..\scripts\dem_2_streams.py -i temp_zones_parameters.ini 6 | python ..\..\..\scripts\crt_fill_parameters.py -i temp_zones_parameters.ini 7 | python ..\..\..\scripts\dem_2_streams.py -i temp_zones_parameters.ini 8 | python ..\..\..\scripts\crt_fill_parameters.py -i temp_zones_parameters.ini 9 | python ..\..\..\scripts\stream_parameters.py -i temp_zones_parameters.ini 10 | python ..\..\..\scripts\veg_parameters.py -i temp_zones_parameters.ini 11 | python ..\..\..\scripts\soil_raster_prep.py -i temp_zones_parameters.ini 12 | python ..\..\..\scripts\soil_parameters.py -i temp_zones_parameters.ini 13 | python ..\..\..\scripts\prism_800m_normals.py -i temp_zones_parameters.ini 14 | python ..\..\..\scripts\temp_adjust_parameters.py -i temp_zones_parameters.ini 15 | python ..\..\..\scripts\impervious_parameters.py -i temp_zones_parameters.ini 16 | python ..\..\..\scripts\prms_template_fill.py -i temp_zones_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | va=TYPEC2ZONE_VALUEN 2 | OUTLET 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.dbf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.dbf -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/temp_zones.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | v 2 | A HRU_TYPEN 3 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/temperature_zones/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/temp_1sta_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs - Sagehen test model - Temperature zones 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | 38 | 39 | ## DEM Parameters 40 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 41 | dem_orig_path = .\dem\ned90m_nad83z11.img 42 | dem_units = meters 43 | # Resampling method: BILINEAR, CUBIC, NEAREST 44 | dem_projection_method = NEAREST 45 | dem_cellsize = 90 46 | # Calculate flow accumulation weighted elevation 47 | calc_flow_acc_dem_flag = False 48 | # flow_acc_dem_factor = 0.001 49 | # Field to initially set DEM_ADJ 50 | dem_adj_copy_field = DEM_MEAN 51 | reset_dem_adj_flag = False 52 | 53 | 54 | ## DEM 2 Streams 55 | # Cells with flow accumulations >= threshold will be designated as stream cells 56 | flow_acc_threshold = 100 57 | # All 1st order streams with a length below threshold will be removed 58 | flow_length_threshold = 3 59 | 60 | calc_flow_dir_points_flag = True 61 | 62 | 63 | ## CRT Parameters 64 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 65 | crt_hruflg = 0 66 | crt_flowflg = 3 67 | crt_dpit = 0.01 68 | crt_outitmax = 100000 69 | 70 | 71 | ## CRT Fill Parameters 72 | use_crt_fill_flag = True 73 | 74 | 75 | ## Vegetation Parameters 76 | # Vegetation Type 77 | # Assume NEAREST Resampling 78 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 79 | veg_type_field = VALUE 80 | veg_type_cellsize = 10 81 | 82 | # Vegetation Cover 83 | # Assume NEAREST Resampling 84 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 85 | veg_cover_cellsize = 10 86 | 87 | 88 | ## Soils Parameters 89 | soil_orig_folder = ..\..\sagehen\soils 90 | ksat_name = ksat.img 91 | awc_name = awc.img 92 | clay_pct_name = clay.img 93 | sand_pct_name = sand.img 94 | soil_cellsize = 50 95 | # Soil rasters are percent (i.e. 25%) 96 | # If false interpret as decimals (i.e. 0.25) 97 | soil_pct_flag = True 98 | # Fill nodata in soil rasters using nibble technique 99 | fill_soil_nodata_flag = True 100 | # Initial soil moisture and recharge are calculated as a fraction of the max 101 | moist_init_ratio = 0.1 102 | rechr_init_ratio = 0.1 103 | ssr2gw_k_default = 0.001 104 | # Setting true will clip the root depth to the soil depth 105 | clip_root_depth_flag = False 106 | # Soil depth will only be read if clip_root_depth_flag = True 107 | # Soil depth must be in inches (?) 108 | soil_depth_name = avg_depth 109 | 110 | 111 | ## Impervious Cover Parameters 112 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 113 | # Resampling method: BILINEAR, CUBIC, NEAREST 114 | impervious_projection_method = BILINEAR 115 | impervious_cellsize = 10 116 | # Impervious cover rasters are percent (i.e. 25%) 117 | # If false interpret as decimals (i.e. 0.25) 118 | impervious_pct_flag = True 119 | 120 | 121 | ## Remap Files 122 | remap_folder = ..\..\..\remaps 123 | aspect_remap = aspect.rmp 124 | # Output values are 10 * value 125 | temp_adj_remap = temp_adj_x10.rmp 126 | # Output values are floats 127 | cov_type_remap = covtype.rmp 128 | covden_sum_remap = covdensum.rmp 129 | covden_win_remap = covdenwin.rmp 130 | snow_intcp_remap = snow_intcp.rmp 131 | srain_intcp_remap = srain_intcp.rmp 132 | wrain_intcp_remap = wrain_intcp.rmp 133 | root_depth_remap = rtdepth.rmp 134 | 135 | 136 | ## PRISM Parameters 137 | prism_folder = ..\..\sagehen\prism 138 | # Resampling method: BILINEAR, CUBIC, NEAREST 139 | prism_projection_method = BILINEAR 140 | # Output projected cellsize, not PRISM input cellsize 141 | prism_cellsize = 90 142 | # Recalculate JH coefficient with PRISM temperature values 143 | calc_prism_jh_coef_flag = True 144 | 145 | 146 | # Temperature Calculation (ZONES, 1STA, LAPSE) 147 | temperature_calc_method = 1STA 148 | temp_obs_units = C 149 | # If temperature_calc_method is 1STA, mean monthly temperature must be set manually 150 | # NOTE - These are made up values! 151 | tmax_obs_list = 5, 5, 7, 10, 15, 20, 25, 25, 20, 15, 10, 5 152 | tmin_obs_list = -7, -8, -6, -4, 0, 3, 6, 6, 3, 0, -4, -7 153 | # If a valid non-zero HRU_ID is set, temperature ratios will be adjusted to be 1 at that HRU_ID 154 | temp_hru_id = 0 155 | 156 | 157 | ## PRMS Parameter 158 | prms_parameter_folder = .\hru_params 159 | 160 | single_param_file_flag = False 161 | # single_param_file_name = prms_inputs.param 162 | param_column_flag = True 163 | 164 | # Default/template values 165 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 166 | prms_param_csv_path = .\hru_params\prms_parameters.csv 167 | 168 | # Parameters passed through to PRMS param file(s) 169 | # Number of air temperature measurement stations in the data file 170 | ntemp = 1 171 | # Elevation units (0=feet, 1=meter) 172 | elev_units = 0 173 | # Index of the temperature station used to compute basin temperature values (1's based) 174 | basin_tsta = 1 175 | # Index of the base temperature station used for lapse rate calculations (1's based) 176 | hru_tsta = 1 177 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 178 | hru_tlaps = 1 179 | # Elevation of each air temperature measurment station 180 | tsta_elev = 6340 181 | 182 | # # If using temp_1sta module, manually define monthly lapse rates 183 | # tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 184 | # tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 185 | 186 | 187 | ## Miscellaneous 188 | # Scratch Workspace ('in_memory' or 'scratch') 189 | # in_memory may not work for computers with limited RAM 190 | # scratch will write all intermediate files to a scratch folder on the disk 191 | scratch_name = in_memory 192 | # scratch_name = scratch 193 | 194 | # Scale floating point values before converting to Int and calculating Median 195 | int_factor = 1 196 | 197 | # Don't change this unless it conflicts with an existing field 198 | orig_fid_field = ORIG_FID 199 | -------------------------------------------------------------------------------- /examples/sagehen_tests/temperature_zones/temp_zones_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs - Sagehen test model - Temperature zones 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | 38 | 39 | ## DEM Parameters 40 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 41 | dem_orig_path = .\dem\ned90m_nad83z11.img 42 | dem_units = meters 43 | # Resampling method: BILINEAR, CUBIC, NEAREST 44 | dem_projection_method = NEAREST 45 | dem_cellsize = 90 46 | # Calculate flow accumulation weighted elevation 47 | calc_flow_acc_dem_flag = False 48 | # flow_acc_dem_factor = 0.001 49 | # Field to initially set DEM_ADJ 50 | dem_adj_copy_field = DEM_MEAN 51 | reset_dem_adj_flag = False 52 | 53 | 54 | ## DEM 2 Streams 55 | # Cells with flow accumulations >= threshold will be designated as stream cells 56 | flow_acc_threshold = 100 57 | # All 1st order streams with a length below threshold will be removed 58 | flow_length_threshold = 3 59 | 60 | calc_flow_dir_points_flag = True 61 | 62 | 63 | ## CRT Parameters 64 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 65 | crt_hruflg = 0 66 | crt_flowflg = 3 67 | crt_dpit = 0.01 68 | crt_outitmax = 100000 69 | 70 | 71 | ## CRT Fill Parameters 72 | use_crt_fill_flag = True 73 | 74 | 75 | ## Vegetation Parameters 76 | # Vegetation Type 77 | # Assume NEAREST Resampling 78 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 79 | veg_type_field = VALUE 80 | veg_type_cellsize = 10 81 | 82 | # Vegetation Cover 83 | # Assume NEAREST Resampling 84 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 85 | veg_cover_cellsize = 10 86 | 87 | 88 | ## Soils Parameters 89 | soil_orig_folder = ..\..\sagehen\soils 90 | ksat_name = ksat.img 91 | awc_name = awc.img 92 | clay_pct_name = clay.img 93 | sand_pct_name = sand.img 94 | soil_cellsize = 50 95 | # Soil rasters are percent (i.e. 25%) 96 | # If false interpret as decimals (i.e. 0.25) 97 | soil_pct_flag = True 98 | # Fill nodata in soil rasters using nibble technique 99 | fill_soil_nodata_flag = True 100 | # Initial soil moisture and recharge are calculated as a fraction of the max 101 | moist_init_ratio = 0.1 102 | rechr_init_ratio = 0.1 103 | ssr2gw_k_default = 0.001 104 | # Setting true will clip the root depth to the soil depth 105 | clip_root_depth_flag = False 106 | # Soil depth will only be read if clip_root_depth_flag = True 107 | # Soil depth must be in inches (?) 108 | soil_depth_name = avg_depth 109 | 110 | 111 | ## Impervious Cover Parameters 112 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 113 | # Resampling method: BILINEAR, CUBIC, NEAREST 114 | impervious_projection_method = BILINEAR 115 | impervious_cellsize = 10 116 | # Impervious cover rasters are percent (i.e. 25%) 117 | # If false interpret as decimals (i.e. 0.25) 118 | impervious_pct_flag = True 119 | 120 | 121 | ## Remap Files 122 | remap_folder = ..\..\..\remaps 123 | aspect_remap = aspect.rmp 124 | # Output values are 10 * value 125 | temp_adj_remap = temp_adj_x10.rmp 126 | # Output values are floats 127 | cov_type_remap = covtype.rmp 128 | covden_sum_remap = covdensum.rmp 129 | covden_win_remap = covdenwin.rmp 130 | snow_intcp_remap = snow_intcp.rmp 131 | srain_intcp_remap = srain_intcp.rmp 132 | wrain_intcp_remap = wrain_intcp.rmp 133 | root_depth_remap = rtdepth.rmp 134 | 135 | 136 | ## PRISM Parameters 137 | prism_folder = ..\..\sagehen\prism 138 | # Resampling method: BILINEAR, CUBIC, NEAREST 139 | prism_projection_method = BILINEAR 140 | # Output projected cellsize, not PRISM input cellsize 141 | prism_cellsize = 90 142 | # Recalculate JH coefficient with PRISM temperature values 143 | calc_prism_jh_coef_flag = True 144 | 145 | 146 | # Temperature Calculation (ZONES, 1STA, LAPSE) 147 | temperature_calc_method = ZONES 148 | temp_obs_units = C 149 | # If temperature_calc_method is ZONES, mean monthly tmax/tmin will be read from a 150 | # temperature zones shapefile. 151 | # Temperature ratios will be adjusted to be 1 at each zone HRU_ID 152 | # (set in temp_hru_id_field) 153 | temp_zone_path = .\shapefiles\temp_zones.shp 154 | temp_zone_id_field = TEMP_ZONE 155 | temp_hru_id_field = TMP_HRU_ID 156 | tmax_obs_field_format = TMAX_{:02d} 157 | tmin_obs_field_format = TMIN_{:02d} 158 | 159 | 160 | ## PRMS Parameter 161 | prms_parameter_folder = .\hru_params 162 | 163 | single_param_file_flag = False 164 | # single_param_file_name = prms_inputs.param 165 | param_column_flag = True 166 | 167 | # Default/template values 168 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 169 | prms_param_csv_path = .\hru_params\prms_parameters.csv 170 | 171 | # Parameters passed through to PRMS param file(s) 172 | # Number of air temperature measurement stations in the data file 173 | ntemp = 1 174 | # Elevation units (0=feet, 1=meter) 175 | elev_units = 0 176 | # Index of the temperature station used to compute basin temperature values (1's based) 177 | basin_tsta = 1 178 | # Index of the base temperature station used for lapse rate calculations (1's based) 179 | # hru_tsta = 1 180 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 181 | hru_tlaps = 1 182 | # Elevation of each air temperature measurment station 183 | tsta_elev = 6340 184 | 185 | 186 | ## Miscellaneous 187 | # Scratch Workspace ('in_memory' or 'scratch') 188 | # in_memory may not work for computers with limited RAM 189 | # scratch will write all intermediate files to a scratch folder on the disk 190 | scratch_name = in_memory 191 | # scratch_name = scratch 192 | 193 | # Scale floating point values before converting to Int and calculating Median 194 | int_factor = 1 195 | 196 | # Don't change this unless it conflicts with an existing field 197 | orig_fid_field = ORIG_FID 198 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/dem/ned90m_nad83z11_terminal_basin.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/dem/ned90m_nad83z11_terminal_basin.img -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/hru_params/prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,1 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,calculated 18 | nlake_hrus,calibration_parameters,calculated 19 | nsnow,calibration_parameters,0 20 | nsol,calibration_parameters,0 21 | nssr,calibration_parameters,calculated 22 | nsub,calibration_parameters,calculated 23 | ntemp,calibration_parameters,config_file 24 | one,calibration_parameters,1 25 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\..\scripts\fishnet_generator.py -i terminal_basin_parameters.ini --overwrite 2 | 3 | python ..\..\..\scripts\hru_parameters.py -i terminal_basin_parameters.ini 4 | python ..\..\..\scripts\dem_parameters.py -i terminal_basin_parameters.ini 5 | python ..\..\..\scripts\dem_2_streams.py -i terminal_basin_parameters.ini 6 | python ..\..\..\scripts\crt_fill_parameters.py -i terminal_basin_parameters.ini 7 | python ..\..\..\scripts\dem_2_streams.py -i terminal_basin_parameters.ini 8 | python ..\..\..\scripts\crt_fill_parameters.py -i terminal_basin_parameters.ini 9 | python ..\..\..\scripts\stream_parameters.py -i terminal_basin_parameters.ini 10 | python ..\..\..\scripts\veg_parameters.py -i terminal_basin_parameters.ini 11 | python ..\..\..\scripts\soil_raster_prep.py -i terminal_basin_parameters.ini 12 | python ..\..\..\scripts\soil_parameters.py -i terminal_basin_parameters.ini 13 | python ..\..\..\scripts\prism_800m_normals.py -i terminal_basin_parameters.ini 14 | python ..\..\..\scripts\ppt_ratio_parameters.py -i terminal_basin_parameters.ini 15 | python ..\..\..\scripts\impervious_parameters.py -i terminal_basin_parameters.ini 16 | python ..\..\..\scripts\prms_template_fill.py -i terminal_basin_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | v 2 | a=TYPEC2ZONE_VALUEN 3 | SWALE 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | v 2 | A HRU_TYPEN 3 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_basin/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_basin/terminal_basin_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | 38 | 39 | ## DEM Parameters 40 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 41 | dem_orig_path = .\dem\ned90m_nad83z11_terminal_basin.img 42 | dem_units = meters 43 | # Resampling method: BILINEAR, CUBIC, NEAREST 44 | dem_projection_method = NEAREST 45 | dem_cellsize = 90 46 | # Calculate flow accumulation weighted elevation 47 | calc_flow_acc_dem_flag = False 48 | # flow_acc_dem_factor = 0.001 49 | # Field to initially set DEM_ADJ 50 | dem_adj_copy_field = DEM_MEAN 51 | reset_dem_adj_flag = False 52 | 53 | 54 | ## DEM 2 Streams 55 | # Cells with flow accumulations >= threshold will be designated as stream cells 56 | flow_acc_threshold = 100 57 | # All 1st order streams with a length below threshold will be removed 58 | flow_length_threshold = 2 59 | 60 | calc_flow_dir_points_flag = True 61 | 62 | 63 | ## CRT Parameters 64 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 65 | crt_hruflg = 0 66 | crt_flowflg = 3 67 | crt_dpit = 0.01 68 | crt_outitmax = 100000 69 | 70 | 71 | ## CRT Fill Parameters 72 | use_crt_fill_flag = True 73 | 74 | 75 | ## Vegetation Parameters 76 | # Vegetation Type 77 | # Assume NEAREST Resampling 78 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 79 | veg_type_field = VALUE 80 | veg_type_cellsize = 10 81 | 82 | # Vegetation Cover 83 | # Assume NEAREST Resampling 84 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 85 | veg_cover_cellsize = 10 86 | 87 | 88 | ## Soils Parameters 89 | soil_orig_folder = ..\..\sagehen\soils 90 | ksat_name = ksat.img 91 | awc_name = awc.img 92 | clay_pct_name = clay.img 93 | sand_pct_name = sand.img 94 | soil_cellsize = 50 95 | # Soil rasters are percent (i.e. 25%) 96 | # If false interpret as decimals (i.e. 0.25) 97 | soil_pct_flag = True 98 | # Fill nodata in soil rasters using nibble technique 99 | fill_soil_nodata_flag = True 100 | # Initial soil moisture and recharge are calculated as a fraction of the max 101 | moist_init_ratio = 0.1 102 | rechr_init_ratio = 0.1 103 | ssr2gw_k_default = 0.001 104 | # Setting true will clip the root depth to the soil depth 105 | clip_root_depth_flag = False 106 | # Soil depth will only be read if clip_root_depth_flag = True 107 | # Soil depth must be in inches (?) 108 | soil_depth_name = avg_depth 109 | 110 | 111 | ## Impervious Cover Parameters 112 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 113 | # Resampling method: BILINEAR, CUBIC, NEAREST 114 | impervious_projection_method = BILINEAR 115 | impervious_cellsize = 10 116 | # Impervious cover rasters are percent (i.e. 25%) 117 | # If false interpret as decimals (i.e. 0.25) 118 | impervious_pct_flag = True 119 | 120 | 121 | ## Remap Files 122 | remap_folder = ..\..\..\remaps 123 | aspect_remap = aspect.rmp 124 | # Output values are 10 * value 125 | temp_adj_remap = temp_adj_x10.rmp 126 | # Output values are floats 127 | cov_type_remap = covtype.rmp 128 | covden_sum_remap = covdensum.rmp 129 | covden_win_remap = covdenwin.rmp 130 | snow_intcp_remap = covtype_to_snow_intcp.rmp 131 | srain_intcp_remap = covtype_to_srain_intcp.rmp 132 | wrain_intcp_remap = covtype_to_wrain_intcp.rmp 133 | root_depth_remap = rtdepth.rmp 134 | 135 | 136 | ## PRISM Parameters 137 | prism_folder = ..\..\sagehen\prism 138 | # Resampling method: BILINEAR, CUBIC, NEAREST 139 | prism_projection_method = BILINEAR 140 | # Output projected cellsize, not PRISM input cellsize 141 | prism_cellsize = 90 142 | # Recalculate JH coefficient with PRISM temperature values 143 | calc_prism_jh_coef_flag = True 144 | 145 | # PPT Ratios 146 | set_ppt_zones_flag = False 147 | ppt_obs_units = mm 148 | 149 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 150 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 151 | 152 | 153 | ## PRMS Parameter 154 | prms_parameter_folder = .\hru_params 155 | 156 | single_param_file_flag = False 157 | # single_param_file_name = prms_inputs.param 158 | param_column_flag = True 159 | 160 | # Default/template values 161 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 162 | prms_param_csv_path = .\hru_params\prms_parameters.csv 163 | 164 | # Parameters passed through to PRMS param file(s) 165 | # Number of air temperature measurement stations in the data file 166 | ntemp = 1 167 | # Elevation units (0=feet, 1=meter) 168 | elev_units = 0 169 | # Index of the temperature station used to compute basin temperature values (1's based) 170 | basin_tsta = 1 171 | # Index of the base temperature station used for lapse rate calculations (1's based) 172 | hru_tsta = 1 173 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 174 | hru_tlaps = 1 175 | # Elevation of each air temperature measurment station 176 | tsta_elev = 6340 177 | 178 | # If using temp_1sta module, manually define monthly lapse rates 179 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 180 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 181 | 182 | 183 | ## Miscellaneous 184 | # Scratch Workspace ('in_memory' or 'scratch') 185 | # in_memory may not work for computers with limited RAM 186 | # scratch will write all intermediate files to a scratch folder on the disk 187 | scratch_name = in_memory 188 | # scratch_name = scratch 189 | 190 | # Scale floating point values before converting to Int and calculating Median 191 | int_factor = 1 192 | 193 | # Don't change this unless it conflicts with an existing field 194 | orig_fid_field = ORIG_FID 195 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/dem/ned90m_nad83z11_terminal_lake.img: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/dem/ned90m_nad83z11_terminal_lake.img -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/hru_params/prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,1 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,calculated 18 | nlake_hrus,calibration_parameters,calculated 19 | nsnow,calibration_parameters,0 20 | nsol,calibration_parameters,0 21 | nssr,calibration_parameters,calculated 22 | nsub,calibration_parameters,calculated 23 | ntemp,calibration_parameters,config_file 24 | one,calibration_parameters,1 25 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/hru_params/prms_parameters.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,DIMENSION_NAMES,TYPE,DEFAULT_VALUE,,,,,,,,,,,,,,,,,,,,, 2 | adjmix_rain,calibration_parameters,nmonths,2,1,1,1,0.01,0.2,1.2,1,1,1,1.1,0.5,1.2,,,,,,,,,, 3 | albset_rna,default_values,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 4 | albset_rnm,default_values,one,2,0.6,,,,,,,,,,,,,,,,,,,,, 5 | albset_sna,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 6 | albset_snm,default_values,one,2,0.2,,,,,,,,,,,,,,,,,,,,, 7 | basin_tsta,calibration_parameters,one,1,config_file,,,,,,,,,,,,,,,,,,,,, 8 | carea_max,calibration_parameters,nhru,2,CAREA_MAX,,,,,,,,,,,,,,,,,,,,, 9 | cascade_flg,cascade,one,1,0,,,,,,,,,,,,,,,,,,,,, 10 | cascade_tol,cascade,one,2,5,,,,,,,,,,,,,,,,,,,,, 11 | cecn_coef,default_values,nmonths,2,5,,,,,,,,,,,,,,,,,,,,, 12 | circle_switch,cascade,one,1,0,,,,,,,,,,,,,,,,,,,,, 13 | cov_type,gis_derived_parameters,nhru,1,COV_TYPE,,,,,,,,,,,,,,,,,,,,, 14 | covden_sum,gis_derived_parameters,nhru,2,COVDEN_SUM,,,,,,,,,,,,,,,,,,,,, 15 | covden_win,gis_derived_parameters,nhru,2,COVDEN_WIN,,,,,,,,,,,,,,,,,,,,, 16 | dday_intcp,calibration_parameters,nmonths,2,-9,-9,-9,-7.8,-10,-22,-38,-36,-15,-15,-15,-10,,,,,,,,,, 17 | dday_slope,calibration_parameters,nmonths,2,3.10E-01,3.10E-01,3.90E-01,3.40E-01,4.00E-01,5.40E-01,6.50E-01,6.50E-01,4.00E-01,3.40E-01,3.90E-01,3.10E-01,,,,,,,,,, 18 | den_init,default_values,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 19 | den_max,default_values,one,2,0.6,,,,,,,,,,,,,,,,,,,,, 20 | elev_units,default_values,one,1,config_file,,,,,,,,,,,,,,,,,,,,, 21 | emis_noppt,default_values,one,2,0.757,,,,,,,,,,,,,,,,,,,,, 22 | epan_coef,default_values,nmonths,2,1,,,,,,,,,,,,,,,,,,,,, 23 | fastcoef_lin,calibration_parameters,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 24 | fastcoef_sq,calibration_parameters,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 25 | freeh2o_cap,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 26 | gwflow_coef,calibration_parameters,one,2,0.014,,,,,,,,,,,,,,,,,,,,, 27 | gwsink_coef,default_values,one,2,0.01,,,,,,,,,,,,,,,,,,,,, 28 | gwstor_init,default_values,one,2,2,,,,,,,,,,,,,,,,,,,,, 29 | gwstor_min,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 30 | gw_up_id,cascade,ncascdgw,1,crt_file,,,,,,,,,,,,,,,,,,,,, 31 | gw_down_id,cascade,ncascdgw,1,crt_file,,,,,,,,,,,,,,,,,,,,, 32 | gvr_cell_id,cascade,nhrucell,1,HRU_ID,,,,,,,,,,,,,,,,,,,,, 33 | gw_pct_up,cascade,ncascdgw,2,crt_file,,,,,,,,,,,,,,,,,,,,, 34 | gw_strmseg_down_id,cascade,ncascdgw,1,crt_file,,,,,,,,,,,,,,,,,,,,, 35 | hru_area,gis_derived_parameters,nhru,2,HRU_AREA,,,,,,,,,,,,,,,,,,,,, 36 | hru_aspect,gis_derived_parameters,nhru,2,DEM_ASPECT,,,,,,,,,,,,,,,,,,,,, 37 | hru_deplcrv,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 38 | hru_down_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 39 | hru_elev,gis_derived_parameters,nhru,2,DEM_ADJ,,,,,,,,,,,,,,,,,,,,, 40 | hru_lat,gis_derived_parameters,nhru,2,HRU_LAT,,,,,,,,,,,,,,,,,,,,, 41 | hru_pct_up,cascade,ncascade,2,crt_file,,,,,,,,,,,,,,,,,,,,, 42 | hru_percent_imperv,calibration_parameters,nhru,2,IMPERV_PCT,,,,,,,,,,,,,,,,,,,,, 43 | hru_psta,gis_derived_parameters,nhru,1,HRU_PSTA,,,,,,,,,,,,,,,,,,,,, 44 | hru_slope,gis_derived_parameters,nhru,2,DEM_SLP_P,,,,,,,,,,,,,,,,,,,,, 45 | hru_strmseg_down_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 46 | hru_subbasin,calibration_parameters,nhru,1,SUB_BASIN,,,,,,,,,,,,,,,,,,,,, 47 | hru_tlaps,calibration_parameters,nhru,1,config_file,,,,,,,,,,,,,,,,,,,,, 48 | hru_tsta,calibration_parameters,nhru,1,config_file,,,,,,,,,,,,,,,,,,,,, 49 | hru_type,gis_derived_parameters,nhru,1,HRU_TYPE,,,,,,,,,,,,,,,,,,,,, 50 | hru_up_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 51 | imperv_stor_max,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 52 | jh_coef,calibration_parameters,nmonths,2,0.02200635,0.02200635,0.02200635,0.02200635,0.01560465,0.01560465,0.0148041,0.01380375,0.012204,0.01300455,0.02000565,0.02000565,,,,,,,,,, 53 | jh_coef_hru,gis_derived_parameters,nhru,2,JH_COEF,,,,,,,,,,,,,,,,,,,,, 54 | lake_hru_id,gis_derived_parameters,nhru,1,LAKE_ID,,,,,,,,,,,,,,,,,,,,, 55 | mapvars_freq,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 56 | mapvars_units,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 57 | max_missing,default_values,one,1,3,,,,,,,,,,,,,,,,,,,,, 58 | melt_force,default_values,one,1,140,,,,,,,,,,,,,,,,,,,,, 59 | melt_look,default_values,one,1,90,,,,,,,,,,,,,,,,,,,,, 60 | outlet_sta,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 61 | potet_sublim,default_values,one,2,0.5,,,,,,,,,,,,,,,,,,,,, 62 | ppt_rad_adj,default_values,nmonths,2,0.02,,,,,,,,,,,,,,,,,,,,, 63 | precip_units,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 64 | pref_flow_den,calibration_parameters,one,2,0.2,,,,,,,,,,,,,,,,,,,,, 65 | print_freq,default_values,one,1,3,,,,,,,,,,,,,,,,,,,,, 66 | print_type,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 67 | prms_warmup,default_values,one,1,2,,,,,,,,,,,,,,,,,,,,, 68 | rad_trncf,gis_derived_parameters,nhru,2,RAD_TRNCF,,,,,,,,,,,,,,,,,,,,, 69 | radadj_intcp,default_values,one,2,1,,,,,,,,,,,,,,,,,,,,, 70 | radadj_slope,default_values,one,2,0.02,,,,,,,,,,,,,,,,,,,,, 71 | radj_sppt,default_values,one,2,0.44,,,,,,,,,,,,,,,,,,,,, 72 | radj_wppt,default_values,one,2,0.5,,,,,,,,,,,,,,,,,,,,, 73 | radmax,default_values,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 74 | rain_adj,calibration_parameters,nhru;nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 75 | runoff_units,default_values,one,1,0,,,,,,,,,,,,,,,,,,,,, 76 | sat_threshold,calibration_parameters,nhru,2,MOIST_MAX,,,,,,,,,,,,,,,,,,,,, 77 | settle_const,default_values,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 78 | slowcoef_lin,calibration_parameters,nhru,2,SLOWCOEF_L,,,,,,,,,,,,,,,,,,,,, 79 | slowcoef_sq,calibration_parameters,nhru,2,SLOWCOEF_S,,,,,,,,,,,,,,,,,,,,, 80 | smidx_coef,calibration_parameters,one,2,0.01,,,,,,,,,,,,,,,,,,,,, 81 | smidx_exp,default_values,one,2,0.3,,,,,,,,,,,,,,,,,,,,, 82 | snarea_curve,default_values,ndeplval,2,5.00E-02,2.40E-01,4.00E-01,5.30E-01,6.50E-01,7.50E-01,8.20E-01,8.80E-01,9.30E-01,9.90E-01,1.00E+00,5.00E-02,2.50E-01,4.00E-01,4.80E-01,5.40E-01,5.80E-01,6.10E-01,6.40E-01,6.60E-01,6.80E-01,7.00E-01 83 | snarea_thresh,default_values,nhru,2,SNAREA_THR,,,,,,,,,,,,,,,,,,,,, 84 | snow_adj,calibration_parameters,nhru;nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 85 | snow_intcp,gis_derived_parameters,nhru,2,SNOW_INTCP,,,,,,,,,,,,,,,,,,,,, 86 | snowinfil_max,default_values,one,2,2,,,,,,,,,,,,,,,,,,,,, 87 | soil2gw_max,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 88 | soil_moist_init,gis_derived_parameters,nhru,2,MOIST_INIT,,,,,,,,,,,,,,,,,,,,, 89 | soil_moist_max,calibration_parameters,nhru,2,MOIST_MAX,,,,,,,,,,,,,,,,,,,,, 90 | soil_rechr_init,gis_derived_parameters,nhru,2,RECHR_INIT,,,,,,,,,,,,,,,,,,,,, 91 | soil_rechr_max,calibration_parameters,nhru,2,RECHR_MAX,,,,,,,,,,,,,,,,,,,,, 92 | soil_type,gis_derived_parameters,nhru,1,SOIL_TYPE,,,,,,,,,,,,,,,,,,,,, 93 | srain_intcp,gis_derived_parameters,nhru,2,SRAIN_INTC,,,,,,,,,,,,,,,,,,,,, 94 | ssr2gw_exp,default_values,one,2,1.2,,,,,,,,,,,,,,,,,,,,, 95 | ssr2gw_rate,calibration_parameters,nssr,2,SSR2G_RATE,,,,,,,,,,,,,,,,,,,,, 96 | ssstor_init,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 97 | subbasin_down,calibration_parameters,nsub,1,calculated,,,,,,,,,,,,,,,,,,,,, 98 | temp_units,default_values,one,1,0,,,,,,,,,,,,,,,,,,,,, 99 | tmax_adj,gis_derived_parameters,nhru,2,TMAX_ADJ,,,,,,,,,,,,,,,,,,,,, 100 | tmax_allrain,calibration_parameters,nmonths,2,6.00E+01,6.00E+01,6.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,6.00E+01,6.00E+01,6.00E+01,,,,,,,,,, 101 | tmax_allsnow,calibration_parameters,one,2,32,,,,,,,,,,,,,,,,,,,,, 102 | tmax_index,gis_derived_parameters,nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 103 | tmin_adj,gis_derived_parameters,nhru,2,TMIN_ADJ,,,,,,,,,,,,,,,,,,,,, 104 | transp_beg,default_values,one,1,4,,,,,,,,,,,,,,,,,,,,, 105 | transp_end,default_values,one,1,10,,,,,,,,,,,,,,,,,,,,, 106 | transp_tmax,default_values,one,2,500,,,,,,,,,,,,,,,,,,,,, 107 | tsta_elev,calibration_parameters,ntemp,2,config_file,,,,,,,,,,,,,,,,,,,,, 108 | tstorm_mo,default_values,nmonths,1,0,,,,,,,,,,,,,,,,,,,,, 109 | wrain_intcp,gis_derived_parameters,nhru,2,WRAIN_INTC,,,,,,,,,,,,,,,,,,,,, 110 | tmax_lapse,calibration_parameters,nmonths,2,config_file,,,,,,,,,,,,,,,,,,,,, 111 | tmin_lapse,calibration_parameters,nmonths,2,config_file,,,,,,,,,,,,,,,,,,,,, 112 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/runscripts.bat: -------------------------------------------------------------------------------- 1 | python ..\..\..\scripts\fishnet_generator.py -i terminal_lake_parameters.ini --overwrite 2 | 3 | python ..\..\..\scripts\hru_parameters.py -i terminal_lake_parameters.ini 4 | python ..\..\..\scripts\dem_parameters.py -i terminal_lake_parameters.ini 5 | python ..\..\..\scripts\dem_2_streams.py -i terminal_lake_parameters.ini 6 | python ..\..\..\scripts\crt_fill_parameters.py -i terminal_lake_parameters.ini 7 | python ..\..\..\scripts\dem_2_streams.py -i terminal_lake_parameters.ini 8 | python ..\..\..\scripts\crt_fill_parameters.py -i terminal_lake_parameters.ini 9 | python ..\..\..\scripts\stream_parameters.py -i terminal_lake_parameters.ini 10 | python ..\..\..\scripts\veg_parameters.py -i terminal_lake_parameters.ini 11 | python ..\..\..\scripts\soil_raster_prep.py -i terminal_lake_parameters.ini 12 | python ..\..\..\scripts\soil_parameters.py -i terminal_lake_parameters.ini 13 | python ..\..\..\scripts\prism_800m_normals.py -i terminal_lake_parameters.ini 14 | python ..\..\..\scripts\ppt_ratio_parameters.py -i terminal_lake_parameters.ini 15 | python ..\..\..\scripts\impervious_parameters.py -i terminal_lake_parameters.ini 16 | python ..\..\..\scripts\prms_template_fill.py -i terminal_lake_parameters.ini 17 | -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.dbf: -------------------------------------------------------------------------------- 1 | v 2 | ANAMEC Fake Lake  -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/lakes.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/lakes.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/lakes.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/lakes.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/lakes.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.dbf: -------------------------------------------------------------------------------- 1 | v 2 | a=TYPEC2ZONE_VALUEN 3 | SWALE 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/model_points.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/model_points.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/model_points.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/model_points.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/model_points.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.dbf: -------------------------------------------------------------------------------- 1 | v 2 | A HRU_TYPEN 3 | 1 -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.prj: -------------------------------------------------------------------------------- 1 | PROJCS["NAD_1983_UTM_Zone_11N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-117.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]] -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.sbn: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/watershed.sbn -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.sbx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/watershed.sbx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.shp: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/watershed.shp -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/shapefiles/watershed.shx: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/examples/sagehen_tests/terminal_lake/shapefiles/watershed.shx -------------------------------------------------------------------------------- /examples/sagehen_tests/terminal_lake/terminal_lake_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 90 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = True 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned90m_nad83z11_terminal_lake.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = NEAREST 48 | dem_cellsize = 90 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = False 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_MEAN 54 | reset_dem_adj_flag = False 55 | 56 | 57 | ## DEM 2 Streams 58 | # Cells with flow accumulations >= threshold will be designated as stream cells 59 | flow_acc_threshold = 100 60 | # All 1st order streams with a length below threshold will be removed 61 | flow_length_threshold = 2 62 | 63 | calc_flow_dir_points_flag = True 64 | 65 | 66 | ## CRT Parameters 67 | crt_exe_path = ..\..\..\crt\CRT_1.3.1.exe 68 | crt_hruflg = 0 69 | crt_flowflg = 3 70 | crt_dpit = 0.01 71 | crt_outitmax = 100000 72 | 73 | 74 | ## CRT Fill Parameters 75 | use_crt_fill_flag = True 76 | 77 | 78 | ## Vegetation Parameters 79 | # Vegetation Type 80 | # Assume NEAREST Resampling 81 | veg_type_orig_path = ..\..\sagehen\veg\us_140evt.img 82 | veg_type_field = VALUE 83 | veg_type_cellsize = 10 84 | 85 | # Vegetation Cover 86 | # Assume NEAREST Resampling 87 | veg_cover_orig_path = ..\..\sagehen\veg\us_140evc.img 88 | veg_cover_cellsize = 10 89 | 90 | 91 | ## Soils Parameters 92 | soil_orig_folder = ..\..\sagehen\soils 93 | ksat_name = ksat.img 94 | awc_name = awc.img 95 | clay_pct_name = clay.img 96 | sand_pct_name = sand.img 97 | soil_cellsize = 50 98 | # Soil rasters are percent (i.e. 25%) 99 | # If false interpret as decimals (i.e. 0.25) 100 | soil_pct_flag = True 101 | # Fill nodata in soil rasters using nibble technique 102 | fill_soil_nodata_flag = True 103 | # Initial soil moisture and recharge are calculated as a fraction of the max 104 | moist_init_ratio = 0.1 105 | rechr_init_ratio = 0.1 106 | ssr2gw_k_default = 0.001 107 | # Setting true will clip the root depth to the soil depth 108 | clip_root_depth_flag = False 109 | # Soil depth will only be read if clip_root_depth_flag = True 110 | # Soil depth must be in inches (?) 111 | soil_depth_name = avg_depth 112 | 113 | 114 | ## Impervious Cover Parameters 115 | impervious_orig_path = ..\..\sagehen\impervious\nlcd2011_imp.img 116 | # Resampling method: BILINEAR, CUBIC, NEAREST 117 | impervious_projection_method = BILINEAR 118 | impervious_cellsize = 10 119 | # Impervious cover rasters are percent (i.e. 25%) 120 | # If false interpret as decimals (i.e. 0.25) 121 | impervious_pct_flag = True 122 | 123 | 124 | ## Remap Files 125 | remap_folder = ..\..\..\remaps 126 | aspect_remap = aspect.rmp 127 | # Output values are 10 * value 128 | temp_adj_remap = temp_adj_x10.rmp 129 | # Output values are floats 130 | cov_type_remap = covtype.rmp 131 | covden_sum_remap = covdensum.rmp 132 | covden_win_remap = covdenwin.rmp 133 | snow_intcp_remap = covtype_to_snow_intcp.rmp 134 | srain_intcp_remap = covtype_to_srain_intcp.rmp 135 | wrain_intcp_remap = covtype_to_wrain_intcp.rmp 136 | root_depth_remap = rtdepth.rmp 137 | 138 | 139 | ## PRISM Parameters 140 | prism_folder = ..\..\sagehen\prism 141 | # Resampling method: BILINEAR, CUBIC, NEAREST 142 | prism_projection_method = BILINEAR 143 | # Output projected cellsize, not PRISM input cellsize 144 | prism_cellsize = 90 145 | # Recalculate JH coefficient with PRISM temperature values 146 | calc_prism_jh_coef_flag = True 147 | 148 | # PPT Ratios 149 | set_ppt_zones_flag = False 150 | ppt_obs_units = mm 151 | 152 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 153 | ppt_obs_list = 140, 133, 119, 54, 38, 18, 11, 15, 22, 52, 103, 142 154 | 155 | 156 | ## PRMS Parameter 157 | prms_parameter_folder = .\hru_params 158 | 159 | single_param_file_flag = False 160 | # single_param_file_name = prms_inputs.param 161 | param_column_flag = True 162 | 163 | # Default/template values 164 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 165 | prms_param_csv_path = .\hru_params\prms_parameters.csv 166 | 167 | # Parameters passed through to PRMS param file(s) 168 | # Number of air temperature measurement stations in the data file 169 | ntemp = 1 170 | # Elevation units (0=feet, 1=meter) 171 | elev_units = 0 172 | # Index of the temperature station used to compute basin temperature values (1's based) 173 | basin_tsta = 1 174 | # Index of the base temperature station used for lapse rate calculations (1's based) 175 | hru_tsta = 1 176 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 177 | hru_tlaps = 1 178 | # Elevation of each air temperature measurment station 179 | tsta_elev = 6340 180 | 181 | # If using temp_1sta module, manually define monthly lapse rates 182 | tmax_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 183 | tmin_lapse = 9, 9, 7, 7, 6, 5, 5, 4, 6, 7, 8, 9 184 | 185 | 186 | ## Miscellaneous 187 | # Scratch Workspace ('in_memory' or 'scratch') 188 | # in_memory may not work for computers with limited RAM 189 | # scratch will write all intermediate files to a scratch folder on the disk 190 | scratch_name = in_memory 191 | # scratch_name = scratch 192 | 193 | # Scale floating point values before converting to Int and calculating Median 194 | int_factor = 1 195 | 196 | # Don't change this unless it conflicts with an existing field 197 | orig_fid_field = ORIG_FID 198 | -------------------------------------------------------------------------------- /examples/template/template_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 300 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is needed for extracting fishnet cell raster values 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, ) 44 | dem_orig_path = .\dem\ned10m_nad83.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = BILINEAR 48 | dem_cellsize = 10 49 | # Calculate flow accumulation weighted elevation 50 | calc_flow_acc_dem_flag = True 51 | # flow_acc_dem_factor = 0.001 52 | # Field to initially set DEM_ADJ 53 | dem_adj_copy_field = DEM_FLOWAC 54 | reset_dem_adj_flag = False 55 | # Round DEM_ADJ to the given number of decimals 56 | dem_adj_decimals = 2 57 | 58 | 59 | ## DEM 2 Streams 60 | # Cells with flow accumulations >= threshold will be designated as stream cells 61 | flow_acc_threshold = 30 62 | # All 1st order streams with a length below threshold will be removed 63 | flow_length_threshold = 3 64 | 65 | calc_flow_dir_points_flag = True 66 | 67 | 68 | ## CRT Parameters 69 | crt_exe_path = ..\..\crt\CRT_1.3.1.exe 70 | crt_hruflg = 0 71 | crt_flowflg = 3 72 | crt_dpit = 0.01 73 | crt_outitmax = 100000 74 | 75 | 76 | ## CRT Fill Parameters 77 | use_crt_fill_flag = False 78 | 79 | 80 | ## Vegetation Parameters 81 | # Vegetation Type 82 | # Assume NEAREST Resampling 83 | veg_type_orig_path = .\landfire\us_130evt.img 84 | veg_type_field = VALUE 85 | veg_type_cellsize = 10 86 | 87 | # Vegetation Cover 88 | # Assume NEAREST Resampling 89 | veg_cover_orig_path = .\landfire\us_130evc.img 90 | veg_cover_cellsize = 10 91 | 92 | 93 | ## Soil Parameters 94 | soil_orig_folder = .\soils 95 | ksat_name = ksat.img 96 | awc_name = awc.img 97 | clay_pct_name = clay.img 98 | sand_pct_name = sand.img 99 | soil_cellsize = 50 100 | # Soil rasters are percent (i.e. 25%) 101 | # If false interpret as decimals (i.e. 0.25) 102 | soil_pct_flag = True 103 | # Fill nodata in soil rasters using nibble technique 104 | fill_soil_nodata_flag = True 105 | # Initial soil moisture and recharge are calculated as a fraction of the max 106 | moist_init_ratio = 0.1 107 | rechr_init_ratio = 0.1 108 | # Geology based multiplier raster for ssr2gw rate 109 | ssr2gw_mult_flag = True 110 | ssr2gw_mult_name = ssr2g_mult.img 111 | # If no raster provided, a default value must be set 112 | ssr2gw_k_default = 0.001 113 | # Read and apply soil depth raster 114 | # Otherwise soil depth will only be derived from rooting depth 115 | soil_depth_flag = False 116 | # Soil depth will only be read if soil_depth_flag = True 117 | # Soil depth must be in inches 118 | soil_depth_name = soil_depth.img 119 | 120 | 121 | ## Impervious Cover Parameters 122 | impervious_orig_path = .\impervious\nlcd2011_imp.img 123 | # Resampling method: BILINEAR, CUBIC, NEAREST 124 | impervious_projection_method = BILINEAR 125 | impervious_cellsize = 10 126 | # Impervious cover rasters are percent (i.e. 25%) 127 | # If false interpret as decimals (i.e. 0.25) 128 | impervious_pct_flag = True 129 | 130 | 131 | ## Remap Files 132 | remap_folder = ..\..\remaps 133 | aspect_remap = aspect.rmp 134 | # Output values are 10 * value 135 | temp_adj_remap = temp_adj_x10.rmp 136 | # Output values are floats 137 | # temp_adj_remap = temp_adj_new.rmp 138 | # cov_type_remap = covtype_systemgrp.rmp 139 | cov_type_remap = covtype.rmp 140 | covden_sum_remap = covdensum.rmp 141 | covden_win_remap = covdenwin.rmp 142 | snow_intcp_remap = covtype_to_snow_intcp.rmp 143 | srain_intcp_remap = covtype_to_srain_intcp.rmp 144 | wrain_intcp_remap = covtype_to_wrain_intcp.rmp 145 | root_depth_remap = rtdepth.rmp 146 | 147 | 148 | ## PRISM Parameters 149 | prism_folder = .\prism_800m_normals 150 | # prism_folder = .\prism_4km_normals 151 | # Resampling method: BILINEAR, CUBIC, NEAREST 152 | prism_projection_method = BILINEAR 153 | # Output projected cellsize, not PRISM input cellsize 154 | prism_cellsize = 300 155 | # Recalculate JH coefficient with PRISM temperature values 156 | calc_prism_jh_coef_flag = True 157 | 158 | # PPT Ratios 159 | set_ppt_zones_flag = False 160 | ppt_obs_units = mm 161 | 162 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 163 | ppt_obs_list = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 164 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 165 | ppt_hru_id = 0 166 | 167 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 168 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 169 | ppt_zone_path = .\shapefiles\ppt_zones.shp 170 | ppt_zone_id_field = PPT_ZONE 171 | ppt_hru_id_field = PPT_HRU_ID 172 | ppt_obs_field_format = PPT_{:02d} 173 | 174 | 175 | ## PRMS Parameters 176 | prms_parameter_folder = .\hru_params 177 | 178 | single_param_file_flag = False 179 | # single_param_file_name = prms_inputs.param 180 | param_column_flag = False 181 | 182 | # Default/template values 183 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 184 | prms_param_csv_path = .\hru_params\prms_parameters.csv 185 | 186 | # Parameters passed through to PRMS param file(s) 187 | # Number of air temperature measurement stations in the data file 188 | ntemp = 2 189 | # Elevation units (0=feet, 1=meter) 190 | elev_units = 0 191 | # Index of the temperature station used to compute basin temperature values (1's based) 192 | basin_tsta = 1 193 | # Index of the base temperature station used for lapse rate calculations (1's based) 194 | hru_tsta = 1 195 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 196 | hru_tlaps = 2 197 | # Elevation of each air temperature measurment station 198 | tsta_elev = 0, 0 199 | 200 | # If using temp_1sta module, manually define monthly lapse rates 201 | # tmax_lapse = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 202 | # tmin_lapse = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 203 | 204 | 205 | ## Miscellaneous 206 | # Scratch Workspace ('in_memory' or 'scratch') 207 | # in_memory may not work for computers with limited RAM 208 | # scratch will write all intermediate files to a scratch folder on the disk 209 | scratch_name = in_memory 210 | # scratch_name = scratch 211 | 212 | # Scale floating point values before converting to Int and calculating Median 213 | int_factor = 1 214 | 215 | # Don't change this unless it conflicts with an existing field 216 | orig_fid_field = ORIG_FID 217 | -------------------------------------------------------------------------------- /gsflow_arcpy_tutorial.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/gsflow_arcpy_tutorial.pdf -------------------------------------------------------------------------------- /misc/MISC.md: -------------------------------------------------------------------------------- 1 | # Miscellaneous Files 2 | 3 | ## Layers 4 | 5 | dem_adj.lyr 6 | flow_dir_points.lyr -------------------------------------------------------------------------------- /misc/dem_adj.lyr: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/misc/dem_adj.lyr -------------------------------------------------------------------------------- /misc/flow_dir_points.lyr: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/gsflow/gsflow-arcpy/0ee50ba27eacac850a603ca4d8c70308952c8e5b/misc/flow_dir_points.lyr -------------------------------------------------------------------------------- /remaps/REMAPS.md: -------------------------------------------------------------------------------- 1 | # ArcGIS 10.2+ Remap Files -------------------------------------------------------------------------------- /remaps/arc10p1/aspect.rmp: -------------------------------------------------------------------------------- 1 | -1 -1:0 2 | 0 23.0:0 3 | 23.00 68.0:45 4 | 68.00 113.0:90 5 | 113.00 158.0:135 6 | 158.00 203.0:180 7 | 203.00 248.0:225 8 | 248.00 293.0:270 9 | 293.00 338.0:315 10 | 338.00 360.0:0 -------------------------------------------------------------------------------- /remaps/arc10p1/convert_remap_arc10p2.py: -------------------------------------------------------------------------------- 1 | #-------------------------------- 2 | # Name: convert_remap_arc10p2.py 3 | # Purpose: PRMS Remap Modfiy 4 | # Author: Charles Morton 5 | # Created 2015-05-18 6 | # Python: 2.7 7 | #-------------------------------- 8 | 9 | import datetime as dt 10 | import logging 11 | import os 12 | import re 13 | import sys 14 | 15 | ################################################################################ 16 | 17 | def prms_remap_modify(remap_folder): 18 | """Calculate PRMS Remap Modify 19 | 20 | Args: 21 | remap_folder: folder path 22 | Returns: 23 | None 24 | """ 25 | 26 | try: 27 | logging.info('\nPRMS Remap Modify') 28 | 29 | ## Build output folder if necessary 30 | input_ws = os.path.join(remap_folder, 'arc10p1') 31 | output_ws = os.path.join(remap_folder, 'arc10p2') 32 | if not os.path.isdir(output_ws): 33 | os.mkdir(output_ws) 34 | 35 | ## Remove comments from ASCII remap files 36 | for remap_name in os.listdir(input_ws): 37 | remap_input_path = os.path.join(input_ws, remap_name) 38 | remap_output_path = os.path.join(output_ws, remap_name) 39 | if not os.path.isfile(remap_input_path): 40 | continue 41 | if not remap_input_path.lower().endswith('.rmp'): 42 | continue 43 | logging.debug(' Modifying: {0}'.format(remap_name)) 44 | with open(remap_input_path, 'r') as remap_f: 45 | input_lines = remap_f.readlines() 46 | 47 | ## Separate remap values and comments 48 | remap_lines, comment_lines = [], [] 49 | for i, line in enumerate(input_lines): 50 | ## Remove newline characters and extra leading/trailing whitespace 51 | line_split = [x.strip() for x in line.strip().split('/*')] 52 | ## Comments will result in line_split have two items 53 | if len(line_split) == 2: 54 | comment_lines.append('# {0} - {1}'.format( 55 | line_split[0].split(':')[0].strip(), line_split[1])) 56 | if line_split: 57 | remap_lines.append(' : '.join([x.strip() for x in line_split[0].split(':')])) 58 | 59 | ## First write the remap values, 60 | ## Then write the comments (if there are any) 61 | with open(remap_output_path, 'w') as remap_f: 62 | for i, line in enumerate(remap_lines): 63 | remap_f.write(line) 64 | ## Don't write newline character on last line 65 | ## This causes an error in ArcGIS 10.2.2 66 | if (i+1) == len(remap_lines) and not comment_lines: 67 | break 68 | remap_f.write('\n') 69 | 70 | for i, line in enumerate(comment_lines): 71 | ## Replace California with abbreviation 72 | ## Limit comments to 80 characters 73 | remap_f.write(line.replace('California', 'CA')[:79]) 74 | ## Don't write newline character on last line 75 | ## This causes an error in ArcGIS 10.2.2 76 | if (i+1) == len(comment_lines): 77 | break 78 | remap_f.write('\n') 79 | 80 | except: 81 | logging.exception('Unhandled Exception Error\n\n') 82 | raw_input('ENTER') 83 | 84 | ################################################################################ 85 | 86 | if __name__ == '__main__': 87 | workspace = os.getcwd() 88 | 89 | ## Create Basic Logger 90 | logging.basicConfig(level=logging.DEBUG, format='%(message)s') 91 | 92 | ## Run Information 93 | logging.info('\n{0}'.format('#'*80)) 94 | log_f = '{0:<20s} {1}' 95 | logging.info(log_f.format( 96 | 'Run Time Stamp:', dt.datetime.now().isoformat(' '))) 97 | logging.info(log_f.format('Current Directory:', workspace)) 98 | logging.info(log_f.format('Script:', os.path.basename(sys.argv[0]))) 99 | 100 | ## Calculate PRMS Remap Modify 101 | prms_remap_modify(workspace) 102 | -------------------------------------------------------------------------------- /remaps/arc10p1/covdensum.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 /* No data 2 | 11 : 0 /* Open Water 3 | 12 : 0 /* Snow/Ice 4 | 13 : 50 /* Developed-Upland Deciduous Forest 5 | 14 : 60 /* Developed-Upland Evergreen Forest 6 | 15 : 55 /* Developed-Upland Mixed Forest 7 | 16 : 45 /* Developed-Upland Herbaceous 8 | 17 : 35 /* Developed-Upland Shrubland 9 | 22 : 10 /* Developed - Low Intensity 10 | 23 : 5 /* Developed - Medium Intensity 11 | 24 : 0 /* Developed - High Intensity 12 | 25 : 0 /* Developed-Roads 13 | 31 : 0 /* Barren 14 | 32 : 0 /* Quarries-Strip Mines-Gravel Pits 15 | 60 : 80 /* NASS-Orchard 16 | 61 : 80 /* NASS-Vineyard 17 | 63 : 80 /* NASS-Row Crop-Close Grown Crop 18 | 64 : 80 /* NASS-Row Crop 19 | 65 : 80 /* NASS-Close Grown Crop 20 | 66 : 80 /* NASS-Pasture and Hayland 21 | 67 : 80 /* NASS-Pasture and Hayland 22 | 68 : 80 /* NASS-Wheat 23 | 69 : 70 /* NASS-Aquaculture 24 | 75 : 60 /* Herbaceous Semi-dry 25 | 76 : 60 /* Herbaceous Semi-wet 26 | 81 : 80 /* Pasture/Hay 27 | 82 : 80 /* Cultivated Crops 28 | 95 : 80 /* Herbaceous Wetlands 29 | 100 : 10 /* Sparse Vegetation Canopy 30 | 101 : 15 /* Tree Cover >= 10 and < 20% 31 | 102 : 25 /* Tree Cover >= 20 and < 30% 32 | 103 : 35 /* Tree Cover >= 30 and < 40% 33 | 104 : 45 /* Tree Cover >= 40 and < 50% 34 | 105 : 55 /* Tree Cover >= 50 and < 60% 35 | 106 : 65 /* Tree Cover >= 60 and < 70% 36 | 107 : 75 /* Tree Cover >= 70 and < 80% 37 | 108 : 85 /* Tree Cover >= 80 and < 90% 38 | 109 : 95 /* Tree Cover >= 90 and <= 100% 39 | 111 : 15 /* Shrub Cover >= 10 and < 20% 40 | 112 : 25 /* Shrub Cover >= 20 and < 30% 41 | 113 : 35 /* Shrub Cover >= 30 and < 40% 42 | 114 : 45 /* Shrub Cover >= 40 and < 50% 43 | 115 : 55 /* Shrub Cover >= 50 and < 60% 44 | 116 : 65 /* Shrub Cover >= 60 and < 70% 45 | 117 : 75 /* Shrub Cover >= 70 and < 80% 46 | 118 : 85 /* Shrub Cover >= 80 and < 90% 47 | 119 : 95 /* Shrub Cover >= 90 and <= 100% 48 | 121 : 15 /* Herb Cover >= 10 and < 20% 49 | 122 : 25 /* Herb Cover >= 20 and < 30% 50 | 123 : 35 /* Herb Cover >= 30 and < 40% 51 | 124 : 45 /* Herb Cover >= 40 and < 50% 52 | 125 : 55 /* Herb Cover >= 50 and < 60% 53 | 126 : 65 /* Herb Cover >= 60 and < 70% 54 | 127 : 75 /* Herb Cover >= 70 and < 80% 55 | 128 : 85 /* Herb Cover >= 80 and < 90% 56 | 129 : 95 /* Herb Cover >= 90 and <= 100% -------------------------------------------------------------------------------- /remaps/arc10p1/covdenwin.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 /* Bare 2 | 1 : 60 /* Grass 3 | 2 : 40 /* Shrub 4 | 3 : 50 /* Deciduous Trees 5 | 4 : 100 /* Coniferous Trees 6 | -------------------------------------------------------------------------------- /remaps/arc10p1/temp_adj_x10.rmp: -------------------------------------------------------------------------------- 1 | -1 : 0 /* Flat 2 | 0 : -18 /* North 3 | 45 : -10 /* North-East 4 | 90 : 0 /* East 5 | 135 : 10 /* South-East 6 | 180 : 12 /* South 7 | 225 : 18 /* South-West 8 | 270 : 0 /* West 9 | 315 : -10 /* North-West 10 | -------------------------------------------------------------------------------- /remaps/aspect.rmp: -------------------------------------------------------------------------------- 1 | -1 -1 : 0 2 | 0 23.0 : 0 3 | 23.00 68.0 : 45 4 | 68.00 113.0 : 90 5 | 113.00 158.0 : 135 6 | 158.00 203.0 : 180 7 | 203.00 248.0 : 225 8 | 248.00 293.0 : 270 9 | 293.00 338.0 : 315 10 | 338.00 360.0 : 0 -------------------------------------------------------------------------------- /remaps/covdensum.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 2 | 11 : 0 3 | 12 : 0 4 | 13 : 50 5 | 14 : 60 6 | 15 : 55 7 | 16 : 45 8 | 17 : 35 9 | 22 : 10 10 | 23 : 5 11 | 24 : 0 12 | 25 : 0 13 | 31 : 0 14 | 32 : 0 15 | 60 : 80 16 | 61 : 80 17 | 63 : 80 18 | 64 : 80 19 | 65 : 80 20 | 66 : 80 21 | 67 : 80 22 | 68 : 80 23 | 69 : 70 24 | 75 : 60 25 | 76 : 60 26 | 81 : 80 27 | 82 : 80 28 | 95 : 80 29 | 100 : 10 30 | 101 : 15 31 | 102 : 25 32 | 103 : 35 33 | 104 : 45 34 | 105 : 55 35 | 106 : 65 36 | 107 : 75 37 | 108 : 85 38 | 109 : 95 39 | 111 : 15 40 | 112 : 25 41 | 113 : 35 42 | 114 : 45 43 | 115 : 55 44 | 116 : 65 45 | 117 : 75 46 | 118 : 85 47 | 119 : 95 48 | 121 : 15 49 | 122 : 25 50 | 123 : 35 51 | 124 : 45 52 | 125 : 55 53 | 126 : 65 54 | 127 : 75 55 | 128 : 85 56 | 129 : 95 57 | # 0 - No data 58 | # 11 - Open Water 59 | # 12 - Snow/Ice 60 | # 13 - Developed-Upland Deciduous Forest 61 | # 14 - Developed-Upland Evergreen Forest 62 | # 15 - Developed-Upland Mixed Forest 63 | # 16 - Developed-Upland Herbaceous 64 | # 17 - Developed-Upland Shrubland 65 | # 22 - Developed - Low Intensity 66 | # 23 - Developed - Medium Intensity 67 | # 24 - Developed - High Intensity 68 | # 25 - Developed-Roads 69 | # 31 - Barren 70 | # 32 - Quarries-Strip Mines-Gravel Pits 71 | # 60 - NASS-Orchard 72 | # 61 - NASS-Vineyard 73 | # 63 - NASS-Row Crop-Close Grown Crop 74 | # 64 - NASS-Row Crop 75 | # 65 - NASS-Close Grown Crop 76 | # 66 - NASS-Pasture and Hayland 77 | # 67 - NASS-Pasture and Hayland 78 | # 68 - NASS-Wheat 79 | # 69 - NASS-Aquaculture 80 | # 75 - Herbaceous Semi-dry 81 | # 76 - Herbaceous Semi-wet 82 | # 81 - Pasture/Hay 83 | # 82 - Cultivated Crops 84 | # 95 - Herbaceous Wetlands 85 | # 100 - Sparse Vegetation Canopy 86 | # 101 - Tree Cover >= 10 and < 20% 87 | # 102 - Tree Cover >= 20 and < 30% 88 | # 103 - Tree Cover >= 30 and < 40% 89 | # 104 - Tree Cover >= 40 and < 50% 90 | # 105 - Tree Cover >= 50 and < 60% 91 | # 106 - Tree Cover >= 60 and < 70% 92 | # 107 - Tree Cover >= 70 and < 80% 93 | # 108 - Tree Cover >= 80 and < 90% 94 | # 109 - Tree Cover >= 90 and <= 100% 95 | # 111 - Shrub Cover >= 10 and < 20% 96 | # 112 - Shrub Cover >= 20 and < 30% 97 | # 113 - Shrub Cover >= 30 and < 40% 98 | # 114 - Shrub Cover >= 40 and < 50% 99 | # 115 - Shrub Cover >= 50 and < 60% 100 | # 116 - Shrub Cover >= 60 and < 70% 101 | # 117 - Shrub Cover >= 70 and < 80% 102 | # 118 - Shrub Cover >= 80 and < 90% 103 | # 119 - Shrub Cover >= 90 and <= 100% 104 | # 121 - Herb Cover >= 10 and < 20% 105 | # 122 - Herb Cover >= 20 and < 30% 106 | # 123 - Herb Cover >= 30 and < 40% 107 | # 124 - Herb Cover >= 40 and < 50% 108 | # 125 - Herb Cover >= 50 and < 60% 109 | # 126 - Herb Cover >= 60 and < 70% 110 | # 127 - Herb Cover >= 70 and < 80% 111 | # 128 - Herb Cover >= 80 and < 90% 112 | # 129 - Herb Cover >= 90 and <= 100% -------------------------------------------------------------------------------- /remaps/covdenwin.rmp: -------------------------------------------------------------------------------- 1 | 0 : 00 2 | 1 : 60 3 | 2 : 40 4 | 3 : 50 5 | 4 : 100 6 | # 0 - Bare 7 | # 1 - Grass 8 | # 2 - Shrub 9 | # 3 - Deciduous Trees 10 | # 4 - Coniferous Trees -------------------------------------------------------------------------------- /remaps/covtype_to_snow_intcp.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 2 | 1 : 0 3 | 2 : 2 4 | 3 : 3 5 | 4 : 5 6 | # This remap uses coverage type to designate summer rain 7 | # interception values (inches) based on the GIS Weasel user's manual 8 | # These values need to be scaled by 0.01 for use in PRMS 9 | # 0 - Bare 10 | # 1 - Grass 11 | # 2 - Shrub 12 | # 3 - Tree (deciduous) 13 | # 4 - Tree (evergreen) -------------------------------------------------------------------------------- /remaps/covtype_to_srain_intcp.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 2 | 1 : 0 3 | 2 : 2 4 | 3 : 3 5 | 4 : 5 6 | # This remap uses coverage type to designate summer rain 7 | # interception values (inches) based on the GIS Weasel user's manual 8 | # These values need to be scaled by 0.01 for use in PRMS 9 | # 0 - Bare 10 | # 1 - Grass 11 | # 2 - Shrub 12 | # 3 - Tree (deciduous) 13 | # 4 - Tree (evergreen) -------------------------------------------------------------------------------- /remaps/covtype_to_wrain_intcp.rmp: -------------------------------------------------------------------------------- 1 | 0 : 0 2 | 1 : 0 3 | 2 : 2 4 | 3 : 3 5 | 4 : 5 -------------------------------------------------------------------------------- /remaps/temp_adj_x10.rmp: -------------------------------------------------------------------------------- 1 | -999 : 0 2 | -1 : 0 3 | 0 44 : -18 4 | 45 89 : -10 5 | 90 134 : 0 6 | 135 179 : 10 7 | 180 224 : 12 8 | 225 269 : 18 9 | 270 314 : 0 10 | 315 360 : -10 11 | # -1 - Flat 12 | # 0 - North 13 | # 45 - North-East 14 | # 90 - East 15 | # 135 - South-East 16 | # 180 - South 17 | # 225 - South-West 18 | # 270 - West 19 | # 315 - North-West -------------------------------------------------------------------------------- /scripts/field_list.ini: -------------------------------------------------------------------------------- 1 | ## GSFLOW Parameter Shapefile Fields 2 | [FIELDS] 3 | 4 | ## HRU Fields 5 | id_field = HRU_ID 6 | type_field = HRU_TYPE 7 | area_field = HRU_AREA 8 | row_field = HRU_ROW 9 | col_field = HRU_COL 10 | x_field = HRU_X 11 | y_field = HRU_Y 12 | lat_field = HRU_LAT 13 | lon_field = HRU_LON 14 | hru_sink_field = HRU_SINK 15 | 16 | ## DEM Fields 17 | dem_mean_field = DEM_MEAN 18 | dem_min_field = DEM_MIN 19 | dem_max_field = DEM_MAX 20 | dem_adj_field = DEM_ADJ 21 | dem_flowacc_field = DEM_FLOWAC 22 | dem_aspect_field = DEM_ASPECT 23 | dem_slope_deg_field = DEM_SLP_D 24 | dem_slope_rad_field = DEM_SLP_R 25 | dem_slope_pct_field = DEM_SLP_P 26 | dem_sink_field = DEM_SINK 27 | crt_elev_field = CRT_ELEV 28 | crt_fill_field = CRT_FILL 29 | 30 | ## Flow accumulation weighted DEM fields 31 | dem_sum_field = DEM_SUM 32 | dem_count_field = DEM_COUNT 33 | dem_flowacc_field = DEM_FLOWAC 34 | 35 | ## Lake Fields 36 | lake_id_field = LAKE_ID 37 | lake_area_field = LAKE_AREA 38 | 39 | ## Precipitation Zone Fields 40 | ppt_zone_id_field = PPT_ZONE 41 | hru_psta_field = HRU_PSTA 42 | 43 | ## Temperature Zone Fields 44 | temp_zone_id_field = TEMP_ZONE 45 | hru_tsta_field = HRU_TSTA 46 | 47 | ## DEM Based Fields 48 | topo_index_field = TOPO_INDEX 49 | deplcrv_field = DEPLCRV 50 | jh_tmin_field = JH_TMIN 51 | jh_tmax_field = JH_TMAX 52 | jh_coef_field = JH_COEF 53 | snarea_thresh_field = SNAREA_THR 54 | tmax_adj_field = TMAX_ADJ 55 | tmin_adj_field = TMIN_ADJ 56 | 57 | ## Impervious Cover Fields 58 | imperv_pct_field = IMPERV_PCT 59 | carea_max_field = CAREA_MAX 60 | 61 | ## Vegetation Fields 62 | cov_type_field = COV_TYPE 63 | covden_sum_field = COVDEN_SUM 64 | covden_win_field = COVDEN_WIN 65 | rad_trncf_field = RAD_TRNCF 66 | snow_intcp_field = SNOW_INTCP 67 | srain_intcp_field = SRAIN_INTC 68 | wrain_intcp_field = WRAIN_INTC 69 | 70 | ## Soil Fields 71 | awc_field = AWC 72 | clay_pct_field = CLAY_PCT 73 | sand_pct_field = SAND_PCT 74 | ksat_field = KSAT 75 | soil_type_field = SOIL_TYPE 76 | soil_root_max_field = SOIL_RT_MX 77 | moist_init_field = MOIST_INIT 78 | moist_max_field = MOIST_MAX 79 | rechr_init_field = RECHR_INIT 80 | rechr_max_field = RECHR_MAX 81 | ssr2gw_rate_field = SSR2G_RATE 82 | ssr2gw_k_field = SSR2G_KFAC 83 | pref_flow_den_field = PREF_FLOW_DN 84 | slowcoef_lin_field = SLOWCOEF_L 85 | slowcoef_sq_field = SLOWCOEF_S 86 | 87 | ## Stream Fields 88 | irunbound_field = IRUNBOUND 89 | iseg_field = ISEG 90 | flow_dir_field = FLOW_DIR 91 | krch_field = KRCH 92 | irch_field = IRCH 93 | jrch_field = JRCH 94 | reach_field = IREACH 95 | rchlen_field = RCHLEN 96 | maxreach_field = MAXREACH 97 | outseg_field = OUTSEG 98 | iupseg_field = IUPSEG 99 | strm_top_field = STRM_TOP 100 | strm_slope_field = STRM_SLOPE 101 | subbasin_field = SUB_BASIN 102 | segbasin_field = SEG_BASIN 103 | outflow_field = OUTFLOWHRU 104 | -------------------------------------------------------------------------------- /scripts/template_parameters.ini: -------------------------------------------------------------------------------- 1 | # GSFLOW Parameter Generator Inputs 2 | [INPUTS] 3 | 4 | # Save temporary rasters & shapefiles in parameter folder 5 | parameter_folder = .\hru_params 6 | 7 | ## HRU Fishnet Parameters 8 | # Fishnet cell size 9 | hru_cellsize = 300 10 | # Main shapefile with all parameters 11 | hru_fishnet_path = .\hru_params\hru_params.shp 12 | # Point shapefile that is used for getting raster values from the rasters 13 | hru_centroid_path = .\hru_params\hru_params_label.shp 14 | 15 | # Parameters needed for generating a new fishnet 16 | # ref_x and ref_y can be used to snap or align the fishnet to an existing point 17 | # If they are not, the lower left corner of the study area will be used 18 | hru_ref_x = 0 19 | hru_ref_y = 0 20 | # Method for adjusting the extent to the snap point: EXPAND, SHRINK, ROUND 21 | hru_param_snap_method = EXPAND 22 | # Number of extra cells to include around study area 23 | hru_buffer_cells = 2 24 | 25 | # Study Area 26 | study_area_path = .\shapefiles\watershed.shp 27 | 28 | # All model outlet, subbasin, and swale points are defined in a single shapefile 29 | # Value in "type" field must be: OUTLET, SUBBASIN, or SWALE 30 | # Subbasin numbering will be defined using value in "zone" field 31 | model_points_path = .\shapefiles\model_points.shp 32 | model_points_zone_field = FID 33 | model_points_type_field = TYPE 34 | 35 | # Lake Parameters 36 | set_lake_flag = False 37 | lake_path = .\shapefiles\lakes.shp 38 | lake_zone_field = FID 39 | lake_area_pct = 40 40 | 41 | 42 | ## DEM Parameters 43 | # Generate all DEM related parameters (mean, min, max, slope, aspect, dem_adj, dem_flowac) 44 | dem_orig_path = .\dem\ned10m_nad83.img 45 | dem_units = meters 46 | # Resampling method: BILINEAR, CUBIC, NEAREST 47 | dem_projection_method = BILINEAR 48 | dem_cellsize = 10 49 | # Calculate topographic index 50 | calc_topo_index_flag = False 51 | # Calculate flow accumulation weighted elevation 52 | calc_flow_acc_dem_flag = True 53 | # flow_acc_dem_factor = 0.001 54 | # Field to initially set DEM_ADJ 55 | dem_adj_copy_field = DEM_FLOWAC 56 | reset_dem_adj_flag = False 57 | # Round DEM_ADJ to the given number of decimals 58 | dem_adj_decimals = 2 59 | 60 | 61 | ## DEM 2 Streams 62 | # Cells with flow accumulations >= threshold will be designated as stream cells 63 | flow_acc_threshold = 30 64 | # All 1st order streams with a length below threshold will be removed 65 | flow_length_threshold = 3 66 | 67 | calc_flow_dir_points_flag = True 68 | 69 | # Inactive water (ocean, not currently supported) 70 | # set_inactive_water_flag = False 71 | # inactive_water_path = .\shapefiles\ocean_inactive.shp 72 | 73 | 74 | ## CRT Parameters 75 | crt_exe_path = ..\..\crt\CRT1.4_beta.exe 76 | crt_hruflg = 0 77 | crt_flowflg = 3 78 | crt_dpit = 0.01 79 | crt_outitmax = 100000 80 | 81 | 82 | ## CRT Fill Parameters 83 | use_crt_fill_flag = False 84 | 85 | 86 | ## Vegetation Parameters 87 | # Vegetation Type 88 | # Assume NEAREST Resampling 89 | veg_type_orig_path = .\landfire\us_130evt.img 90 | veg_type_field = VALUE 91 | veg_type_cellsize = 10 92 | 93 | # Vegetation Cover 94 | # Assume NEAREST Resampling 95 | veg_cover_orig_path = .\landfire\us_130evc.img 96 | veg_cover_cellsize = 10 97 | 98 | 99 | ## Soil Parameters 100 | soil_orig_folder = .\soils 101 | ksat_name = ksat.img 102 | awc_name = awc.img 103 | clay_pct_name = clay.img 104 | sand_pct_name = sand.img 105 | soil_cellsize = 50 106 | # Soil rasters are percent (i.e. 25%) 107 | # If false interpret as decimals (i.e. 0.25) 108 | soil_pct_flag = True 109 | # Fill nodata in soil rasters using nibble technique 110 | fill_soil_nodata_flag = True 111 | # Initial soil moisture and recharge are calculated as a fraction of the max 112 | moist_init_ratio = 0.1 113 | rechr_init_ratio = 0.1 114 | ssr2gw_k_default = 0.001 115 | # Read and apply soil depth raster 116 | # Otherwise soil depth will only be derived from rooting depth 117 | soil_depth_flag = False 118 | # Soil depth will only be read if soil_depth_flag = True 119 | # Soil depth must be in inches 120 | soil_depth_name = soil_depth.img 121 | 122 | 123 | ## Impervious Cover Parameters 124 | impervious_orig_path = .\impervious\nlcd2011_imp.img 125 | # Resampling method: BILINEAR, CUBIC, NEAREST 126 | impervious_projection_method = BILINEAR 127 | impervious_cellsize = 10 128 | # Impervious cover rasters are percent (i.e. 25%) 129 | # If false interpret as decimals (i.e. 0.25) 130 | impervious_pct_flag = True 131 | 132 | 133 | ## Remap Files 134 | remap_folder = ..\..\remaps 135 | aspect_remap = aspect.rmp 136 | # Output values are 10 * value 137 | temp_adj_remap = temp_adj_x10.rmp 138 | # Output values are floats 139 | cov_type_remap = covtype.rmp 140 | covden_sum_remap = covdensum.rmp 141 | covden_win_remap = covdenwin.rmp 142 | snow_intcp_remap = covtype_to_snow_intcp.rmp 143 | srain_intcp_remap = covtype_to_srain_intcp.rmp 144 | wrain_intcp_remap = covtype_to_wrain_intcp.rmp 145 | root_depth_remap = rtdepth.rmp 146 | 147 | 148 | ## PRISM Parameters 149 | prism_folder = .\prism_800m_normals 150 | # prism_folder = .\prism_4km_normals 151 | # Resampling method: BILINEAR, CUBIC, NEAREST 152 | prism_projection_method = BILINEAR 153 | # Output projected cellsize, not PRISM input cellsize 154 | prism_cellsize = 300 155 | # Recalculate JH coefficient with PRISM temperature values 156 | calc_prism_jh_coef_flag = True 157 | 158 | # PPT Ratios 159 | set_ppt_zones_flag = False 160 | ppt_obs_units = mm 161 | 162 | # If set_ppt_zones_flag is false, mean monthly ppt must be set manually 163 | ppt_obs_list = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 164 | # If a valid HRU_ID is set, PPT ratios will be adjusted to be 1 at that HRU_ID 165 | ppt_hru_id = 0 166 | 167 | # If set_ppt_zones_flag is True, mean monthly ppt will be read from a PPT zones shapefile 168 | # PPT ratios will be adjusted to be 1 at each zone HRU_ID (set in ppt_hru_id_field) 169 | ppt_zone_path = .\shapefiles\ppt_zones.shp 170 | ppt_zone_id_field = PPT_ZONE 171 | ppt_hru_id_field = PPT_HRU_ID 172 | ppt_obs_field_format = PPT_{:02d} 173 | 174 | 175 | ## PRMS Parameters 176 | prms_parameter_folder = .\hru_params 177 | 178 | single_param_file_flag = False 179 | # single_param_file_name = prms_inputs.param 180 | param_column_flag = False 181 | 182 | # Default/template values 183 | prms_dimen_csv_path = .\hru_params\prms_dimensions.csv 184 | prms_param_csv_path = .\hru_params\prms_parameters.csv 185 | 186 | # Parameters passed through to PRMS param file(s) 187 | # Number of air temperature measurement stations in the data file 188 | ntemp = 2 189 | # Elevation units (0=feet, 1=meter) 190 | elev_units = 0 191 | # Index of the temperature station used to compute basin temperature values (1's based) 192 | basin_tsta = 1 193 | # Index of the base temperature station used for lapse rate calculations (1's based) 194 | hru_tsta = 1 195 | # Index of the lapse temperature station used for lapse rate calculations (1's based) 196 | hru_tlaps = 2 197 | # Elevation of each air temperature measurment station 198 | tsta_elev = 0, 0 199 | 200 | # If using temp_1sta module, manually define monthly lapse rates 201 | # tmax_lapse = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 202 | # tmin_lapse = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 203 | 204 | 205 | ## Miscellaneous 206 | # Scratch Workspace ('in_memory' or 'scratch') 207 | # in_memory may not work for computers with limited RAM 208 | # scratch will write all intermediate files to a scratch folder on the disk 209 | scratch_name = in_memory 210 | # scratch_name = scratch 211 | 212 | # Scale floating point values before converting to Int and calculating Median 213 | int_factor = 1 214 | 215 | # Don't change this unless it conflicts with an existing field 216 | orig_fid_field = ORIG_FID 217 | -------------------------------------------------------------------------------- /scripts/template_prms_dimensions.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,SIZE 2 | ncascade,calibration_parameters,calculated 3 | ncascdgw,calibration_parameters,calculated 4 | ndays,calibration_parameters,366 5 | ndepl,calibration_parameters,2 6 | ndeplval,calibration_parameters,22 7 | nevap,calibration_parameters,0 8 | ngw,calibration_parameters,calculated 9 | ngwcell,calibration_parameters,calculated 10 | nhru,calibration_parameters,calculated 11 | nhrucell,calibration_parameters,calculated 12 | nmonths,calibration_parameters,12 13 | nobs,calibration_parameters,2 14 | nrain,calibration_parameters,1 15 | nreach,calibration_parameters,calculated 16 | nsegment,calibration_parameters,calculated 17 | nlake,calibration_parameters,calculated 18 | nsnow,calibration_parameters,0 19 | nsol,calibration_parameters,0 20 | nssr,calibration_parameters,calculated 21 | nsub,calibration_parameters,calculated 22 | ntemp,calibration_parameters,config_file 23 | one,calibration_parameters,1 24 | nlake_hrus,calibration_parameters,calculated 25 | -------------------------------------------------------------------------------- /scripts/template_prms_parameters.csv: -------------------------------------------------------------------------------- 1 | NAME,PARAM_FILE,DIMENSION_NAMES,TYPE,DEFAULT_VALUE,,,,,,,,,,,,,,,,,,,,, 2 | adjmix_rain,calibration_parameters,nmonths,2,1,1,1,0.01,0.2,1.2,1,1,1,1.1,0.5,1.2,,,,,,,,,, 3 | albset_rna,default_values,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 4 | albset_rnm,default_values,one,2,0.6,,,,,,,,,,,,,,,,,,,,, 5 | albset_sna,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 6 | albset_snm,default_values,one,2,0.2,,,,,,,,,,,,,,,,,,,,, 7 | basin_tsta,calibration_parameters,one,1,config_file,,,,,,,,,,,,,,,,,,,,, 8 | carea_max,calibration_parameters,nhru,2,CAREA_MAX,,,,,,,,,,,,,,,,,,,,, 9 | cascade_flg,cascade,one,1,0,,,,,,,,,,,,,,,,,,,,, 10 | cascade_tol,cascade,one,2,5,,,,,,,,,,,,,,,,,,,,, 11 | cecn_coef,default_values,nmonths,2,5,,,,,,,,,,,,,,,,,,,,, 12 | circle_switch,cascade,one,1,0,,,,,,,,,,,,,,,,,,,,, 13 | cov_type,gis_derived_parameters,nhru,1,COV_TYPE,,,,,,,,,,,,,,,,,,,,, 14 | covden_sum,gis_derived_parameters,nhru,2,COVDEN_SUM,,,,,,,,,,,,,,,,,,,,, 15 | covden_win,gis_derived_parameters,nhru,2,COVDEN_WIN,,,,,,,,,,,,,,,,,,,,, 16 | dday_intcp,calibration_parameters,nmonths,2,-9,-9,-9,-7.8,-10,-22,-38,-36,-15,-15,-15,-10,,,,,,,,,, 17 | dday_slope,calibration_parameters,nmonths,2,3.10E-01,3.10E-01,3.90E-01,3.40E-01,4.00E-01,5.40E-01,6.50E-01,6.50E-01,4.00E-01,3.40E-01,3.90E-01,3.10E-01,,,,,,,,,, 18 | den_init,default_values,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 19 | den_max,default_values,one,2,0.6,,,,,,,,,,,,,,,,,,,,, 20 | elev_units,default_values,one,1,config_file,,,,,,,,,,,,,,,,,,,,, 21 | emis_noppt,default_values,one,2,0.757,,,,,,,,,,,,,,,,,,,,, 22 | epan_coef,default_values,nmonths,2,1,,,,,,,,,,,,,,,,,,,,, 23 | fastcoef_lin,calibration_parameters,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 24 | fastcoef_sq,calibration_parameters,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 25 | freeh2o_cap,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 26 | gwflow_coef,calibration_parameters,one,2,0.014,,,,,,,,,,,,,,,,,,,,, 27 | gwsink_coef,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 28 | gwstor_init,default_values,one,2,2,,,,,,,,,,,,,,,,,,,,, 29 | gwstor_min,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 30 | gw_up_id,cascade,ncascdgw,1,crt_file,,,,,,,,,,,,,,,,,,,,, 31 | gw_down_id,cascade,ncascdgw,1,crt_file,,,,,,,,,,,,,,,,,,,,, 32 | hru_area,gis_derived_parameters,nhru,2,HRU_AREA,,,,,,,,,,,,,,,,,,,,, 33 | hru_aspect,gis_derived_parameters,nhru,2,DEM_ASPECT,,,,,,,,,,,,,,,,,,,,, 34 | hru_deplcrv,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 35 | hru_down_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 36 | hru_elev,gis_derived_parameters,nhru,2,DEM_ADJ,,,,,,,,,,,,,,,,,,,,, 37 | hru_lat,gis_derived_parameters,nhru,2,HRU_LAT,,,,,,,,,,,,,,,,,,,,, 38 | hru_pct_up,cascade,ncascade,2,crt_file,,,,,,,,,,,,,,,,,,,,, 39 | hru_percent_imperv,calibration_parameters,nhru,2,IMPERV_PCT,,,,,,,,,,,,,,,,,,,,, 40 | hru_psta,gis_derived_parameters,nhru,1,HRU_PSTA,,,,,,,,,,,,,,,,,,,,, 41 | hru_slope,gis_derived_parameters,nhru,2,DEM_SLP_P,,,,,,,,,,,,,,,,,,,,, 42 | hru_strmseg_down_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 43 | hru_subbasin,calibration_parameters,nhru,1,SUB_BASIN,,,,,,,,,,,,,,,,,,,,, 44 | hru_tlaps,calibration_parameters,nhru,1,config_file,,,,,,,,,,,,,,,,,,,,, 45 | hru_tsta,calibration_parameters,nhru,1,config_file,,,,,,,,,,,,,,,,,,,,, 46 | hru_type,gis_derived_parameters,nhru,1,HRU_TYPE,,,,,,,,,,,,,,,,,,,,, 47 | hru_up_id,cascade,ncascade,1,crt_file,,,,,,,,,,,,,,,,,,,,, 48 | imperv_stor_max,default_values,one,2,0.05,,,,,,,,,,,,,,,,,,,,, 49 | jh_coef,calibration_parameters,nmonths,2,0.016301,0.016301,0.016301,0.016301,0.011559,0.011559,0.010966,0.010225,0.00904,0.009633,0.014819,0.014819,,,,,,,,,, 50 | jh_coef_hru,gis_derived_parameters,nhru,2,JH_COEF,,,,,,,,,,,,,,,,,,,,, 51 | lake_hru_id,gis_derived_parameters,nhru,1,LAKE_ID,,,,,,,,,,,,,,,,,,,,, 52 | mapvars_freq,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 53 | mapvars_units,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 54 | max_missing,default_values,one,1,3,,,,,,,,,,,,,,,,,,,,, 55 | melt_force,default_values,one,1,140,,,,,,,,,,,,,,,,,,,,, 56 | melt_look,default_values,one,1,90,,,,,,,,,,,,,,,,,,,,, 57 | outlet_sta,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 58 | potet_sublim,default_values,one,2,0.5,,,,,,,,,,,,,,,,,,,,, 59 | ppt_rad_adj,default_values,nmonths,2,0.02,,,,,,,,,,,,,,,,,,,,, 60 | precip_units,default_values,one,1,0,,,,,,,,,,,,,,,,,,,,, 61 | pref_flow_den,calibration_parameters,one,2,0.2,,,,,,,,,,,,,,,,,,,,, 62 | print_freq,default_values,one,1,3,,,,,,,,,,,,,,,,,,,,, 63 | print_type,default_values,one,1,1,,,,,,,,,,,,,,,,,,,,, 64 | prms_warmup,default_values,one,1,2,,,,,,,,,,,,,,,,,,,,, 65 | rad_trncf,gis_derived_parameters,nhru,2,RAD_TRNCF,,,,,,,,,,,,,,,,,,,,, 66 | radadj_intcp,default_values,one,2,1,,,,,,,,,,,,,,,,,,,,, 67 | radadj_slope,default_values,one,2,0.02,,,,,,,,,,,,,,,,,,,,, 68 | radj_sppt,default_values,one,2,0.44,,,,,,,,,,,,,,,,,,,,, 69 | radj_wppt,default_values,one,2,0.5,,,,,,,,,,,,,,,,,,,,, 70 | radmax,default_values,one,2,0.8,,,,,,,,,,,,,,,,,,,,, 71 | rain_adj,calibration_parameters,nhru;nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 72 | runoff_units,default_values,one,1,0,,,,,,,,,,,,,,,,,,,,, 73 | sat_threshold,calibration_parameters,nhru,2,MOIST_MAX,,,,,,,,,,,,,,,,,,,,, 74 | settle_const,default_values,one,2,0.1,,,,,,,,,,,,,,,,,,,,, 75 | slowcoef_lin,calibration_parameters,nhru,2,SLOWCOEF_L,,,,,,,,,,,,,,,,,,,,, 76 | slowcoef_sq,calibration_parameters,nhru,2,SLOWCOEF_S,,,,,,,,,,,,,,,,,,,,, 77 | smidx_coef,calibration_parameters,one,2,0.01,,,,,,,,,,,,,,,,,,,,, 78 | smidx_exp,default_values,one,2,0.3,,,,,,,,,,,,,,,,,,,,, 79 | snarea_curve,default_values,ndeplval,2,5.00E-02,2.40E-01,4.00E-01,5.30E-01,6.50E-01,7.50E-01,8.20E-01,8.80E-01,9.30E-01,9.90E-01,1.00E+00,5.00E-02,2.50E-01,4.00E-01,4.80E-01,5.40E-01,5.80E-01,6.10E-01,6.40E-01,6.60E-01,6.80E-01,7.00E-01 80 | snarea_thresh,default_values,nhru,2,SNAREA_THR,,,,,,,,,,,,,,,,,,,,, 81 | snow_adj,calibration_parameters,nhru;nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 82 | snow_intcp,gis_derived_parameters,nhru,2,SNOW_INTCP,,,,,,,,,,,,,,,,,,,,, 83 | snowinfil_max,default_values,one,2,2,,,,,,,,,,,,,,,,,,,,, 84 | soil2gw_max,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 85 | soil_moist_init,gis_derived_parameters,nhru,2,MOIST_INIT,,,,,,,,,,,,,,,,,,,,, 86 | soil_moist_max,calibration_parameters,nhru,2,MOIST_MAX,,,,,,,,,,,,,,,,,,,,, 87 | soil_rechr_init,gis_derived_parameters,nhru,2,RECHR_INIT,,,,,,,,,,,,,,,,,,,,, 88 | soil_rechr_max,calibration_parameters,nhru,2,RECHR_MAX,,,,,,,,,,,,,,,,,,,,, 89 | soil_type,gis_derived_parameters,nhru,1,SOIL_TYPE,,,,,,,,,,,,,,,,,,,,, 90 | srain_intcp,gis_derived_parameters,nhru,2,SRAIN_INTC,,,,,,,,,,,,,,,,,,,,, 91 | ssr2gw_exp,default_values,one,2,1.2,,,,,,,,,,,,,,,,,,,,, 92 | ssr2gw_rate,calibration_parameters,nssr,2,SSR2G_RATE,,,,,,,,,,,,,,,,,,,,, 93 | ssstor_init,default_values,one,2,0,,,,,,,,,,,,,,,,,,,,, 94 | subbasin_down,calibration_parameters,nsub,1,calculated,,,,,,,,,,,,,,,,,,,,, 95 | temp_units,default_values,one,1,0,,,,,,,,,,,,,,,,,,,,, 96 | tmax_adj,gis_derived_parameters,nhru,2,TMAX_ADJ,,,,,,,,,,,,,,,,,,,,, 97 | tmax_allrain,calibration_parameters,nmonths,2,6.00E+01,6.00E+01,6.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,7.00E+01,6.00E+01,6.00E+01,6.00E+01,,,,,,,,,, 98 | tmax_allsnow,calibration_parameters,one,2,32,,,,,,,,,,,,,,,,,,,,, 99 | tmax_index,gis_derived_parameters,nmonths,2,calculated,,,,,,,,,,,,,,,,,,,,, 100 | tmin_adj,gis_derived_parameters,nhru,2,TMIN_ADJ,,,,,,,,,,,,,,,,,,,,, 101 | transp_beg,default_values,one,1,4,,,,,,,,,,,,,,,,,,,,, 102 | transp_end,default_values,one,1,10,,,,,,,,,,,,,,,,,,,,, 103 | transp_tmax,default_values,one,2,500,,,,,,,,,,,,,,,,,,,,, 104 | tsta_elev,calibration_parameters,ntemp,2,config_file,,,,,,,,,,,,,,,,,,,,, 105 | tstorm_mo,default_values,nmonths,1,0,,,,,,,,,,,,,,,,,,,,, 106 | wrain_intcp,gis_derived_parameters,nhru,2,WRAIN_INTC,,,,,,,,,,,,,,,,,,,,, 107 | --------------------------------------------------------------------------------