For help: get output variables form Rayleigh

Rayleigh
1

Hi all,

I am new to Rayleigh and I am trying to simulate magnetic filed by running it. To prepare the input file and get results, I copied the &output_namelist from benchmark_diagnostics_input file to the main_input file I was using. However, I found that there are no magnetic variables in the*&output_namelist* (please see below). I was wondering should I add the magnetic variables in the main_input file?
Also I did find the magnetic variables listed in Rayleigh_Output_Variables.pdf. But I don’t know how should I write the code to make it work. Appreciate if you can help or have any suggestions!

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The out_namelist from benchmark_diagnostics_input:

&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.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, 1438, 1468
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
/

%%%%%%%%%%%%%%%%
snap shot of the variable list:

image
image.png1266Ă—936 89.1 KB

2

Hi Zhou,
Glad to hear of someone new using Rayleigh! The short answer to your question is “yes,” in general you can just add the magnetic variables to the output list BUT your input file needs to be set up to use magnetism. An example input file is the c2001_case1_input file, and details are discussed in Rayleigh_user_Guide.pdf. Note that you may also find the Diagnostic_Plotting.html/pdf files useful as well (html looks a lot better than pdf as far as formatting goes).
All the Best,
Nick

1 Like
3

Thanks Nick. With your reply I reread the user guide, and was able to add magnetic parameters into the output. It is easy by simply adding the the variable codes into the output lines. Those new-added variables were also able to be visualized by adjusting the Diagnostic_plotting code.

Thanks for your answer.