00001 subroutine compute_fnc_division(pot1,pot2,pot3)
00002   use phys_constant,  only : long
00003   use grid_parameter, only : nrg, ntg, npg
00004   implicit none
00005   real(long), pointer    :: pot1(:,:,:), pot2(:,:,:), pot3(:,:,:)
00006 
00007   pot3(0:nrg,0:ntg,0:npg) = pot1(0:nrg,0:ntg,0:npg)/pot2(0:nrg,0:ntg,0:npg)
00008 end subroutine compute_fnc_division