Jump to content
GMS, SMS, and WMS User Forum


  • Content count

  • Joined

  • Last visited

Everything posted by cansudemir

  1. cansudemir

    ERROR: transient modflow forward run

    Hi, I have been using Modflow-usg and made successful runs. However, there are still many bugs in USG since it's still in development. Therefore I wanted to switch to Modflow with the same data set. I have 20336 cells with cell size of 150m in x and y direction, HK and starting head values were assigned to the model from 2-D scatter data- Interpolate to Modflow. When all the conceptual model data is converted to modflow, all the things look fine (as they were before in the modflow-usg, everything is same). Simulation check is not giving any errors. However, when I run the model I continuously receive the same error and MODFLOW terminates with error without doing nothing. I have tried all the solvers. I have decreased the number of cells. Instead of starting heads I also tried checking steady state in stress period dialog. The error is: ERROR READING ARRAY. Error opening data set: Pilot Points/Index1/start vals. Error opening data set: Pilot Points/Index1/weights. MODFLOW terminated with error! Error reading ......hed. File formatted incorrectly. Other output encountered: ---Premature end of file file_name.ccf Could you please help me solving this problem. Thanks! Cansu
  2. cansudemir

    Observation Wells in Transient Model

    Hello, I am having a serious problem with the observations in my transient model. The observation wells have measurement data from 2009, 2010, 2011, 2012 and 2013. Each well has different measurement dates. Some wells have measurement of all years but some of them have less. When I run the model and plotted computed vs observed graph, in the view values part I see observations of different months and years. For example, in the plot of 16.03.2013 there are some observation values of 2009 (a well whose last observation is 2009). What is the problem? How can I erase these wells from the times that they don't belong to? Need help ASAP. Thank you!! Cansu
  3. cansudemir

    Transient PEST Run with Pilot Points

    Hello dear users and the developers, I am having a problem with Transient PEST calibration (GMS version 10.2.3). My model consists of 1 layer UGrid and have boundary conditions of specified head (along the coast, 0) and GHB. I have done several PEST runs with 90 pilot points in steady state condition. It worked and got results. In the calibration window, it was showing the flow of parameters changing, the graph was being drawn. Then, I changed it to 7 years long transient model. The forward transient run with MODFLOW-USG was successfully done. But when it comes to transient PEST calibration with 90 pilot points, there is nothing going on in the window. No parameters showing up, just the time passes and it says parameter estimation finished. But no results. When I close it with "read the solution", it gives an error saying that: boost::filesystem::remove: The process cannot access the file because it is being used by another process: "C:\Users\user\AppData\Local\Temp\GMS_10980\ppest\slave1\GMS FILE NAME.ba6" During the Parallel PEST run, several cmd.exe are working besides pslave.exe. In the cmd.exe, there is something written by the developer: "Note: This program reads a MODFLOW-USG unformatted file. Sometimes there are problems in reading files of this type due to incompatibilities between different FORTRAN compliers. If there are any such problems please contact ... and I will send you an alternative copy of this program." Can this be the problem of mine in my case? Should I contact the person to ask for the alternative copy? Or what do you recommend? Thank you! Best, Cansu
  4. Hello, For a very long time, I have been trying to create a numerical model of an aquifer with complex hydrogeology (by using borehole-to-Ugrid option). I had several layers of 4 different materials. But all my efforts were trash. Because, I have never obtain a successful Modflow-usg run. They never converged. Although wetting option was on, the aquifer had so many dry cells. And these dry cells were always at same locations (not near/at the location of wells). It was meaningless. I couldn't run the PEST. So I decided to make my system more basic and created a one-layer model with HK zones. It worked well without the wetting. I have applied pilot points approach for the calibration. And pest is working because the forward run has a successful result. However, this time the calculated residual is too high. I have done many things to decrease the error. But the best result until now is about 7000 (Sum of Squared Weighted Residual Head+Flow) and for Root mean squared residual it is 16. These values are so big. Interestingly, the computed vs. observed head plot (attached) looks not so bad that there should be a way to fit them to the line. I only have head observations in my hand and I don't have any flow observations data. May this be the problem? Or what may be the problem? Could you please help me on this? I would be so happy to hear your recommendations to lower my calibration error to an acceptable value. Note: The computed heads are higher than the measured heads. The boundary conditions are ; specified head along the coast and GHB condition was assigned to some part of the boundary. The model is steady-state. Thank you! Best, Cansu
  5. cansudemir

    LPF package

    Dear @Alan Lemon, Regarding your above reply; In my case similar to above, when I assign negative values for HK zones, the run results completely wrong. With the negative values I do can create key parameters by using "initialize from model". But the model doesn't run properly. On the other hand, with positive HK values I got a successfully finished run. However in this case, I couldn't initialize the key parameters from the model and don't know how to create parameters by hand. What should I do? Should I just add new parameters? If so, how does the program understands that these parameters indicate my assigned HKs? Best, Cansu
  6. cansudemir

    Tracking particles in Modflow-USG doesn't works

    Hello Dear @Alan Lemon, Has the bug regarding the contour labeling for UGrid been solved? I am still using GMS 10.1.3. and have a strange kind of fear that if I download the updated versions, I won't be able to open the older versions on which my whole study depends. But this contour labeling is a serious problem for my studies right now so if this error was solved in one of the following versions of GMS I will download that. I need your help. Thank you. Best, Cansu
  7. Hi there, I'm having a problem related with assigning specified flow boundary conditions to my UGrid / Modflow-USG model. The error appears when I do the '' Map to MODFLOW'' option while converting my conceptual model to UGrid numerical model. I have read the helping document saying that If you want to use ''Specified flow boundary condition'' the UGrid must be a stacked grid. A stacked grid is one in which there is no vertical discritization of layers and the horizontal discritization of all layers is the same. (i.e. the Index of Vertical Sub-Discretization (IVSD) in the UnStructured Grid Discretization (DISU) package could be set to -1) (See the Figure attached). Firstly, I did not understand what this means and how to find the IVSD in the DISU package and set it to -1 to be able to use specified flow boundary condition. In my model there are 2 groundwater divide boundaries which must be zero flow and 2 formation boundaries where the aquifer feeds from neighboring formations with a specified flow value. If I ll have to change my specified flow bound. cond. with something else or If there is a way to use them, could you please help me to find the solution? Thank You!