updating values of an array in a subroutine in fortran

Question

i have the following module:

module prestel
implicit none
contains
  subroutine gradient(rho,N,dr,gradiant_rho)
    implicit none
    real(kind=dp), allocatable