I run the model under steady state condition by separating the model in to 15 layers and 16 slices, thickness is about 300m. I assigned initial groundwater level, K, n, Ss and boundary condition. Now problem settings, Saturated condition, top slice as free & movable and below layers are phreatic and bottom most layer is under fixed condition. After simulation the head values for every layer is different (not uniform). I dont know how to fix this problem and where i made mistake. I would happy if some one solve this problem.