00001 subroutine update_grfield(pot,grfield,convf)
00002   use phys_constant,  only : long
00003   implicit none
00004   real(long), pointer    :: pot(:,:,:)
00005   real(long), pointer    :: grfield(:,:,:)
00006   real(long), intent(in) :: convf
00007 
00008   grfield = convf*pot + (1.0d0-convf)*grfield
00009 
00010 end subroutine update_grfield