mirror has sign consistent with being on RHS of GKE; it is the factor multiplying dg/dvpa in the mirror term mirror_sign set to +/- 1 depending on the sign of the mirror term. NB: mirror_sign = -1 corresponds to positive advection velocity set up the tridiagonal matrix that must be inverted for the implicit treatment of the mirror operator