(^ for the simulation in this figure I had changed the pRef value to = 0 for debugging)