iframe-proxy | Sunbelt Computer Software

119 lines (96 loc) · 3.46 KB
# coding: utf-8
# Copyright 2020-2022 IFREMER (Brest, FRANCE), all rights reserved.
# contact -- mailto:nicolas.raillard@ifremer.fr
# This file is part of Resourcecode.
# Written by Logilab SA (contact@logilab.fr)
# Resourcecode is free software: you can redistribute it and/or modify it under
# the terms of the GNU General Public License as published by the Free
# Software Foundation, either version 3.0 of the License, or any later version.
# Resourcecode is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
# FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
# You should have received a copy of the GNU General Public License along
# with Resourcecode. If not, see <https://www.gnu.org/licenses/>.
import configparser
import logging
from pathlib import Path
from typing import Union, Tuple
import numpy as np
from numpy import triu_indices, tril_indices
CONFIG_FILEPATHS = [
    os.environ.get("RESOURCECODE_CONFIG_FILEPATH"),
    "./resourcecode.ini",
    "~/.config/resourcecode.ini",
    f"{sys.prefix}/etc/resourcecode/config.ini",
    Path(__file__).parent / "config" / "config.ini",
LOGGER = logging.getLogger("resourcecode.default")
LOGGER.addHandler(logging.StreamHandler())
LOGGER.setLevel(os.environ.get("RESOURCECODE_LOG_THRESHOLD", "WARNING"))
def get_config():
    config = configparser.ConfigParser()
    for config_filepath in CONFIG_FILEPATHS:
        LOGGER.debug("try reading config from %s", config_filepath)
        if not config_filepath:
            continue
        config_filepath = Path(config_filepath).expanduser()
        if not config_filepath.exists():
            continue
        config.read(config_filepath)
        LOGGER.info("config loaded from %s", config_filepath)
        break
        raise FileNotFoundError("no config file was found")
    return config
def set_trig(m, values, part="upper"):
    """Set the `values` upper/lower `part` of the square matrix `m`"""
    part = part.lower()
    if part == "upper":
        tri_indices = triu_indices(len(values), k=1)
    elif part == "lower":
        tri_indices = tril_indices(len(values), k=-1)
        raise ValueError(f"part must be 'lower' or 'upper', '{part}' given")
    for n, (i, j) in enumerate(zip(*tri_indices)):
        m[i, j] = values[n]
def haversine(lon1, lat1, lon2, lat2):
    Calculate the great circle distance between two points
    on the earth (specified in decimal degrees)
    EARTH_RADIUS_METER = 6367e3
    lon1, lat1, lon2, lat2 = map(np.radians, [lon1, lat1, lon2, lat2])
    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = np.sin(dlat / 2.0) ** 2 + np.cos(lat1) * np.cos(lat2) * np.sin(dlon / 2.0) ** 2
    c = 2 * np.arcsin(np.sqrt(a))
    return c * EARTH_RADIUS_METER
def zmcomp2metconv(
    u: Union[float, np.ndarray], v: Union[float, np.ndarray]
) -> Union[Tuple[float, float], Tuple[np.ndarray, np.ndarray]]:
    Converts wind or current zonal and meridional velocity components to
    magnitude and direction according to meteorological convention.
    Parameters
    ----------
        zonal velocity (ISU)
        meridional velocity (ISU)
    Results
    -------
        magnitude of speed (ISU)
        direction from which flow comes (<B0>)
    V = np.sqrt(u**2 + v**2)
    D = (270 - np.arctan2(v, u) * 180 / np.pi) % 360
    return (V, D)
Sunbelt Computer Software

PL/B Language Development and Support

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Uh oh!

Sunbelt Computer Software

PL/B Language Development and Support

Uh oh!

FilesExpand file tree

utils.py

Latest commit

History

utils.py

File metadata and controls
