Sylvain 0 Posted February 6, 2018 Report Share Posted February 6, 2018 Hi, I would like to know if there is another way than manually assign transient recharge values from 2d array in the recharge package window. I have 120 raster (10 years of monthly recharge) of spatially distributed recharge to load in GMS. I'm thinking of going into MODFLOW input files. Any ideas? Regards Quote Link to post Share on other sites
Michal 0 Posted February 7, 2018 Report Share Posted February 7, 2018 Hi, I think it could be possible by interpolating transient TIN data (or maybe even rasters) into Recharge package using Interpolate to Modflow layers command. But I have not done it before. Let us know how you manage to handle it. Quote Link to post Share on other sites
Bruce Campbell 2 Posted February 7, 2018 Report Share Posted February 7, 2018 Hi Sylvain -- I've had good results using a 2-D grid. I) create a 2-D grid from the 3-D grid 2) do a linear (or other type) interpolation of your raster recharge data to the to 2-D grid. These will stack up for your 120 rasters under the 2-D grid 3) Go to MODFLOW -> Optional Packages -> Recharge and the use the 2D Dataset - > Array to put the 120 datasets into the model at the appropriate stress period Quote Link to post Share on other sites
Sylvain 0 Posted February 12, 2018 Author Report Share Posted February 12, 2018 Hi, Thanks for the reply. But i was looking for a more automated process than manually import the interpolated data into the Recharge window (2D Dataset -> Array) and assign them to each stress period. Lets say it takes 1 min to import each 2d grid for each stress period and i have 120! Thanks! Quote Link to post Share on other sites
Sylvain 0 Posted February 12, 2018 Author Report Share Posted February 12, 2018 I have found out how do directly modify the recharge values in the raw MODFLOW .h5 file. You can use R package: https://www.bioconductor.org/packages/release/bioc/vignettes/rhdf5/inst/doc/rhdf5.pdf Good luck and thanks for you ideas! Quote Link to post Share on other sites
Michal 0 Posted February 16, 2018 Report Share Posted February 16, 2018 Thank you for sharing the solution. Nice one Have you used some sort of interpolation to obtain modflow compatible dataset from the raster data? Or is the grid and rasters identical in number of rows and columns? Quote Link to post Share on other sites
Alan Lemon 3 Posted March 6, 2018 Report Share Posted March 6, 2018 We will be adding a feature to the next version of GMS to apply rasters to various stress periods in the Areal packages (RCH, EVT, ETS). There will not be a need to import these into your project they can be processed directly from the RCH package. Quote Link to post Share on other sites
Jef D. 0 Posted May 28, 2018 Report Share Posted May 28, 2018 Hi Sylvain, Alan, I also have >100 recharge rasters for different stress periods that I would like to import in my MODFLOW model. The recharge raster grid is identical to my MODFLOW grid. Sylvain, would it be possible to give some insight in your R code to directly write the data to the MODFLOW GMS hdf5 file? I managed to read my current recharge in the .h5 file using the package you mentioned but I don't know how to correctly write new recharge values for multiple stress periods to the .h5 file? I have my recharge rasters ready as R data frames. Alan, is this new feature already available in version 10.3.6? Quote Link to post Share on other sites
Alan Lemon 3 Posted June 5, 2018 Report Share Posted June 5, 2018 The new feature will be in the next version of GMS (most likely 10.4). Quote Link to post Share on other sites
Petra Furychova 0 Posted November 20, 2018 Report Share Posted November 20, 2018 Hi Alan, could I ask if this feature is now included in the recently upgraded version 10.4? Thanks Quote Link to post Share on other sites
Alan Lemon 3 Posted January 11, 2019 Report Share Posted January 11, 2019 Sorry for the late reply. This feature is in GMS 10.4 in the recharge package. There is a button "Rasters -> Array". You choose all the raster files that you would like to apply to your grid and then there is a spreadsheet where you can specify the stress period for each raster. Quote Link to post Share on other sites
Recommended Posts
Join the conversation
You can post now and register later. If you have an account, sign in now to post with your account.