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:

where and 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(:,:)