00001 subroutine update_matter(potf,mtfield,convf) 00002 use phys_constant, only : long 00003 use grid_parameter, only : nrf, ntf, npf 00004 implicit none 00005 real(long), pointer :: potf(:,:,:) 00006 real(long), pointer :: mtfield(:,:,:) 00007 real(long), intent(in) :: convf 00008 mtfield(0:nrf,0:ntf,0:npf) = convf * potf(0:nrf,0:ntf,0:npf) & 00009 & + (1.d0-convf)*mtfield(0:nrf,0:ntf,0:npf) 00010 end subroutine update_matter