"""
Parameter file for use in the Dedalus 2D anelastic convection script.
"""

import numpy as np

Lx, Lz = 2, 1                       # Domain size
Nx, Nz = 256, 128                    # Number of
Pr = 1.                             # Prandtl number
Ra = 1.8e8                          # Rayleigh number
Ta = 1e7
Np = 0.5                              # Number of density scale heights
m = 1.5                             # Polytropic index
Phi = np.pi/4
theta = 1 - np.exp(-Np/m)           # Dimensionaless inverse T scale height

initial_timestep = 1.5e-6           # Initial timestep
max_dt = 1e-5                       # max dt

snapshot_freq = 1.5e-3              # Frequency snapshot files are outputted
analysis_freq = 1.5e-4              # Frequency analysis files are outputted

end_sim_time = 2.                   # Stop time in simulations units
end_wall_time = np.inf              # Stop time in wall time
end_iterations = np.inf             # Stop time in iterations
