00001 subroutine msec_store_f_vector(f_vector) 00002 use phys_constant, only : long 00003 use def_bh_parameter 00004 use def_quantities 00005 use def_iter_quantities 00006 implicit none 00007 real(long), pointer :: f_vector(:) 00008 integer :: i 00009 ! 00010 i=0 00011 i=i+1; f_vector(i) = (admmass - komarmass)/admmass 00012 i=i+1; f_vector(i) = qua_loc_spin/irredmass/irredmass - 0.0d0 00013 ! 00014 end subroutine msec_store_f_vector