Jump to content
GMS, SMS, and WMS User Forum

Jonas

Members
  • Content count

    26
  • Joined

  • Last visited

Everything posted by Jonas

  1. I have problems running PEST NSMC with fixed pilot points. Both in my own model and in the NSMC tutorial files, I get the same error message "Cannot open parameter value file nsmcII_200_random_nul.par." and the simulation finishes without completing any iteration. Is this a bug i in GMS? Kind regards, Jonas
  2. Fixed pilot points in PEST NSMC

    An update, my own model was possible to run after using the SUBREG1 utility in PEST. A problem is that my own model with a fixed pilot point is only possible to run for log-uniform distributions with the RANDPAR utility. When I try to create log-normal distributions in RANDPAR i get the error message seen in the screenshot below. Ideas on how to solve this issue are appreciated.
  3. Fixed pilot points in PEST NSMC

    I used Kriging for interpolation of the pilot points. Yes, that was the procedure for the tutorial file. When I run my own model something seems to go wrong with svd-assist since the sensitivity file for svd-assist is empty.
  4. PEST, SVD and SVD-assist

    I often experience problems when running a PEST parameter estimation model with SVD-assist. I find it quite arbitrary when it work and when it does not work. For some occasions it seems to work just by moving the modelling files to another computer and run the model on that workstation. For my current model, SVD works but SVD-assist does not. To reach reasonable results with NSMC is even more problematic, in particular for transient models. I run the latest nightly build and NWT as solver. Someone with similar experiences or suggestions to a more streamlined and stable model procedure?
  5. Fixed pilot points in PEST NSMC

    An update on this matter. As suggested by Michal, I was able to run the PARREP and RANDPAR utility in PEST to create the *.par files for the different runs. Before running PARREP, I ran a parameter estimation run in GMS with NOPTMAX = -1 to create a *.bpa file. After running RANDPAR I was able to perform the normal procedure for a NSMC inverse calibration run in GMS. This worked out well for the NSMC tutorial with "sc1v56" as a fixed pilot points (keeping value 400 in screenshot below). I cross my fingers that it also works for my own model. However, I would appreciate if this procedure could work directly in GMS.
  6. PEST, SVD and SVD-assist

    Thank you again Michal. Now I was able to follow your suggested procedure with a fixed pilot point in the NSMC tutorial and it worked out quite good. Please refer to my other thread on this matter.
  7. PEST, SVD and SVD-assist

    Now I try to run the RANDPAR utility in the command prompt but get the error message "Cannot open covariance matrix file hk_12.matrix."
  8. PEST, SVD and SVD-assist

    The exact error message while running PEST with NSMC is "cannot open parameter value file m20s20000_random_nul1.par". If someone has an idea on this error message it would be very appreciated.
  9. PEST, SVD and SVD-assist

    Thanks Michal and Julian for your replies. Eventually I might end up running PEST outside of GMS if it does not work. I managed to get a parameter estimation run with svd-assist to work. I also managed to get a PEST NSMC run to work. But, when I want to repeat the procedure with fixed pilot points I get a message: "Cannot find sm20s20000_random_nul1.par". I cannot find the corresponding file in the folder of the model file neither. As I wrote in my first post I find it quite arbitrary when PEST works together with GMS. I would be happy with a response from someone from GMS's developing team on this issue.
  10. PEST, SVD and SVD-assist

    Thanks for your quick reply Michal, I agree on that it seems like PEST is doing its job but GMS does not read in the result properly. But I can't find the reason for this since the same modelling procedure has worked for other models. Quite frustrating. How do use the PARREP utility? Should I use it outside of GMS interface? Kind regards, Jonas
  11. PEST, SVD and SVD-assist

    This is what typically happens when running with SVD-assist. As can be seen from the blue marked line, PEST is able to reduce the error to 98.885 from the initial error of 155.76. But, then something goes wrong since this solution never comes up as a iteration before the parameter estimation is finished.
  12. I use Kriging and variograms to describe spatial variability of pilot points for HK and RCH in a PEST parameter estimation simulation. I use the log values for both HK and RCH. For one pilot point set it works perfectly fine to use a logarithmic semivariogram, which also looks reasonable. These HK values are in the range between 1e-4 to 1e-6 m/sek. For pilot point sets with lower values (in the range between 1e-8 and 1e-10) I get the following message when I try to build a logarithmic semivariogram: Warning. There were no points calculated; bad lag parameters! A semivariogram (that is not logarithmic) does however allow me to build a variogram. I have two questions: 1. What is the reason for this warning? 2. Would it be a significant difference in interpolation between a semivariogram and a logarithmic semivariogram? I am happy to receive any comments and ideas on these issues. Kind regards, Jonas
  13. Kriging with logarithmic semivariogram

    Hi Julian, I decided to use another variogram method than the semi-logarithmic for the cases when I got the error. If you have problems with the Kriging procedure, have you tried one of the recent nightly builds?
  14. Hi, I experience the tool for defining stress periods as dates/times as quite buggy. Please refer to screenshot below. First, there is a different format between the first row and the second two rows (month as number "10" instead of letters "Oct"). But the most significant issue is that the "length" changes to weird numbers when the date and time is changed. In the example below the first stress period has a negative number but I have also experienced zero as length. My workaround has been to not make changes as date/times by unmarking the "use dates/times" box, calculate the corresponding length between and then change back to "use dates/times". This works but not so convenient. Another issue is that I have seconds as time format in my model but to get the length right I have to divide the seconds by 2 in the "length" box to get the date/times right. Finally, the model terminates with an error when I use date/times but works when I am un-tickning the box mentioned above (not using date/times). Any ideas on how to get these issues sorted out? Kind regards, Jonas
  15. Stress periods defined as dates/times

    Unfortunately I still experience the stress-periods as quite buggy. The images below show stress periods as date/times and as seconds respectively. The result is from a PEST NSMC simulation with a transient pump-test as calibration target. During the initial forward runs and PEST calibration without NSMC I did not experience the same issue with the date-time series. The simulation itself seems to work fine but the problem arises when creating graphs and comparing the simulated results with my observations. Since I think there is a bug behind the issue I would appreciate if someone from the GMS development team could look into it. Kind regards Jonas
  16. Stress periods defined as dates/times

    Thank's Dirk. You advice worked for me as well.
  17. PEST and Kriging

    Hi, When i customize the variogram for the pilot points of the HK-field in the nsmc I and II tutorials, and also in my own models, I encounter the problem showed in the attached figure. It seems like the whole field is generated from only one point and not all available ones. I also edited the already set up variogram in the tutorials, and then changed back to its original form and still got the same problem. Any ideas on how to solve this issue? Kind regards, Jonas
  18. PEST and Kriging

    Thanks Alan, the nightly build solved the issue.
  19. Null Space Monte Carlo

    Hi again! I attempt to use Null Space Monte Carlo in a similar manner as presented in PEST - Null Space Monte Carlo I tutorial (http://gmstutorials-10.1.aquaveo.com/MODFLOW-StochasticModeling-NullSpaceMonteCarloI.pdf) As presented in the attached text file I have some model that converge although most models does not converge for some reason (Ideas on how to adress this are also appreciated). As seen in the text file, there are different calculated values for the pilot points between the four converged solutions. The bpa files for the different solutions are also different. All good so far. My issue is that there is no difference between the grid files of the HK parameter for the four different solutions. Any ideas on how to solve this issue? Could it have something to do with the Kriging process? Best regards, Jonas T40_nsmc.txt
  20. Null Space Monte Carlo

    For others interested in this issue, the support team has now fixed the bug that caused the issue of PEST and null space MC with the NWT solver. Download the latest nightly build (Jan 25th) from: http://www.aquaveo.com/downloads-gms Woodward, it was a pleasure to read your article, "Uncertainty in the modelling of spatial and temporal patterns of shallow groundwater flow paths: The role of geological and hydrological site information" in Journal of Hydrology where you have applied this method.
  21. Null Space Monte Carlo

    Thank you Alan, I have done so. Looking forward for their response. Woodward, the more simplified converges in almost all cases. The problem is that non-unique grids are created.
  22. Null Space Monte Carlo

    Now I have simplified the model a lot to only one layer and use the LPF package. It seems that the only solver that creates unique grid files is 2000. For 2005 and NWT it is the same issue as previous specified. Others with similar experience? Kind regards, Jonas
  23. Null Space Monte Carlo

    I use the NWT solver with UPW in my model. When I changed to the same solver in the NSMC I tutorial, the same thing happens. Different values in the bpa files but non-unique grid files with HK values between the different solutions. Ideas on this issue are appreciated!
  24. Hi! I have set up a model for a municipal water supply in an glacifluvial esker. I want to calibrate the hydraulic conductivity (HK) against steady-state heads in the esker material with PEST and pilot points. I also have other materials in the model but with a fixed HK. The recharge is also calibrated with PEST in three different polygons. The picture below shows layer four (of five). Left shows materials (only bedrock and esker in this layer) together with locations of the pilot points. The middle part shows calibrated heads after PEST has been running. Right shows HK after calibration with PEST. My issue is that the area with esker material (left) does not overlap the area that is calibrated with pilot points in the rightmost picture (non red area). Something that looks strange in particular is the straight paths from the middle to the south of the model that get the same HK as the bedrock although the material is set to esker. It is the same issue in the other layers of the model. Ideas on how to solve this issue are appreciated. Kind regards, Jonas
  25. I created the zones based on solids. For some reason the materials got the right values according to the zones but some cells got wrong values of HK. I changed this manually by selecting these cells. This worked. Could this be a bug in GMS maybe?
×