Changed primitive variable gradient calculation from analytic to WLS or GG, enabled axisymmetric simulation for y > or < 0