Hi,
I am running Rayleigh without mentioning any symmetry in main input file, but it is expressing fourfold symmetry.
Initially, I run this example:
Contents — Rayleigh documentation and get the plots of various outputs described in the link. The Equatorial_Slices obtained are fourfold symmetric.
When I run the code with higher Rayleigh number this symmetry is still imposed to the solution. I have tried various inputs changes such as omitting benchmark at all, but nothing changed.
Hope this email address my question. I am attaching one of the main input files (which I used to run with different Rayleigh numbers) as your review.
Please investigate the issue.
&problemsize_namelist
n_r = 48
n_theta = 64
nprow = 4
npcol = 4
aspect_ratio = 0.35d0
shell_depth = 1.0d0
/
&numerical_controls_namelist
/
&physical_controls_namelist
rotation = .True.
magnetism = .false.
viscous_heating = .false.
ohmic_heating = .false.
advect_reference_state = .false.
/
&temporal_controls_namelist
max_time_step = 1.0d-4
max_iterations = 40000
checkpoint_interval = 5000
cflmin = 0.4d0
cflmax = 0.6d0
/
&io_controls_namelist
/
&output_namelist
! New equatorial slices output
equatorial_values = 1,2,3 ! velocity components
equatorial_frequency = 5000
equatorial_nrec = 2
!New meridional slices output. You must also specify the number of phi-indices
!Valid values range from 1 through nphi = 2*n_theta
!The indices below range from phi = 0 up to phi = pi
meridional_values = 1,2,3 ! radial and phi components of velocity; temperature
meridional_frequency = 5000
meridional_nrec = 2
meridional_indices_nrm = 0.7
shellslice_levels_nrm =0.5, 0.9
shellslice_values = 1,2,3 ! velocity components
shellslice_frequency = 10000
shellslice_nrec = 2
shellspectra_levels_nrm = 0.5 , 0.9
shellspectra_values = 1,2,3 ! velocity and temperature
shellspectra_frequency = 10000
shellspectra_nrec = 2
azavg_values = 1,2,3,201,202,501 ! same as above + r- and theta- mass flux
azavg_frequency = 1000
azavg_nrec = 10
point_probe_values = 1,2,3, 1201,1202,1203,1216, 1219,1220, 1221, 1228,1229,1230
point_probe_r_nrm = 0.2, 0.6, 0.9
point_probe_theta_nrm= 0.35
point_probe_phi_nrm = 0.21, 0.72
point_probe_frequency=10
point_probe_nrec=500
point_probe_cache_size=100
! velocity, temperature, energy fluxes, and Kinetic Energy
shellavg_values = 1,2,3,501, 1440, 1470
shellavg_frequency = 200
shellavg_nrec = 50
! Kinetic energy, Mean KE, Diff-Rot KE, and Convective KE
globalavg_values = 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, 411, 412
globalavg_frequency = 100
globalavg_nrec = 100
sph_mode_ell = 2,4,8
sph_mode_levels_nrm = 0.4, 0.7, 0.9
sph_mode_values = 1,2,3
sph_mode_frequency = 100
sph_mode_nrec = 50
full3d_values = 64 ! temperature
full3d_frequency = 9000000
/
&Boundary_Conditions_Namelist
no_slip_boundaries = .true.
strict_L_Conservation = .false.
dtdr_bottom = 0.0d0
T_Top = 0.0d0
T_Bottom = 1.0d0
fix_tvar_top = .true.
fix_tvar_bottom = .true.
/
&Initial_Conditions_Namelist
init_type = -1 ! Benchmark init
temp_amp = 1.0d1
temp_w = 0.01d4
restart_iter = 00025000
/
&Test_Namelist
/
&Reference_Namelist
Ekman_Number = 1.0d-3
Rayleigh_Number = 1.0d5
Prandtl_Number = 1.0d0
Magnetic_Prandtl_Number = 5.0d0
reference_type = 1
heating_type = 0 ! No heating
gravity_power = 1.0d0 ! g ~ radius
/
&Transport_Namelist
/