from MODFLOWtxt import *

#HighPlains Model

#Enter cell sizes and path to working folder and path to output folder
DEM_cellsize = 250.0 #Cellsize of DEM being used for data creation, must be less than or equal to smallest cell size entered below
Recharge_cellsize = 500.0
dem_minus = 3 #Value to subtract from DEM to create initial head values
cellsizes = [5000] #Cellsizes to use for raster creation
workingFolder = "C:\\Tutorial_Data\\From_Andy\\HighPlainsGIS\\" #Folder containing primary GIS files
outputFolder = "C:\\Tutorial_Data\\From_Andy\\Outputs\\"  #Folder to write output to

#Enter the names of your working files
SPYD = "spyd.shp"
Hyd_K = "hyd_k.shp"
Bedrock = "bedrock_elev.shp"
Recharge = "rech_in_yr_ss"
PredevWaterLv = "predevwlv_elev_Project.shp"
HighPlainsRivers = "HighPlainsMajorRiversNHD.shp"
DEM = "dem_1000"
AquiferCoverage = "hp_extents.shp"

#Enter the following folder locations
ToolboxFolder = "C:/Program Files (x86)/ArcGis/Desktop10.0/ArcToolbox/Toolboxes/"
Spatial_Ref_Folder = "C:/Program Files (x86)/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems/UTM/NAD 1983/"


names=[SPYD,Hyd_K,Bedrock,Recharge,PredevWaterLv,HighPlainsRivers,DEM,AquiferCoverage,ToolboxFolder,Spatial_Ref_Folder]

    
for cellsize in cellsizes:
    outDirec = outputFolder + str(cellsize)
    MODFLOWtxt(cellsize,DEM_cellsize,Recharge_cellsize,dem_minus,workingFolder,outDirec,names)