Commit 4d91eee3 authored by Veit Lüschow's avatar Veit Lüschow

restoring temp indstead of set fixed value at north and south boundary

parent a7aac2a4
......@@ -37,7 +37,6 @@ class DWBC2Setup(VerosSetup):
vs.dwbc2 = True
vs.eq_of_state_type = 1
vs.enable_tempsalt_sources = True
# vs.restart_input_filename='merjet_9637.restart.h5'
......@@ -238,17 +237,19 @@ class DWBC2Setup(VerosSetup):
# vs.t_restoring[-3:,:,:] = 1 / (1e1 * vs.dt_mom)
# vs.t_restoring[-5:-3,:,:] = 1 / (2e1 * vs.dt_mom)
vs.tscl_tr = 2e-5
vs.t_restoring_tr[:,:,-1] = 2e-5
vs.t_restoring_tr[:,:,-2] = 1e-5
vs.t_restoring_tr[:,:,-1] = vs.tscl_tr
vs.t_restoring_tr[:,:,-2] = vs.tscl_tr / 2
@veros_method
def set_forcing(self, vs):
pass
# vs.dv_rest[...] = vs.maskV * vs.t_restoring * (vs.v_ini[...] - vs.v[...,vs.tau])
# vs.temp_source[...] = vs.maskT * vs.t_restoring_tr * (vs.temp_surface - vs.temp[... , vs.tau])
# vs.salt_source[...] = vs.maskT * vs.t_restoring_tr * (35 - vs.salt[... , vs.tau])
vs.temp_source[...] = vs.maskT * vs.t_restoring_tr * (vs.temp_surface - vs.temp[... , vs.tau])
vs.salt_source[...] = vs.maskT * vs.t_restoring_tr * (35 - vs.salt[... , vs.tau])
@veros_method
......
......@@ -61,7 +61,7 @@ p = section('dwbc2.snapshot.nc','v')
high=0.2
plt.figure()
p.select(350,'zt',15).plot.contourf(vmin=-high,vmax=high,levels=30,cmap='RdBu_r')
p.select(360,'zt',15).plot.contourf(vmin=-high,vmax=high,levels=30,cmap='RdBu_r')
#
# plt.colorbar()
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment