00001 subroutine allocate_iter_quantities(niq)
00002   use phys_constant, only : long
00003   use def_iter_quantities
00004   use def_bh_parameter, only : ome_bh, spin_bh
00005   use def_quantities
00006   use make_array_1d
00007   implicit none
00008   integer, intent(in) :: niq
00009   integer :: i
00010 
00011   call alloc_array1d(msec_x_oold,1,niq)
00012   call alloc_array1d(msec_x_old ,1,niq)
00013   call alloc_array1d(msec_f_oold,1,niq)
00014   call alloc_array1d(msec_f_old ,1,niq)
00015 
00016   i=0
00017   i=i+1; msec_x_oold(i) = ome_bh
00018   i=i+1; msec_x_oold(i) = spin_bh
00019 
00020   i=0
00021   i=i+1; msec_f_oold(i) = 0.0d0
00022   i=i+1; msec_f_oold(i) = 0.0d0
00023 
00024 end subroutine allocate_iter_quantities