Source code for opendrift.models.openoil.adios.dirjs

# This file is part of OpenDrift.
#
# OpenDrift 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, version 2
#
# OpenDrift 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 details.
#
# You should have received a copy of the GNU General Public License
# along with OpenDrift.  If not, see <https://www.gnu.org/licenses/>.
#
# Copyright 2021, Gaute Hope, MET Norway

from importlib import resources
from pathlib import Path
import logging
import json
import itertools
from functools import lru_cache

from .oil import OpendriftOil

logger = logging.getLogger(__name__)


[docs] @lru_cache(maxsize=None) def __get_archive__(): import lzma oils = [] with resources.open_binary('opendrift.models.openoil.adios', 'oils.xz') as archive: with lzma.open(archive, 'rt') as c: oils = json.load(c) # Add additional oils for f in resources.contents('opendrift.models.openoil.adios.extra_oils'): if Path(f).suffix == '.json': o = json.loads(resources.read_text('opendrift.models.openoil.adios.extra_oils', f)) o = { 'data': { 'attributes' : o } } o['data']['_id'] = o['data']['attributes']['oil_id'] o['data']['attributes']['metadata']['location'] = 'NORWAY' logger.debug(f"Adding additional oil: {f}..: {o['data']['_id']}, {o['data']['attributes']['metadata']['name']}") oils.append(o) return oils
[docs] def get_oil_names(location=None): a = __get_archive__() if location is not None: a = filter( lambda o: o['data']['attributes']['metadata'].get( 'location', None) == location, a) return [o['data']['attributes']['metadata']['name'] for o in a]
[docs] def oils(limit=50, query=''): oils = filter( lambda o: query in o['data']['attributes']['metadata']['name'], __get_archive__()) return list(OpendriftOil(o) for o in itertools.islice(oils, limit))
[docs] def find_full_oil_from_name(name) -> 'OpendriftOil': logger.info(f'Querying ADIOS database for oil: {name}') o = oils(query=name) if len(o) > 1: logger.warning( f"Several oils found with name: {name}: {[oo.id for oo in o]}, using first." ) elif len(o) < 1: raise ValueError(f"No oil found with name: {name}") return o[0]
[docs] def get_full_oil_from_id(_id) -> 'OpendriftOil': logger.debug(f"Fetching full oil: {_id}") oils = filter(lambda o: _id == o['data']['_id'], __get_archive__()) oil = next(OpendriftOil(o) for o in oils) return oil