Giter VIP home page Giter VIP logo

massbalance-sandbox's Introduction

docs/_static/logo.png


OGGM is a modular open source model for glacier dynamics

OGGM is able to simulate past and future mass balance, volume and geometry of (almost) any glacier in the world, in a fully automated and extensible workflow.

The model accounts for glacier geometry (including contributory branches) and includes an explicit ice dynamics module. We rely exclusively on publicly available data for calibration and validation. OGGM is modular and supports novel modelling workflows: it LOVES to be remixed and reused!

docs/_static/ex_tasman.jpg

Installation, documentation

The documentation is hosted on ReadTheDocs: http://docs.oggm.org

Get in touch

About

Version:Pypi version Supported python versions
Citation:GMD Paper Zenodo
Tests:Code coverage Linux build status Mass balance cross validation Documentation status Benchmark status
License:BSD-3-Clause License
Authors:

See the version history for a list of all contributors.

massbalance-sandbox's People

Contributors

fmaussion avatar lilianschuster avatar sarah-hanus avatar

Stargazers

 avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar  avatar

massbalance-sandbox's Issues

Issues with processing climate data

Hi Lilian,

First of all, thanks a lot for all these tools. They are really exhaustive, I'm sure I'll be able to find a lot of useful stuff for my project. I'm particularly interested in the daily climate processing, temperature scaling and initialization of snow/ice fractions.

I've been following the how_to_daily_input_daily_output.ipynb notebook in order to try a few things. I'm encountering a couple of issues when trying to process daily climate data for a gdir. When I do (bear in mind that the source code here is in Julia, but all the libraries called are in Python), I get the following error:

MBsandbox = pyimport("MBsandbox.mbmod_daily_oneflowline")

PARAMS["hydro_month_nh"]=1
climate="W5E5"
MBsandbox.process_w5e5_data(gdir, climate_type=climate, temporal_resol="daily") 


OSError(-101, 'NetCDF: HDF error')
  File "/Users/Bolib001/Python/oggm/oggm/utils/_workflow.py", line 490, in _entity_task
    out = task_func(gdir, **kwargs)
  File "/Users/Bolib001/Python/massbalance-sandbox/MBsandbox/mbmod_daily_oneflowline.py", line 435, in process_w5e5_data
    with xr.open_dataset(path_tmp) as ds:
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/api.py", line 497, in open_dataset
    backend_ds = backend.open_dataset(
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/netCDF4_.py", line 551, in open_dataset
    store = NetCDF4DataStore.open(
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/netCDF4_.py", line 380, in open
    return cls(manager, group=group, mode=mode, lock=lock, autoclose=autoclose)
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/netCDF4_.py", line 328, in __init__
    self.format = self.ds.data_model
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/netCDF4_.py", line 389, in ds
    return self._acquire()
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/netCDF4_.py", line 383, in _acquire
    with self._manager.acquire_context(needs_lock) as root:
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/contextlib.py", line 119, in __enter__
    return next(self.gen)
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/file_manager.py", line 187, in acquire_context
    file, cached = self._acquire_with_cache_info(needs_lock)
  File "/Users/Bolib001/miniconda3/envs/oggm_env/lib/python3.9/site-packages/xarray/backends/file_manager.py", line 205, in _acquire_with_cache_info
    file = self._opener(*self._args, **kwargs)
  File "src/netCDF4/_netCDF4.pyx", line 2307, in netCDF4._netCDF4.Dataset.__init__
  File "src/netCDF4/_netCDF4.pyx", line 1925, in netCDF4._netCDF4._ensure_nc_success
PyError ($(Expr(:escape, :(ccall(#= /Users/Bolib001/.julia/packages/PyCall/BD546/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'OSError'>
OSError(-101, 'NetCDF: HDF error')

After looking around a little bit, this seems to be quite a common bug, which can be due to different things: either corrupted .nc files or problems with the order of installation of the h5py and netCDF4 libraries. I have tried multiple things with no success so far. I was just wondering if you have already encountered this problem and if you know of any easy fix :)

I also tested MBsandbox.process_era5_daily_data(gdir) as suggested in the notebook, but when it just seems to run forever (I stop it after several minutes running). Is it OK to use this function or am I missing something?

Thanks a lot again!

dev-pointMB

Hi OGGM developpers!

In the dev-pointMB branch i'm going to develop a point mass balance, that is, add a class in here named TIModel_Sfc_Type_point, very similar to TIModel_Sfc_Type . This is going to be applied to a future 2-D mb over glacier surface.

FRA.

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.