コード例 #1
0
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
sim.g = 9.81  # Gravity                     (m/sec^2)
sim.f0 = 1.e-4  # Coriolis                    (1/sec)
sim.beta = 0e-10  # Coriolis beta               (1/m/sec)
sim.Hs = [500.]  # Vector of mean layer depths (m)
sim.rho = [1025.]  # Vector of layer densities   (kg/m^3)
sim.end_time = 20 * 24. * hour  # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 32

# Plotting parameters
sim.plott = 12. * hour  # Period of plots
sim.animate = 'Save'  # 'Save' to create video frames,
# 'Anim' to animate,
# 'None' otherwise
sim.plot_vars = ['vort', 'v', 'u', 'h']
sim.clims = [[-0.8, 0.8], [-0.5, 0.5], [], []]

# Output parameters
sim.output = True  # True or False
sim.savet = 10. * day  # Time between saves

# Diagnostics parameters
sim.diagt = 2. * minute  # Time for output
sim.diagnose = False  # True or False

# Initialize the grid and zero solutions
コード例 #2
0
ファイル: example_2D_geoadjust.py プロジェクト: PyRsw/PyRsw
sim.Nx  = 128               # 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)
sim.f0  = 1.e-4           # Coriolis                    (1/sec)
sim.beta = 0e-11          # Coriolis beta parameter     (1/m/sec)
sim.cfl = 0.1             # CFL coefficient             (m)
sim.Hs  = [100.]          # Vector of mean layer depths (m)
sim.rho = [1025.]         # Vector of layer densities   (kg/m^3)
sim.end_time = 2.*24.*hour   # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 4

# Plotting parameters
sim.plott   = 15.*minute  # Period of plots
sim.animate = 'Anim'      # 'Save' to create video frames,
                          # 'Anim' to animate,
                          # 'None' otherwise
sim.plot_vars = ['h']
#sim.plot_vars = ['vort','div']
#sim.clims = [[-0.015, 0.015],[-0.001, 0.001]]
                         
# Output parameters
sim.output = False        # True or False
sim.savet  = 1.*hour      # Time between saves

# Diagnostics parameters
sim.diagt    = 2.*minute  # Time for output
sim.diagnose = False      # True or False
コード例 #3
0
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)
sim.beta = 0e-10          # Coriolis beta parameter     (1/m/sec)
sim.Hs  = [100.]          # Vector of mean layer depths (m)
sim.rho = [1025.]         # Vector of layer densities   (kg/m^3)
sim.end_time = 2*24.*hour   # End Time                    (sec)

# Parallel: Only applies to the FFTWs
sim.num_threads = 4

# Plotting parameters
sim.plott   = 15.*minute  # Period of plots
sim.animate = 'Anim'      # 'Save' to create video frames,
                          # 'Anim' to animate,
                          # 'None' otherwise

sim.plot_vars = ['u','v','h']   # Specify which variables to plot
                                # Specify manual ylimits if desired
                                # An empty list uses default limits
sim.ylims=[[-0.18,0.18],[-0.18,0.18],[-0.5,1.0]] 

# Output parameters
sim.output = False        # True or False
sim.savet  = 1.*hour      # Time between saves

# Diagnostics parameters
sim.diagt    = 2.*minute  # Time for output
コード例 #4
0
ファイル: example_1D_geoadjust2.py プロジェクト: PyRsw/PyRsw
sim.Nx = 1  # Grid points in x
sim.Ny = 512  # Grid points in y
sim.Nz = 1  # Number of layers
sim.g = 9.81  # Gravity                     (m/sec^2)
sim.f0 = 1.0e-4  # Coriolis                    (1/sec)
sim.beta = 0e-10  # Coriolis beta parameter     (1/m/sec)
sim.cfl = 0.15  # CFL coefficient             (m)
sim.Hs = [100.0]  # Vector of mean layer depths (m)
sim.rho = [1025.0]  # Vector of layer densities   (kg/m^3)
sim.end_time = 6.0 * 24.0 * hour  # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 4

