def c_gridrec(tomo, center, recon, theta, **kwargs):

    if len(tomo.shape) == 2:
        # no y-axis (only one slice)
        dy = 1