# plot param.plot_var = 'pv' param.freq_plot = 10 a = 4. param.cax = [-a, a] param.plot_interactive = True param.plot_psi = True param.plot_pvback = True param.colorscheme = 'symmetric' param.imshow_interpolation = 'bilinear' param.generate_mp4 = True # physics param.beta = 1. param.Rd = 0.1 # 2*grid.dx param.bottom_torque = True param.forcing = False param.forcing_module = 'forcing_dipoles' param.noslip = False param.diffusion = False grid = Grid(param) param.Kdiff = 0.1e-3 * grid.dx # add an island # grid.msk[28:32,34:38]=0 # grid.finalize() f2d = Fluid2d(param, grid) model = f2d.model
param.list_diag = ['pv', 'pv2']
param.freq_his = 1
param.freq_diag = 1
# plot
param.plot_var = 'pv'
param.freq_plot = 10
a = 0.2
param.cax = [-a, a]
param.plot_interactive = False
param.colorscheme = 'symmetric'
param.generate_mp4 = False
# physics
param.beta = 0.
param.Rd = 10. # 2*grid.dx
param.forcing = False
param.noslip = False
param.diffusion = False
grid = Grid(param)
param.Kdiff = 0.1e-3 * grid.dx
f2d = Fluid2d(param, grid)
model = f2d.model
xr, yr = grid.xr, grid.yr
vor = model.var.get('pv')
def vortex(param, grid, x0, y0, sigma, vortex_type, ratio=1):
param.list_diag = ['ke', 'pv', 'pv2'] param.freq_his = 50 param.freq_diag = 10 # plot param.plot_var = 'pv' param.freq_plot = 10 a = 0.5 param.cax = [-a, a] param.plot_interactive = True param.colorscheme = 'imposed' param.generate_mp4 = False # physics param.beta = 1. param.Rd = 0.01 # 2*grid.dx param.forcing = False param.noslip = False param.diffusion = False grid = Grid(param) param.Kdiff = 0.1e-3 * grid.dx # add an island # grid.msk[28:32,34:38]=0 # grid.finalize() f2d = Fluid2d(param, grid) model = f2d.model xr, yr = grid.xr, grid.yr
# *** plot ***
param.plot_var='pv' # variable to plot (depends on the model)
param.freq_plot=10 # number of time step between two plots
param.cax=[-1,1] # colorscale if colorscheme='imposed'
param.plot_interactive=False # activate the interactive plot
param.colorscheme='imposed' # 'imposed', 'minmax', 'symmetric'
param.plotting_module='plotting_adv' # name of the plotting
# script. The default one is
# core/plotting.py
# *** physics ***
param.beta=1. # beta parameter (for QG model only)
param.Rd=1. # Rossby deformation radius (for QG only)
param.gravity=1. # for Boussinesq model
param.forcing = True # activate the forcing
param.forcing_module='forcing_dbl_gyre' #python module name where the forcing is defined
param.noslip = False # active the no-slip boundary condition (causes energy dissipation)
param.diffusion=False # activate a diffusion on tracer.
param.additional_tracer=['tracer'] # you may add an additional passive
# tracer, by simply adding this
# line. Example in inverse_cascade.py
grid = Grid(param) # define the grid coordinates grid.xr, grid.yr and
# the mask, grid.msk. You may modify the mask just below
param.Kdiff=0.5e-4*grid.dx # diffusion coefficient (same for all
# tracers), there is a possibility to
param.list_diag = 'all' param.freq_his = 50 param.freq_diag = 10 # plot param.plot_var = 'pv' param.freq_plot = 10 a = 0.5 param.cax = [-a, a] param.plot_interactive = True param.colorscheme = 'imposed' param.generate_mp4 = True # physics param.beta = 1. param.Rd = 1. # 2*grid.dx param.forcing = True param.noslip = False param.diffusion = False grid = Grid(param) param.Kdiff = 0.5e-4 * grid.dx # add an island # grid.msk[28:32,34:38]=0 # grid.finalize_msk() f2d = Fluid2d(param, grid) model = f2d.model # set the forcing
param.list_diag = 'all'
param.freq_his = 20
param.freq_diag = 5.
# plot
param.plot_var = 'pv'
param.freq_plot = 10
a = 0.5
param.cax = [-a, a]
param.plot_interactive = True
param.colorscheme = 'imposed'
param.generate_mp4 = False
# physics
param.beta = 1.
param.Rd = .1
param.forcing = False
param.forcing_module = 'forcing' # not yet implemented
param.noslip = False
param.diffusion = False
param.isisland = True
psi0 = -5e-4 # this sets psi on the Northern wall (psi=0 on the Southern wall)
grid = Grid(param)
nh = grid.nh
def disc(param, grid, x0, y0, sigma):
""" function to set up a circular island
param.list_diag = ['ke', 'pv', 'pv2'] param.freq_plot = 5 param.freq_his = 10 param.freq_diag = 1 # plot param.plot_interactive = True param.plot_psi = True param.cmap = 'inferno' param.plot_var = 'pv' param.cax = np.array([-1, 1]) * .5 param.colorscheme = 'imposed' # physics param.beta = 1. param.Rd = 50. # 2*grid.dx param.forcing = False param.noslip = False param.diffusion = False # # you may activate the forcing and use the forcing below # it's a white noise forcing with some time correlation # such forcing is often used in turbulence studies # param.forcing_module='forcing_euler' grid = Grid(param) param.Kdiff = 5e-4 * grid.dx f2d = Fluid2d(param, grid) model = f2d.model
