Looking for help with a school project

42 views
Skip to first unread message

Matt Q

unread,
Jan 31, 2018, 9:45:43 PM1/31/18
to pitlakq-users
Hello,

I am an undergrad student at Lakehead University in Thunder Bay Ontario in Canada. My thesis advisory has given me a project that requires that I use pitlakq because it is free and open source due to the lack of budget available for undergraduate students theses. My thesis project is to model the Steep Rock Iron Mine pit lake with pitlakq using old input data that is available from a previous modeling project that used DYRESM (Dynamic Reservoir Simulation Model). The purpose is to compare and contrast the results of the 2 models and put forward the advantages and disadvantages of them, one being that pitlakq is free and DYRESM though once free is now licensed and must be purchased.

Going into this I had no programing or computer software experience so I have been struggling and relying on friends to assist me as best as they can. I have been trying for a long time with no success at getting pitlakq to run with the data that I have been provided. Initially I thought it could be due to my lack of knowledge in programming but I am starting to wonder if it's due to the data itself that I'm trying to work with. DYRESM used a very different spacial system and file type and it could be that what I thought was converted to something pitlakq could use may not be the case.

The completion of my degree is dependent on completing a thesis project and it must be done by end of April. If it turns out that the data that I have been provided with cannot be ran in pitlakq than I must come up with a new topic soon.

I have attached the files that I have been trying to run in pitlakq. I am happy to provide as much additional information as possible because I am unsure what may be needed in order to help. Any guidance, clarity or assistance would be greatly appreciated and alleviate some stress that I have been having.

Thank you,
Matt
bath_asci.grd
preprocessing.yaml
reservoir.txt

Dr. Mike Mueller

unread,
Feb 5, 2018, 3:24:51 PM2/5/18
to pitlak...@googlegroups.com
Hi Matt,

I can pre-process the files you sent:

setting new values for bathymetry
setting new values for path
setting new values for bathymetry
preprocessing is running ...
total_volume n = 4242.782339752240 Mio m3
total_volume entire lake= 4242.782339752240 Mio m3
set bottom_elevation to: 170.00
surface_area 13941619.1124
runtime was 390.6 seconds

More than 4.000 million m3. This is a pretty big pit lake. ;)

I think the problem is your lake contour.
Looks like my version of Surfer cannot import its own Surfer 6 text format
anymore. I drew the contour map with matplotlib. See the enclosed picture.
You need to have a neat closed contour. Otherwise all the areas outside of
the lake with a level below the max. lake water surface will become part of
the lake. Please work on the contour in Surfer and export only the parts
that are actually part of the lake.

It should look the lake in the tutorial:
http://pitlakq.readthedocs.io/en/latest/_images/tut_lake.png
All the orange area above the maximum water level of 200 is excluded.
Please modify your contour and set all unused areas to levels higher than
411, which is your maximum water surface.

Best,
Mike


Am 01.02.18 um 03:45 schrieb Matt Q:
> --
> You received this message because you are subscribed to the Google Groups
> "pitlakq-users" group.
> To unsubscribe from this group and stop receiving emails from it, send an email
> to pitlakq-user...@googlegroups.com
> <mailto:pitlakq-user...@googlegroups.com>.
> For more options, visit https://groups.google.com/d/optout.

lake_contour.png

Matt Q

unread,
Feb 8, 2018, 3:55:58 PM2/8/18
to pitlakq-users
Hi Mike,

Thank you for your help! I was able to edit the contours of my map data and produced a very promising image with do_show_bath. I am now encountering another issue with the next step when trying to run the project (named 'attempt_four'). I'm not sure what changes I need to make to the input files I'm using. I have used copies of the tutorial file but with the information changed to suite my project, but I have also looked through the template files and am unsure how to use them if necessary for what I'm trying to process.

I have attached what I think would be needed to understand where my error is coming from.

Thank you!
Matt


Error Message Feb 8th.txt
ncdump bath.txt
bath_asci.grd
successful map figure.png
pitlakq.yaml
w2.yaml

Dr. Mike Mueller

unread,
Feb 22, 2018, 6:57:00 AM2/22/18
to pitlak...@googlegroups.com
Hi Matt,

I've been traveling the last weeks. This is took up all my time.
I will have a look at your model and will answer your questions
early next week.

Best,
Mike

Am 08.02.18 um 21:55 schrieb Matt Q:

Matt Q

unread,
Feb 22, 2018, 11:10:51 AM2/22/18
to pitlakq-users
Thank you Mike,

I have been working on it myself but I haven't solved what it is that I'm missing or setting incorrectly. I have added initial concentrations as the minimum detection limits since posting the last message here.

Thank you,
Matt

Dr. Mike Mueller

unread,
Mar 2, 2018, 12:05:57 PM3/2/18
to pitlak...@googlegroups.com
Hi Matt,

Your bathymetry makes no sense.

I generated the enclosed picture with
(see
http://pitlakq.readthedocs.io/en/latest/tut.html#preprocess-the-bathymetry-data
for how do do this):

from pitlakq.postprocessing import show_bath
show_bath.main("output/bath.nc", bottom=360)

As you can see, there are essentially two lakes, i.e. white areas.
This doses not work. There needs be to only one lake body.
Also, your segments have a width of 1 m meter. Typical segments should
have a with of 50 to 200 m (or more if the lake is large).
Furthermore, you have many inactive, i.e. gray cells. Looks like
your bottom value is way to deep.

Once you fixed your bathymetry you need to echo the dimensions
form bath.nc file in your w2.yaml. This would be the relevant parts:

bounds:
number_of_segments:
value: 174
number_of_layers:
value: 52
number_of_constituents_segments:
value: 174
number_of_constituents_layers:
value: 52

bathymetry:
starting_water_level:
value: 360

Best,
Mike

Am 22.02.18 um 17:10 schrieb Matt Q:
bath.png
w2.yaml

Dr. Mike Mueller

unread,
Mar 6, 2018, 7:21:04 AM3/6/18
to pitlak...@googlegroups.com
Hi Matt,

I just realized that I overlooked your successful "map figure.png".
I could not reproduce this with the bathymetry you attached.
Can you send me the bath.nc and the code you use to produce this
figure.

Best,
Mike


Am 02.03.18 um 18:05 schrieb Dr. Mike Mueller:
Reply all
Reply to author
Forward
0 new messages