Exemplo n.º 1
0
param = Param('default.xml')
param.modelname = 'euler'
param.expname = 'freedecay'

# domain and resolution
param.nx = 64 * 2
param.ny = param.nx
param.Ly = param.Lx
param.npx = 1
param.npy = 1
param.geometry = 'perio'

# time
param.tend = 1200.
param.cfl = 1.2
param.adaptable_dt = True
param.dt = .05
param.dtmax = 10.

# discretization
param.order = 5

# output
param.var_to_save = ['vorticity', 'psi', 'tracer']
param.list_diag = 'all'
param.freq_plot = 5
param.freq_his = 10
param.freq_diag = 1
param.exacthistime = False

# plot
Exemplo n.º 2
0
# domain and resolution
ratio = 1
param.ny = 64 * 2
param.nx = param.ny * ratio
param.Ly = 1.
param.Lx = param.Ly * ratio
param.npx = 1
param.npy = 1
param.geometry = 'closed'

param.hydroepsilon = 1.

# time
param.tend = 50.
param.cfl = 0.1
param.adaptable_dt = False
param.dt = 1e-2
param.dtmax = .1

# discretization
param.order = 5
param.timestepping = 'RK3_SSP'  # 'LFAM3'

# output
param.var_to_save = ['vorticity', 'psi', 'u', 'v', 'buoyancy', 'tracer']
param.list_diag = ['ke', 'pe', 'energy', 'brms', 'vorticity', 'enstrophy']
param.freq_plot = 20
param.freq_his = .2
param.freq_diag = .1

# plot