Commit bcd30dce authored by Veit Lüschow's avatar Veit Lüschow

topo version with gap ready for test on mistral

parent 5e9cd543
......@@ -7,6 +7,7 @@
#SBATCH -t 12:00:00 # Set a limit on the total run time
#SBATCH -A mh0256 # Charge resources on this project account
#SBATCH -o my_job.o%j # File name for standard and error output
#SBATCH --mail-type=FAIL # receive mail only when job fails
# Execute a serial program, e.g.
source activate veros
......
......@@ -2,8 +2,8 @@
#from IPython import get_ipython
#get_ipython().magic('reset -sf')
# from side_tools import *
# clean()
from side_tools import *
clean()
from veros import VerosSetup, veros_method
from veros.tools import cli
......@@ -126,10 +126,12 @@ class DWBC2Setup(VerosSetup):
vs.kbot[3,:] = np.int(22) # slope
vs.kbot[4,:] = np.int(16)
vs.kbot[5,:] = np.int(10)
vs.kbot[6,:] = np.int(7)
vs.kbot[7,:] = np.int(5)
vs.kbot[8,:] = np.int(3)
vs.kbot[9,:] = np.int(2)
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[4:,30:70] = np.int(1) # steeper slope
......@@ -145,7 +147,7 @@ class DWBC2Setup(VerosSetup):
Flow initial and restoring profile
"""
zcenter = 1500
zvariance = 2e5
zvariance = 3e5
xvariance = 1.1
xcenter = 6.5 # normally 16
xx,zz = np.meshgrid(vs.xt,vs.zt,indexing='ij')
......
......@@ -63,7 +63,7 @@ p = section('dwbc2.snapshot.nc','v')
high=0.2
plt.figure()
p.select(0,'yu',50).plot.contourf(levels=30,cmap='RdBu_r')
p.select(10,'zt',5).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