crustal_deformation/crustal_model_2D.prm fails

[tjhei]

cd cookbooks; make -f check.mk BUILD=$BUILDDIR crustal_deformation/ fails with

    The iterative Stokes solver did not converge. It reported the                                                   
    following error:  
                  Additional information:                                                                                         
    Iterative method reported convergence failure in step 5000. The                                                 
    residual in the last step was 1.69553e+15.                                                                      
                                                                                                                    
    The required residual for convergence is: 50735275.339261.
    See output-crustal_model_2D/solver_history.txt for convergence
    history.
                                                                                

Originally posted by @tjhei in #4069 (comment)

[tjhei]

@naliboff can you take a look please? You can just run crustal_model_2D.prm to see the problem.