1 SUBROUTINE denmf(trigsu, trigsv, ifaxu, ifaxv, irho)
15 INTEGER,
DIMENSION(13) :: ifaxu, ifaxv
16 REAL(rprec),
DIMENSION(3*nthetah/2 + 1) :: trigsu
17 REAL(rprec),
DIMENSION(2*nzetah) :: trigsv
21 REAL(rprec),
PARAMETER :: zero = 0, d18 = 1.0e-18_dp,
26 REAL(rprec),
DIMENSION(nthetah + 2,nzetah) :: a11
35 IF (isymm0 .eq. 0)
THEN
36 a11(:nthetah,:nzetah) =
37 1 (abs(one - bfield(:nthetah,:nzetah)) + d18)**1.5_dp
38 2 *b2avg(irho)/(bmax1(irho)*bfield(:nthetah,:nzetah))**2
39 a11(nthetah+1,:nzetah) = zero
40 a11(nthetah+2,:nzetah) = zero
42 CALL do_fft (a11, dmn, trigsu, trigsv, ifaxu, ifaxv, nthetah,
43 1 nzetah, mbuse, nbuse)
47 avgbobm2 = b2avg(irho)/bmax1(irho)**2