Please help with updating parameter wave in IntegrateODE function. solving PDE through NUMOL

Hello,

I am using the function integrateODE as the integrator to solve a system of PDEs through the numerical method of lines. I need the integrator to update the parameter wave (pw) every time it calls my derivative function; all the elements in PW are a function of either yw (the wave to integrate) or dydx (its derivative, computed by the derivative function). Is there any way to do this? Another option I have tried is to make the derivative function calculate dydx as a function of Yw and the current dydx guess, but this does not seem to be possible either (Im not sure it is possible to make the function calculate dydx as a function of the current yw and dydx guess). Is there any example project where somebody has used NUMOL in igor with integrateODE?

Thank you for your help!