fixed rk for X point etc., fixed poisson, however rk slow and not nice programmed and dz_geom_t has to be further optimzed