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
# 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