anchor.diff:        dRhou(:,i)=0d0
         dRhoE(i)=0d0
       end forall
  -    call ode%init(nEq,rhs)
  +    call ode%init(nEq,rhs,pSet=pset,pSol=psol)
       call ode%setTol(RTOL,ATOL)
       call ode%setMaxSteps(-1)
       call ode%setIV(t,y)
  
positive.code:     foo=dat
    end function
    
  end module
  
  !> Euler solver
  program foeuler
    use modEuler
    character(20)::tmpStr
    integer::ier
    
    call init()
    write(tmpStr,*)iOut
  
negative.code:   end do
    call clear()
  end program
  