In [1]:
%reset
from netCDF4 import Dataset
import numpy as np
import matplotlib
import matplotlib.pyplot as plt
import xarray as xr
from xgcm import Grid
import numpy as np
import cmocean
import warnings
warnings.filterwarnings("ignore")
In [2]:
#%cd /glade/scratch/eyankovsky/Backscatter_tests_with_Scott/Half_degree/Bachman_EBT_test4/

%cd /glade/scratch/eyankovsky/Backscatter_tests_with_Scott/Half_degree/Jansen_EBT_test2/

av = xr.open_dataset('averages_00030002.nc', decode_times=False)  
MEKEds = xr.open_dataset('MEKEaverages_00031502.nc', decode_times=False)


fs = xr.open_dataset('static.nc', decode_times=False)

/glade/scratch/eyankovsky/Backscatter_tests_with_Scott/Half_degree/Jansen_EBT_test2
In [ ]:
In [3]:
#Read in various grid variables from the static file:

depth = np.array(fs['depth_ocean'])
area_t=np.array(fs['area_t'])

xh = fs.xh.values
yh = fs.yh.values


MEKE_src_adv      = MEKEds['MEKE_src_adv'].mean('time')
MEKE_src_mom_K4   = MEKEds['MEKE_src_mom_K4'].mean('time')
MEKE_src_btm_drag = MEKEds['MEKE_src_btm_drag'].mean('time')
MEKE_src_GM       = MEKEds['MEKE_src_GM'].mean('time')
MEKE_src_mom_lp   = MEKEds['MEKE_src_mom_lp'].mean('time')
MEKE_src_mom_bh   = MEKEds['MEKE_src_mom_bh'].mean('time')

MEKE = MEKEds['MEKE'].mean('time')
Sum = MEKE_src_adv + MEKE_src_mom_K4 + MEKE_src_btm_drag + MEKE_src_GM + MEKE_src_mom_lp + MEKE_src_mom_bh

print(Sum.shape)

np.nanmean(MEKE_src_mom_lp)

(280, 120)
Out[3]:
-1.0985548e-08
In [ ]:
In [4]:
import matplotlib.colors as colors

class MidpointNormalize(colors.Normalize):
    def __init__(self, vmin=None, vmax=None, midpoint=None, clip=False):
        self.midpoint = midpoint
        colors.Normalize.__init__(self, vmin, vmax, clip)

    def __call__(self, value, clip=None):
        # I'm ignoring masked values and all kinds of edge cases to make a
        # simple example...
        x, y = [self.vmin, self.midpoint, self.vmax], [0, 0.5, 1]
        return np.ma.masked_array(np.interp(value, x, y))
In [5]:
plt.figure(figsize=(14,10),dpi=100)


cmin=-1e-7; cmax=1e-7;


plt.subplot(2,4,1)
plt.pcolormesh(xh,yh,Sum,cmap='RdBu_r',vmin=cmin/100,vmax=cmax/100)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('Error')
cbar=plt.colorbar()

