forked from danshapero/greenland_inversions
-
Notifications
You must be signed in to change notification settings - Fork 0
/
grid_dependence.py
49 lines (37 loc) · 1.1 KB
/
grid_dependence.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
#!/usr/bin/env python
import sys
import os
import glob
import shutil
from meshes import make_meshes
from elmer import make_elmer_meshes
import run
# ------------
def main(argv):
"""
This script runs inversions for Helheim at successively higher
mesh resolution.
"""
glacier = "helheim"
dxs = [300, 150, 75]
for dx in dxs:
os.chdir("meshes")
make_meshes.main(["-l", str(dx)])
os.chdir("../")
os.chdir("elmer")
make_elmer_meshes.main(argv)
os.chdir("../")
directory = glacier + "_grid_dependence"
os.mkdir(directory)
log_file_name = directory + "/log_dx-" + str(dx) + ".txt"
run.main(["-g", glacier,
"-i", "20",
"-o", log_file_name])
for filename in glob.glob("meshes/" + glacier + "/"
+ glacier + ".[1-2]*"):
shutil.move(filename, directory)
shutil.move("elmer/" + glacier + "3d", directory)
shutil.rmtree("elmer/" + glacier)
# -----------------------
if __name__ == "__main__":
main(sys.argv[1:])