by mewo2 » 01 May 2015, 14:23
Initially, using the actual ice shelf geometry, I thought that the GCR was failing to converge. However, I've been trying to reduce things to a simpler test case, and when I use a simple square geometry, I get (fairly slow) convergence of the GCR, but to a nonsensical result.
I haven't been using any grounding line dynamics  the domain is fully floating, with velocities specified along the upstream and lateral boundaries, and pressure boundary conditions (using the Buoyancy USF) on the downstream and lower boundaries.
I've attached a test case you can try for yourself. The "normal" Stokes solution is generated by seq.sif and the ParStokes by par.sif.
 Attachments

 parstokesshelf.tgz
 Test case .sif and .geo
 (2.47 KiB) Downloaded 428 times