plt.subplot(2,4,2)
plt.pcolormesh(xh,yh,MEKE_src_GM,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_GM')
cbar=plt.colorbar()

plt.subplot(2,4,3)
plt.pcolormesh(xh,yh,MEKE_src_mom_lp,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_mom_lp')
cbar=plt.colorbar()

plt.subplot(2,4,4)
plt.pcolormesh(xh,yh,MEKE_src_mom_bh,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_mom_bh')
cbar=plt.colorbar()


plt.subplot(2,4,5)
plt.pcolormesh(xh,yh,MEKE,cmap='magma',vmin=0,vmax=0.005)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE')
cbar=plt.colorbar()

plt.subplot(2,4,6)
plt.pcolormesh(xh,yh,MEKE_src_btm_drag,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_btm_drag')
cbar=plt.colorbar()
plt.subplot(2,4,7)

plt.pcolormesh(xh,yh,MEKE_src_mom_K4,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_mom_K4')
cbar=plt.colorbar()

plt.subplot(2,4,8)
plt.pcolormesh(xh,yh,MEKE_src_adv,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE_src_adv')
cbar=plt.colorbar()
No description has been provided for this image
In [21]:
fig, ax = plt.subplots(nrows=3, ncols=3, sharex=True, sharey=True, figsize=(8, 10),dpi=150)

cmin=-5e-8; cmax=5e-8;

ax[0,0]=plt.subplot(2,3,1)
plt.pcolormesh(xh,yh,MEKE,cmap='Reds',vmin=0,vmax=0.025)
plt.gca().set_aspect('equal', adjustable='box')
plt.title('MEKE [m$^2$/s$^2$]')
plt.ylabel('Latitude')
#plt.set.yticklabels([])
cbar=plt.colorbar()
cbar.formatter.set_powerlimits((0, 0))


# plt.subplot(2,3,2)
# plt.pcolormesh(xh,yh,MEKE_src_GM,cmap='RdBu_r',vmin=cmin,vmax=cmax)
# plt.gca().set_aspect('equal', adjustable='box')
# plt.title('MEKE_src_GM')
# cbar=plt.colorbar()

ax[0,1]=plt.subplot(2,3,2)
plt.pcolormesh(xh,yh,MEKE_src_mom_lp,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title(r'$\.e_{\rm BScat}$ [m$^2$/s$^3$]')
cbar=plt.colorbar()

ax[0,2]=plt.subplot(2,3,3)
plt.pcolormesh(xh,yh,MEKE_src_mom_bh,cmap='RdBu_r',vmin=cmin,vmax=cmax)
plt.gca().set_aspect('equal', adjustable='box')
plt.title(r'$\.e_{\rm Smag}$ [m$^2$/s$^3$]')

cbar=plt.colorbar()
ax[1,0]=plt.subplot(2,3,4)
plt.pcolormesh(xh,yh,MEKE_src_adv,cmap='RdBu_r',vmin=cmin/20,vmax=cmax/20)
plt.gca().set_aspect('equal', adjustable='box')
plt.title(r'$\.e_{\rm adv}$ [m$^2$/s$^3$]')
plt.xlabel('Longitude')
plt.ylabel('Latitude')

cbar=plt.colorbar()

ax[1,1]=plt.subplot(2,3,5)
plt.pcolormesh(xh,yh,MEKE_src_btm_drag,cmap='RdBu_r',vmin=cmin/20,vmax=cmax/20)
plt.gca().set_aspect('equal', adjustable='box')
plt.title(r'$\.e_{\rm diss}$ [m$^2$/s$^3$]')
plt.xlabel('Longitude')

cbar=plt.colorbar()



ax[1,2]=plt.subplot(2,3,6)
plt.pcolormesh(xh,yh,Sum,cmap='RdBu_r',vmin=cmin/20,vmax=cmax/20)
plt.gca().set_aspect('equal', adjustable='box')
plt.title(r'$\.e_{\rm residual}$ [m$^2$/s$^3$]')
cbar=plt.colorbar()
plt.xlabel('Longitude')

ax[0,0].tick_params(labelbottom=False)
ax[0,1].tick_params(labelbottom=False)
ax[0,2].tick_params(labelbottom=False)
ax[0,1].tick_params(labelleft=False)
ax[0,2].tick_params(labelleft=False)
ax[1,1].tick_params(labelleft=False)
ax[1,2].tick_params(labelleft=False)
No description has been provided for this image
In [10]:
MEKEds
Out[10]:
/* CSS stylesheet for displaying xarray objects in jupyterlab. * */ :root { --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1)); --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54)); --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38)); --xr-border-color: var(--jp-border-color2, #e0e0e0); --xr-disabled-color: var(--jp-layout-color3, #bdbdbd); --xr-background-color: var(--jp-layout-color0, white); --xr-background-color-row-even: var(--jp-layout-color1, white); --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee); } html[theme=dark], body.vscode-dark { --xr-font-color0: rgba(255, 255, 255, 1); --xr-font-color2: rgba(255, 255, 255, 0.54); --xr-font-color3: rgba(255, 255, 255, 0.38); --xr-border-color: #1F1F1F; --xr-disabled-color: #515151; --xr-background-color: #111111; --xr-background-color-row-even: #111111; --xr-background-color-row-odd: #313131; } .xr-wrap { display: block; min-width: 300px; max-width: 700px; } .xr-text-repr-fallback { /* fallback to plain text repr when CSS is not injected (untrusted notebook) */ display: none; } .xr-header { padding-top: 6px; padding-bottom: 6px; margin-bottom: 4px; border-bottom: solid 1px var(--xr-border-color); } .xr-header > div, .xr-header > ul { display: inline; margin-top: 0; margin-bottom: 0; } .xr-obj-type, .xr-array-name { margin-left: 2px; margin-right: 10px; } .xr-obj-type { color: var(--xr-font-color2); } .xr-sections { padding-left: 0 !important; display: grid; grid-template-columns: 150px auto auto 1fr 20px 20px; } .xr-section-item { display: contents; } .xr-section-item input { display: none; } .xr-section-item input + label { color: var(--xr-disabled-color); } .xr-section-item input:enabled + label { cursor: pointer; color: var(--xr-font-color2); } .xr-section-item input:enabled + label:hover { color: var(--xr-font-color0); } .xr-section-summary { grid-column: 1; color: var(--xr-font-color2); font-weight: 500; } .xr-section-summary > span { display: inline-block; padding-left: 0.5em; } .xr-section-summary-in:disabled + label { color: var(--xr-font-color2); } .xr-section-summary-in + label:before { display: inline-block; content: '►'; font-size: 11px; width: 15px; text-align: center; } .xr-section-summary-in:disabled + label:before { color: var(--xr-disabled-color); } .xr-section-summary-in:checked + label:before { content: '▼'; } .xr-section-summary-in:checked + label > span { display: none; } .xr-section-summary, .xr-section-inline-details { padding-top: 4px; padding-bottom: 4px; } .xr-section-inline-details { grid-column: 2 / -1; } .xr-section-details { display: none; grid-column: 1 / -1; margin-bottom: 5px; } .xr-section-summary-in:checked ~ .xr-section-details { display: contents; } .xr-array-wrap { grid-column: 1 / -1; display: grid; grid-template-columns: 20px auto; } .xr-array-wrap > label { grid-column: 1; vertical-align: top; } .xr-preview { color: var(--xr-font-color3); } .xr-array-preview, .xr-array-data { padding: 0 5px !important; grid-column: 2; } .xr-array-data, .xr-array-in:checked ~ .xr-array-preview { display: none; } .xr-array-in:checked ~ .xr-array-data, .xr-array-preview { display: inline-block; } .xr-dim-list { display: inline-block !important; list-style: none; padding: 0 !important; margin: 0; } .xr-dim-list li { display: inline-block; padding: 0; margin: 0; } .xr-dim-list:before { content: '('; } .xr-dim-list:after { content: ')'; } .xr-dim-list li:not(:last-child):after { content: ','; padding-right: 5px; } .xr-has-index { font-weight: bold; } .xr-var-list, .xr-var-item { display: contents; } .xr-var-item > div, .xr-var-item label, .xr-var-item > .xr-var-name span { background-color: var(--xr-background-color-row-even); margin-bottom: 0; } .xr-var-item > .xr-var-name:hover span { padding-right: 5px; } .xr-var-list > li:nth-child(odd) > div, .xr-var-list > li:nth-child(odd) > label, .xr-var-list > li:nth-child(odd) > .xr-var-name span { background-color: var(--xr-background-color-row-odd); } .xr-var-name { grid-column: 1; } .xr-var-dims { grid-column: 2; } .xr-var-dtype { grid-column: 3; text-align: right; color: var(--xr-font-color2); } .xr-var-preview { grid-column: 4; } .xr-var-name, .xr-var-dims, .xr-var-dtype, .xr-preview, .xr-attrs dt { white-space: nowrap; overflow: hidden; text-overflow: ellipsis; padding-right: 10px; } .xr-var-name:hover, .xr-var-dims:hover, .xr-var-dtype:hover, .xr-attrs dt:hover { overflow: visible; width: auto; z-index: 1; } .xr-var-attrs, .xr-var-data { display: none; background-color: var(--xr-background-color) !important; padding-bottom: 5px !important; } .xr-var-attrs-in:checked ~ .xr-var-attrs, .xr-var-data-in:checked ~ .xr-var-data { display: block; } .xr-var-data > table { float: right; } .xr-var-name span, .xr-var-data, .xr-attrs { padding-left: 25px !important; } .xr-attrs, .xr-var-attrs, .xr-var-data { grid-column: 1 / -1; } dl.xr-attrs { padding: 0; margin: 0; display: grid; grid-template-columns: 125px auto; } .xr-attrs dt, .xr-attrs dd { padding: 0; margin: 0; float: left; padding-right: 10px; width: auto; } .xr-attrs dt { font-weight: normal; grid-column: 1; } .xr-attrs dt:hover span { display: inline-block; background: var(--xr-background-color); padding-right: 10px; } .xr-attrs dd { grid-column: 2; white-space: pre-wrap; word-break: break-all; } .xr-icon-database, .xr-icon-file-text2 { display: inline-block; vertical-align: middle; width: 1em; height: 1.5em !important; stroke-width: 0; stroke: currentColor; fill: currentColor; } 
<xarray.Dataset>
Dimensions:            (nv: 2, time: 100, xh: 120, xq: 121, yh: 280, yq: 281)
Coordinates:
  * xh                 (xh) float64 0.25 0.75 1.25 1.75 ... 58.75 59.25 59.75
  * yh                 (yh) float64 -69.75 -69.25 -68.75 ... 68.75 69.25 69.75
  * time               (time) float64 3.15e+04 3.151e+04 ... 3.199e+04 3.2e+04
  * nv                 (nv) float64 1.0 2.0
  * xq                 (xq) float64 0.0 0.5 1.0 1.5 2.0 ... 58.5 59.0 59.5 60.0
  * yq                 (yq) float64 -70.0 -69.5 -69.0 -68.5 ... 69.0 69.5 70.0
Data variables:
    MEKE               (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_KH            (time, yh, xh) float32 ...
    MEKE_KU            (time, yh, xh) float32 ...
    MEKE_AU            (time, yh, xh) float32 ...
    MEKE_Ue            (time, yh, xh) float32 ...
    MEKE_Ub            (time, yh, xh) float32 ...
    MEKE_Ut            (time, yh, xh) float32 ...
    MEKE_src           (time, yh, xh) float32 ...
    MEKE_decay         (time, yh, xh) float32 ...
    MEKE_GM_src        (time, yh, xh) float32 ...
    MEKE_mom_src       (time, yh, xh) float32 ...
    MEKE_GME_snk       (time, yh, xh) float32 ...
    MEKE_Le            (time, yh, xh) float32 ...
    MEKE_Lrhines       (time, yh, xh) float32 ...
    MEKE_Leady         (time, yh, xh) float32 ...
    MEKE_gamma_b       (time, yh, xh) float32 ...
    MEKE_gamma_t       (time, yh, xh) float32 ...
    KHMEKE_u           (time, yh, xq) float32 ...
    KHMEKE_v           (time, yq, xh) float32 ...
    MEKE_src_adv       (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_src_mom_K4    (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_src_btm_drag  (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_src_GM        (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_src_mom_lp    (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    MEKE_src_mom_bh    (time, yh, xh) float32 nan nan nan nan ... nan nan nan
    average_T1         (time) float64 3.15e+04 3.150e+04 ... 3.199e+04 3.2e+04
    average_T2         (time) float64 3.150e+04 3.151e+04 ... 3.2e+04 3.2e+04
    average_DT         (time) float64 5.0 5.0 5.0 5.0 5.0 ... 5.0 5.0 5.0 5.0
    time_bnds          (time, nv) float64 3.15e+04 3.150e+04 ... 3.2e+04 3.2e+04
Attributes:
    filename:          MEKEaverages_00031502.nc
    title:             NeverWorld2
    associated_files:  area_t: static.nc
    grid_type:         regular
    grid_tile:         N/A
xarray.Dataset
  • Dimensions:
    • nv: 2
    • time: 100
    • xh: 120
    • xq: 121
    • yh: 280
    • yq: 281
  • Coordinates: (6)
    • xh
      (xh)
      float64
      0.25 0.75 1.25 ... 59.25 59.75
      long_name : h point nominal longitude units : degrees_east cartesian_axis : X
      array([ 0.25,  0.75,  1.25,  1.75,  2.25,  2.75,  3.25,  3.75,  4.25,  4.75,
        5.25,  5.75,  6.25,  6.75,  7.25,  7.75,  8.25,  8.75,  9.25,  9.75,
       10.25, 10.75, 11.25, 11.75, 12.25, 12.75, 13.25, 13.75, 14.25, 14.75,
       15.25, 15.75, 16.25, 16.75, 17.25, 17.75, 18.25, 18.75, 19.25, 19.75,
       20.25, 20.75, 21.25, 21.75, 22.25, 22.75, 23.25, 23.75, 24.25, 24.75,
       25.25, 25.75, 26.25, 26.75, 27.25, 27.75, 28.25, 28.75, 29.25, 29.75,
       30.25, 30.75, 31.25, 31.75, 32.25, 32.75, 33.25, 33.75, 34.25, 34.75,
       35.25, 35.75, 36.25, 36.75, 37.25, 37.75, 38.25, 38.75, 39.25, 39.75,
       40.25, 40.75, 41.25, 41.75, 42.25, 42.75, 43.25, 43.75, 44.25, 44.75,
       45.25, 45.75, 46.25, 46.75, 47.25, 47.75, 48.25, 48.75, 49.25, 49.75,
       50.25, 50.75, 51.25, 51.75, 52.25, 52.75, 53.25, 53.75, 54.25, 54.75,
       55.25, 55.75, 56.25, 56.75, 57.25, 57.75, 58.25, 58.75, 59.25, 59.75])
    • yh
      (yh)
      float64
      -69.75 -69.25 ... 69.25 69.75
      long_name : h point nominal latitude units : degrees_north cartesian_axis : Y
      array([-69.75, -69.25, -68.75, ...,  68.75,  69.25,  69.75])
    • time
      (time)
      float64
      3.15e+04 3.151e+04 ... 3.2e+04
      long_name : time units : days since 0001-01-01 00:00:00 cartesian_axis : T calendar_type : THIRTY_DAY_MONTHS calendar : THIRTY_DAY_MONTHS bounds : time_bnds
      array([31502.5, 31507.5, 31512.5, 31517.5, 31522.5, 31527.5, 31532.5, 31537.5,
       31542.5, 31547.5, 31552.5, 31557.5, 31562.5, 31567.5, 31572.5, 31577.5,
       31582.5, 31587.5, 31592.5, 31597.5, 31602.5, 31607.5, 31612.5, 31617.5,
       31622.5, 31627.5, 31632.5, 31637.5, 31642.5, 31647.5, 31652.5, 31657.5,
       31662.5, 31667.5, 31672.5, 31677.5, 31682.5, 31687.5, 31692.5, 31697.5,
       31702.5, 31707.5, 31712.5, 31717.5, 31722.5, 31727.5, 31732.5, 31737.5,
       31742.5, 31747.5, 31752.5, 31757.5, 31762.5, 31767.5, 31772.5, 31777.5,
       31782.5, 31787.5, 31792.5, 31797.5, 31802.5, 31807.5, 31812.5, 31817.5,
       31822.5, 31827.5, 31832.5, 31837.5, 31842.5, 31847.5, 31852.5, 31857.5,
       31862.5, 31867.5, 31872.5, 31877.5, 31882.5, 31887.5, 31892.5, 31897.5,
       31902.5, 31907.5, 31912.5, 31917.5, 31922.5, 31927.5, 31932.5, 31937.5,
       31942.5, 31947.5, 31952.5, 31957.5, 31962.5, 31967.5, 31972.5, 31977.5,
       31982.5, 31987.5, 31992.5, 31997.5])
    • nv
      (nv)
      float64
      1.0 2.0
      long_name : vertex number units : none cartesian_axis : N
      array([1., 2.])
    • xq
      (xq)
      float64
      0.0 0.5 1.0 1.5 ... 59.0 59.5 60.0
      long_name : q point nominal longitude units : degrees_east cartesian_axis : X
      array([ 0. ,  0.5,  1. ,  1.5,  2. ,  2.5,  3. ,  3.5,  4. ,  4.5,  5. ,  5.5,
        6. ,  6.5,  7. ,  7.5,  8. ,  8.5,  9. ,  9.5, 10. , 10.5, 11. , 11.5,
       12. , 12.5, 13. , 13.5, 14. , 14.5, 15. , 15.5, 16. , 16.5, 17. , 17.5,
       18. , 18.5, 19. , 19.5, 20. , 20.5, 21. , 21.5, 22. , 22.5, 23. , 23.5,
       24. , 24.5, 25. , 25.5, 26. , 26.5, 27. , 27.5, 28. , 28.5, 29. , 29.5,
       30. , 30.5, 31. , 31.5, 32. , 32.5, 33. , 33.5, 34. , 34.5, 35. , 35.5,
       36. , 36.5, 37. , 37.5, 38. , 38.5, 39. , 39.5, 40. , 40.5, 41. , 41.5,
       42. , 42.5, 43. , 43.5, 44. , 44.5, 45. , 45.5, 46. , 46.5, 47. , 47.5,
       48. , 48.5, 49. , 49.5, 50. , 50.5, 51. , 51.5, 52. , 52.5, 53. , 53.5,
       54. , 54.5, 55. , 55.5, 56. , 56.5, 57. , 57.5, 58. , 58.5, 59. , 59.5,
       60. ])
    • yq
      (yq)
      float64
      -70.0 -69.5 -69.0 ... 69.5 70.0
      long_name : q point nominal latitude units : degrees_north cartesian_axis : Y
      array([-70. , -69.5, -69. , ...,  69. ,  69.5,  70. ])
  • Data variables: (29)
    • MEKE
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : Mesoscale Eddy Kinetic Energy units : m2 s-2 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[         nan,          nan, ...,          nan,          nan],
        [         nan, 1.325917e-05, ..., 4.736797e-06,          nan],
        ...,
        [         nan, 9.418671e-06, ..., 1.076484e-06,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 5.507454e-05, ..., 5.342876e-06,          nan],
        ...,
        [         nan, 1.961855e-05, ..., 1.685546e-06,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       ...,

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 4.079768e-04, ..., 3.010977e-06,          nan],
        ...,
        [         nan, 6.825547e-04, ..., 2.096902e-06,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 2.347056e-04, ..., 3.283013e-06,          nan],
        ...,
        [         nan, 3.990509e-04, ..., 2.079619e-06,          nan],
        [         nan,          nan, ...,          nan,          nan]]],
      dtype=float32)
    • MEKE_KH
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived diffusivity units : m2 s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_KU
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived lateral viscosity units : m2 s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_AU
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived lateral biharmonic viscosity units : m4 s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Ue
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived eddy-velocity scale units : m s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Ub
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived bottom eddy-velocity scale units : m s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Ut
      (time, yh, xh)
      float32
      ...
      long_name : MEKE derived barotropic eddy-velocity scale units : m s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_src
      (time, yh, xh)
      float32
      ...
      long_name : MEKE energy source units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_decay
      (time, yh, xh)
      float32
      ...
      long_name : MEKE decay rate units : s-1 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_GM_src
      (time, yh, xh)
      float32
      ...
      long_name : MEKE energy available from thickness mixing units : W m-2 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_mom_src
      (time, yh, xh)
      float32
      ...
      long_name : MEKE energy available from momentum units : W m-2 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_GME_snk
      (time, yh, xh)
      float32
      ...
      long_name : MEKE energy lost to GME backscatter units : W m-2 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Le
      (time, yh, xh)
      float32
      ...
      long_name : Eddy mixing length used in the MEKE derived eddy diffusivity units : m cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Lrhines
      (time, yh, xh)
      float32
      ...
      long_name : Rhines length scale used in the MEKE derived eddy diffusivity units : m cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_Leady
      (time, yh, xh)
      float32
      ...
      long_name : Eady length scale used in the MEKE derived eddy diffusivity units : m cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_gamma_b
      (time, yh, xh)
      float32
      ...
      long_name : Ratio of bottom-projected eddy velocity to column-mean eddy velocity units : nondim cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • MEKE_gamma_t
      (time, yh, xh)
      float32
      ...
      long_name : Ratio of barotropic eddy velocity to column-mean eddy velocity units : nondim cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      [3360000 values with dtype=float32]
    • KHMEKE_u
      (time, yh, xq)
      float32
      ...
      long_name : Zonal diffusivity of MEKE units : m2 s-1 cell_methods : yh:mean xq:point time: mean time_avg_info : average_T1,average_T2,average_DT interp_method : none
      [3388000 values with dtype=float32]
    • KHMEKE_v
      (time, yq, xh)
      float32
      ...
      long_name : Meridional diffusivity of MEKE units : m2 s-1 cell_methods : yq:point xh:mean time: mean time_avg_info : average_T1,average_T2,average_DT interp_method : none
      [3372000 values with dtype=float32]
    • MEKE_src_adv
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the horizontal advection of MEKE units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[          nan,           nan, ...,           nan,           nan],
        [          nan,  2.850493e-11, ...,  2.928214e-12,           nan],
        ...,
        [          nan, -1.055516e-12, ...,  2.737292e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan,  2.326759e-11, ...,  2.771116e-12,           nan],
        ...,
        [          nan, -3.440152e-12, ...,  2.704325e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       ...,

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -1.699153e-10, ...,  1.381902e-12,           nan],
        ...,
        [          nan, -1.749502e-10, ...,  6.997423e-14,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -4.606194e-11, ...,  1.980042e-12,           nan],
        ...,
        [          nan,  8.697668e-12, ..., -1.904941e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]]],
      dtype=float32)
    • MEKE_src_mom_K4
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the biharmonic of MEKE units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       ...,

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]]], dtype=float32)
    • MEKE_src_btm_drag
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the bottom drag acting on MEKE units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[          nan,           nan, ...,           nan,           nan],
        [          nan, -5.955678e-11, ..., -1.899339e-11,           nan],
        ...,
        [          nan, -9.961407e-11, ..., -2.525525e-11,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -2.711997e-10, ..., -2.144404e-11,           nan],
        ...,
        [          nan, -2.054272e-10, ..., -3.957972e-11,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       ...,

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -2.412892e-09, ..., -1.260603e-11,           nan],
        ...,
        [          nan, -5.296749e-09, ..., -4.903627e-11,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -1.340101e-09, ..., -1.321974e-11,           nan],
        ...,
        [          nan, -3.624543e-09, ..., -4.886812e-11,           nan],
        [          nan,           nan, ...,           nan,           nan]]],
      dtype=float32)
    • MEKE_src_GM
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the thickness mixing (GM scheme) units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       ...,

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]],

       [[nan, nan, ..., nan, nan],
        [nan,  0., ...,  0., nan],
        ...,
        [nan,  0., ...,  0., nan],
        [nan, nan, ..., nan, nan]]], dtype=float32)
    • MEKE_src_mom_lp
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the Laplacian of resolved flows units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[          nan,           nan, ...,           nan,           nan],
        [          nan, -3.570098e-12, ..., -1.334514e-12,           nan],
        ...,
        [          nan, -5.332010e-12, ..., -3.634904e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -2.041340e-11, ..., -8.013071e-13,           nan],
        ...,
        [          nan, -3.927317e-12, ..., -3.705581e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       ...,

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -1.361619e-10, ..., -5.275114e-13,           nan],
        ...,
        [          nan, -4.991446e-10, ..., -5.904024e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]],

       [[          nan,           nan, ...,           nan,           nan],
        [          nan, -3.815689e-11, ..., -5.214029e-13,           nan],
        ...,
        [          nan, -3.053067e-10, ..., -4.545817e-13,           nan],
        [          nan,           nan, ...,           nan,           nan]]],
      dtype=float32)
    • MEKE_src_mom_bh
      (time, yh, xh)
      float32
      nan nan nan nan ... nan nan nan nan
      long_name : MEKE energy source from the biharmonic of resolved flows units : m2 s-3 cell_methods : area:mean yh:mean xh:mean time: mean cell_measures : area: area_t time_avg_info : average_T1,average_T2,average_DT
      array([[[         nan,          nan, ...,          nan,          nan],
        [         nan, 6.543775e-11, ..., 2.655131e-11,          nan],
        ...,
        [         nan, 1.789059e-10, ..., 2.694154e-11,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 2.609439e-10, ..., 1.605143e-11,          nan],
        ...,
        [         nan, 1.884744e-10, ..., 4.102047e-11,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       ...,

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 2.703765e-09, ..., 1.177210e-11,          nan],
        ...,
        [         nan, 6.955740e-09, ..., 5.410191e-11,          nan],
        [         nan,          nan, ...,          nan,          nan]],

       [[         nan,          nan, ...,          nan,          nan],
        [         nan, 7.969787e-10, ..., 1.277304e-11,          nan],
        ...,
        [         nan, 2.514419e-09, ..., 4.511322e-11,          nan],
        [         nan,          nan, ...,          nan,          nan]]],
      dtype=float32)
    • average_T1
      (time)
      float64
      3.15e+04 3.150e+04 ... 3.2e+04
      long_name : Start time for average period units : days since 0001-01-01 00:00:00
      array([31500., 31505., 31510., 31515., 31520., 31525., 31530., 31535., 31540.,
       31545., 31550., 31555., 31560., 31565., 31570., 31575., 31580., 31585.,
       31590., 31595., 31600., 31605., 31610., 31615., 31620., 31625., 31630.,
       31635., 31640., 31645., 31650., 31655., 31660., 31665., 31670., 31675.,
       31680., 31685., 31690., 31695., 31700., 31705., 31710., 31715., 31720.,
       31725., 31730., 31735., 31740., 31745., 31750., 31755., 31760., 31765.,
       31770., 31775., 31780., 31785., 31790., 31795., 31800., 31805., 31810.,
       31815., 31820., 31825., 31830., 31835., 31840., 31845., 31850., 31855.,
       31860., 31865., 31870., 31875., 31880., 31885., 31890., 31895., 31900.,
       31905., 31910., 31915., 31920., 31925., 31930., 31935., 31940., 31945.,
       31950., 31955., 31960., 31965., 31970., 31975., 31980., 31985., 31990.,
       31995.])
    • average_T2
      (time)
      float64
      3.150e+04 3.151e+04 ... 3.2e+04
      long_name : End time for average period units : days since 0001-01-01 00:00:00
      array([31505., 31510., 31515., 31520., 31525., 31530., 31535., 31540., 31545.,
       31550., 31555., 31560., 31565., 31570., 31575., 31580., 31585., 31590.,
       31595., 31600., 31605., 31610., 31615., 31620., 31625., 31630., 31635.,
       31640., 31645., 31650., 31655., 31660., 31665., 31670., 31675., 31680.,
       31685., 31690., 31695., 31700., 31705., 31710., 31715., 31720., 31725.,
       31730., 31735., 31740., 31745., 31750., 31755., 31760., 31765., 31770.,
       31775., 31780., 31785., 31790., 31795., 31800., 31805., 31810., 31815.,
       31820., 31825., 31830., 31835., 31840., 31845., 31850., 31855., 31860.,
       31865., 31870., 31875., 31880., 31885., 31890., 31895., 31900., 31905.,
       31910., 31915., 31920., 31925., 31930., 31935., 31940., 31945., 31950.,
       31955., 31960., 31965., 31970., 31975., 31980., 31985., 31990., 31995.,
       32000.])
    • average_DT
      (time)
      float64
      5.0 5.0 5.0 5.0 ... 5.0 5.0 5.0 5.0
      long_name : Length of average period units : days
      array([5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5.,
       5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5.,
       5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5.,
       5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5.,
       5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5., 5.,
       5., 5., 5., 5., 5., 5., 5., 5., 5., 5.])
    • time_bnds
      (time, nv)
      float64
      3.15e+04 3.150e+04 ... 3.2e+04
      long_name : time axis boundaries units : days
      array([[31500., 31505.],
       [31505., 31510.],
       [31510., 31515.],
       [31515., 31520.],
       [31520., 31525.],
       [31525., 31530.],
       [31530., 31535.],
       [31535., 31540.],
       [31540., 31545.],
       [31545., 31550.],
       [31550., 31555.],
       [31555., 31560.],
       [31560., 31565.],
       [31565., 31570.],
       [31570., 31575.],
       [31575., 31580.],
       [31580., 31585.],
       [31585., 31590.],
       [31590., 31595.],
       [31595., 31600.],
       [31600., 31605.],
       [31605., 31610.],
       [31610., 31615.],
       [31615., 31620.],
       [31620., 31625.],
       [31625., 31630.],
       [31630., 31635.],
       [31635., 31640.],
       [31640., 31645.],
       [31645., 31650.],
       [31650., 31655.],
       [31655., 31660.],
       [31660., 31665.],
       [31665., 31670.],
       [31670., 31675.],
       [31675., 31680.],
       [31680., 31685.],
       [31685., 31690.],
       [31690., 31695.],
       [31695., 31700.],
       [31700., 31705.],
       [31705., 31710.],
       [31710., 31715.],
       [31715., 31720.],
       [31720., 31725.],
       [31725., 31730.],
       [31730., 31735.],
       [31735., 31740.],
       [31740., 31745.],
       [31745., 31750.],
       [31750., 31755.],
       [31755., 31760.],
       [31760., 31765.],
       [31765., 31770.],
       [31770., 31775.],
       [31775., 31780.],
       [31780., 31785.],
       [31785., 31790.],
       [31790., 31795.],
       [31795., 31800.],
       [31800., 31805.],
       [31805., 31810.],
       [31810., 31815.],
       [31815., 31820.],
       [31820., 31825.],
       [31825., 31830.],
       [31830., 31835.],
       [31835., 31840.],
       [31840., 31845.],
       [31845., 31850.],
       [31850., 31855.],
       [31855., 31860.],
       [31860., 31865.],
       [31865., 31870.],
       [31870., 31875.],
       [31875., 31880.],
       [31880., 31885.],
       [31885., 31890.],
       [31890., 31895.],
       [31895., 31900.],
       [31900., 31905.],
       [31905., 31910.],
       [31910., 31915.],
       [31915., 31920.],
       [31920., 31925.],
       [31925., 31930.],
       [31930., 31935.],
       [31935., 31940.],
       [31940., 31945.],
       [31945., 31950.],
       [31950., 31955.],
       [31955., 31960.],
       [31960., 31965.],
       [31965., 31970.],
       [31970., 31975.],
       [31975., 31980.],
       [31980., 31985.],
       [31985., 31990.],
       [31990., 31995.],
       [31995., 32000.]])
  • Attributes: (5)
    filename : MEKEaverages_00031502.nc title : NeverWorld2 associated_files : area_t: static.nc grid_type : regular grid_tile : N/A
In [ ]: