Table of Contents
Crashing on first time step on the western boundary at 2 465 14 in fortran.
Going to try smoothing the restart file.
In [1]:
import matplotlib.pyplot as plt
import numpy as np
import xarray as xr
%matplotlib inline
In [2]:
restart_file = xr.open_dataset('/results/SalishSea/nowcast-agrif.201702/13nov20/SalishSea_04914000_restart.nc')
restart_file
Out[2]:
xarray.Dataset
- t: 1
- x: 398
- y: 898
- z: 40
- nav_lon(y, x)float32...
[357404 values with dtype=float32]
- nav_lat(y, x)float32...
[357404 values with dtype=float32]
- nav_lev(z)float32...
array([ 0.5 , 1.500003, 2.500011, 3.500031, 4.500071, 5.500151, 6.50031 , 7.500623, 8.501236, 9.502433, 10.504766, 11.509312, 12.518167, 13.535412, 14.568982, 15.634288, 16.761173, 18.007135, 19.481785, 21.389978, 24.100256, 28.229916, 34.685757, 44.517723, 58.484333, 76.58559 , 98.06296 , 121.866516, 147.08946 , 173.11449 , 199.57304 , 226.2603 , 253.06664 , 279.93454 , 306.8342 , 333.75018 , 360.67453 , 387.6032 , 414.5341 , 441.4661 ], dtype=float32) - time_counter(t)float64...
array([4914000.])
- kt()float64...
array(4914000.)
- ndastp()float64...
array(20201113.)
- adatrj()float64...
array(2275.)
- ssh_ibb(t, y, x)float64...
[357404 values with dtype=float64]
- rnf_b(t, y, x)float64...
[357404 values with dtype=float64]
- rnf_hc_b(t, y, x)float64...
[357404 values with dtype=float64]
- rnf_sc_b(t, y, x)float64...
[357404 values with dtype=float64]
- utau_b(t, y, x)float64...
[357404 values with dtype=float64]
- vtau_b(t, y, x)float64...
[357404 values with dtype=float64]
- qns_b(t, y, x)float64...
[357404 values with dtype=float64]
- emp_b(t, y, x)float64...
[357404 values with dtype=float64]
- sfx_b(t, y, x)float64...
[357404 values with dtype=float64]
- en(t, z, y, x)float64...
[14296160 values with dtype=float64]
- avt(t, z, y, x)float64...
[14296160 values with dtype=float64]
- avm(t, z, y, x)float64...
[14296160 values with dtype=float64]
- avmu(t, z, y, x)float64...
[14296160 values with dtype=float64]
- avmv(t, z, y, x)float64...
[14296160 values with dtype=float64]
- mxln(t, z, y, x)float64...
[14296160 values with dtype=float64]
- ub2_b(t, y, x)float64...
[357404 values with dtype=float64]
- vb2_b(t, y, x)float64...
[357404 values with dtype=float64]
- sbc_hc_b(t, y, x)float64...
[357404 values with dtype=float64]
- sbc_sc_b(t, y, x)float64...
[357404 values with dtype=float64]
- qsr_hc_b(t, z, y, x)float64...
[14296160 values with dtype=float64]
- fraqsr_1lev(t, y, x)float64...
[357404 values with dtype=float64]
- fse3t_b(t, z, y, x)float64...
[14296160 values with dtype=float64]
- fse3t_n(t, z, y, x)float64...
[14296160 values with dtype=float64]
- rdt()float64...
array(40.)
- rdttra1()float64...
array(40.)
- ub(t, z, y, x)float64...
[14296160 values with dtype=float64]
- vb(t, z, y, x)float64...
[14296160 values with dtype=float64]
- tb(t, z, y, x)float64...
[14296160 values with dtype=float64]
- sb(t, z, y, x)float64...
[14296160 values with dtype=float64]
- rotb(t, z, y, x)float64...
[14296160 values with dtype=float64]
- hdivb(t, z, y, x)float64...
[14296160 values with dtype=float64]
- sshb(t, y, x)float64...
[357404 values with dtype=float64]
- un(t, z, y, x)float64...
[14296160 values with dtype=float64]
- vn(t, z, y, x)float64...
[14296160 values with dtype=float64]
- tn(t, z, y, x)float64...
[14296160 values with dtype=float64]
- sn(t, z, y, x)float64...
[14296160 values with dtype=float64]
- rotn(t, z, y, x)float64...
[14296160 values with dtype=float64]
- hdivn(t, z, y, x)float64...
[14296160 values with dtype=float64]
- sshn(t, y, x)float64...
[357404 values with dtype=float64]
- rhop(t, z, y, x)float64...
[14296160 values with dtype=float64]
- file_name :
- SalishSea_04914000_restart.nc
- TimeStamp :
- 13/11/2020 10:56:10 -0800
In [7]:
prob = 465-1
pz = 14-1
Look at velocity
In [8]:
ib, ie = 385, 467
uveln19 = restart_file.un[0, pz]
uvelb19 = restart_file.ub[0, pz]
uvelb19[ib:ie, 1].plot(label="Before 1", linewidth=5)
uveln19[ib:ie, 1].plot(label="Next 1")
uvelb19[ib:ie, 3].plot(label="Before 3", linewidth=5)
uveln19[ib:ie, 3].plot(label="Next 3")
plt.plot([prob-ib, prob-ib], [-0.5, 2], color='k');
plt.legend();
In [9]:
ib, ie = 385, 467
vveln19 = restart_file.vn[0, pz]
vvelb19 = restart_file.vb[0, pz]
vvelb19[ib:ie, 1].plot(label="Before 1", linewidth=5)
vveln19[ib:ie, 1].plot(label="Next 1")
vvelb19[ib:ie, 3].plot(label="Before 3", linewidth=5)
vveln19[ib:ie, 3].plot(label="Next 3")
plt.plot([prob-ib, prob-ib], [-0.5, 2], color='k');
plt.legend();
In [10]:
restart_file.un[0, :, :, 1] = restart_file.un[0, :, :, 2]
restart_file.ub[0, :, :, 1] = restart_file.ub[0, :, :, 2]
restart_file.vn[0, :, :, 1] = restart_file.vn[0, :, :, 2]
restart_file.vb[0, :, :, 1] = restart_file.vb[0, :, :, 2]
In [11]:
restart_file.attrs['Modified'] = 'copied velocity one away from west boundary to boundary'
restart_file.attrs
Out[11]:
{'file_name': 'SalishSea_04914000_restart.nc',
'TimeStamp': '13/11/2020 10:56:10 -0800',
'Modified': 'copied velocity one away from west boundary to boundary'}
In [12]:
restart_file.to_netcdf('SalishSea_04914000_restart.nc')
In [ ]: