anchor.diff:          iwork = 0
           work = 0.0
   
  -        ipar => ipar_vals
  -        rpar => rpar_vals
  +        allocate(ipar(num_deriv_ipar),rpar(num_deriv_rpar))
   
           ipar(i_id) = NScool_id
           ipar(i_num_terminal_writes) = 0
  
positive.code:     integer, parameter :: i_id = 1
      integer, parameter :: i_num_terminal_writes = 2
      integer, parameter :: num_deriv_ipar = 2
      integer, parameter :: num_deriv_rpar = 1    ! not used
      
  contains
      subroutine do_integrate_crust(NScool_id,ierr)
  
negative.code:         integer :: lwork, liwork, lrd, lid
          integer :: imas,mumas,mlmas
          real(dp), dimension(1) :: rtol, atol
          real(dp), pointer, dimension(:) :: rpar
          integer, pointer, dimension(:) :: ipar
          real(dp), pointer, dimension(:) :: rpar_decsol, work
  