/
KHM_Friedrichshagen_Plot_Cl-Konzentation_Szenarios_kombiniert.py
75 lines (59 loc) · 2.27 KB
/
KHM_Friedrichshagen_Plot_Cl-Konzentation_Szenarios_kombiniert.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
import os
import matplotlib.pyplot as plt
import flopy
ws1 = './Szenario 1 FOK_flach -15.0 FOK_tief -21.0' # Pfad zum Modellordner des Szenarios
ws2 = './Szenario 0 Basismodell' # Pfad zum Modellordner des
sz_nr = 1 # Nummerierung des Szenarios für die Legende des Plots
f1 = os.path.join(ws1, 'transport.nam')
mt1 = flopy.mt3d.mt.Mt3dms.load(f1)
f2 = os.path.join(ws2, 'transport.nam')
mt2 = flopy.mt3d.mt.Mt3dms.load(f1)
fname = os.path.join(ws1, "MT3D001.UCN")
ucnobj = flopy.utils.UcnFile(fname)
times = ucnobj.get_times()
conc = ucnobj.get_alldata()
fname = os.path.join(ws1, "MT3D001.OBS")
if os.path.isfile(fname):
cvt = mt1.load_obs(fname)
else:
cvt = None
fname = os.path.join(ws1, "MT3D001.MAS")
mvt = mt1.load_mas(fname)
fname2 = os.path.join(ws2, "MT3D001.UCN")
ucnobj2 = flopy.utils.UcnFile(fname2)
times2 = ucnobj2.get_times()
conc2 = ucnobj2.get_alldata()
fname2 = os.path.join(ws2, "MT3D001.OBS")
if os.path.isfile(fname2):
cvt2 = mt2.load_obs(fname2)
else:
cvt2 = None
fname2 = os.path.join(ws2, "MT3D001.MAS")
mvt2 = mt2.load_mas(fname2)
plt.rcParams.update({'font.size': 18})
x = cvt["time"]
y = cvt["(3, 37, 25)"] * 1000 + 50
plt.plot(x, y, label="Flache Brunnen Szenario {}".format(sz_nr), color='blue')
x = cvt2["time"]
y = cvt2["(3, 37, 25)"] * 1000 + 50
plt.plot(x, y, label="Flache Brunnen Basismodell", color='blue', linestyle='dashed')
x = cvt["time"]
y = cvt["(5, 25, 25)"] * 1000 + 50
plt.plot(x, y, label="Tiefer Brunnen Szenario {}".format(sz_nr), color='red')
x = cvt2["time"]
y = cvt2["(5, 25, 25)"] * 1000 + 50
plt.plot(x, y, label="Tiefer Brunnen Basismodell", color='red', linestyle='dashed')
plt.xticks(fontsize=18)
plt.yticks(fontsize=18)
# plt.yticks([50, 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000, 1100, 1200, 1300, 1400,
# 1500, 1600, 1700, 1800, 1900, 2000])
plt.xlim(0, 3500)
plt.ylim(0, 2.0 * 1000)
plt.xlabel("ZEIT IN TAGEN", fontsize=18)
plt.ylabel(("Chlorid-Konzentration in mg/L"), fontsize=18)
# plt.rcParams.update({'font.size': 18})
plt.grid(visible=True)
# title1 = "Änderung der Chlorid-Konzentrationen im KHM Friedrichshagen " \
# "\nabhängig von der Tiefe der Filteroberkante (FL) der flachen Brunnen"
plt.legend()
plt.show()