Fortran Wrotsplat

Design of the program wrotsplat.f90 (09 Dec 2004)

• Using the rotate-and-splat method, since it is simpler.
• we want to calculate the velocity moments of different emission lines.

Data structures

Numerical grid is (r,z,θ). Each cell makes contributions to the velocity moments:

\[ M_n = \int \eta \, v_z^n \, dz \]

where $v_z = \mathbf{v} \cdot \mathbf{n}$ and $\mathbf{n}$ is the LOS vector.

  ! velocity moments of the emission along the LOS
  ! for each of the different line types
  type moment_t 
    real(DP), dimension(3) :: er, ei, en
  end type moment_t
  ! big map of all moments of all lines
  type(moment_t), allocatable :: momentmap(:,:)