comSol["indFinal"] = int(np.round((xCenter - xllc) / csz) + 1)