27 lines
638 B
Python
27 lines
638 B
Python
|
import configparser
|
||
|
|
import pathlib
|
||
|
|
|
||
|
|
import numpy as np
|
||
|
|
import pandas as pd
|
||
|
|
import matplotlib.pyplot as plt
|
||
|
|
|
||
|
|
from .read_swash import *
|
||
|
|
|
||
|
|
|
||
|
|
config = configparser.ConfigParser()
|
||
|
|
config.read("config.ini")
|
||
|
|
|
||
|
|
cache = pathlib.Path(config.get("out", "root"))
|
||
|
|
root = pathlib.Path(config.get("swash", "out"))
|
||
|
|
|
||
|
|
bathy = pd.read_hdf(cache.joinpath("bathy.h5"), "bathy")
|
||
|
|
n = bathy.index.size
|
||
|
|
|
||
|
|
botl = read_nohead_scalar(root.joinpath("botl.dat"), n)
|
||
|
|
dep = read_nohead_scalar(root.joinpath("dep.dat"), n)
|
||
|
|
watl = read_nohead_scalar(root.joinpath("watl.dat"), n)
|
||
|
|
|
||
|
|
plt.plot(dep.T-botl.T, label="dep")
|
||
|
|
plt.plot(-botl.T[:, 0], label="botl")
|
||
|
|
plt.show(block=True)
|