Contact US

Log In

Come Join Us!

Are you a
Computer / IT professional?
Join Tek-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!

*Tek-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Need help in Fortran programming - MODFLOW

Need help in Fortran programming - MODFLOW

Need help in Fortran programming - MODFLOW



I need someone who is familiar with MODFLOW code, to help me about my research project.
For sensitivity analysis in this specific project some kind of changes in MODFLOW code is required; I think this change should be occurred in the SIP module. This method is called adjoint problem. Based on some search, I find the following description:

"The adjoint problem can be readily solved using MODFLOW with only minor changes to the original code. The solution is obtained by (1) running MODFLOW (forward problem) to obtain the A matrix, (2) setting the RHS vector equal to zero, except for the grid cell containing the head of interest, which is set to one, and (3) rerunning MODFLOW using the new RHS to obtain the importance function" (M. I. Jyrkama and J. F. Sykes, 2006)

After obtaining the "importance function", sensitivity analysis would be possible.
Now I need to put RHS equal to zero. Except for the target cell; the value for the target cell would be 1.
Note that RHS means is Right Hand Side of flow equation and this is a vector for finite difference grid... and all of this process is part of the reliability analysis method.
These changes are big challenges for me. Where and how these changes should be done in the module?

I have attached the code module here in several formats. In addition, the theory of SIP computation is attached (download pdf) read only pages 423-453.
In order to run and compile the MODFLOW code, description that has been listed in this link is important. These changes should be applied on the Aquaveo MODFLOW-2005 code.

Hoping to compensate for your kindness

Behzad Sarhadi


Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Tek-Tips Forums free from inappropriate posts.
The Tek-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Tek-Tips forums is a member-only feature.

Click Here to join Tek-Tips and talk with other members! Already a Member? Login

Close Box

Join Tek-Tips® Today!

Join your peers on the Internet's largest technical computer professional community.
It's easy to join and it's free.

Here's Why Members Love Tek-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close