# Parallel? Only applies to the FFTWs sim.num_threads = 4 # Plotting parameters sim.plott = 15.*minute # Period of plots sim.animate = 'Anim' # 'Save' to create video frames, # 'Anim' to animate, # 'None' otherwise sim.plot_vars = ['h'] #sim.plot_vars = ['vort','div'] #sim.clims = [[-0.015, 0.015],[-0.001, 0.001]] # Output parameters sim.output = False # True or False sim.savet = 1.*hour # Time between saves # Diagnostics parameters sim.diagt = 2.*minute # Time for output sim.diagnose = False # True or False # Initialize the grid and zero solutions sim.initialize() for ii in range(sim.Nz): # Set mean depths sim.soln.h[:,:,ii] = sim.Hs[ii] # Gaussian initial conditions W = 200.e3 # Width amp = 1. # Amplitude sim.soln.h[:,:,0] += amp*np.exp(-(sim.X/W)**2 - (sim.Y/W)**2)
sim.num_threads = 4 # Plotting parameters sim.plott = 15.*minute # Period of plots sim.animate = 'Anim' # 'Save' to create video frames, # 'Anim' to animate, # 'None' otherwise sim.plot_vars = ['u','v','h'] # Specify which variables to plot # Specify manual ylimits if desired # An empty list uses default limits sim.ylims=[[-0.18,0.18],[-0.18,0.18],[-0.5,1.0]] # Output parameters sim.output = False # True or False sim.savet = 1.*hour # Time between saves # Diagnostics parameters sim.diagt = 2.*minute # Time for output sim.diagnose = False # True or False # Initialize the grid and zero solutions sim.initialize() for ii in range(sim.Nz): # Set mean depths sim.soln.h[:,:,ii] = sim.Hs[ii] # Gaussian initial conditions x0 = 1.*sim.Lx/2. # Centre W = 200.e3 # Width amp = 1. # Amplitude
sim.end_time = 20 * 24. * hour # End Time (sec) # Parallel? Only applies to the FFTWs sim.num_threads = 32 # Plotting parameters sim.plott = 12. * hour # Period of plots sim.animate = 'Save' # 'Save' to create video frames, # 'Anim' to animate, # 'None' otherwise sim.plot_vars = ['vort', 'v', 'u', 'h'] sim.clims = [[-0.8, 0.8], [-0.5, 0.5], [], []] # Output parameters sim.output = True # True or False sim.savet = 10. * day # Time between saves # Diagnostics parameters sim.diagt = 2. * minute # Time for output sim.diagnose = False # True or False # Initialize the grid and zero solutions sim.initialize() for ii in range(sim.Nz): # Set mean depths sim.soln.h[:, :, ii] = sim.Hs[ii] # Bickley Jet initial conditions # First we define the jet Ljet = 10e3 # Jet width amp = 0.1 # Elevation of free-surface in basic state
sim.end_time = 20*24.*hour # End Time (sec) # Parallel? Only applies to the FFTWs sim.num_threads = 32 # Plotting parameters sim.plott = 1.*hour # Period of plots sim.animate = 'Save' # 'Save' to create video frames, # 'Anim' to animate, # 'None' otherwise sim.plot_vars = ['vort', 'v', 'u', 'h'] sim.clims = [ [-0.8, 0.8], [-0.5,0.5], [], []] # Output parameters sim.output = True # True or False sim.savet = 10.*day # Time between saves # Diagnostics parameters sim.diagt = 2.*minute # Time for output sim.diagnose = False # True or False # Initialize the grid and zero solutions sim.initialize() for ii in range(sim.Nz): # Set mean depths sim.soln.h[:,:,ii] = sim.Hs[ii] # Bickley Jet initial conditions # First we define the jet Ljet = 10e3 # Jet width amp = 0.1 # Elevation of free-surface in basic state