00001 subroutine compute_fnc_inversion_index(pot1,pot2,index) 00002 use phys_constant, only : long 00003 use grid_parameter, only : nrg, ntg, npg 00004 implicit none 00005 real(long), pointer :: pot1(:,:,:), pot2(:,:,:) 00006 real(long) :: index 00007 ! 00008 pot2(0:nrg,0:ntg,0:npg) = pot1(0:nrg,0:ntg,0:npg)**index 00009 end subroutine compute_fnc_inversion_index