In [1]:
import netCDF4 as nc
import datetime as dt
import subprocess
import requests
import matplotlib.pyplot as plt
import cmocean
import numpy as np
import os
import glob
import dateutil as dutil
from salishsea_tools import viz_tools
%matplotlib inline
In [2]:
with nc.Dataset('/ocean/eolson/MEOPAR/NEMO-forcing/grid/mesh_mask201702_noLPE.nc') as fm:
tmask=np.copy(fm.variables['tmask'])
umask=np.copy(fm.variables['umask'])
vmask=np.copy(fm.variables['vmask'])
navlon=np.copy(fm.variables['nav_lon'])
navlat=np.copy(fm.variables['nav_lat'])
e3t_0=np.copy(fm.variables['e3t_0'])
e3u_0=np.copy(fm.variables['e3u_0'])
e3v_0=np.copy(fm.variables['e3v_0'])
e1t=np.copy(fm.variables['e1t'])
e2t=np.copy(fm.variables['e2t'])
e1v=np.copy(fm.variables['e1v'])
e2u=np.copy(fm.variables['e2u'])
A=fm.variables['e1t'][0,:,:]*fm.variables['e2t'][0,:,:]*tmask[0,0,:,:]
In [3]:
t0=dt.datetime(2015,2,6) # 1st start date of run
#te=dt.datetime(2016,12,1)# last start date of runfnum=18
stm=np.shape(tmask)
SiN=2.0
nlen=28
#nlen=37
dlist=[t0+dt.timedelta(days=ii*10) for ii in range(0,nlen)]
In [4]:
sdir0='/results/SalishSea/nowcast-green/'
sdir1='/results/SalishSea/hindcast_v1/'
sdir3='/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_Z3/'
In [5]:
tmaskC=np.copy(tmask)
tmaskC[:,:,370:490,:12]=0
tmaskC[:,:,887:,30:70]=0
In [6]:
tlist=dlist
SiGlobalTot=dict()
SiTot=dict()
BSiTot=dict()
DiatTot=dict()
changeSiGlobalTot=dict()
for idir in (sdir0,sdir1,sdir3):
fformat1='%d%b%y/'
if idir.startswith('/data/eolson/MEOPAR/SS36runs/CedarRuns/'):
fformatT='SalishSea_1h_*_ptrc_T_%Y%m%d-*.nc'
fformatP='SalishSea_1h_*_ptrc_T_%Y%m%d-*.nc'
elif idir==sdir0:
fformatT='SalishSea_1h_%Y%m%d_%Y%m%d_ptrc_T.nc'
fformatP='SalishSea_1h_%Y%m%d_%Y%m%d_grid_T.nc'
elif idir==sdir1:
fformatT='SalishSea_1h_%Y%m%d_%Y%m%d_ptrc_T.nc'
fformatP='SalishSea_1h_%Y%m%d_%Y%m%d_aux_T_.nc'
sumSi=np.zeros((len(tlist),stm[2],stm[3]))
sumBSi=np.zeros((len(tlist),stm[2],stm[3]))
sumDiat=np.zeros((len(tlist),stm[2],stm[3]))
ind=-1
for idt0 in tlist:
ind=ind+1
cdir=idt0.strftime(fformat1).lower()
iffT=idt0.strftime(fformatT)
iffP=idt0.strftime(fformatP)
if idir.startswith('/data/eolson/MEOPAR/SS36runs/CedarRuns/'):
sffT=idir+iffT
sffP=idir+iffP
elif idir.startswith('/results/'):
sffT=idir+cdir+iffT
sffP=idir+cdir+iffP
f=nc.Dataset(glob.glob(sffT)[0])
print(sffT)
fP=nc.Dataset(glob.glob(sffP)[0])
if idir==sdir0:
e3t=np.expand_dims((1+fP.variables['sossheig'][0,:,:]/np.sum(e3t_0*tmask,1)),0)*e3t_0
elif idir==sdir1:
e3t=fP.variables['e3t'][:,:,:,:]
Vol=A*e3t
sumSi[ind,:,:]=1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['silicon'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol
sumBSi[ind,:,:]=1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['biogenic_silicon'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol
sumDiat[ind,:,:]=SiN*1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['diatoms'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol
SiGlobalTot[idir]=np.sum(np.sum(sumSi+sumBSi+sumDiat,2),1)
SiTot[idir]=np.sum(np.sum(sumSi,2),1)
BSiTot[idir]=np.sum(np.sum(sumBSi,2),1)
DiatTot[idir]=np.sum(np.sum(sumDiat,2),1)
changeSiGlobalTot[idir]=[SiGlobalTot[idir][ii+1]-SiGlobalTot[idir][ii] for ii in range(0,len(tlist)-1)]
/results/SalishSea/nowcast-green/06feb15/SalishSea_1h_20150206_20150206_ptrc_T.nc /results/SalishSea/nowcast-green/16feb15/SalishSea_1h_20150216_20150216_ptrc_T.nc /results/SalishSea/nowcast-green/26feb15/SalishSea_1h_20150226_20150226_ptrc_T.nc /results/SalishSea/nowcast-green/08mar15/SalishSea_1h_20150308_20150308_ptrc_T.nc /results/SalishSea/nowcast-green/18mar15/SalishSea_1h_20150318_20150318_ptrc_T.nc /results/SalishSea/nowcast-green/28mar15/SalishSea_1h_20150328_20150328_ptrc_T.nc /results/SalishSea/nowcast-green/07apr15/SalishSea_1h_20150407_20150407_ptrc_T.nc /results/SalishSea/nowcast-green/17apr15/SalishSea_1h_20150417_20150417_ptrc_T.nc /results/SalishSea/nowcast-green/27apr15/SalishSea_1h_20150427_20150427_ptrc_T.nc /results/SalishSea/nowcast-green/07may15/SalishSea_1h_20150507_20150507_ptrc_T.nc /results/SalishSea/nowcast-green/17may15/SalishSea_1h_20150517_20150517_ptrc_T.nc /results/SalishSea/nowcast-green/27may15/SalishSea_1h_20150527_20150527_ptrc_T.nc /results/SalishSea/nowcast-green/06jun15/SalishSea_1h_20150606_20150606_ptrc_T.nc /results/SalishSea/nowcast-green/16jun15/SalishSea_1h_20150616_20150616_ptrc_T.nc /results/SalishSea/nowcast-green/26jun15/SalishSea_1h_20150626_20150626_ptrc_T.nc /results/SalishSea/nowcast-green/06jul15/SalishSea_1h_20150706_20150706_ptrc_T.nc /results/SalishSea/nowcast-green/16jul15/SalishSea_1h_20150716_20150716_ptrc_T.nc /results/SalishSea/nowcast-green/26jul15/SalishSea_1h_20150726_20150726_ptrc_T.nc /results/SalishSea/nowcast-green/05aug15/SalishSea_1h_20150805_20150805_ptrc_T.nc /results/SalishSea/nowcast-green/15aug15/SalishSea_1h_20150815_20150815_ptrc_T.nc /results/SalishSea/nowcast-green/25aug15/SalishSea_1h_20150825_20150825_ptrc_T.nc /results/SalishSea/nowcast-green/04sep15/SalishSea_1h_20150904_20150904_ptrc_T.nc /results/SalishSea/nowcast-green/14sep15/SalishSea_1h_20150914_20150914_ptrc_T.nc /results/SalishSea/nowcast-green/24sep15/SalishSea_1h_20150924_20150924_ptrc_T.nc /results/SalishSea/nowcast-green/04oct15/SalishSea_1h_20151004_20151004_ptrc_T.nc /results/SalishSea/nowcast-green/14oct15/SalishSea_1h_20151014_20151014_ptrc_T.nc /results/SalishSea/nowcast-green/24oct15/SalishSea_1h_20151024_20151024_ptrc_T.nc /results/SalishSea/nowcast-green/03nov15/SalishSea_1h_20151103_20151103_ptrc_T.nc /results/SalishSea/hindcast_v1/06feb15/SalishSea_1h_20150206_20150206_ptrc_T.nc /results/SalishSea/hindcast_v1/16feb15/SalishSea_1h_20150216_20150216_ptrc_T.nc /results/SalishSea/hindcast_v1/26feb15/SalishSea_1h_20150226_20150226_ptrc_T.nc /results/SalishSea/hindcast_v1/08mar15/SalishSea_1h_20150308_20150308_ptrc_T.nc /results/SalishSea/hindcast_v1/18mar15/SalishSea_1h_20150318_20150318_ptrc_T.nc /results/SalishSea/hindcast_v1/28mar15/SalishSea_1h_20150328_20150328_ptrc_T.nc /results/SalishSea/hindcast_v1/07apr15/SalishSea_1h_20150407_20150407_ptrc_T.nc /results/SalishSea/hindcast_v1/17apr15/SalishSea_1h_20150417_20150417_ptrc_T.nc /results/SalishSea/hindcast_v1/27apr15/SalishSea_1h_20150427_20150427_ptrc_T.nc /results/SalishSea/hindcast_v1/07may15/SalishSea_1h_20150507_20150507_ptrc_T.nc /results/SalishSea/hindcast_v1/17may15/SalishSea_1h_20150517_20150517_ptrc_T.nc /results/SalishSea/hindcast_v1/27may15/SalishSea_1h_20150527_20150527_ptrc_T.nc /results/SalishSea/hindcast_v1/06jun15/SalishSea_1h_20150606_20150606_ptrc_T.nc /results/SalishSea/hindcast_v1/16jun15/SalishSea_1h_20150616_20150616_ptrc_T.nc /results/SalishSea/hindcast_v1/26jun15/SalishSea_1h_20150626_20150626_ptrc_T.nc /results/SalishSea/hindcast_v1/06jul15/SalishSea_1h_20150706_20150706_ptrc_T.nc /results/SalishSea/hindcast_v1/16jul15/SalishSea_1h_20150716_20150716_ptrc_T.nc /results/SalishSea/hindcast_v1/26jul15/SalishSea_1h_20150726_20150726_ptrc_T.nc /results/SalishSea/hindcast_v1/05aug15/SalishSea_1h_20150805_20150805_ptrc_T.nc /results/SalishSea/hindcast_v1/15aug15/SalishSea_1h_20150815_20150815_ptrc_T.nc /results/SalishSea/hindcast_v1/25aug15/SalishSea_1h_20150825_20150825_ptrc_T.nc /results/SalishSea/hindcast_v1/04sep15/SalishSea_1h_20150904_20150904_ptrc_T.nc /results/SalishSea/hindcast_v1/14sep15/SalishSea_1h_20150914_20150914_ptrc_T.nc /results/SalishSea/hindcast_v1/24sep15/SalishSea_1h_20150924_20150924_ptrc_T.nc /results/SalishSea/hindcast_v1/04oct15/SalishSea_1h_20151004_20151004_ptrc_T.nc /results/SalishSea/hindcast_v1/14oct15/SalishSea_1h_20151014_20151014_ptrc_T.nc /results/SalishSea/hindcast_v1/24oct15/SalishSea_1h_20151024_20151024_ptrc_T.nc /results/SalishSea/hindcast_v1/03nov15/SalishSea_1h_20151103_20151103_ptrc_T.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150206-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150216-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150226-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150308-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150318-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150328-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150407-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150417-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150427-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150507-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150517-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150527-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150606-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150616-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150626-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150706-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150716-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150726-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150805-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150815-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150825-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150904-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150914-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150924-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151004-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151014-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151024-*.nc /data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151103-*.nc
In [7]:
plt.plot(SiGlobalTot[sdir0]-SiGlobalTot[sdir0][0],'r-')
plt.plot(SiGlobalTot[sdir1]-SiGlobalTot[sdir1][0],'b-')
plt.plot(SiGlobalTot[sdir3]-SiGlobalTot[sdir3][0],'g-')
Out[7]:
[<matplotlib.lines.Line2D at 0x7f8bd4591b70>]
In [ ]: