https://github.com/perrette/pyglacier

Utilities to help running a glacier model written in fortran (not included)
https://github.com/perrette/pyglacier

Last synced: 12 months ago
JSON representation

Utilities to help running a glacier model written in fortran (not included)

• Host: GitHub
• URL: https://github.com/perrette/pyglacier
• Owner: perrette
• License: other
• Created: 2015-04-09T14:33:58.000Z (about 11 years ago)
• Default Branch: master
• Last Pushed: 2015-04-24T11:51:24.000Z (about 11 years ago)
• Last Synced: 2025-02-05T08:51:34.596Z (over 1 year ago)
• Language: Python
• Size: 195 KB
• Stars: 0
• Watchers: 2
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.md
  • License: LICENSE.md

Awesome Lists containing this project

README

          
# pyglacier

Wrapper to run and tune the glacier model

Install

-------

python setup.py install

Getting started

---------------

    import pyglacier

    pyglacier.setup('/path/to/glacier/code')

The following output will be printed to screen:

    ===============

    Setup pyglacier

    ===============

    ...pyglacier version:  0.0.0.dev-15faa90

    ...fortran code directory:  ../glaciercode//

    ...fortran code version:  73b7f72

    ...default executable:  $CODEDIR/main.exe

    ...default glacier lib:  $CODEDIR/.obj/wrapper.so

    ...default namelist:  $CODEDIR/params.nml

    ...default output directory:  ./out/

You can then define a glacier class from restart file 

(or whatever netCDF file with the expected glacier fields):

    gl = pyglacier.Glacier.read('path/to/restart.nc')

    stress = gl.compute_stress()  # return a dimarray.Dataset of stress calculations

    gl.update_velocity()  # update velocity so that stress residual is zero

    stress2 = gl.compute_stress() # check the new stress balance

The glacier can also be defined more simply with four fields:

    gl = pyglacier.Glacier(x, zb, w, h, u)

and parameters can be changed:

    gl.set_param(name, value, group=group) # group is only needed if param name is ambiguous

A more exhaustive documentation will be made later. For now just check inline help for more information 

about glacier methods:

    help(pyglacier.Glacier)