Search found 64 matches

by rgladstone
07 Aug 2014, 08:31
Forum: Elmer/Ice
Topic: non-reproducible machine-specific convergence failure
Replies: 3
Views: 8336

Re: non-reproducible machine-specific convergence failure

Thanks Thomas. I am guessing either something to do with MUMPS or MPI libraries, but I'll check after installing the same Elmer version first as that is the easiest thing to do. I guess I'll try installing MUMPS next rather than using the default MUMPS installation we have here. Which libraries does...
by rgladstone
07 Aug 2014, 04:48
Forum: Elmer/Ice
Topic: non-reproducible machine-specific convergence failure
Replies: 3
Views: 8336

non-reproducible machine-specific convergence failure

Hi all, not sure whether to put this in the Elmer/Ice section or in the main Elmer section as it is probably not ice-specific, though I am not yet certain! I am running flowline marine ice sheet simulations. I am running on Taito (CSC machine in Finland) and on Raijin (NCI machine in Australia). Som...
by rgladstone
24 Jun 2014, 12:43
Forum: Elmer/Ice
Topic: Trouble with lateral friction
Replies: 1
Views: 5545

Re: Trouble with lateral friction

I've now tried this with the Weertman and also the Coulomb/Cavity sliding laws and get the same problem: non-convergence in the nonlinear Stokes iterations when a significant lateral drag is applied.

I've varied the thickness of the initial slab, still the same problem.
by rgladstone
24 Jun 2014, 02:34
Forum: Elmer/Ice
Topic: Trouble with lateral friction
Replies: 1
Views: 5545

Trouble with lateral friction

Hi all, I am currently trying to set up a flowline marine ice sheet simulation with parameterised lateral drag. The lateral drag code is described here: http://elmerice.elmerfem.org/wiki/doku.php?id=userfunctions:lateralfriction The basal drag law is a function of overburden pressure. I am trying to...
by rgladstone
16 Jun 2014, 06:53
Forum: Bug reports
Topic: ElmerGrid partition option
Replies: 0
Views: 3058

ElmerGrid partition option

Hi, when using ElmerGrid with the -partition option for less than 3 dimensions (my case is 2 dimensions) a dummy argument is needed in the place of the third dimension in order for the partition strategy argument to be recognised. I need to set this argument to 2 otherwise the structured mesh mapper...
by rgladstone
05 Jun 2014, 11:55
Forum: Bug reports
Topic: problem when using "Calculate Velocity = Logical True"
Replies: 3
Views: 6623

Re: problem when using "Calculate Velocity = Logical True"

Hi Peter and co, I wonder if anyone has had a chance to look into this? I have a more-or-less identical problem, though the setup posted above is a smaller one, and so is perhaps more convenient for debugging purposes. I've run the above setup with a recent Elmer checkout, compiled with debug flags ...
by rgladstone
26 Feb 2014, 06:45
Forum: Installation & compilation
Topic: ParallelUtils error
Replies: 3
Views: 7028

Re: ParallelUtils error

Juha asked about the validity of this code. There is a reply here: http://gcc.gnu.org/ml/fortran/2012-05/msg00051.html The bottom line seems to be that the intel compilers seem to be over-cautious, but perhaps in a sensible way. The intel compilers throw an error if a pointer is assigned to a target...
by rgladstone
05 Dec 2013, 02:13
Forum: Elmer/Ice
Topic: Grid2DInterpolator updated
Replies: 7
Views: 9904

Re: Grid2DInterpolator updated

Thanks Olivier. I'll check in the bug fix. For interpolation in presence of noData, I calculate weights for normal bilinear interpolation, set the weights to zero for noData points, then scale the remaining weights so that they total one. You are right, this isn't obvious/trivial and can lead to pro...
by rgladstone
30 Nov 2013, 00:50
Forum: Elmer/Ice
Topic: Troubles with grounding line simulations
Replies: 4
Views: 7019

Re: Troubles with grounding line simulations

Thanks so much Olivier, it would have taken me a looooong time to figure that out on my own.
by rgladstone
28 Nov 2013, 05:11
Forum: Elmer/Ice
Topic: Grid2DInterpolator updated
Replies: 7
Views: 9904

Re: Grid2DInterpolator updated

The default in the 2d grid interpolator in the event of one to three neighbouring points having the noData value is to take the mean of valid points. I have written a few lines to instead carry out linear interpolation. Shall I set this as the default (which makes more sense to me) or as an option (...