# Plotting parameters
sim.plott = 20.0 * minute  # Period of plots
sim.animate = "Anim"  # 'Save' to create video frames,
# 'Anim' to animate,
# 'None' otherwise

sim.plot_vars = ["vort", "div", "h"]  # Specify which variables to plot
# sim.plot_vars = ['u','v','h']   # Specify which variables to plot
# Specify manual ylimits if desired
# An empty list uses default limits
sim.ylims = [[-0.01, 0.01], [-2.0, 2.0], [-1.0, 1.0]]
# sim.ylims=[[-0.3,0.3],[-0.2,0.2],[-1.0,1.0]]

# Output parameters
sim.output = False  # True or False
sim.savet = 1.0 * hour  # Time between saves
コード例 #5
0
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)
sim.f0 = 1.e-4  # Coriolis                    (1/sec)
sim.beta = 0e-10  # Coriolis beta               (1/m/sec)
sim.cfl = 0.1  # CFL coefficient             (m)
sim.Hs = [100.]  # Vector of mean layer depths (m)
sim.rho = [1025.]  # Vector of layer densities   (kg/m^3)
sim.end_time = 14 * 24. * hour  # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 4

# Plotting parameters
sim.plott = 30. * minute  # Period of plots
sim.animate = 'Save'  # 'Save' to create video frames,
# 'Anim' to animate,
# 'None' otherwise
sim.plot_vars = ['h']
sim.ylims = [[-0.1, 0.1]]
#sim.plot_vars = ['vort','div']
#sim.clims = [ [-0.8, 0.8],[-0.1, 0.1]]

# Output parameters
sim.output = False  # True or False
sim.savet = 1. * hour  # Time between saves

# Diagnostics parameters
sim.diagt = 2. * minute  # Time for output
sim.diagnose = False  # True or False
コード例 #6
0
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)
sim.f0  = 1.e-4           # Coriolis                    (1/sec)
sim.beta = 0e-10          # Coriolis beta               (1/m/sec)
sim.Hs  = [500.]          # Vector of mean layer depths (m)
sim.rho = [1025.]         # Vector of layer densities   (kg/m^3)
sim.end_time = 20*24.*hour   # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 32

# Plotting parameters
sim.plott   = 1.*hour  # Period of plots
sim.animate = 'Save'      # 'Save' to create video frames,
                          # 'Anim' to animate,
                          # 'None' otherwise
sim.plot_vars = ['vort', 'v', 'u', 'h']
sim.clims = [ [-0.8, 0.8], [-0.5,0.5], [], []]                         

# Output parameters
sim.output = True        # True or False
sim.savet  = 10.*day      # Time between saves

# Diagnostics parameters
sim.diagt    = 2.*minute  # Time for output
sim.diagnose = False      # True or False

# Initialize the grid and zero solutions
コード例 #7
0
ファイル: example_1D_BickleyJet.py プロジェクト: PyRsw/PyRsw
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)
sim.f0  = 1.e-4           # Coriolis                    (1/sec)
sim.beta = 0e-10          # Coriolis beta               (1/m/sec)
sim.cfl = 0.1             # CFL coefficient             (m)
sim.Hs  = [100.]          # Vector of mean layer depths (m)
sim.rho = [1025.]         # Vector of layer densities   (kg/m^3)
sim.end_time = 14*24.*hour   # End Time                    (sec)

# Parallel? Only applies to the FFTWs
sim.num_threads = 4

# Plotting parameters
sim.plott   = 30.*minute  # Period of plots
sim.animate = 'Save'      # 'Save' to create video frames,
                          # 'Anim' to animate,
                          # 'None' otherwise
sim.plot_vars = ['h']
sim.ylims=[[-0.1,0.1]] 
#sim.plot_vars = ['vort','div']
#sim.clims = [ [-0.8, 0.8],[-0.1, 0.1]]                         

# Output parameters
sim.output = False        # True or False
sim.savet  = 1.*hour      # Time between saves

# Diagnostics parameters
sim.diagt    = 2.*minute  # Time for output
sim.diagnose = False      # True or False