00001 subroutine update_surface(potrs,rsnew,convf) 00002 use phys_constant, only : long 00003 use grid_parameter, only : ntg, npg 00004 implicit none 00005 real(long), pointer :: potrs(:,:) 00006 real(long), pointer :: rsnew(:,:) 00007 real(long), intent(in) :: convf 00008 rsnew(0:ntg,0:npg) = convf *potrs(0:ntg,0:npg) & 00009 & + (1.d0-convf)*rsnew(0:ntg,0:npg) 00010 end subroutine update_surface