Zooplankton biomass
Contents
Zooplankton biomass#
%load_ext autoreload
%autoreload 2
import xarray as xr
from dask.distributed import Client
import numpy as np
import catalog
import util
import pandas as pd
xr.set_options(keep_attrs=True)
import matplotlib.pyplot as plt
import cartopy
import cartopy.crs as ccrs
import pop_tools
/glade/work/kristenk/miniconda/envs/cesm-exp/lib/python3.7/site-packages/dask_jobqueue/core.py:20: FutureWarning: tmpfile is deprecated and will be removed in a future release. Please use dask.utils.tmpfile instead.
from distributed.utils import tmpfile
ds_grid = pop_tools.get_grid('POP_gx1v7')
lons=ds_grid.TLONG
lats=ds_grid.TLAT
Parameters#
casename = 'g.e22a06.G1850ECOIAF_JRA_PHYS_DEV.TL319_g17.scope_v1'
component = 'pop'
stream = 'h'
cluster_scheduler_address = None
# Parameters
casename = "g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001branch"
component = "pop"
stream = "h"
cluster_scheduler_address = "tcp://10.12.206.3:35506"
assert component in ['pop']
assert stream in ['h', 'h.ecosys.nday1']
Connect to cluster#
if cluster_scheduler_address is None:
cluster, client = util.get_ClusterClient()
cluster.scale(12)
else:
client = Client(cluster_scheduler_address)
client
Client
Client-5b853e36-9359-11ed-aaa7-3cecef1acc46
| Connection method: Direct | |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/8787/status |
Scheduler Info
Scheduler
Scheduler-7b28f886-17c1-40f8-ae24-3360603786e8
| Comm: tcp://10.12.206.3:35506 | Workers: 32 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/8787/status | Total threads: 32 |
| Started: 31 minutes ago | Total memory: 744.96 GiB |
Workers
Worker: PBSCluster-0
| Comm: tcp://10.12.206.60:40381 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/37473/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.60:43642 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-she593m8 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 515.69 MiB | Spilled bytes: 0 B |
| Read bytes: 310.15 MiB | Write bytes: 246.93 kiB |
Worker: PBSCluster-1
| Comm: tcp://10.12.206.57:34913 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/43168/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.57:40520 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-ais4w5nl | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 555.05 MiB | Spilled bytes: 0 B |
| Read bytes: 20.45 MiB | Write bytes: 11.19 kiB |
Worker: PBSCluster-10
| Comm: tcp://10.12.206.55:40610 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/41580/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:41770 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-24q12p4a | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 4.0% | Last seen: Just now |
| Memory usage: 547.65 MiB | Spilled bytes: 0 B |
| Read bytes: 21.91 kiB | Write bytes: 14.17 kiB |
Worker: PBSCluster-11
| Comm: tcp://10.12.206.35:36760 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/34231/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.35:39163 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-tm589a2g | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 608.83 MiB | Spilled bytes: 0 B |
| Read bytes: 21.61 kiB | Write bytes: 12.73 kiB |
Worker: PBSCluster-12
| Comm: tcp://10.12.206.42:46667 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/44609/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:43503 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-63bl396p | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 512.89 MiB | Spilled bytes: 0 B |
| Read bytes: 19.54 kiB | Write bytes: 8.24 kiB |
Worker: PBSCluster-13
| Comm: tcp://10.12.206.55:34741 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/41987/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:39933 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-3taktpky | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 519.87 MiB | Spilled bytes: 0 B |
| Read bytes: 20.31 kiB | Write bytes: 8.74 kiB |
Worker: PBSCluster-14
| Comm: tcp://10.12.206.42:44723 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/38884/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:37139 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-btgwv71o | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 563.05 MiB | Spilled bytes: 0 B |
| Read bytes: 19.55 kiB | Write bytes: 12.28 kiB |
Worker: PBSCluster-15
| Comm: tcp://10.12.206.60:34840 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/36356/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.60:36517 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-s50upym9 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 541.54 MiB | Spilled bytes: 0 B |
| Read bytes: 187.61 MiB | Write bytes: 83.42 kiB |
Worker: PBSCluster-16
| Comm: tcp://10.12.206.55:33480 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/44604/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:34515 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-20gqvxq8 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 521.15 MiB | Spilled bytes: 0 B |
| Read bytes: 20.27 kiB | Write bytes: 11.14 kiB |
Worker: PBSCluster-17
| Comm: tcp://10.12.206.42:46005 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/44117/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:32947 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-18iilri4 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 551.21 MiB | Spilled bytes: 0 B |
| Read bytes: 19.56 kiB | Write bytes: 9.86 kiB |
Worker: PBSCluster-18
| Comm: tcp://10.12.206.60:33226 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/42773/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.60:46337 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-8qq6z3y3 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 4.0% | Last seen: Just now |
| Memory usage: 524.71 MiB | Spilled bytes: 0 B |
| Read bytes: 187.59 MiB | Write bytes: 84.45 kiB |
Worker: PBSCluster-19
| Comm: tcp://10.12.206.36:36934 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/38605/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.36:40345 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-0wibvgh_ | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 549.18 MiB | Spilled bytes: 0 B |
| Read bytes: 96.65 kiB | Write bytes: 40.51 kiB |
Worker: PBSCluster-2
| Comm: tcp://10.12.206.42:35532 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/37233/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:38029 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-tvn8qi89 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 506.30 MiB | Spilled bytes: 0 B |
| Read bytes: 19.53 kiB | Write bytes: 9.85 kiB |
Worker: PBSCluster-20
| Comm: tcp://10.12.206.57:44618 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/35828/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.57:39590 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-j87b3i6z | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 526.30 MiB | Spilled bytes: 0 B |
| Read bytes: 20.50 MiB | Write bytes: 9.60 kiB |
Worker: PBSCluster-21
| Comm: tcp://10.12.206.42:36146 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/42013/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:42567 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-te65hpsa | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 573.22 MiB | Spilled bytes: 0 B |
| Read bytes: 19.73 kiB | Write bytes: 14.00 kiB |
Worker: PBSCluster-22
| Comm: tcp://10.12.206.37:45266 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/41617/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.37:46778 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-6zjl_rf2 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 524.25 MiB | Spilled bytes: 0 B |
| Read bytes: 5.36 MiB | Write bytes: 140.84 kiB |
Worker: PBSCluster-23
| Comm: tcp://10.12.206.57:36646 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/33282/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.57:39220 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-tedt0oqb | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 562.36 MiB | Spilled bytes: 0 B |
| Read bytes: 20.52 MiB | Write bytes: 10.42 kiB |
Worker: PBSCluster-24
| Comm: tcp://10.12.206.35:43781 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/41280/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.35:41682 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-oenmf7ye | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 493.55 MiB | Spilled bytes: 0 B |
| Read bytes: 16.61 kiB | Write bytes: 4.67 kiB |
Worker: PBSCluster-25
| Comm: tcp://10.12.206.42:38958 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/45762/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:41096 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-r999w54s | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 502.65 MiB | Spilled bytes: 0 B |
| Read bytes: 19.59 kiB | Write bytes: 9.88 kiB |
Worker: PBSCluster-26
| Comm: tcp://10.12.206.35:41019 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/46197/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.35:41093 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-n8z_ugpj | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 582.85 MiB | Spilled bytes: 0 B |
| Read bytes: 21.66 kiB | Write bytes: 15.19 kiB |
Worker: PBSCluster-27
| Comm: tcp://10.12.206.40:37756 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/34868/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.40:46061 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-ah4h4uwy | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 520.84 MiB | Spilled bytes: 0 B |
| Read bytes: 13.97 kiB | Write bytes: 10.07 kiB |
Worker: PBSCluster-28
| Comm: tcp://10.12.206.55:40536 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/40763/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:43688 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-mqh1_swd | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 516.66 MiB | Spilled bytes: 0 B |
| Read bytes: 20.28 kiB | Write bytes: 10.34 kiB |
Worker: PBSCluster-29
| Comm: tcp://10.12.206.35:43010 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/45869/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.35:38155 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-0zp2yifw | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 550.70 MiB | Spilled bytes: 0 B |
| Read bytes: 21.56 kiB | Write bytes: 11.82 kiB |
Worker: PBSCluster-3
| Comm: tcp://10.12.206.55:43243 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/45154/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:40506 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-d77lpe1m | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 511.52 MiB | Spilled bytes: 0 B |
| Read bytes: 20.29 kiB | Write bytes: 7.93 kiB |
Worker: PBSCluster-30
| Comm: tcp://10.12.206.42:40059 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/38562/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:35150 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-s1j98ibm | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 534.24 MiB | Spilled bytes: 0 B |
| Read bytes: 19.54 kiB | Write bytes: 8.24 kiB |
Worker: PBSCluster-31
| Comm: tcp://10.12.206.55:40509 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/46123/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:42505 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-0193oowl | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 491.86 MiB | Spilled bytes: 0 B |
| Read bytes: 20.29 kiB | Write bytes: 11.15 kiB |
Worker: PBSCluster-4
| Comm: tcp://10.12.206.60:33761 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/38547/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.60:43086 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-go_20xsh | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 553.15 MiB | Spilled bytes: 0 B |
| Read bytes: 310.89 MiB | Write bytes: 247.07 kiB |
Worker: PBSCluster-5
| Comm: tcp://10.12.206.56:36183 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/39960/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.56:36523 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-ak_et0qx | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 528.58 MiB | Spilled bytes: 0 B |
| Read bytes: 296.40 MiB | Write bytes: 145.73 MiB |
Worker: PBSCluster-6
| Comm: tcp://10.12.206.55:42418 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/43316/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:35755 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-w069mt42 | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 499.32 MiB | Spilled bytes: 0 B |
| Read bytes: 25.86 kiB | Write bytes: 177.44 kiB |
Worker: PBSCluster-7
| Comm: tcp://10.12.206.42:38064 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/33959/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.42:37237 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-8g92ok8a | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 514.12 MiB | Spilled bytes: 0 B |
| Read bytes: 19.71 kiB | Write bytes: 13.18 kiB |
Worker: PBSCluster-8
| Comm: tcp://10.12.206.35:46716 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/38950/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.35:46823 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-t0bj5mka | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 0.0% | Last seen: Just now |
| Memory usage: 520.34 MiB | Spilled bytes: 0 B |
| Read bytes: 21.61 kiB | Write bytes: 12.73 kiB |
Worker: PBSCluster-9
| Comm: tcp://10.12.206.55:36074 | Total threads: 1 |
| Dashboard: https://jupyterhub.hpc.ucar.edu/stable/user/kristenk/proxy/41381/status | Memory: 23.28 GiB |
| Nanny: tcp://10.12.206.55:35112 | |
| Local directory: /glade/scratch/kristenk/dask-workers/dask-worker-space/worker-vzrt3wyj | |
| Tasks executing: 0 | Tasks in memory: 0 |
| Tasks ready: 0 | Tasks in flight: 0 |
| CPU usage: 2.0% | Last seen: Just now |
| Memory usage: 490.55 MiB | Spilled bytes: 0 B |
| Read bytes: 20.28 kiB | Write bytes: 7.93 kiB |
Load the data#
dsets = catalog.to_dataset_dict(
case=casename,
component=component,
stream=stream,
)
dsets.keys()
g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001branch pop h /glade/scratch/kristenk/archive
dict_keys(['g.e22.GOMIPECOIAF_JRA-1p4-2018.TL319_g17.4p2z.001branch.pop.h'])
ds = dsets[f'{casename}.{component}.{stream}']
Variable list#
variables =['mesozooC', 'microzooC']
keep_vars=['z_t','z_t_150m','dz','time_bound','time','TAREA','TLAT','TLONG'] + variables
ds = ds.drop([v for v in ds.variables if v not in keep_vars])
Make a monthly climatology dataset#
time = pd.date_range(start="2009-01-01", end="2019-01-01",freq="M")
mon_ds = ds.copy()
mon_ds['time'] = time
mon_ds = ds.groupby('time.month').mean('time')
Compute long-term mean#
ds = ds.mean(dim='time')
Plot mesozooplankton and microzooplankton biomass at the surface#
######
fig = plt.figure(figsize=(8,5))
ax = fig.add_subplot(2,1,1, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title('microzooC at surface', fontsize=12)
lon, lat, field = util.adjust_pop_grid(lons, lats, ds.microzooC.isel(z_t_150m=0))
pc=ax.pcolormesh(lon, lat, field, cmap='jet',vmin=0,vmax=2,transform=ccrs.PlateCarree())
cbar1 = fig.colorbar(pc, ax=ax,extend='max',label='microzooC (mmol m$^{-3}$)')
ax = fig.add_subplot(2,1,2, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title('mesozooC at surface', fontsize=12)
lon, lat, field = util.adjust_pop_grid(lons, lats, ds.mesozooC.isel(z_t_150m=0))
pc=ax.pcolormesh(lon, lat, field, cmap='jet',vmin=0,vmax=4,transform=ccrs.PlateCarree())
cbar1 = fig.colorbar(pc, ax=ax,extend='max',label='mesozooC (mmol m$^{-3}$)')
Compare mesozooplankton biomass to COPEPOD database#
file = '/glade/u/home/kristenk/work/copepod-biomass/data/POP_gx1v7/copepod-2012__cmass-m00-qtr.zarr'
ds_copepod = xr.open_zarr(file, consolidated=False)
ds_copepod['copepod_C']=ds_copepod.copepod_C/12.011
fig = plt.figure(figsize=(12,3))
ax = fig.add_subplot(1,2,1, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title('COPEPOD dataset', fontsize=12)
lon, lat, field = util.adjust_pop_grid(lons, lats, ds_copepod.copepod_C)
pc=ax.pcolormesh(lon, lat, field, cmap='Reds',vmin=0,vmax=2,transform=ccrs.PlateCarree())
ax.coastlines('110m',linewidth=0.5)
ax = fig.add_subplot(1,2,2, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title('CESM ${\it Mesozooplankton}$ biomass', fontsize=12)
lon, lat, field = util.adjust_pop_grid(lons, lats, ds.mesozooC.mean(dim='z_t_150m'))
pc=ax.pcolormesh(lon, lat, field, cmap='Reds',vmin=0,vmax=2,transform=ccrs.PlateCarree())
ax.coastlines('110m',linewidth=0.5)
fig.subplots_adjust(right=0.8)
cbar_ax = fig.add_axes([0.85, 0.15, 0.02, 0.7])
fig.colorbar(pc, cax=cbar_ax,extend='max', label='top 150m/200m mean (mmol m$^{-3}$)');
Compare to observation-based GLMM of global mesozooplankton biomass climatology#
Petrik et al., submitted to Global Biogeochemical Cycles
In the comparison below the model seems to have somewhat too high biomass, but in this recent paper by Yang, Atkinson, et al., estimates of zooplankton biomass in the Southern Ocean are higher and more in the line with model results.
file = '/glade/u/home/kristenk/work/copepod-biomass/obsglmm_zmeso_vint_200m_monthly_climatology.nc'
ds_copepod_clim = xr.open_dataset(file)
ds_copepod_clim.zmeso200.attrs['units'] = 'mgC m-2'
months = ['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec']
### depth integrate and convert model to mol/m2
mon_ds['mesozooC_zint'] = ((mon_ds.mesozooC) * 10.).sum(dim='z_t_150m') #in mmol/m2
mon_ds['mesozooC_zint'] = mon_ds['mesozooC_zint'] * 12.011 #convert to mgC/m2
mon_ds['mesozooC_zint'].attrs['units'] = 'mgC m-2'
fig = plt.figure(figsize=(5,18))
for row in np.arange(1,13):
ts=row-1
plot = row*2 - 1
ax = fig.add_subplot(12,2,plot, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title(months[ts]+' obs', fontsize=12)
pc=ax.pcolormesh(ds_copepod_clim.Lon, ds_copepod_clim.Lat, ds_copepod_clim.zmeso200.isel(month=ts),
cmap='Reds',vmin=0,vmax=4000,transform=ccrs.PlateCarree())
ax.coastlines('110m',linewidth=0.5)
plot = row*2
ax = fig.add_subplot(12,2,plot, projection=ccrs.Robinson(central_longitude=305.0))
ax.set_title(months[ts]+' CESM', fontsize=12)
tmp = mon_ds.mesozooC_zint.isel(month=ts)
lon, lat, field = util.adjust_pop_grid(lons, lats, tmp)
pc=ax.pcolormesh(lon, lat, field, cmap='Reds',vmin=0,vmax=4000,transform=ccrs.PlateCarree())
ax.coastlines('110m',linewidth=0.5)
cbar_ax = fig.add_axes([0.92, 0.15, 0.03, 0.7])
fig.colorbar(pc, cax=cbar_ax,extend='max', label='Depth-integrated copepod biomass (mg m$^{-2}$)');
TO DO#
• Additional comparisons to make (Colleen’s suggestions): look at the seasonal cycle in different regions and also the relationship between zooplankton biomass and chlorophyll.
• Add production fields
• z-ratio comparison
• Comparisons to MAREDAT