from param import Param from grid import Grid from fluid2d import Fluid2d import numpy as np param = Param('default.xml') param.modelname = 'quasigeostrophic' param.expname = 'vortex_topo_0' # domain and resolution ratio = 2 param.nx = 64 * 4 param.ny = param.nx // ratio param.npx = 1 param.npy = 1 param.Lx = 2. param.Ly = param.Lx / ratio param.geometry = 'perio' # time param.tend = 50. param.cfl = 0.8 param.adaptable_dt = True param.dt = 1. param.dtmax = 100. # discretization param.order = 5 param.flux_splitting_method = 'minmax' # output
from param import Param from grid import Grid from fluid2d import Fluid2d import numpy as np param = Param("toto") param.modelname = 'sqg' param.expname = 'sqg_vortex' # domain and resolution param.nx = 64 * 2 param.ny = param.nx param.npy = 1 param.Lx = 1. param.Ly = 1. param.geometry = 'perio' # time param.tend = 40 param.cfl = 0.8 param.adaptable_dt = True param.dt = 1. param.dtmax = 100. # discretization param.order = 5 # output param.ageostrophic = False param.var_to_save = ['pv', 'psi', 'u', 'v', 'vorticity'] param.list_diag = ['pv', 'pv2']
from numpy import exp,sqrt,pi,cos,random,shape param = Param('default.xml') # the xml file contains all the defaults values + a doc on variable param.modelname = 'euler' # advection, euler, boussinesq, quasigeostrophic param.expname = 'myfirstexp' # string used to name the netcdf files: # expname_his.nc and expname_diag.nc # *** domain and resolution *** param.nx = 32*2 # nb of grid points, should of the form 2**p, 2**p*3 or 2**p*5 param.ny = 32 # idem param.npy=4 # nb of processors in y: should be a power of 2 param.Lx = 2. # domain size in x [Lx can be dimensional if the user # wishes] param.Ly = param.Lx/2 # domain size in y. The model imposes that dx=dy # => Lx/nx = Ly/ny param.geometry='square' # square, disc, xchannel, ychannel, # perio. 'square' is a closed domain. For a # more complex domain, adjust the grid.msk # once grid is available # *** time *** param.tend=4e4 # time of integration [could be dimensional of not] param.adaptable_dt=True
from grid import Grid from fluid2d import Fluid2d from numpy import exp,sqrt,pi,cos,where,random,shape from restart import Restart from island import Island param = Param('default.xml') param.modelname = 'euler' param.expname = 'karman_0' # domain and resolution ratio = 2 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='xchannel' # time param.tend=600 param.cfl=1. param.adaptable_dt=True param.dt = 1e-2 param.dtmax=1. # discretization param.order=5 param.timestepping='RK3_SSP'
from param import Param from grid import Grid from fluid2d import Fluid2d import numpy as np param = Param('default.xml') param.modelname = 'boussinesq' param.expname = 'RB_00' # domain and resolution param.nx = 64 * 2 param.ny = param.nx / 4 param.npx = 1 param.Lx = 4. param.Ly = 1. param.geometry = 'xchannel' # time param.tend = 40. param.cfl = 1. param.adaptable_dt = True param.dt = .1 param.dtmax = .1 param.exacthistime = False # discretization param.order = 5 param.aparab = 0.02 # output param.plot_var = 'buoyancy'
from fluid2d import Fluid2d from numpy import exp, sqrt, pi, cos, sin, where, random, shape, tanh, cumsum, cosh import numpy as np from restart import Restart from island import Island param = Param('default.xml') param.modelname = 'boussinesq' param.expname = 'Wave_NH' # domain and resolution ratio = 2 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 = 10 param.cfl = 0.3 param.adaptable_dt = False param.dt = 2e-2 param.dtmax = 1. # discretization param.order = 5
from param import Param from grid import Grid from fluid2d import Fluid2d import numpy as np param = Param('default.xml') param.modelname = 'boussinesq' param.expname = 'khi_0' # domain and resolution ratio = 2 param.ny = 64*2 param.nx = param.ny*ratio param.Ly = 1. param.Lx = 1*ratio param.npx = 1 param.npy = 1 param.geometry = 'xchannel' # time param.tend = 5. param.cfl = 1.5 param.adaptable_dt = True param.dt = 0.01 param.dtmax = 0.02 # discretization param.order = 5 param.timestepping = 'RK3_SSP' # output