1 SUBROUTINE amplitud(rcenter, zcenter, angin, r0c, z0c, rhoc, rhos,
12 REAL(rprec) rcenter, zcenter, r0c, z0c
13 REAL(rprec),
DIMENSION(ntheta) :: angin, xpts, xin, yin
14 REAL(rprec),
DIMENSION(0:mrho-1) :: rhoc, rhos
20 REAL(rprec) :: xmult, arg, xi, yi, t1, t2,
21 1 r1c(mu), r1s(mu), z1c(mu), z1s(mu), tnorm
30 xpts(:ntheta) = angin(:ntheta)
45 xi = xmult*(xin(j)-rcenter)
46 yi = xmult*(yin(j)-zcenter)
47 r1c(m) = r1c(m) + cos(m*arg)*xi
48 r1s(m) = r1s(m) + sin(m*arg)*xi
49 z1c(m) = z1c(m) + cos(m*arg)*yi
50 z1s(m) = z1s(m) + sin(m*arg)*yi
54 r10 = sqrt( r1c(1)**2 + r1s(1)**2 + z1c(1)**2 + z1s(1)**2 )
55 WRITE(3,
'(/,3(a,1pe10.3))')
56 1
' RAXIS = ', rcenter,
' ZAXIS = ', zcenter,
' R10 = ',r10
57 WRITE(*,
'(/,3(a,1pe10.3))')
58 1
' RAXIS = ', rcenter,
' ZAXIS = ', zcenter,
' R10 = ',r10
67 rhoc(m) = 0.5_dp*(r1c(m+1) + z1s(m+1))/t1
68 rhos(m) = 0.5_dp*(r1s(m+1) - z1c(m+1))/t1
73 IF (tnorm .EQ. 0) cycle
75 rhoc(m) = tnorm*( (r1c(m+1) + z1s(m+1))*t1
76 1 + (r1c(m-1) - z1s(m-1))*t2 )
77 rhos(m) = tnorm*( (r1s(m+1) - z1c(m+1))*t1
78 1 + (r1s(m-1) + z1c(m-1))*t2 )
82 END SUBROUTINE amplitud