Added average constraint. Also, free_surface_fault is troublesome
Work on free surface fault. Also added a displacement disk compression
problem. Doesn't work...