Commit 065cf981 authored by Veit Lüschow's avatar Veit Lüschow

topo version vor 0.1deg ready

parent 8f4a942e
......@@ -43,10 +43,10 @@ class DWBC2Setup(VerosSetup):
"""
GRID, TIMESTEP,...
"""
vs.nx, vs.ny, vs.nz = 45, 100, 30
vs.dt_mom = 300 # normally 2400
vs.dt_tracer = 86400 / 16
vs.runlen = 86400 * 365 * 10
vs.nx, vs.ny, vs.nz = 100, 150, 30
vs.dt_mom = 150 # normally 2400
vs.dt_tracer = 86400 / 32
vs.runlen = 86400 * 365 * 2
# vs.runlen = 86400 * 100
# vs.restart_frequency = 86400 * 100
# vs.runlen = 100
......@@ -71,9 +71,9 @@ class DWBC2Setup(VerosSetup):
'HARMONIC'
vs.enable_hor_friction = True
vs.A_h = 1000
vs.K_h = 1000
vs.A_h = 10
vs.K_h = 10
'VERTICAL'
vs.enable_implicit_vert_friction = True
vs.kappaM = 1e-4
......@@ -90,7 +90,7 @@ class DWBC2Setup(VerosSetup):
'BOTTOM'
vs.enable_bottom_friction = True
vs.r_bot = 1e-3
vs.r_bot = 3e-3
......@@ -100,8 +100,8 @@ class DWBC2Setup(VerosSetup):
vs.dzt[...] = np.array([240.,200.,160,140,120,100,90,80,70,60,50,40,40,40,40, \
40,40,40,40,50,60,70,80,90,100,120,140,160,200,240])
vs.dxt[...] = 0.25
vs.dyt[...] = 0.25
vs.dxt[...] = 0.1
vs.dyt[...] = 0.1
vs.x_origin = 5.0
vs.y_origin = 10
......@@ -123,21 +123,22 @@ class DWBC2Setup(VerosSetup):
vs.kbot[:,:] = np.int(1)
vs.kbot[:3,:] = np.int(0) # western boundary
vs.kbot[3,:] = np.int(22) # slope
vs.kbot[4,:] = np.int(16)
vs.kbot[5,:] = np.int(10)
vs.kbot[6:8,:] = np.int(7)
vs.kbot[8,:] = np.int(5)
vs.kbot[9,:] = np.int(3)
vs.kbot[9:,:] = np.int(2)
vs.kbot[3:4,:] = np.int(22) # slope
vs.kbot[4,:] = np.int(18) # slope
vs.kbot[5,:] = np.int(16)
vs.kbot[6,:] = np.int(10)
vs.kbot[7:9,:] = np.int(7)
vs.kbot[9:11,:] = np.int(5)
vs.kbot[12:14,:] = np.int(3)
vs.kbot[14:,:] = np.int(2)
# vs.kbot[4:,30:70] = np.int(1) # steeper slope
vs.kbot[4:10,50:55] = np.int(16) # seamount
vs.kbot[4:9,48:50] = np.int(16)
vs.kbot[4:9,55:57] = np.int(16)
vs.kbot[4:7,46:48] = np.int(16)
vs.kbot[4:7,57:59] = np.int(16)
# vs.kbot[4:10,50:55] = np.int(16) # seamount
# vs.kbot[4:9,48:50] = np.int(16)
# vs.kbot[4:9,55:57] = np.int(16)
# vs.kbot[4:7,46:48] = np.int(16)
# vs.kbot[4:7,57:59] = np.int(16)
# import matplotlib.pyplot as plt
# plt.contourf(x,y,vs.kbot)
......@@ -152,8 +153,8 @@ class DWBC2Setup(VerosSetup):
"""
zcenter = 1500
zvariance = 3e5
xvariance = 1.1
xcenter = 6.5 # normally 16
xvariance = 0.8
xcenter = 5.5 # normally 16
xx,zz = np.meshgrid(vs.xt,vs.zt,indexing='ij')
uz = np.exp(- ( (zz+(zcenter) )**2 / zvariance ) ) # exponential decay in z
ux = np.exp(- ( (xx-(xcenter) )**2 / xvariance ) ) # gaussian in x
......
......@@ -63,7 +63,7 @@ p = section('dwbc2.snapshot.nc','v')
high=0.2
plt.figure()
p.select(10,'zt',15).plot.contourf(levels=30,cmap='RdBu_r')
p.select(2,'zt',15).plot.contourf(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