import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step import Fluxes as Flux from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '2D_Bickley_Jet_H500' # Geometry and Model Equations sim.geomx = 'periodic' # Geometry Types: 'periodic' or 'walls' sim.geomy = 'walls' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.method = 'Spectral' # Numerical method: 'Spectral' sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' sim.flux_method = Flux.spectral_sw # Flux method: spectral_sw is only option currently # Specify paramters sim.Lx = 200e3 # Domain extent (m) sim.Ly = 200e3 # Domain extent (m) sim.Nx = 128 # Grid points in x sim.Ny = 128 # Grid points in y sim.Nz = 1 # Number of layers
import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step import Fluxes as Flux from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '2D GeoAdjust' # Geometry and Model Equations sim.geomx = 'periodic' # Geometry Types: 'periodic' or 'walls' sim.geomy = 'periodic' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.method = 'Spectral' # Numerical method: 'Spectral' sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' sim.flux_method = Flux.spectral_sw # Flux method: spectral_sw is only option currently # Specify paramters sim.Lx = 4000e3 # Domain extent (m) sim.Ly = 3000e3 # Domain extent (m) sim.Nx = 128 # Grid points in x sim.Ny = 128 # Grid points in y sim.Nz = 1 # Number of layers
import numpy as np import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '1D_Sadourny' # Geometry and Model Equations sim.geomy = 'walls' # Geometry Types: 'periodic' or 'walls' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' sim.method = 'Sadourny' # Numerical method: 'Sadourny' sim.method = 'Spectral' # Specify paramters sim.Lx = 4000e3 # Domain extent (m) sim.Ly = 4000e3 # Domain extent (m) sim.Nx = 1 # Grid points in x sim.Ny = 128 # Grid points in y sim.Nz = 1 # Number of layers sim.g = 9.81 # Gravity (m/sec^2)
import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step import Fluxes as Flux from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '1D GeoAdjust' # Geometry and Model Equations sim.geomy = 'periodic' # Geometry Types: 'periodic' or 'walls' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.method = 'Spectral' # Numerical method: 'Spectral' sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' sim.flux_method = Flux.spectral_sw # Flux method: spectral_sw is only option currently # Specify paramters sim.Ly = 4000e3 # Domain extent (m) sim.Nx = 1 # Grid points in x sim.Ny = 1024 # Grid points in y sim.Nz = 1 # Number of layers sim.g = 9.81 # Gravity (m/sec^2) sim.f0 = 1.e-4 # Coriolis (1/sec)
import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step #import Fluxes as Flux from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '2D_Bickley_Jet_Sadourny_N256' # Geometry and Model Equations sim.geomx = 'periodic' # Geometry Types: 'periodic' or 'walls' sim.geomy = 'walls' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.method = 'Sadourny' # Numerical method: 'Spectral' sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' # Specify paramters sim.Lx = 150e3 # Domain extent (m) sim.Ly = 100e3 # Domain extent (m) sim.Nx = 256 # Grid points in x sim.Ny = 256 # Grid points in y sim.Nz = 1 # Number of layers sim.g = 9.81 # Gravity (m/sec^2)
import matplotlib.pyplot as plt import sys # Add the PyRsw tools to the path # At the moment it is given explicitely. # In the future, it could also be added to the # pythonpath environment variable sys.path.append('../src') import Steppers as Step import Fluxes as Flux from PyRsw import Simulation from constants import minute, hour, day sim = Simulation() # Create a simulation object sim.run_name = '1D Bickley Jet' # Geometry and Model Equations sim.geomx = 'periodic' # Geometry Types: 'periodic' or 'walls' sim.geomy = 'walls' sim.stepper = Step.AB3 # Time-stepping algorithm: Euler, AB2, RK4 sim.method = 'Spectral' # Numerical method: 'Spectral' sim.dynamics = 'Nonlinear' # Dynamics: 'Nonlinear' or 'Linear' sim.flux_method = Flux.spectral_sw # Flux method: spectral_sw is only option currently # Specify paramters sim.Lx = 200e3 # Domain extent (m) sim.Ly = 200e3 # Domain extent (m) sim.Nx = 1 # Grid points in x sim.Ny = 128 # Grid points in y sim.Nz = 1 # Number of layers