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