import argparse import configparser import logging from multiprocessing import ThreadPool import pathlib import numpy as np import pandas as pd from .read_swash import ReadSwash parser = argparse.ArgumentParser(description="Convert swash output to numpy") parser.add_argument("-v", "--verbose", action="count", default=0) args = parser.parse_args() logging.basicConfig(level=max((10, 20 - 10 * args.verbose))) log = logging.getLogger("sws_npz") log.info("Starting sws -> npz converter") config = configparser.ConfigParser() config.read("config.ini") sws_out = pathlib.Path(config.get("swash", "out")) inp = pathlib.Path(config.get("post", "inp")) log.info(f"Reading swash output from '{sws_out}'") rsws = ReadSwash() rsws.read_time(sws_out.joinpath("tsec.dat")) rsws.read_x(sws_out.joinpath("xp.dat")) var = { "dep": rsws.read_scalar, "botl": rsws.read_const, "watl": rsws.read_scalar, "vel": rsws.read_vector, "press": rsws.read_scalar, "zk": rsws.read_scalar_lay, "velk": rsws.read_vector_lay, } inp.mkdir(exist_ok=True) with ThreadPool() as pool: pool.map( lambda name, f: np.save( inp.joinpath(name), f(sws_out.joinpath(name).with_suffix(".dat")) ), var.items(), )