GMS is a pure groundwater modeling software which technically based on the transient numerical solution of 3D Laplase's equation. Mathlab is a mathematical software for the solution of complex mathematical forms.

Interesting.

Hi Azade,

I suggest you to use Python for the job. Python has built up modules and packages for groundwater modelling (MODFLOW). You can use the FLOPY package which is very convenient. Once you write the model input files (.nam and other files) with GMS then you can load the model and run with FLOPY (python). I am doing the same for my research so far I have had great experiences with FLOPY (using python jupyter) so far. It's easy to learn and has lots of potentials (there are many Flopy examples you can help yourself with). Some links if you're interested:

https://water.usgs.gov/ogw/flopy/

https://github.com/modflowpy/flopy

@Fahad Khan Khadim

Thanks a lot for your suggestion, It is very nice that you are working on such project too.

I will try to earn information about linking GMS in FLOPY(python) from your suggestion links.

But since i need to link MIKE11 with GMS in continuing my work, maybe, do you know that MIKE11 software is compatible with python or not?

Regards,

AZADE

Hi Azade, Unfortunately I have no experience working with MIKE and I don't think there is any Python package to couple Mike with Python. I am not sure though. Thanks!

Thanks,

Good luck in your research dear Fahad!

Hello Fahad

I have a question, if is possible.

can be implemented optimization approach with a linked GMS - Python program too?

Regards,

Azade

Hi Azade,

You can do parameter optimization using PEST (Python has packages to support PEST). Unfortunately, I am doing my research for about 6/7 months and haven't yet gone that far into optimization. For my model, I am still manually changing parameters (HK, riverbed conductance, interflow) etc. to see how much the RMSE is. Once I have a clear understanding of ranges of HK's to calibrate my model I will run a sensitivity test to see for which HK the RMSE is the lowest. As of now I don't have plans to use PEST (which if you don't know from before might be difficult to master).

Fahad

Azadeh,

All you need for your project is a Modflow executable (which is publicly available) and Python (and its Flopy package).

I would suggest to pause using Matlab and GMS for a while and spend some time with Flopy (using the links provided by Fahad).

One thing you can do. First create the model withe hep of GMS. Then in the installation directory see some source codes are provided. Also you will find some files created where your model is saved. Now you can copy the required source code in the model directory. Then with help of nam file (which is already created by GMS in the model directory) you can run the desire model. The next step is to call this thing in MATLAB environment. use following command:

command='mf2k_h5.exe your_model_name.mfn'

status=dos(command)

Please note that the exe and the associates files should be in the same directory.

Hi,

Recently, the GMS-MODFLOW code has been written and implemented in MATLAB.

Ali Jamalifard

You can consider comment of Subhajit Dey and Vasileios Christelis in this discuss.

Hi

I know this topic is old, but I would like to share my experience.

It is possible to link Modflow and Matlab easily, but to link GMS and Matlab you need to have a huge experience in advanced programming.

I developed the Modflow-Matlab interface in such a way that I take advantage of GMS pre and post-processing capabilities as well as Matlab powerful numerical capabilities.

All you need to learn is how to write the Modflow inputs packages and how to retrieve the Modflow outputs.

I suggest, first, see maflab code. It provided a very good description of the reading/writing of Modflow/MT3DMS models.

DO NOT FORGET to read user manuals of Modflow and MT3DMS.

I suggest seeing C/C++ documents in the reading/writing of files topics.

As mentioned earlier, using Flopy is another way to programmatically control the Modflow model(optimization, sensitivity and so on). both Matlab and Python method has a drawback in parallelism as long as they use the reading/writing of files to interface.

I hope this helpful.