The primary unknowns are the wave_scalar and its first time derivative wave_fscalar ( as TOCHNOG only solves first order in time equations, the first time derivative of also becomes primary unknown in order to turn this second order in time equation into a set of first order in time equations). Further notation: space coordinate, time and speed of sound.