ready version for mistral

veros/setup/dwbc2/dwbc2.py

@@ -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
@@ -44,9 +44,9 @@ class DWBC2Setup(VerosSetup):
         """
         GRID, TIMESTEP,...
         """
-        vs.nx, vs.ny, vs.nz = 20, 35, 30
-        vs.dt_mom = 1200 # normally 2400
-        vs.dt_tracer = 86400 / 4
+        vs.nx, vs.ny, vs.nz = 30, 60, 30
+        vs.dt_mom = 300 # normally 2400
+        vs.dt_tracer = 86400 / 16
         vs.runlen = 86400 * 365 * 1
         # vs.runlen = 86400 * 100
         # vs.restart_frequency = 86400 * 100
@@ -101,8 +101,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[...] = 1
-        vs.dyt[...] = 1
+        vs.dxt[...] = 0.25
+        vs.dyt[...] = 0.25
         vs.x_origin = 5.0
         vs.y_origin = 10
         
@@ -137,8 +137,8 @@ class DWBC2Setup(VerosSetup):
         """
         zcenter = 1500
         zvariance = 3e5
-        xvariance = 2.
-        xcenter =  7.5 # normally 16
+        xvariance = 1.2
+        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

veros/setup/dwbc2/plot.py

@@ -61,7 +61,7 @@ p = section('dwbc2.snapshot.nc','v')
 
 high=0.2
 plt.figure()
-p.select(360,'zt',15).plot.contourf(vmin=-high,vmax=high,levels=30,cmap='RdBu_r')
+p.select(350,'zt',15).plot.contourf(vmin=-high,vmax=high,levels=30,cmap='RdBu_r')
 # 
 
 # plt.colorbar()