Graphic files G_#.TAG and G_ave.TAG¶
These are fortran unformatted files containing 3D data (in the form
vector_array(phi, theta, r) ) which can be used to visualize the solution.
They are written after a fixed number of time steps as specified by the user
in the Output Control namelist using the parameters
listed in the section on output of graphic files. In case
l_average is set to .true., then an average graphic file,
named G_ave.TAG, containing time averaged values of 3D data, is also written
at the end of the simulation.
These files are written in chunks of latitude for one radial level at a time by the
subroutine graphOut or by
graphOut_mpi depending on whether
USE_MPI is set to Yes or No in the Makefile. The structure of the
file looks like below:
!------------- ! Line 1 !------------- version !Graphout_version_9 (using MPI without comp. without pressure) !Graphout_version_10 (using MPI, without comp. with pressure) !Graphout_version_11 (using MPI, with comp. without pressure) !Graphout_version_12 (using MPI, with comp. with pressure) !Graphout_version_5 (without MPI, with pressure and comp.) !Graphout_version_6 (without MPI, with comp. without pressure) !Graphout_version_7 (without MPI, without comp. without pressure) !Graphout_version_8 (without MPI, without comp. with pressure) !---------- ! Line 2 !---------- runid !---------- ! Line 3 !---------- time, n_r_max, n_theta_max, n_phi_tot, !time = Time of writing n_r_ic_max-1, minc, nThetasBs, !(Simulation time), ra, ek, pr, prmag, !nThetasBs = no. of radratio, sigma_ratio !theta blocks !---------- ! Line 4 !---------- theta(1:n_theta_max) !----------------------------------------------------------------------- !----------------------------- !Graphout_version_[9/10/11/12] !----------------------------- ! These versions are written when the code uses MPI (USE_MPI=yes). Parallel ! chunks of fields are written for different radial levels. Chunks in theta ! are written in parallel using OpenMP !--------------- ! Data !--------------- !------------- ! Block N !------------- !------------ ! Line 4 + N !------------ n_r-1, r(n_r)/r(1), n_theta_start, n_theta_stop !Radial index, radius in terms !of r_cmb, start and stop of !the theta block !--------------- ! Line 4 + (N+1) !--------------- sr(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Entropy !--------------- ! Line 4 + (N+2) !--------------- vr(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Radial velocity !--------------- ! Line 4 + (N+3) !--------------- vt(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Theta component of velocity !--------------- ! Line 4 + (N+4) !--------------- vp(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Zonal (phi component) of !velocity if (l_chemical_conv): !If composition is stored !--------------- ! Line 4 + (N+5) !--------------- xir(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !composition if (l_PressGraph): !If pressure is stored !----------------- ! Line 4 + (N+6/7) !----------------- pr(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !pressure if (l_mag): !For a magnetic run !------------------- ! Line 4 + (N+5/6/7) !------------------- br(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Radial magnetic field !------------------- ! Line 4 + (N+6/7/8) !------------------- bt(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Theta component of !magnetic field !------------------- ! Line 4 + (N+7/8/9) !------------------- bp(1:n_phi_tot, n_theta_start:n_theta_stop, n_r) !Zonal (phi component) !of magnetic field !-------------------------- !Graphout_version_[5/6/7/8] !-------------------------- !This version is written when the code does not use MPI (USE_MPI=no). !Chunks in theta are written in parallel with OpenMP. !--------------- ! Data !--------------- !------------- ! Block N !------------- !--------------- ! Line 4 + (N+1) !--------------- n_r-1, r(n_r)/r(1), n_theta_start, n_theta_stop !---------------------------------------------------------- ! Each of the following data point is written in a new line !---------------------------------------------------------- !----------------- ! Entropy !----------------- sr(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r sr(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... sr(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r sr(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... sr(n_phi_tot,n_theta_start+1,n_r) ... sr(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r sr(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... sr(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !----------------- ! Radial velocity !----------------- vr(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r vr(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... vr(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r vr(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... vr(n_phi_tot,n_theta_start+1,n_r) ... vr(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r vr(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... vr(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !----------------------------- ! Theta component of velocity !----------------------------- vt(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r vt(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... vt(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r vt(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... vt(n_phi_tot,n_theta_start+1,n_r) ... vt(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r vt(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... vt(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !---------------------------------- ! Zonal (phi component) of velocity !---------------------------------- vp(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r vp(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... vp(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r vp(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... vp(n_phi_tot,n_theta_start+1,n_r) ... vp(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r vp(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... vp(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r if (l_chemical_conv): !If chemical composition is stored !------------- ! Composition !------------- xi(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r xi(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... xi(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r xi(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... xi(n_phi_tot,n_theta_start+1,n_r) ... xi(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r xi(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... xi(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r if (l_PressGraph): !If pressure is stored !----------- ! Pressure !----------- pr(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r pr(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... pr(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r pr(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... pr(n_phi_tot,n_theta_start+1,n_r) ... pr(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r pr(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... pr(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r if (l_mag): !Only if it is a magnetic case !---------------------- ! Radial magnetic field !---------------------- br(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r br(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... br(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r br(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... br(n_phi_tot,n_theta_start+1,n_r) ... br(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r br(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... br(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !---------------------------------- ! Theta component of magnetic field !---------------------------------- bt(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r bt(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... bt(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r bt(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... bt(n_phi_tot,n_theta_start+1,n_r) ... bt(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r bt(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... bt(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !---------------------------------------- ! Zonal (phi component) of magnetic field !---------------------------------------- bp(1,n_theta_start,n_r) !n_phi = 1, n_theta = n_theta_start, n_r bp(2,n_theta_start,n_r) !n_phi = 2, n_theta = n_theta_start, n_r ... bp(n_phi_tot,n_theta_start,n_r) !n_phi = n_phi_tot, n_theta = n_theta_start, n_r bp(1,n_theta_start+1,n_r) !n_phi = 1, n_theta = n_theta_start+1, n_r ... bp(n_phi_tot,n_theta_start+1,n_r) ... bp(1,n_theta_stop,n_r) !n_phi = 1, n_theta = n_theta_stop, n_r bp(2,n_theta_stop,n_r) !n_phi = 2, n_theta = n_theta_stop, n_r ... bp(n_phi_tot,n_theta_stop,n_r) !n_phi = n_phi_tot, n_theta = n_theta_stop, n_r !----------------- !Subsequent blocks !----------------- !Block N+1 in both cases have data at the same radial level but the next !theta chunk (n_theta_start + nThetaB, n_theta_stop + n_thetaB) !After data for all the theta blocks have been written for one radial !level, everything above is repeated for the next radial level
The graphic files can be read using the python class MagicGraph.
>>> gr = MagicGraph(ivar = 1, tag='TAG')
>>> # print radial velocity
>>> print(gr.vr)
They can be visualized using the Surf class:
>>> s = Surf(tag='TAG')
>>> # Surface map of radial velocity:
>>> s.surf(field = 'vr', r = 0.5, cm = 'jet', levels = 50)
>>> s.slice(field = 'br', lon_0 = [0]) # Longitudinal Slice of radial magnetic field
>>> s.equat(field = 'entropy') # Equatorial slice of entropy