Calculate the parallel spectrum
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex, | intent(in), | dimension(:,:,:) | :: | an | ||
complex, | intent(out), | dimension(:,:,:) | :: | an2 |
subroutine par_spectrum(an, an2)
use gs2_transforms, only: kz_spectrum
use theta_grid, only: ntgrid
implicit none
complex, dimension(:,:,:), intent(in) :: an
complex, dimension(:,:,:), intent(out) :: an2
real :: scale
call kz_spectrum (an, an2)
scale = 1./real(4*ntgrid**2)
an2 = an2*scale
end subroutine par_spectrum