416   static PetscScalar cl, cr, wl, zl, wr, zr, pst, durl, skpr1, skpr2;  in riem1mdt()  local
442   pst     = *pl / *pr;  in riem1mdt()
443   skpr1   = cr * (pst - 1.) * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst)));  in riem1mdt()
445   rarepr2 = gascr4 * 2. * cr * (1. - PetscPowScalar(pst, d__1));  in riem1mdt()
446   pst     = *pr / *pl;  in riem1mdt()
447   skpr2   = cl * (pst - 1.) * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)));  in riem1mdt()
449   rarepr1 = gascl4 * 2. * cl * (1. - PetscPowScalar(pst, d__1));  in riem1mdt()
499       pst     = *pstar / *pl;  in riem1mdt()
500 …ustarl  = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst)…  in riem1mdt()
501 …* 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. + (*gaml + 1.) * pst) / (*gaml * 3. - 1. +…  in riem1mdt()
522       pst    = *pstar / *pl;  in riem1mdt()
523 …ustarl = *uxl - (pst - 1.) * cl * PetscSqrtScalar(2. / (*gaml * (*gaml - 1. + (*gaml + 1.) * pst))…  in riem1mdt()
524 …* 2. * (*gaml - 1. + (*gaml + 1.) * pst)) * (*gaml - 1. + (*gaml + 1.) * pst) / (*gaml * 3. - 1. +…  in riem1mdt()
525       pst    = *pstar / *pr;  in riem1mdt()
526 …ustarr = *uxr + (pst - 1.) * cr * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst))…  in riem1mdt()
527 …* 2. * (*gamr - 1. + (*gamr + 1.) * pst)) * (*gamr - 1. + (*gamr + 1.) * pst) / (*gamr * 3. - 1. +…  in riem1mdt()
548       pst     = *pstar / *pr;  in riem1mdt()
549 …ustarr  = *uxr + (pst - 1.) * cr * PetscSqrtScalar(2. / (*gamr * (*gamr - 1. + (*gamr + 1.) * pst)…  in riem1mdt()
550 …* 2. * (*gamr - 1. + (*gamr + 1.) * pst)) * (*gamr - 1. + (*gamr + 1.) * pst) / (*gamr * 3. - 1. +…  in riem1mdt()
564     pst     = *pstar / *pl;  in riem1mdt()
565     *rstarl = ((*gaml + 1.) * pst + *gaml - 1.) / ((*gaml - 1.) * pst + *gaml + 1.) * *rl;  in riem1mdt()
568     pst     = *pstar / *pr;  in riem1mdt()
569     *rstarr = ((*gamr + 1.) * pst + *gamr - 1.) / ((*gamr - 1.) * pst + *gamr + 1.) * *rr;  in riem1mdt()