Search found 36 matches

by joeatodd
27 Mar 2020, 17:05
Forum: Elmer/Ice
Topic: Zero velocities in a finer mesh resolution
Replies: 3
Views: 6934

Re: Zero velocities in a finer mesh resolution

Hi Martin, I had a quick look at your case. I think what's actually happening is that MUMPS is throwing an allocation error and giving up, resulting in no velocity computation. Basically, your mesh is too fine for whatever system you are running on (and indeed too large for my own desktop computer)....
by joeatodd
03 Feb 2020, 14:50
Forum: Elmer/Ice
Topic: Time step dependency for Stokes grounding line
Replies: 5
Views: 8431

Re: Time step dependency for Stokes grounding line

Hi all, I had previously done some investigations with a basal channel mesh. I have been rather busy recently and haven't had time to continue investigations into this, but Brandon provided me with the mesh they use in their investigations (subaerial calving at the right margin). I can provide this ...
by joeatodd
21 Dec 2015, 16:14
Forum: General
Topic: Coding style in Fortran code of Elmer
Replies: 1
Views: 9141

Re: Coding style in Fortran code of Elmer

Small piece of info regarding emacs, the default indentation level is 3, and this can be changed with: (setq f90-do-indent 2) (setq f90-if-indent 2) (setq f90-structure-indent 2) in the .emacs config file. Also, the following will take care of syntax highlighting and should cause keywords to appear ...
by joeatodd
15 Oct 2015, 14:12
Forum: Elmer/Ice
Topic: the temperature of the bottom of the bed of the ice
Replies: 1
Views: 6237

Re: the temperature of the bottom of the bed of the ice

Hey, For the temperature problem, did you set the following in the temperature solver section? Loop While Unconstrained Nodes = Logical True For the values of the parameters you mention, reasonable values are listed on this page: http://elmerice.elmerfem.org/wiki/doku.php?id=problems:rheology&s[...
by joeatodd
02 Jul 2015, 13:49
Forum: Elmer/Ice
Topic: Install Elmer/ice in Redhat
Replies: 4
Views: 8424

Re: Install Elmer/ice in Redhat

Hi Saisai, It'd be more useful to see the logs of your installation to see why the elmer/ice installation failed. When you do bash elmer_install.sh you could pipe the errors to a log file like this: bash elmer_install.sh 2>&1 | tee LOG and then either upload the whole LOG file somewhere, or have...
by joeatodd
25 Jun 2015, 11:17
Forum: Elmer/Ice
Topic: Install Elmer/ice in Redhat
Replies: 4
Views: 8424

Re: Install Elmer/ice in Redhat

Hi Saisai,

What errors did you encounter when installing elmer/ice?

Joe
by joeatodd
10 Jun 2015, 13:36
Forum: Elmer/Ice
Topic: TemperateIce stability issues
Replies: 2
Views: 6481

Re: TemperateIce stability issues

Hey Peter, Thanks for getting back to me. I tried modifying the temperature solver to set MeshVelocity = 0.0 internally when running in steady state but this didn't make a difference. I've tracked the problem down to the recomputation (I think) of the stabilization parameters within DisplaceMesh. I ...
by joeatodd
09 Jun 2015, 16:19
Forum: Elmer/Ice
Topic: TemperateIce stability issues
Replies: 2
Views: 6481

TemperateIce stability issues

Hi all, I'm encountering some strange convergence issues in the temperature solver in a coupled NS/Heat equation + mesh update spinup run. Although my simulation is transient, I'm attempting to solve the steady state temperature equation by specifying: Steady State Condition = Real 1.0 in the Temper...
by joeatodd
09 Jun 2015, 15:58
Forum: Elmer/Ice
Topic: Glacier width in FreeSurfaceSolver
Replies: 1
Views: 5494

Re: Glacier width in FreeSurfaceSolver

Hi Sharon, It's possible to do this as an accumulation term in the freesurface equation, but I think that it's more accurate to directly account for the addition of mass in the stokes equations. I wrote some code to do this in a flowline case by adding a "Pressure Load" following Peter's a...
by joeatodd
09 Apr 2015, 11:41
Forum: ElmerSolver
Topic: GlobalToLocal: did not converge
Replies: 5
Views: 3736

Re: GlobalToLocal: did not converge

Hey Peter,

In all the cases where this error occurs, it seems like the point is not within the element, so I suppose within the context of interpolation, it's probably not an important error.

Cheers,

Joe