A relevant routine for the helical forcing is:

  st_calcPhases.F90

First it takes the solenoidal component and then it chooses real
and imaginary components appropriate for the helical forcing.
