right state
rhRR = uR[i,1]
uuRR = uR[i,2]/rhRR
eeRR = uR[i,3]/rhRR
ppRR = gm*(eeRR*rhRR - 0.5*rhRR*(uuRR*uuRR))
aaRR = sqrt(abs(gamma*ppRR/rhRR))