Skip to content

Instantly share code, notes, and snippets.

@abkfenris

abkfenris/run3.ipynb

Last active October 18, 2019 19:05
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save abkfenris/ab2415124c4c6e3fbbeccaf890a311f5 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save abkfenris/ab2415124c4c6e3fbbeccaf890a311f5 to your computer and use it in GitHub Desktop.
Download ZIP
ROMS exploration
Display the source blob
Display the rendered blob
Raw
run3.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import xarray as xr"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/usr/local/miniconda3/envs/IOOS/lib/python3.7/site-packages/ipykernel_launcher.py:1: FutureWarning: In xarray version 0.14 the default behaviour of `open_mfdataset`\n",
"will change. To retain the existing behavior, pass\n",
"combine='nested'. To use future default behavior, pass\n",
"combine='by_coords'. See\n",
"http://xarray.pydata.org/en/stable/combining.html#combining-multi\n",
"\n",
" \"\"\"Entry point for launching an IPython kernel.\n",
"/usr/local/miniconda3/envs/IOOS/lib/python3.7/site-packages/xarray/backends/api.py:934: FutureWarning: The datasets supplied have global dimension coordinates. You may want\n",
"to use the new `combine_by_coords` function (or the\n",
"`combine='by_coords'` option to `open_mfdataset`) to order the datasets\n",
"before concatenation. Alternatively, to continue concatenating based\n",
"on the order the datasets are supplied in future, please use the new\n",
"`combine_nested` function (or the `combine='nested'` option to\n",
"open_mfdataset).\n",
" from_openmfds=True,\n"
]
},
{
"data": {
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (boundary: 4, eta_rho: 362, eta_u: 362, eta_v: 361, ocean_time: 5, s_rho: 40, s_w: 41, tracer: 26, xi_rho: 722, xi_u: 721, xi_v: 722)\n",
"Coordinates:\n",
" * s_w (s_w) float64 -1.0 -0.975 -0.95 ... -0.05 -0.025 0.0\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 ... -0.0375 -0.0125\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Dimensions without coordinates: boundary, eta_rho, eta_u, eta_v, tracer, xi_rho, xi_u, xi_v\n",
"Data variables:\n",
" AKt (ocean_time, s_w, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 41, 362, 722), meta=np.ndarray>\n",
" AKv (ocean_time, s_w, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 41, 362, 722), meta=np.ndarray>\n",
" Akt_bak (ocean_time, tracer) float64 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" Akv_bak (ocean_time) float64 1e-05 1e-05 1e-05 1e-05 1e-05\n",
" BioIter (ocean_time) int32 1 1 1 1 1\n",
" Cs_r (ocean_time, s_rho) float64 dask.array<chunksize=(1, 40), meta=np.ndarray>\n",
" Cs_w (ocean_time, s_w) float64 dask.array<chunksize=(1, 41), meta=np.ndarray>\n",
" FSobc_in (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" FSobc_out (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" Falpha (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" Fbeta (ocean_time) float64 4.0 4.0 4.0 4.0 4.0\n",
" Fgamma (ocean_time) float64 0.284 0.284 0.284 0.284 0.284\n",
" Hsbl (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" Lm2CLM (ocean_time) int32 1 1 1 1 1\n",
" Lm3CLM (ocean_time) int32 1 1 1 1 1\n",
" LnudgeM2CLM (ocean_time) int32 1 1 1 1 1\n",
" LnudgeM3CLM (ocean_time) int32 1 1 1 1 1\n",
" LnudgeTCLM (ocean_time, tracer) int32 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" LsshCLM (ocean_time) int32 1 1 1 1 1\n",
" LtracerCLM (ocean_time, tracer) int32 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" LtracerSponge (ocean_time, tracer) int32 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" LtracerSrc (ocean_time, tracer) int32 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" LuvSponge (ocean_time) int32 0 0 0 0 0\n",
" LuvSrc (ocean_time) int32 0 0 0 0 0\n",
" LwSrc (ocean_time) int32 0 0 0 0 0\n",
" M2nudg (ocean_time) float64 0.002778 0.002778 ... 0.002778\n",
" M2obc_in (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" M2obc_out (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" M3nudg (ocean_time) float64 0.002778 0.002778 ... 0.002778\n",
" M3obc_in (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" M3obc_out (ocean_time, boundary) float64 dask.array<chunksize=(1, 4), meta=np.ndarray>\n",
" NKML (ocean_time) int32 1 1 1 1 1\n",
" P_C_max_Di (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" P_C_max_Lg (ocean_time) float64 1.302e-05 ... 1.302e-05\n",
" P_C_max_Sm (ocean_time) float64 1.302e-05 ... 1.302e-05\n",
" Pair (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" RHO_0 (ocean_time) float64 1.035e+03 ... 1.035e+03\n",
" Runoff (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" Tcline (ocean_time) float64 250.0 250.0 250.0 250.0 250.0\n",
" Tnudg (ocean_time, tracer) float64 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" Tnudg_SSS (ocean_time) float64 0.01111 0.01111 ... 0.01111\n",
" Tobc_in (ocean_time, boundary, tracer) float64 dask.array<chunksize=(1, 4, 26), meta=np.ndarray>\n",
" Tobc_out (ocean_time, boundary, tracer) float64 dask.array<chunksize=(1, 4, 26), meta=np.ndarray>\n",
" Uwind (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" Vstretching (ocean_time) int32 4 4 4 4 4\n",
" Vtransform (ocean_time) int32 2 2 2 2 2\n",
" Vwind (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" Znudg (ocean_time) float64 0.002778 0.002778 ... 0.002778\n",
" Zob (ocean_time) float64 0.02 0.02 0.02 0.02 0.02\n",
" Zos (ocean_time) float64 0.02 0.02 0.02 0.02 0.02\n",
" a1_co2 (ocean_time) float64 2.069e+03 ... 2.069e+03\n",
" a1_o2 (ocean_time) float64 1.93e+03 1.93e+03 ... 1.93e+03\n",
" a2_co2 (ocean_time) float64 -118.6 -118.6 ... -118.6 -118.6\n",
" a2_o2 (ocean_time) float64 -117.5 -117.5 ... -117.5 -117.5\n",
" a3_co2 (ocean_time) float64 2.931 2.931 2.931 2.931 2.931\n",
" a3_o2 (ocean_time) float64 3.116 3.116 3.116 3.116 3.116\n",
" a4_co2 (ocean_time) float64 -0.027 -0.027 ... -0.027 -0.027\n",
" a4_o2 (ocean_time) float64 -0.0306 -0.0306 ... -0.0306\n",
" a_0 (ocean_time) float64 2.009 2.009 2.009 2.009 2.009\n",
" a_1 (ocean_time) float64 3.22 3.22 3.22 3.22 3.22\n",
" a_2 (ocean_time) float64 4.05 4.05 4.05 4.05 4.05\n",
" a_3 (ocean_time) float64 4.945 4.945 4.945 4.945 4.945\n",
" a_4 (ocean_time) float64 -0.2568 -0.2568 ... -0.2568\n",
" a_5 (ocean_time) float64 3.888 3.888 3.888 3.888 3.888\n",
" agg_Di (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" agg_Lg (ocean_time) float64 3.472 3.472 3.472 3.472 3.472\n",
" agg_Sm (ocean_time) float64 1.157 1.157 1.157 1.157 1.157\n",
" alk (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" alk_2_n_denit (ocean_time) float64 1.169 1.169 1.169 1.169 1.169\n",
" alpha_Di (ocean_time) float64 4.6e-05 4.6e-05 ... 4.6e-05\n",
" alpha_Lg (ocean_time) float64 4.6e-05 4.6e-05 ... 4.6e-05\n",
" alpha_Sm (ocean_time) float64 9.2e-05 9.2e-05 ... 9.2e-05\n",
" alpha_fescav (ocean_time) float64 4.756e-07 ... 4.756e-07\n",
" b_0 (ocean_time) float64 -0.006245 ... -0.006245\n",
" b_1 (ocean_time) float64 -0.007376 ... -0.007376\n",
" b_2 (ocean_time) float64 -0.01034 -0.01034 ... -0.01034\n",
" b_3 (ocean_time) float64 -0.008171 ... -0.008171\n",
" beta_fescav (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" bresp_Di (ocean_time) float64 2.894e-07 ... 2.894e-07\n",
" bresp_Lg (ocean_time) float64 2.894e-07 ... 2.894e-07\n",
" bresp_Sm (ocean_time) float64 2.604e-07 ... 2.604e-07\n",
" bresp_bact (ocean_time) float64 8.681e-08 ... 8.681e-08\n",
" bresp_lgz (ocean_time) float64 3.704e-08 ... 3.704e-08\n",
" bresp_mdz (ocean_time) float64 9.259e-08 ... 9.259e-08\n",
" bresp_smz (ocean_time) float64 2.315e-07 ... 2.315e-07\n",
" c_0 (ocean_time) float64 -4.887e-07 ... -4.887e-07\n",
" c_2_n (ocean_time) float64 6.625 6.625 6.625 6.625 6.625\n",
" ca_2_n_arag (ocean_time) float64 0.1325 0.1325 ... 0.1325 0.1325\n",
" ca_2_n_calc (ocean_time) float64 0.06625 0.06625 ... 0.06625\n",
" caco3_sat_max (ocean_time) float64 10.0 10.0 10.0 10.0 10.0\n",
" cadet_arag (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" cadet_calc (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" coef_hp (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" dic (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" dstart (ocean_time) datetime64[ns] 1980-01-28T13:00:00.002880 ... 1980-01-28T13:00:00.002880\n",
" dt (ocean_time) float64 240.0 240.0 240.0 240.0 240.0\n",
" dtfast (ocean_time) float64 8.0 8.0 8.0 8.0 8.0\n",
" el (ocean_time) float64 2.493e+06 ... 2.493e+06\n",
" evaporation (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" exu_Di (ocean_time) float64 0.13 0.13 0.13 0.13 0.13\n",
" exu_Lg (ocean_time) float64 0.13 0.13 0.13 0.13 0.13\n",
" exu_Sm (ocean_time) float64 0.13 0.13 0.13 0.13 0.13\n",
" fe_2_n_max_Di (ocean_time) float64 0.003312 0.003312 ... 0.003312\n",
" fe_2_n_max_Lg (ocean_time) float64 0.003312 0.003312 ... 0.003312\n",
" fe_2_n_max_Sm (ocean_time) float64 0.0003313 ... 0.0003313\n",
" fe_2_n_sed (ocean_time) float64 0.006625 0.006625 ... 0.006625\n",
" fe_2_n_upt_fac (ocean_time) float64 1.5e-05 1.5e-05 ... 1.5e-05\n",
" fe_coast (ocean_time) float64 1e-11 1e-11 1e-11 1e-11 1e-11\n",
" felig_2_don (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" felig_bkg (ocean_time) float64 1e-09 1e-09 1e-09 1e-09 1e-09\n",
" gamma2 (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" gamma_cadet_arag (ocean_time) float64 1.523e-06 ... 1.523e-06\n",
" gamma_cadet_calc (ocean_time) float64 8.618e-07 ... 8.618e-07\n",
" gamma_fescav (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" gamma_irr_mem (ocean_time) float64 1.157e-05 ... 1.157e-05\n",
" gamma_mu_mem (ocean_time) float64 1.157e-05 ... 1.157e-05\n",
" gamma_ndet (ocean_time) float64 6.156e-06 ... 6.156e-06\n",
" gamma_nitrif (ocean_time) float64 3.858e-07 ... 3.858e-07\n",
" gamma_sidet (ocean_time) float64 5.787e-07 ... 5.787e-07\n",
" gamma_sldon (ocean_time) float64 1.286e-07 ... 1.286e-07\n",
" gamma_sldop (ocean_time) float64 1.286e-07 ... 1.286e-07\n",
" gamma_srdon (ocean_time) float64 1.76e-09 1.76e-09 ... 1.76e-09\n",
" gamma_srdop (ocean_time) float64 7.922e-09 ... 7.922e-09\n",
" gge_max_bact (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" gge_max_lgz (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" gge_max_mdz (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" gge_max_smz (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" hc (ocean_time) float64 250.0 250.0 250.0 250.0 250.0\n",
" hp_ipa_bact (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_ipa_det (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_ipa_diaz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_ipa_lgp (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_ipa_lgz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" hp_ipa_mdz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" hp_ipa_smp (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_ipa_smz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_det (ocean_time) float64 0.35 0.35 0.35 0.35 0.35\n",
" hp_phi_ldon (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_ldop (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_nh4 (ocean_time) float64 0.65 0.65 0.65 0.65 0.65\n",
" hp_phi_po4 (ocean_time) float64 0.65 0.65 0.65 0.65 0.65\n",
" hp_phi_sldon (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_sldop (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_srdon (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" hp_phi_srdop (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" htotal_in (ocean_time) float64 1e-08 1e-08 1e-08 1e-08 1e-08\n",
" htotal_scale_hi (ocean_time) float64 100.0 100.0 100.0 100.0 100.0\n",
" htotal_scale_lo (ocean_time) float64 0.01 0.01 0.01 0.01 0.01\n",
" imax_hp (ocean_time) float64 1.042e-06 ... 1.042e-06\n",
" imax_lgz (ocean_time) float64 2.662e-06 ... 2.662e-06\n",
" imax_mdz (ocean_time) float64 6.597e-06 ... 6.597e-06\n",
" imax_smz (ocean_time) float64 1.644e-05 ... 1.644e-05\n",
" io_fescav (ocean_time) float64 10.0 10.0 10.0 10.0 10.0\n",
" irr_inhibit (ocean_time) float64 0.1 0.1 0.1 0.1 0.1\n",
" k_fe_2_n_Di (ocean_time) float64 0.0001656 ... 0.0001656\n",
" k_fe_2_n_Lg (ocean_time) float64 3.975e-05 ... 3.975e-05\n",
" k_fe_2_n_Sm (ocean_time) float64 1.987e-05 ... 1.987e-05\n",
" k_fed_Di (ocean_time) float64 5e-10 5e-10 5e-10 5e-10 5e-10\n",
" k_fed_Lg (ocean_time) float64 5e-10 5e-10 5e-10 5e-10 5e-10\n",
" k_fed_Sm (ocean_time) float64 1e-10 1e-10 1e-10 1e-10 1e-10\n",
" k_ldon_bact (ocean_time) float64 5e-07 5e-07 5e-07 5e-07 5e-07\n",
" k_lith (ocean_time) float64 1.157e-11 ... 1.157e-11\n",
" k_n_inhib_Di (ocean_time) float64 1e-06 1e-06 1e-06 1e-06 1e-06\n",
" k_nh4_Di (ocean_time) float64 5e-07 5e-07 5e-07 5e-07 5e-07\n",
" k_nh4_Lg (ocean_time) float64 5e-07 5e-07 5e-07 5e-07 5e-07\n",
" k_nh4_Sm (ocean_time) float64 1e-07 1e-07 1e-07 1e-07 1e-07\n",
" k_no3_Di (ocean_time) float64 2.5e-06 2.5e-06 ... 2.5e-06\n",
" k_no3_Lg (ocean_time) float64 2.5e-06 2.5e-06 ... 2.5e-06\n",
" k_no3_Sm (ocean_time) float64 5e-07 5e-07 5e-07 5e-07 5e-07\n",
" k_no3_denit (ocean_time) float64 1e-06 1e-06 1e-06 1e-06 1e-06\n",
" k_o2 (ocean_time) float64 2e-05 2e-05 2e-05 2e-05 2e-05\n",
" k_po4_Di (ocean_time) float64 5e-08 5e-08 5e-08 5e-08 5e-08\n",
" k_po4_Lg (ocean_time) float64 5e-08 5e-08 5e-08 5e-08 5e-08\n",
" k_po4_Sm (ocean_time) float64 1e-08 1e-08 1e-08 1e-08 1e-08\n",
" k_sio4_Lg (ocean_time) float64 2e-06 2e-06 2e-06 2e-06 2e-06\n",
" kappa_eppley (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" kfe_eq_lig_hl (ocean_time) float64 1e+08 1e+08 1e+08 1e+08 1e+08\n",
" kfe_eq_lig_ll (ocean_time) float64 1e+12 1e+12 1e+12 1e+12 1e+12\n",
" ki_fescav (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" ki_hp (ocean_time) float64 1.2e-06 1.2e-06 ... 1.2e-06\n",
" ki_lgz (ocean_time) float64 1.25e-06 1.25e-06 ... 1.25e-06\n",
" ki_mdz (ocean_time) float64 1.25e-06 1.25e-06 ... 1.25e-06\n",
" ki_smz (ocean_time) float64 1.25e-06 1.25e-06 ... 1.25e-06\n",
" ktemp_bact (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" ktemp_hp (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" ktemp_lgz (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" ktemp_mdz (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" ktemp_smz (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" ktemp_vir (ocean_time) float64 0.063 0.063 0.063 0.063 0.063\n",
" latent (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" ldon (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" lgz_ipa_bact (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" lgz_ipa_det (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" lgz_ipa_diaz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" lgz_ipa_lgp (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" lgz_ipa_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" lgz_ipa_mdz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" lgz_ipa_smp (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" lgz_ipa_smz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" lith (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" lithdet (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" lwrad (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" mass_2_n (ocean_time) float64 148.7 148.7 148.7 148.7 148.7\n",
" mdz_ipa_bact (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" mdz_ipa_det (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" mdz_ipa_diaz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" mdz_ipa_lgp (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" mdz_ipa_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" mdz_ipa_mdz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" mdz_ipa_smp (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" mdz_ipa_smz (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" mswitch_hp (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" mswitch_lgz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" mswitch_mdz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" mswitch_smz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" mu_max_bact (ocean_time) float64 1.157e-05 ... 1.157e-05\n",
" nAVG (ocean_time) int32 360 360 360 360 360\n",
" nDIA (ocean_time) int32 360 360 360 360 360\n",
" nHIS (ocean_time) int32 360 360 360 360 360\n",
" nRST (ocean_time) int32 5400 5400 5400 5400 5400\n",
" n_2_n_denit (ocean_time) float64 5.9 5.9 5.9 5.9 5.9\n",
" nbact (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" ndefAVG (ocean_time) int32 360 360 360 360 360\n",
" ndefDIA (ocean_time) int32 360 360 360 360 360\n",
" ndefHIS (ocean_time) int32 360 360 360 360 360\n",
" ndet (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" ndi (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" ndtfast (ocean_time) int32 30 30 30 30 30\n",
" nh4 (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nl_tnu2 (ocean_time, tracer) float64 dask.array<chunksize=(1, 26), meta=np.ndarray>\n",
" nl_visc2 (ocean_time) float64 25.0 25.0 25.0 25.0 25.0\n",
" nlg (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nlgz (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nmdz (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" no3 (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nsm (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nsmz (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" nswitch_hp (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" nswitch_lgz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" nswitch_mdz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" nswitch_smz (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" ntimes (ocean_time) int32 3554640 3554640 ... 3554640\n",
" ntsAVG (ocean_time) int32 1 1 1 1 1\n",
" ntsDIA (ocean_time) int32 1 1 1 1 1\n",
" o2 (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" o2_2_c (ocean_time) float64 1.415 1.415 1.415 1.415 1.415\n",
" o2_2_nfix (ocean_time) float64 118.4 118.4 118.4 118.4 118.4\n",
" o2_2_nh4 (ocean_time) float64 7.375 7.375 7.375 7.375 7.375\n",
" o2_2_nitrif (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" o2_2_no3 (ocean_time) float64 9.375 9.375 9.375 9.375 9.375\n",
" o2_inhib_Di_pow (ocean_time) float64 4.0 4.0 4.0 4.0 4.0\n",
" o2_inhib_Di_sat (ocean_time) float64 0.0003 0.0003 ... 0.0003 0.0003\n",
" o2_min (ocean_time) float64 1e-06 1e-06 1e-06 1e-06 1e-06\n",
" omega (ocean_time, s_w, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 41, 362, 722), meta=np.ndarray>\n",
" p_2_n_static (ocean_time) int32 1 1 1 1 1\n",
" p_2_n_static_Di (ocean_time) float64 0.025 0.025 0.025 0.025 0.025\n",
" p_2_n_static_Lg (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" p_2_n_static_Sm (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" phi_det_lgz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_det_mdz (ocean_time) float64 0.2 0.2 0.2 0.2 0.2\n",
" phi_det_smz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_ldon_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_ldon_mdz (ocean_time) float64 0.055 0.055 0.055 0.055 0.055\n",
" phi_ldon_smz (ocean_time) float64 0.165 0.165 0.165 0.165 0.165\n",
" phi_ldon_vir (ocean_time) float64 0.55 0.55 0.55 0.55 0.55\n",
" phi_ldop_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_ldop_mdz (ocean_time) float64 0.045 0.045 0.045 0.045 0.045\n",
" phi_ldop_smz (ocean_time) float64 0.135 0.135 0.135 0.135 0.135\n",
" phi_ldop_vir (ocean_time) float64 0.45 0.45 0.45 0.45 0.45\n",
" phi_lith (ocean_time) float64 0.002 0.002 0.002 0.002 0.002\n",
" phi_nh4_lgz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_nh4_mdz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_nh4_smz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_po4_lgz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_po4_mdz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_po4_smz (ocean_time) float64 0.3 0.3 0.3 0.3 0.3\n",
" phi_sldon_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_sldon_mdz (ocean_time) float64 0.04 0.04 0.04 0.04 0.04\n",
" phi_sldon_smz (ocean_time) float64 0.12 0.12 0.12 0.12 0.12\n",
" phi_sldon_vir (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" phi_sldop_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_sldop_mdz (ocean_time) float64 0.04 0.04 0.04 0.04 0.04\n",
" phi_sldop_smz (ocean_time) float64 0.12 0.12 0.12 0.12 0.12\n",
" phi_sldop_vir (ocean_time) float64 0.4 0.4 0.4 0.4 0.4\n",
" phi_srdon_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_srdon_mdz (ocean_time) float64 0.005 0.005 0.005 0.005 0.005\n",
" phi_srdon_smz (ocean_time) float64 0.015 0.015 0.015 0.015 0.015\n",
" phi_srdon_vir (ocean_time) float64 0.05 0.05 0.05 0.05 0.05\n",
" phi_srdop_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" phi_srdop_mdz (ocean_time) float64 0.015 0.015 0.015 0.015 0.015\n",
" phi_srdop_smz (ocean_time) float64 0.045 0.045 0.045 0.045 0.045\n",
" phi_srdop_vir (ocean_time) float64 0.15 0.15 0.15 0.15 0.15\n",
" po4 (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" q_p_2_n_bact (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" q_p_2_n_lgz (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" q_p_2_n_mdz (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" q_p_2_n_smz (ocean_time) float64 0.0625 0.0625 ... 0.0625 0.0625\n",
" rain (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" rdrg (ocean_time) float64 0.0003 0.0003 ... 0.0003 0.0003\n",
" rdrg2 (ocean_time) float64 0.003 0.003 0.003 0.003 0.003\n",
" remin_eff_fedet (ocean_time) float64 0.1 0.1 0.1 0.1 0.1\n",
" rho (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" rho0 (ocean_time) float64 1.025e+03 ... 1.025e+03\n",
" rpcaco3 (ocean_time) float64 0.08805 0.08805 ... 0.08805\n",
" rplith (ocean_time) float64 0.0008176 ... 0.0008176\n",
" rpsio2 (ocean_time) float64 0.01962 0.01962 ... 0.01962\n",
" salt (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" sensible (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" shflux (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" si_2_n_max_Lg (ocean_time) float64 5.0 5.0 5.0 5.0 5.0\n",
" si_2_n_static_Lg (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" sidet (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" silg (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" sio4 (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" sldon (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" smz_ipa_bact (ocean_time) float64 0.25 0.25 0.25 0.25 0.25\n",
" smz_ipa_det (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" smz_ipa_diaz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" smz_ipa_lgp (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" smz_ipa_lgz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" smz_ipa_mdz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" smz_ipa_smp (ocean_time) float64 1.0 1.0 1.0 1.0 1.0\n",
" smz_ipa_smz (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" spherical (ocean_time) int32 1 1 1 1 1\n",
" srdon (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" ssflux (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" sss_mismatch_threshold (ocean_time) float64 8.64e+03 8.64e+03 ... 8.64e+03\n",
" sustr (ocean_time, eta_u, xi_u) float32 dask.array<chunksize=(1, 362, 721), meta=np.ndarray>\n",
" svstr (ocean_time, eta_v, xi_v) float32 dask.array<chunksize=(1, 361, 722), meta=np.ndarray>\n",
" swrad (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" temp (ocean_time, s_rho, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 40, 362, 722), meta=np.ndarray>\n",
" theta_b (ocean_time) float64 2.0 2.0 2.0 2.0 2.0\n",
" theta_s (ocean_time) float64 7.0 7.0 7.0 7.0 7.0\n",
" thetamax_Di (ocean_time) float64 0.03 0.03 0.03 0.03 0.03\n",
" thetamax_Lg (ocean_time) float64 0.05 0.05 0.05 0.05 0.05\n",
" thetamax_Sm (ocean_time) float64 0.03 0.03 0.03 0.03 0.03\n",
" thetamin (ocean_time) float64 0.002 0.002 0.002 0.002 0.002\n",
" thetamin_nolim (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" u (ocean_time, s_rho, eta_u, xi_u) float32 dask.array<chunksize=(1, 40, 362, 721), meta=np.ndarray>\n",
" ubar (ocean_time, eta_u, xi_u) float32 dask.array<chunksize=(1, 362, 721), meta=np.ndarray>\n",
" v (ocean_time, s_rho, eta_v, xi_v) float32 dask.array<chunksize=(1, 40, 361, 722), meta=np.ndarray>\n",
" vbar (ocean_time, eta_v, xi_v) float32 dask.array<chunksize=(1, 361, 722), meta=np.ndarray>\n",
" vir_Bact (ocean_time) float64 0.3819 0.3819 ... 0.3819 0.3819\n",
" vir_Di (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" vir_Lg (ocean_time) float64 0.0 0.0 0.0 0.0 0.0\n",
" vir_Sm (ocean_time) float64 0.2894 0.2894 ... 0.2894 0.2894\n",
" w (ocean_time, s_w, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 41, 362, 722), meta=np.ndarray>\n",
" wsink (ocean_time) float64 0.001157 0.001157 ... 0.001157\n",
" xl (ocean_time) float64 4.997e+06 ... 4.997e+06\n",
" z_sed (ocean_time) float64 0.1 0.1 0.1 0.1 0.1\n",
" zeta (ocean_time, eta_rho, xi_rho) float32 dask.array<chunksize=(1, 362, 722), meta=np.ndarray>\n",
" zpllgr (ocean_time) float64 0.05 0.05 0.05 0.05 0.05\n",
"Attributes:\n",
" file: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" format: netCDF-3 64bit offset file\n",
" Conventions: CF-1.4\n",
" type: ROMS/TOMS nonlinear model averages file\n",
" title: Northwest Atlantic 7km Grid\n",
" rst_file: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" his_base: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" avg_base: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" dia_base: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" grd_file: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob03T/Inputs/NWA_gr...\n",
" ini_file: /t0/scratch/shuwen/tmpdir_NWA-SZ.HCob05T/NWA-SZ.HCob05...\n",
" frc_file_01: /t0/scratch/shuwen/MERRA_Forcing//MERRA_lwrad_down_199...\n",
" frc_file_02: /t0/scratch/shuwen/MERRA_Forcing//MERRA_Pair_1999-2009...\n",
" frc_file_03: /t0/scratch/shuwen/MERRA_Forcing//MERRA_Qair_1999-2009...\n",
" frc_file_04: /t0/scratch/shuwen/MERRA_Forcing//MERRA_rain_1999-2009...\n",
" frc_file_05: /t0/scratch/shuwen/MERRA_Forcing//MERRA_swrad_1999-200...\n",
" frc_file_06: /t0/scratch/shuwen/MERRA_Forcing//MERRA_Tair_1999-2009...\n",
" frc_file_07: /t0/scratch/shuwen/MERRA_Forcing//MERRA_Uwind_1999-200...\n",
" frc_file_08: /t0/scratch/shuwen/MERRA_Forcing//MERRA_Vwind_1999-200...\n",
" frc_file_09: /t1/scratch/dujuan/NWA/Inputs/Runoff/runoff_NWA_daitre...\n",
" frc_file_10: /t1/scratch/dujuan/NWA/Inputs/Forcing/sss_monthly_clim...\n",
" frc_file_11: /t1/scratch/dujuan/NWA/Inputs/Tides/NWA_tides_otps.nc\n",
" frc_file_12: /t1/scratch/raphael/CCS1/CCS1-I/cobalt/composite_CO2_M...\n",
" frc_file_13: /t1/scratch/raphael/CCS1/CCS1-I/cobalt/Soluble_Fe_Flux...\n",
" frc_file_14: /t1/scratch/shuwen/NWA_COBALT_Inputs/Iron_sediment/iro...\n",
" frc_file_15: /t1/scratch/raphael/CCS1/CCS1-I/cobalt/Mineral_Fe_Flux...\n",
" frc_file_16: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_17: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_18: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_19: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_20: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_21: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_22: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_23: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" frc_file_24: /t1/scratch/shuwen/NWA_COBALT_Inputs/Runoff_bio/NWA_ri...\n",
" bry_file_01: /t0/scratch/shuwen/NWA_COBALT_Inputs/Bdry_SODA3.3.1/NW...\n",
" bry_file_02: /t0/scratch/shuwen/NWA_COBALT_Inputs/Bdry_bio/NWA_bdry...\n",
" clm_file_01: /t0/scratch/shuwen/NWA_COBALT_Inputs/Clim_SODA3.3.1/NW...\n",
" clm_file_02: /t0/scratch/shuwen/NWA_COBALT_Inputs/Clim_bio/NWA_clim...\n",
" nud_file: /t0/scratch/shuwen/NWA_COBALT_Inputs/NWA_nudging_coef.nc\n",
" script_file: \n",
" bpar_file: bio_Cobalt.in.full\n",
" NLM_LBC: \\nEDGE: WEST SOUTH EAST NORTH \\nzeta: ...\n",
" svn_url: https:://myroms.org/svn/src\n",
" svn_rev: exported\n",
" code_dir: /home/shuwen/ROMS\n",
" header_dir: /home/shuwen/RMANAGER/user/shuwen/NWA-SZ.HCob05T/src\n",
" header_file: nwa.h\n",
" os: Linux\n",
" cpu: x86_64\n",
" compiler_system: gfortran\n",
" compiler_command: /opt/modeling/libs/openmpi/1.8.5/gcc/4.4.7/bin/mpif90\n",
" compiler_flags: -frepack-arrays -O3 -ffast-math -I/usr/include -ffree-...\n",
" tiling: 012x016\n",
" history: Fri Feb 3 16:57:34 2017: ncks -4 -L 1 NWA-SZ.HCob05T_...\n",
" ana_file: /home/shuwen/RMANAGER/user/shuwen/NWA-SZ.HCob05T/src/a...\n",
" bio_file: ------------------------------------------------------...\n",
" CPP_options: NWA, ADD_FSOBC, ADD_M2OBC, ALBEDO, ALBEDO_CURVE, ANA_B...\n",
" NCO: \"4.6.1\""
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds = xr.open_mfdataset(\"2007/*.nc\")\n",
"ds"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' (ocean_time: 5, s_rho: 40, eta_rho: 362, xi_rho: 722)>\n",
"dask.array<concatenate, shape=(5, 40, 362, 722), dtype=float32, chunksize=(1, 40, 362, 722), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0375 -0.0125\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Dimensions without coordinates: eta_rho, xi_rho\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds[\"temp\"]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' (ocean_time: 5, s_rho: 40, eta_rho: 362, xi_rho: 722)>\n",
"dask.array<concatenate, shape=(5, 40, 362, 722), dtype=float32, chunksize=(1, 40, 362, 722), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0375 -0.0125\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Dimensions without coordinates: eta_rho, xi_rho\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp = ds[\"temp\"]\n",
"temp"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' (s_rho: 40, eta_rho: 362, xi_rho: 722)>\n",
"dask.array<getitem, shape=(40, 362, 722), dtype=float32, chunksize=(40, 362, 722), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0375 -0.0125\n",
" ocean_time datetime64[ns] 2007-02-01T01:00:00\n",
"Dimensions without coordinates: eta_rho, xi_rho\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp[0]"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' ()>\n",
"dask.array<getitem, shape=(), dtype=float32, chunksize=(), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" s_rho float64 -0.9875\n",
" ocean_time datetime64[ns] 2007-02-01T01:00:00\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"first = temp[0,0,0,0]\n",
"first"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<tr>\n",
"<td>\n",
"<table>\n",
" <thead>\n",
" <tr><td> </td><th> Array </th><th> Chunk </th></tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr><th> Bytes </th><td> 209.09 MB </td> <td> 41.82 MB </td></tr>\n",
" <tr><th> Shape </th><td> (5, 40, 362, 722) </td> <td> (1, 40, 362, 722) </td></tr>\n",
" <tr><th> Count </th><td> 15 Tasks </td><td> 5 Chunks </td></tr>\n",
" <tr><th> Type </th><td> float32 </td><td> numpy.ndarray </td></tr>\n",
" </tbody>\n",
"</table>\n",
"</td>\n",
"<td>\n",
"<svg width=\"380\" height=\"130\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
"\n",
" <!-- Horizontal lines -->\n",
" <line x1=\"0\" y1=\"0\" x2=\"25\" y2=\"0\" style=\"stroke-width:2\" />\n",
" <line x1=\"0\" y1=\"25\" x2=\"25\" y2=\"25\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Vertical lines -->\n",
" <line x1=\"0\" y1=\"0\" x2=\"0\" y2=\"25\" style=\"stroke-width:2\" />\n",
" <line x1=\"5\" y1=\"0\" x2=\"5\" y2=\"25\" />\n",
" <line x1=\"10\" y1=\"0\" x2=\"10\" y2=\"25\" />\n",
" <line x1=\"15\" y1=\"0\" x2=\"15\" y2=\"25\" />\n",
" <line x1=\"20\" y1=\"0\" x2=\"20\" y2=\"25\" />\n",
" <line x1=\"25\" y1=\"0\" x2=\"25\" y2=\"25\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Colored Rectangle -->\n",
" <polygon points=\"0.000000,0.000000 25.412617,0.000000 25.412617,25.412617 0.000000,25.412617\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
"\n",
" <!-- Text -->\n",
" <text x=\"12.706308\" y=\"45.412617\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >5</text>\n",
" <text x=\"45.412617\" y=\"12.706308\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,45.412617,12.706308)\">1</text>\n",
"\n",
"\n",
" <!-- Horizontal lines -->\n",
" <line x1=\"95\" y1=\"0\" x2=\"115\" y2=\"20\" style=\"stroke-width:2\" />\n",
" <line x1=\"95\" y1=\"60\" x2=\"115\" y2=\"80\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Vertical lines -->\n",
" <line x1=\"95\" y1=\"0\" x2=\"95\" y2=\"60\" style=\"stroke-width:2\" />\n",
" <line x1=\"115\" y1=\"20\" x2=\"115\" y2=\"80\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Colored Rectangle -->\n",
" <polygon points=\"95.000000,0.000000 115.663483,20.663483 115.663483,80.829688 95.000000,60.166205\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
"\n",
" <!-- Horizontal lines -->\n",
" <line x1=\"95\" y1=\"0\" x2=\"215\" y2=\"0\" style=\"stroke-width:2\" />\n",
" <line x1=\"115\" y1=\"20\" x2=\"235\" y2=\"20\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Vertical lines -->\n",
" <line x1=\"95\" y1=\"0\" x2=\"115\" y2=\"20\" style=\"stroke-width:2\" />\n",
" <line x1=\"215\" y1=\"0\" x2=\"235\" y2=\"20\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Colored Rectangle -->\n",
" <polygon points=\"95.000000,0.000000 215.000000,0.000000 235.663483,20.663483 115.663483,20.663483\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
"\n",
" <!-- Horizontal lines -->\n",
" <line x1=\"115\" y1=\"20\" x2=\"235\" y2=\"20\" style=\"stroke-width:2\" />\n",
" <line x1=\"115\" y1=\"80\" x2=\"235\" y2=\"80\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Vertical lines -->\n",
" <line x1=\"115\" y1=\"20\" x2=\"115\" y2=\"80\" style=\"stroke-width:2\" />\n",
" <line x1=\"235\" y1=\"20\" x2=\"235\" y2=\"80\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Colored Rectangle -->\n",
" <polygon points=\"115.663483,20.663483 235.663483,20.663483 235.663483,80.829688 115.663483,80.829688\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
"\n",
" <!-- Text -->\n",
" <text x=\"175.663483\" y=\"100.829688\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >722</text>\n",
" <text x=\"255.663483\" y=\"50.746586\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(-90,255.663483,50.746586)\">362</text>\n",
" <text x=\"95.331742\" y=\"90.497947\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(45,95.331742,90.497947)\">40</text>\n",
"</svg>\n",
"</td>\n",
"</tr>\n",
"</table>"
],
"text/plain": [
"dask.array<concatenate, shape=(5, 40, 362, 722), dtype=float32, chunksize=(1, 40, 362, 722), chunktype=numpy.ndarray>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp.data"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' (ocean_time: 5, s_rho: 40, eta_rho: 362, xi_rho: 722)>\n",
"dask.array<concatenate, shape=(5, 40, 362, 722), dtype=float32, chunksize=(1, 40, 362, 722), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0375 -0.0125\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Dimensions without coordinates: eta_rho, xi_rho\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'temp' (ocean_time: 5)>\n",
"dask.array<getitem, shape=(5,), dtype=float32, chunksize=(1,), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" s_rho float64 -0.9875\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Attributes:\n",
" long_name: time-averaged potential temperature\n",
" units: Celsius\n",
" time: ocean_time\n",
" coordinates: lon_rho lat_rho s_rho ocean_time\n",
" field: temperature, scalar, series"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp[:,0,0,0]"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<table>\n",
"<tr>\n",
"<td>\n",
"<table>\n",
" <thead>\n",
" <tr><td> </td><th> Array </th><th> Chunk </th></tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr><th> Bytes </th><td> 20 B </td> <td> 4 B </td></tr>\n",
" <tr><th> Shape </th><td> (5,) </td> <td> (1,) </td></tr>\n",
" <tr><th> Count </th><td> 20 Tasks </td><td> 5 Chunks </td></tr>\n",
" <tr><th> Type </th><td> float32 </td><td> numpy.ndarray </td></tr>\n",
" </tbody>\n",
"</table>\n",
"</td>\n",
"<td>\n",
"<svg width=\"170\" height=\"92\" style=\"stroke:rgb(0,0,0);stroke-width:1\" >\n",
"\n",
" <!-- Horizontal lines -->\n",
" <line x1=\"0\" y1=\"0\" x2=\"120\" y2=\"0\" style=\"stroke-width:2\" />\n",
" <line x1=\"0\" y1=\"42\" x2=\"120\" y2=\"42\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Vertical lines -->\n",
" <line x1=\"0\" y1=\"0\" x2=\"0\" y2=\"42\" style=\"stroke-width:2\" />\n",
" <line x1=\"24\" y1=\"0\" x2=\"24\" y2=\"42\" />\n",
" <line x1=\"48\" y1=\"0\" x2=\"48\" y2=\"42\" />\n",
" <line x1=\"72\" y1=\"0\" x2=\"72\" y2=\"42\" />\n",
" <line x1=\"96\" y1=\"0\" x2=\"96\" y2=\"42\" />\n",
" <line x1=\"120\" y1=\"0\" x2=\"120\" y2=\"42\" style=\"stroke-width:2\" />\n",
"\n",
" <!-- Colored Rectangle -->\n",
" <polygon points=\"0.000000,0.000000 120.000000,0.000000 120.000000,42.009890 0.000000,42.009890\" style=\"fill:#ECB172A0;stroke-width:0\"/>\n",
"\n",
" <!-- Text -->\n",
" <text x=\"60.000000\" y=\"62.009890\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" >5</text>\n",
" <text x=\"140.000000\" y=\"21.004945\" font-size=\"1.0rem\" font-weight=\"100\" text-anchor=\"middle\" transform=\"rotate(0,140.000000,21.004945)\">1</text>\n",
"</svg>\n",
"</td>\n",
"</tr>\n",
"</table>"
],
"text/plain": [
"dask.array<getitem, shape=(5,), dtype=float32, chunksize=(1,), chunktype=numpy.ndarray>"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp[:,0,0,0].data"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 's_rho' (s_rho: 40)>\n",
"array([-0.9875, -0.9625, -0.9375, -0.9125, -0.8875, -0.8625, -0.8375, -0.8125,\n",
" -0.7875, -0.7625, -0.7375, -0.7125, -0.6875, -0.6625, -0.6375, -0.6125,\n",
" -0.5875, -0.5625, -0.5375, -0.5125, -0.4875, -0.4625, -0.4375, -0.4125,\n",
" -0.3875, -0.3625, -0.3375, -0.3125, -0.2875, -0.2625, -0.2375, -0.2125,\n",
" -0.1875, -0.1625, -0.1375, -0.1125, -0.0875, -0.0625, -0.0375, -0.0125])\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0625 -0.0375 -0.0125\n",
"Attributes:\n",
" long_name: S-coordinate at RHO-points\n",
" valid_min: -1.0\n",
" valid_max: 0.0\n",
" positive: up\n",
" standard_name: ocean_s_coordinate_g2\n",
" formula_terms: s: s_rho C: Cs_r eta: zeta depth: h depth_c: hc\n",
" field: s_rho, scalar"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp.coords[\"s_rho\"]"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 'ocean_time' (ocean_time: 5)>\n",
"array(['2007-02-01T01:00:00.000000000', '2007-02-02T01:00:00.000000000',\n",
" '2007-02-03T01:00:00.000000000', '2007-02-04T01:00:00.000000000',\n",
" '2007-02-05T01:00:00.000000000'], dtype='datetime64[ns]')\n",
"Coordinates:\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00\n",
"Attributes:\n",
" long_name: averaged time since initialization\n",
" field: time, scalar, series"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"temp.coords[\"ocean_time\"]"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<xarray.DataArray 's_rho' (s_rho: 40)>\n",
"array([-0.9875, -0.9625, -0.9375, -0.9125, -0.8875, -0.8625, -0.8375, -0.8125,\n",
" -0.7875, -0.7625, -0.7375, -0.7125, -0.6875, -0.6625, -0.6375, -0.6125,\n",
" -0.5875, -0.5625, -0.5375, -0.5125, -0.4875, -0.4625, -0.4375, -0.4125,\n",
" -0.3875, -0.3625, -0.3375, -0.3125, -0.2875, -0.2625, -0.2375, -0.2125,\n",
" -0.1875, -0.1625, -0.1375, -0.1125, -0.0875, -0.0625, -0.0375, -0.0125])\n",
"Coordinates:\n",
" * s_rho (s_rho) float64 -0.9875 -0.9625 -0.9375 ... -0.0625 -0.0375 -0.0125\n",
"Attributes:\n",
" long_name: S-coordinate at RHO-points\n",
" valid_min: -1.0\n",
" valid_max: 0.0\n",
" positive: up\n",
" standard_name: ocean_s_coordinate_g2\n",
" formula_terms: s: s_rho C: Cs_r eta: zeta depth: h depth_c: hc\n",
" field: s_rho, scalar"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds[\"s_rho\"]"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Coordinates:\n",
" * s_w (s_w) float64 -1.0 -0.975 -0.95 -0.925 ... -0.05 -0.025 0.0\n",
" * ocean_time (ocean_time) datetime64[ns] 2007-02-01T01:00:00 ... 2007-02-05T01:00:00"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds[\"AKv\"].coords"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"OrderedDict([('long_name',\n",
" 'zpllgr: Phytoplankton light limitation/growth rate')])"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds[var].attrs"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Units that have lat/lon"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"AKt - time-averaged temperature vertical diffusion coefficient - meter2 second-1\n",
"AKv - time-averaged vertical viscosity coefficient - meter2 second-1\n",
"Hsbl - time-averaged depth of oceanic surface boundary layer - meter\n",
"Pair - time-averaged surface air pressure - Pascal\n",
"Runoff - time-averaged Fresh water runoff from land - kg/s/m^2\n",
"Uwind - time-averaged surface u-wind component - meter second-1\n",
"Vwind - time-averaged surface v-wind component - meter second-1\n",
"alk - time-averaged Alkalinity - mol/kg\n",
"cadet_arag - time-averaged Detrital CaCO3 - mol/kg\n",
"cadet_calc - time-averaged Detrital CaCO3 - mol/kg\n",
"dic - time-averaged Dissolved Inorganic Carbon - mol/kg\n",
"evaporation - time-averaged evaporation rate - kilogram meter-2 second-1\n",
"latent - time-averaged net latent heat flux - watt meter-2\n",
"ldon - time-averaged labile DON - mol/kg\n",
"lith - time-averaged Lithogenic Aluminosilicate - g/kg\n",
"lithdet - time-averaged lithdet - g/kg\n",
"lwrad - time-averaged net longwave radiation flux - watt meter-2\n",
"nbact - time-averaged bacterial - mol/kg\n",
"ndet - time-averaged ndet - mol/kg\n",
"ndi - time-averaged Diazotroph Nitrogen - mol/kg\n",
"nh4 - time-averaged Ammonia - mol/kg\n",
"nlg - time-averaged Large Phytoplankton Nitrogen - mol/kg\n",
"nlgz - time-averaged large Zooplankton Nitrogen - mol/kg\n",
"nmdz - time-averaged Medium-sized zooplankton Nitrogen - mol/kg\n",
"no3 - time-averaged Nitrate - mol/kg\n",
"nsm - time-averaged Small Phytoplankton Nitrogen - mol/kg\n",
"nsmz - time-averaged Small Zooplankton Nitrogen - mol/kg\n",
"o2 - time-averaged Oxygen - mol/kg\n",
"omega - time-averaged S-coordinate vertical momentum component - meter3 second-1\n",
"po4 - time-averaged Phosphate - mol/kg\n",
"rain - time-averaged rain fall rate - kilogram meter-2 second-1\n",
"rho - time-averaged density anomaly - kilogram meter-3\n",
"sensible - time-averaged net sensible heat flux - watt meter-2\n",
"shflux - time-averaged surface net heat flux - watt meter-2\n",
"sidet - time-averaged Detrital Silicon - mol/kg\n",
"silg - time-averaged Large Phytoplankton Silicon - mol/kg\n",
"sio4 - time-averaged Silicate - mol/kg\n",
"sldon - time-averaged Semilabile DON - mol/kg\n",
"srdon - time-averaged Semi-Refractory DON - mol/kg\n",
"ssflux - time-averaged surface net salt flux, (E-P)*SALT - meter second-1\n",
"sustr - time-averaged surface u-momentum stress - newton meter-2\n",
"svstr - time-averaged surface v-momentum stress - newton meter-2\n",
"swrad - time-averaged solar shortwave radiation flux - watt meter-2\n",
"temp - time-averaged potential temperature - Celsius\n",
"u - time-averaged u-momentum component - meter second-1\n",
"ubar - time-averaged vertically integrated u-momentum component - meter second-1\n",
"v - time-averaged v-momentum component - meter second-1\n",
"vbar - time-averaged vertically integrated v-momentum component - meter second-1\n",
"w - time-averaged vertical momentum component - meter second-1\n",
"zeta - time-averaged free-surface - meter\n"
]
}
],
"source": [
"for var in ds.variables:\n",
" subset = ds[var]\n",
" try:\n",
" if \"lat\" in subset.attrs[\"coordinates\"] and \"lon\" in subset.attrs[\"coordinates\"]: \n",
" print(f'{var} - {subset.attrs[\"long_name\"]} - {subset.attrs[\"units\"]}')\n",
" except KeyError:\n",
" pass"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "IOOS Notebooks",
"language": "python",
"name": "ioos"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment