• WARNING
• Installation
• Quickstart
  • Ice Discharge
    • Find the 2020 discharge for the top 3 discharging glaciers
  • Watson River discharge

• This is alpha level software under development.
• Please contact Ken Mankoff <[email protected]> before using.

pip install git+https://gitlab.com/mankoff/PROMICE.git

Can be run in new environment, for example:

Conda:

conda create -n PROMICE python=3.8
conda activate PROMICE
pip install --upgrade git+https://gitlab.com/mankoff/PROMICE.git

Test the installation

conda activate PROMICE
python

Then,

import promice
promice.ice_discharge.discharge(resolution="GIS").to_dataframe().tail()

from promice import ice_discharge as ice
ds = ice.discharge(resolution="GIS")
ds.to_dataframe().tail()

Load an xarray dataset of all discharge

from promice import ice_discharge as ice
ds = ice.discharge(resolution="gate")
print(ds)

<xarray.Dataset>
Dimensions:        (gate: 267, time: 449)
Coordinates:
  * time           (time) datetime64[ns] 1986-04-15 1986-05-15 ... 2020-07-29
  * gate           (gate) int64 1 3 4 5 6 7 8 10 ... 350 352 353 354 355 356 357
Data variables:
    discharge      (gate, time) float32 ...
    err            (gate, time) float32 ...
    coverage       (gate, time) float32 ...
    mean_x         (gate) int32 ...
    mean_y         (gate) int32 ...
    mean_lon       (gate) float32 ...
    mean_lat       (gate) float32 ...
    sector         (gate) int32 ...
    region         (gate) object ...
    name_Bjørk     (gate) object ...
    name_Mouginot  (gate) object ...
Attributes:
    featureType:      timeSeries
    title:            Greenland discharge
    summary:          Greenland discharge per gate
    keywords:         Greenland; Ice Discharge; Calving; Submarine Melt
    source:           git commit: b4f3e88
    creator_name:     Ken Mankoff
    creator_email:    [email protected]
    creator_url:      http://kenmankoff.com
    institution:      GEUS
    references:       10.22008/promice/ice_discharge
    product_version:  [2.]

top3 = ds.resample(time='A')\
         .mean()['discharge']\
         .to_dataframe()['discharge']\
         .unstack()\
         .iloc[-1]\
         .sort_values(ascending=False)[0:3]
top3

gate
184    34.403828
231    33.083332
190    27.318832
Name: 2020-12-31 00:00:00, dtype: float32

df_3 = ds.sel(gate=top3.index)\
         .to_dataframe()['discharge']\
         .unstack()\
         .T

df_3.head()

• Convert the gate ID to named glacier
• Subset to just 2020
• Monthly average

df_3.rename(inplace=True, columns=dict([(i, str(ds.sel({'gate':i})['name_Bjørk'].values)) for i in df_3.columns]))
df_3[df_3.index.year >= 2020].resample('M').mean()

from promice import watson
df = watson.discharge()
df['2012-06-01':].head(3)