Search found 34 matches

by joeatodd
21 Dec 2015, 16:14
Forum: General
Topic: Coding style in Fortran code of Elmer
Replies: 1
Views: 2367

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: 3005

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[]=glen
by joeatodd
02 Jul 2015, 13:49
Forum: Elmer/Ice
Topic: Install Elmer/ice in Redhat
Replies: 4
Views: 3677

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 a l...
by joeatodd
25 Jun 2015, 11:17
Forum: Elmer/Ice
Topic: Install Elmer/ice in Redhat
Replies: 4
Views: 3677

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: 2910

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: 2910

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: 2544

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 advice in t...
by joeatodd
09 Apr 2015, 11:41
Forum: ElmerSolver
Topic: GlobalToLocal: did not converge
Replies: 5
Views: 528

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
by joeatodd
08 Apr 2015, 12:23
Forum: ElmerSolver
Topic: GlobalToLocal: did not converge
Replies: 5
Views: 528

Re: GlobalToLocal: did not converge

Hey Peter, Thanks for the info. In my case, it's failing in 3D elements, and as far as I can tell, they don't appear to be degenerate. I don't know enough about the transformation between global and local coordinate systems, but is it possible that a given element could be sufficiently distorted fro...
by joeatodd
07 Apr 2015, 13:52
Forum: ElmerSolver
Topic: GlobalToLocal: did not converge
Replies: 5
Views: 528

GlobalToLocal: did not converge

Hi all, A quick search of the forums reveals that the warning "GlobalToLocal: did not converge" has come up in a few instances. I'm experiencing this error myself, and I'm wondering if anyone can help me understand the root cause of this? The backtrace for this error is roughly: GlobalToLocal PointI...