FIXME : Add documentation
subroutine allocate_arrays_3d
implicit none
integer :: unit, ntor
! character(200) :: line
call get_unused_unit(unit)
open (unit=unit, file=trim(source), status="old")
read (unit=unit, fmt=*) ntgridin, nthetain, ntor
close (unit=unit)
drhodpsi = 1.0 ! assumes our radial variable matches VVBAL
nthetain = nthetain * iperiod
allocate (thetain(nthetain+1),bmagin(nthetain+1),bmagsm(nthetain+1))
allocate (thetaout(nthetaout),bmagout(nthetaout),alambdaout(nlambdaout))
allocate (thetagrid(-ntgridin:ntgridin))
allocate (gbdrift(-ntgridin:ntgridin))
allocate (gradpar(-ntgridin:ntgridin))
allocate (grho(-ntgridin:ntgridin))
allocate (cvdrift(-ntgridin:ntgridin))
allocate (gds2(-ntgridin:ntgridin))
allocate (bmag(-ntgridin:ntgridin))
allocate (gds21(-ntgridin:ntgridin))
allocate (gds22(-ntgridin:ntgridin))
allocate (cvdrift0(-ntgridin:ntgridin))
allocate (gbdrift0(-ntgridin:ntgridin))
allocate (Rplot(-ntgridin:ntgridin))
allocate (Rprime(-ntgridin:ntgridin))
allocate (Zplot(-ntgridin:ntgridin))
allocate (Zprime(-ntgridin:ntgridin))
allocate (aplot(-ntgridin:ntgridin))
allocate (aprime(-ntgridin:ntgridin))
cvdrift0 = 0. ! assumes theta0 = 0.
gbdrift0 = 0. ! assumes theta0 = 0.
gds21 = 0. ! assumes theta0 = 0.
gds22 = 0. ! assumes theta0 = 0.
grho = 1. ! assumes linear calculation
Rplot = 1. ! pure fiction
Rprime = 1. ! pure fiction
Zplot = 1. ! pure fiction
Zprime = 1. ! pure fiction
aplot = 1. ! pure fiction
aprime = 1. ! pure fiction
end subroutine allocate_arrays_3d