In [2]:
import cmocean.cm as cm
import h5py
import matplotlib as mpl
import matplotlib.colorbar as colorbar
import matplotlib.colors as mplcolours
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
import xarray as xr
%matplotlib inline
In [3]:
mesh = xr.open_dataset('~/MEOPAR/grid/mesh_mask201702.nc')
In [4]:
oil0d = pd.read_csv('/data/sallen/results/MIDOSS/ParticleNoTests/SB_2000_bunkerc/resOilOutput.sro', sep='\s+', skiprows=4)
oil0d = oil0d.drop([0], axis=0)
length =len(oil0d)
oil0d = oil0d.drop([length-3, length-2, length-1, length], axis=0)
In [5]:
fig, ax = plt.subplots(1, 1, figsize=(15,5))
oil0d.Area.plot(ax=ax);
plt.grid();
In [6]:
fig, ax = plt.subplots(1, 1, figsize=(15,5))
oil0d.MDissolved.plot();
oil0d.MEvaporated.plot();
oil0d.MassOil.plot();
oil0d.MDispersed.plot();
oil0d.MBio.plot()
massbeached = (oil0d.VolOilBeached*oil0d.Density/(1-oil0d.VWaterContent)
*(1-oil0d.MWaterContent))
massbeached.plot(label="Beached")
plt.plot(oil0d.MDissolved + oil0d.MEvaporated + oil0d.MassOil
+ massbeached + oil0d.MDispersed);
plt.legend();
plt.grid();
In [7]:
(oil0d.Density/(1-oil0d.VWaterContent)
*(1-oil0d.MWaterContent)).plot();
In [8]:
oil0d.MDispersed.plot();
In [9]:
fig, ax = plt.subplots(1, 1, figsize=(15, 3.5))
ax.plot(oil0d.MDissolved + oil0d.MEvaporated + oil0d.MassOil
+ massbeached + oil0d.MDispersed + oil0d.MBio);
ax.grid();
In [10]:
massbeached.plot(label="Beached");
plt.grid()
In [11]:
(oil0d.MDispersed/(oil0d.MDispersed + oil0d.MassOil))[:-2].plot()
plt.grid()
In [12]:
oil0d.VolumeOil.plot();
plt.legend();
plt.grid();
In [13]:
oil0d.AnalyteMass1.plot()
oil0d.AnalyteMass2.plot()
oil0d.AnalyteMass3.plot()
oil0d.AnalyteMass4.plot()
oil0d.AnalyteMass5.plot();
plt.legend()
Out[13]:
<matplotlib.legend.Legend at 0x7f0295c69710>
In [14]:
oil0d.Thickness.plot();
In [15]:
oilLag = xr.open_dataset('/data/sallen/results/MIDOSS/ParticleNoTests/SB_2000_bunkerc/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc.nc')
oilLag
Out[15]:
xarray.Dataset
- grid_x: 396
- grid_x_latlon: 397
- grid_y: 896
- grid_y_latlon: 897
- grid_z: 40
- time: 191
- time(time)datetime64[ns]2018-01-15T00:30:00 ... 2018-01-22T22:30:00
- standard_name :
- time
- long_name :
- time
- tz_name :
- UTC
array(['2018-01-15T00:30:00.000000000', '2018-01-15T01:30:00.000000000', '2018-01-15T02:30:00.000000000', '2018-01-15T03:30:00.000000000', '2018-01-15T04:30:00.000000000', '2018-01-15T05:30:00.000000000', '2018-01-15T06:30:00.000000000', '2018-01-15T07:30:00.000000000', '2018-01-15T08:30:00.000000000', '2018-01-15T09:30:00.000000000', '2018-01-15T10:30:00.000000000', '2018-01-15T11:30:00.000000000', '2018-01-15T12:30:00.000000000', '2018-01-15T13:30:00.000000000', '2018-01-15T14:30:00.000000000', '2018-01-15T15:30:00.000000000', '2018-01-15T16:30:00.000000000', '2018-01-15T17:30:00.000000000', '2018-01-15T18:30:00.000000000', '2018-01-15T19:30:00.000000000', '2018-01-15T20:30:00.000000000', '2018-01-15T21:30:00.000000000', '2018-01-15T22:30:00.000000000', '2018-01-15T23:30:00.000000000', '2018-01-16T00:30:00.000000000', '2018-01-16T01:30:00.000000000', '2018-01-16T02:30:00.000000000', '2018-01-16T03:30:00.000000000', '2018-01-16T04:30:00.000000000', '2018-01-16T05:30:00.000000000', '2018-01-16T06:30:00.000000000', '2018-01-16T07:30:00.000000000', '2018-01-16T08:30:00.000000000', '2018-01-16T09:30:00.000000000', '2018-01-16T10:30:00.000000000', '2018-01-16T11:30:00.000000000', '2018-01-16T12:30:00.000000000', '2018-01-16T13:30:00.000000000', '2018-01-16T14:30:00.000000000', '2018-01-16T15:30:00.000000000', '2018-01-16T16:30:00.000000000', '2018-01-16T17:30:00.000000000', '2018-01-16T18:30:00.000000000', '2018-01-16T19:30:00.000000000', '2018-01-16T20:30:00.000000000', '2018-01-16T21:30:00.000000000', '2018-01-16T22:30:00.000000000', '2018-01-16T23:30:00.000000000', '2018-01-17T00:30:00.000000000', '2018-01-17T01:30:00.000000000', '2018-01-17T02:30:00.000000000', '2018-01-17T03:30:00.000000000', '2018-01-17T04:30:00.000000000', '2018-01-17T05:30:00.000000000', '2018-01-17T06:30:00.000000000', '2018-01-17T07:30:00.000000000', '2018-01-17T08:30:00.000000000', '2018-01-17T09:30:00.000000000', '2018-01-17T10:30:00.000000000', '2018-01-17T11:30:00.000000000', '2018-01-17T12:30:00.000000000', '2018-01-17T13:30:00.000000000', '2018-01-17T14:30:00.000000000', '2018-01-17T15:30:00.000000000', '2018-01-17T16:30:00.000000000', '2018-01-17T17:30:00.000000000', '2018-01-17T18:30:00.000000000', '2018-01-17T19:30:00.000000000', '2018-01-17T20:30:00.000000000', '2018-01-17T21:30:00.000000000', '2018-01-17T22:30:00.000000000', '2018-01-17T23:30:00.000000000', '2018-01-18T00:30:00.000000000', '2018-01-18T01:30:00.000000000', '2018-01-18T02:30:00.000000000', '2018-01-18T03:30:00.000000000', '2018-01-18T04:30:00.000000000', '2018-01-18T05:30:00.000000000', '2018-01-18T06:30:00.000000000', '2018-01-18T07:30:00.000000000', '2018-01-18T08:30:00.000000000', '2018-01-18T09:30:00.000000000', '2018-01-18T10:30:00.000000000', '2018-01-18T11:30:00.000000000', '2018-01-18T12:30:00.000000000', '2018-01-18T13:30:00.000000000', '2018-01-18T14:30:00.000000000', '2018-01-18T15:30:00.000000000', '2018-01-18T16:30:00.000000000', '2018-01-18T17:30:00.000000000', '2018-01-18T18:30:00.000000000', '2018-01-18T19:30:00.000000000', '2018-01-18T20:30:00.000000000', '2018-01-18T21:30:00.000000000', '2018-01-18T22:30:00.000000000', '2018-01-18T23:30:00.000000000', '2018-01-19T00:30:00.000000000', '2018-01-19T01:30:00.000000000', '2018-01-19T02:30:00.000000000', '2018-01-19T03:30:00.000000000', '2018-01-19T04:30:00.000000000', '2018-01-19T05:30:00.000000000', '2018-01-19T06:30:00.000000000', '2018-01-19T07:30:00.000000000', '2018-01-19T08:30:00.000000000', '2018-01-19T09:30:00.000000000', '2018-01-19T10:30:00.000000000', '2018-01-19T11:30:00.000000000', '2018-01-19T12:30:00.000000000', '2018-01-19T13:30:00.000000000', '2018-01-19T14:30:00.000000000', '2018-01-19T15:30:00.000000000', '2018-01-19T16:30:00.000000000', '2018-01-19T17:30:00.000000000', '2018-01-19T18:30:00.000000000', '2018-01-19T19:30:00.000000000', '2018-01-19T20:30:00.000000000', '2018-01-19T21:30:00.000000000', '2018-01-19T22:30:00.000000000', '2018-01-19T23:30:00.000000000', '2018-01-20T00:30:00.000000000', '2018-01-20T01:30:00.000000000', '2018-01-20T02:30:00.000000000', '2018-01-20T03:30:00.000000000', '2018-01-20T04:30:00.000000000', '2018-01-20T05:30:00.000000000', '2018-01-20T06:30:00.000000000', '2018-01-20T07:30:00.000000000', '2018-01-20T08:30:00.000000000', '2018-01-20T09:30:00.000000000', '2018-01-20T10:30:00.000000000', '2018-01-20T11:30:00.000000000', '2018-01-20T12:30:00.000000000', '2018-01-20T13:30:00.000000000', '2018-01-20T14:30:00.000000000', '2018-01-20T15:30:00.000000000', '2018-01-20T16:30:00.000000000', '2018-01-20T17:30:00.000000000', '2018-01-20T18:30:00.000000000', '2018-01-20T19:30:00.000000000', '2018-01-20T20:30:00.000000000', '2018-01-20T21:30:00.000000000', '2018-01-20T22:30:00.000000000', '2018-01-20T23:30:00.000000000', '2018-01-21T00:30:00.000000000', '2018-01-21T01:30:00.000000000', '2018-01-21T02:30:00.000000000', '2018-01-21T03:30:00.000000000', '2018-01-21T04:30:00.000000000', '2018-01-21T05:30:00.000000000', '2018-01-21T06:30:00.000000000', '2018-01-21T07:30:00.000000000', '2018-01-21T08:30:00.000000000', '2018-01-21T09:30:00.000000000', '2018-01-21T10:30:00.000000000', '2018-01-21T11:30:00.000000000', '2018-01-21T12:30:00.000000000', '2018-01-21T13:30:00.000000000', '2018-01-21T14:30:00.000000000', '2018-01-21T15:30:00.000000000', '2018-01-21T16:30:00.000000000', '2018-01-21T17:30:00.000000000', '2018-01-21T18:30:00.000000000', '2018-01-21T19:30:00.000000000', '2018-01-21T20:30:00.000000000', '2018-01-21T21:30:00.000000000', '2018-01-21T22:30:00.000000000', '2018-01-21T23:30:00.000000000', '2018-01-22T00:30:00.000000000', '2018-01-22T01:30:00.000000000', '2018-01-22T02:30:00.000000000', '2018-01-22T03:30:00.000000000', '2018-01-22T04:30:00.000000000', '2018-01-22T05:30:00.000000000', '2018-01-22T06:30:00.000000000', '2018-01-22T07:30:00.000000000', '2018-01-22T08:30:00.000000000', '2018-01-22T09:30:00.000000000', '2018-01-22T10:30:00.000000000', '2018-01-22T11:30:00.000000000', '2018-01-22T12:30:00.000000000', '2018-01-22T13:30:00.000000000', '2018-01-22T14:30:00.000000000', '2018-01-22T15:30:00.000000000', '2018-01-22T16:30:00.000000000', '2018-01-22T17:30:00.000000000', '2018-01-22T18:30:00.000000000', '2018-01-22T19:30:00.000000000', '2018-01-22T20:30:00.000000000', '2018-01-22T21:30:00.000000000', '2018-01-22T22:30:00.000000000'], dtype='datetime64[ns]') - grid_y(grid_y)int160 1 2 3 4 5 ... 891 892 893 894 895
- standard_name :
- model_y_index
- long_name :
- y index
array([ 0, 1, 2, ..., 893, 894, 895], dtype=int16)
- grid_x(grid_x)int160 1 2 3 4 5 ... 391 392 393 394 395
- standard_name :
- model_x_index
- long_name :
- x index
array([ 0, 1, 2, ..., 393, 394, 395], dtype=int16)
- grid_z(grid_z)int160 1 2 3 4 5 6 ... 34 35 36 37 38 39
- standard_name :
- model_level_index
- long_name :
- depth level
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39], dtype=int16) - grid_y_latlon(grid_y_latlon)float320.0 1.0 2.0 ... 894.0 895.0 896.0
- standard_name :
- model_latlon_y_index
- long_name :
- latlon y index
- units :
- degrees_north
array([ 0., 1., 2., ..., 894., 895., 896.], dtype=float32)
- grid_x_latlon(grid_x_latlon)float320.0 1.0 2.0 ... 394.0 395.0 396.0
- standard_name :
- model_latlon_x_index
- long_name :
- latlon x index
- units :
- degrees_east
array([ 0., 1., 2., ..., 394., 395., 396.], dtype=float32)
- OilConcentration_2D(time, grid_y, grid_x)float64...
- standard_name :
- OilConcentration_2D
- long_name :
- OilConcentration_2D
- units :
- ppm
[67769856 values with dtype=float64]
- Thickness_2D(time, grid_y, grid_x)float64...
- standard_name :
- Thickness_2D
- long_name :
- Thickness_2D
- units :
- microns
[67769856 values with dtype=float64]
- Dissolution_3D(time, grid_z, grid_y, grid_x)float64...
- standard_name :
- Dissolution_3D
- long_name :
- Dissolution_3D
- units :
- Kg/m3
[2710794240 values with dtype=float64]
- OilConcentration_3D(time, grid_z, grid_y, grid_x)float64...
- standard_name :
- OilConcentration_3D
- long_name :
- OilConcentration_3D
- units :
- Kg/m3
[2710794240 values with dtype=float64]
- OilWaterColumnOilVol_3D(time, grid_z, grid_y, grid_x)float64...
- standard_name :
- OilWaterColumnOilVol_3D
- long_name :
- OilWaterColumnOilVol_3D
- units :
- m3/gridcell
[2710794240 values with dtype=float64]
- Latitude(grid_y_latlon, grid_x_latlon)float64...
- standard_name :
- Latitude
- long_name :
- Latitude
- units :
- degrees_north
[356109 values with dtype=float64]
- Longitude(grid_y_latlon, grid_x_latlon)float64...
- standard_name :
- Longitude
- long_name :
- Longitude
- units :
- degrees_east
[356109 values with dtype=float64]
- Beaching_Time(grid_y, grid_x)datetime64[ns]...
- standard_name :
- Beaching_Time
- long_name :
- Beaching Time
[354816 values with dtype=datetime64[ns]]
- Beaching_Volume(grid_y, grid_x)float64...
- standard_name :
- Beaching_Volume
- long_name :
- Beaching Volume
- units :
- m3
[354816 values with dtype=float64]
- Oil_Arrival_Time(grid_y, grid_x)datetime64[ns]...
- standard_name :
- Oil_Arrival_Time
- long_name :
- Oil Arrival Time
[354816 values with dtype=datetime64[ns]]
- nco_openmp_thread_number :
- 1
- history :
- Tue Nov 24 20:44:09 2020: ncks -4 -L4 -A /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_oil_times.nc /localscratch/allen.41131240.0/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc.nc Tue Nov 24 20:38:13 2020: ncrcat -4 -L4 -O -o /localscratch/allen.41131240.0/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T0030.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-15T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-16T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-17T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-18T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-19T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-20T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T22:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-21T23:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T00:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T01:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T02:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T03:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T04:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T05:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T06:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T07:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T08:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T09:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T10:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T11:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T12:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T13:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T14:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T15:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T16:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T17:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T18:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T19:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T20:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T21:30.nc /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_2018-01-22T22:30.nc
- NCO :
- 4.6.6
- history_of_appended_files :
- Tue Nov 24 20:44:09 2020: Appended file /localscratch/allen.41131240.0/tmpqiqkuvx8/Lagrangian_SB_2000_bunkerc_SB_2000_bunkerc_oil_times.nc had no "history" attribute
In [16]:
imin, imax = 200, 350
jmin, jmax = 150, 280
In [17]:
fig, axs = plt.subplots(1, 2, figsize=(15, 5))
it = 20
field = oilLag.OilWaterColumnOilVol_3D[it]
(field[:, imin:imax , jmin:jmax].sum(axis=0)).plot(ax=axs[0], cmap='copper')
print (oilLag.OilWaterColumnOilVol_3D[it, :, imin:imax, jmin:jmax].sum(axis=0).sum(axis=0).sum(axis=0))
axs[1].plot(oilLag.time[0:72], oilLag.OilWaterColumnOilVol_3D[0:72].sum(axis=1).sum(axis=1).sum(axis=1)
)
<xarray.DataArray 'OilWaterColumnOilVol_3D' ()>
array(10000.2349)
Coordinates:
time datetime64[ns] 2018-01-15T20:30:00
Out[17]:
[<matplotlib.lines.Line2D at 0x7f0295ab9650>]
In [18]:
oil0d.VolumeOil[0:72].plot()
Out[18]:
<matplotlib.axes._subplots.AxesSubplot at 0x7f02959b9190>
In [19]:
tmax = 20
fig, axs = plt.subplots(1, 3, figsize=(15, 5))
oilLag.OilWaterColumnOilVol_3D[1:tmax+1, 39, 200:300, 200:300].sum(axis=0).plot(ax=axs[0]);
oilLag.OilWaterColumnOilVol_3D[1:tmax+1, 38, 200:300, 200:300].sum(axis=0).plot(ax=axs[1]);
oilLag.OilWaterColumnOilVol_3D[1:tmax+1, 0:38, 200:300, 200:300].sum(axis=0).sum(axis=0).plot(ax=axs[2]);
In [20]:
fig, ax = plt.subplots(1, 1, figsize=(5, 5))
imin, imax = 250, 290
jmin, jmax = 230, 270
oilLag.OilWaterColumnOilVol_3D[:, 39, imin:imax, jmin:jmax].sum(axis=0).plot(ax=ax, cmap='gist_ncar');
In [21]:
fig, ax = plt.subplots(1, 1, figsize=(5, 5))
imin, imax = 220, 350
jmin, jmax = 150, 280
oilLag.OilWaterColumnOilVol_3D[:, 39, imin:imax, jmin:jmax].sum(axis=0).plot(ax=ax, cmap='gist_ncar');
In [22]:
fig, ax = plt.subplots(1, 1, figsize=(5, 5))
imin, imax = 250, 350
jmin, jmax = 150, 250
oilLag.Beaching_Volume[imin:imax , jmin:jmax].plot(ax=ax, cmap='copper')
print (oilLag.Beaching_Volume.sum(axis=0).sum(axis=0))
<xarray.DataArray 'Beaching_Volume' ()> array(9994.64466358)
In [23]:
bt = np.array(oilLag.Beaching_Time[imin:imax, jmin:jmax] - oilLag.Beaching_Time[imin:imax, jmin:jmax].min())
plt.plot(bt/1e9/3600, 'ro');
In [24]:
watercolour = 'lightskyblue'
landcolour = 'papayawhip'
waterland_cmap = mplcolours.LinearSegmentedColormap.from_list('mycmap', [(0, watercolour), (0.85, watercolour),
(0.850001, landcolour), (1, landcolour)])
xs = range(jmin, jmax)
ys = range(imin, imax)
xx, yy = np.meshgrid(xs, ys)
mymask = oilLag.Beaching_Time[imin:imax, jmin:jmax] == oilLag.Beaching_Time[imin:imax, jmin:jmax].min()
days = np.ma.array(bt/1e9/86400., mask=mymask)
In [27]:
cnorm = mplcolours.Normalize(vmin=0, vmax=7)
mycmap = mpl.cm.copper
fig, ax = plt.subplots(1, 1, figsize=(5, 5))
ax.set_aspect(500/440.)
ax.tick_params(axis='both', which='both', bottom=False, left=False, labelbottom=False, labelleft=False)
ax.pcolormesh(1-mesh.tmask[0, 0, imin:imax, jmin:jmax], cmap=waterland_cmap)
ax.scatter(xx-jmin+1, yy-imin+1, s=10, c=days, cmap='copper', vmin=0, vmax=7, marker='o', edgecolors='face');
ax_cbar = fig.add_axes([0.9, 0.12, 0.05, 0.76])
cb = colorbar.ColorbarBase(ax=ax_cbar, cmap=mycmap, norm=cnorm);
cb.set_label('Days to Beaching')
fig.suptitle('10,000,000 $\ell$ Oil Spilled on Jan 15, 2018 at Salmon Bank');
In [27]:
fig, ax = plt.subplots(1, 1, figsize=(5, 5))
imin, imax = 250, 350
jmin, jmax = 150, 250
np.log(oilLag.Beaching_Volume[imin:imax , jmin:jmax]).plot(ax=ax, cmap='copper_r',
vmax=5, vmin=0);
/home/sallen/anaconda/envs/py37/lib/python3.7/site-packages/xarray/core/computation.py:604: RuntimeWarning: divide by zero encountered in log result_data = func(*input_data)
In [ ]: