1 SUBROUTINE tolicu (torcur)
5 USE parallel_include_module
10 REAL(dp),
INTENT(IN) :: torcur
14 INTEGER :: i, kper, kv
15 REAL(dp) :: current(1), ttolion, ttolioff
16 REAL(dp),
DIMENSION(:,:),
ALLOCATABLE :: xpts
23 ALLOCATE (xpts(3,nvp), stat=i)
24 IF (i .ne. 0) stop
' allocation error in tolicu'
32 xpts(1,i) = raxis_nestor(kv)*(cosper(kper)*cosuv(kv)
33 1 - sinper(kper)*sinuv(kv))
34 xpts(2,i) = raxis_nestor(kv)*(sinper(kper)*cosuv(kv)
35 1 + cosper(kper)*sinuv(kv))
36 xpts(3,i) = zaxis_nestor(kv)
43 CALL initialize_biotsavart (current, xpt=xpts)
45 CALL second0(ttolioff)
46 s_tolicu_time = s_tolicu_time + (ttolioff - ttolion)