[CIG-LONG] Behavior without initial temperature conditions?

[Brian Wilson]

Hi all,

Ok it appears that GALE does set the temperature to 0 if you don't
specify an initial temperature. For some reason it wasn't doing this
the other day but now when I run it without initial conditions it
gives me an error in the nonNewtonian rheology because T cannot be 0.

However, when I tried removing the Stokes equations the code simply
crashed with a memory error, and when I ran it with vx=0 the nonlinear
iterations did not converge for the first time step.

Any ideas?

Thanks,
Brian

On 11/12/11, Patrice Rey <patrice.rey at me.com> wrote:
> Hi Brian
>
>> Is there a capability for GALE to initially compute a steady state
>> temperature solution?
>
> Here is what I suggest.
>
> 1/ Run your model with the boundary conditions driving deformation turned
> off (e.g. vx=0).  You may have to turn off the uzawa solver as well.  In
> fact, turning off the uzawa (Stoke) solver may be the only thing necessary.
>
> 2/ When you are happy that the geotherm is close to steady state, then
> restart your model using
> --restartTimestep="your_timestep_for_steady_state_geotherm" with driving
> boundary conditions and uzawa solver turned on.
>
> Hope this help.
>
> Patrice
>
>
> On 12/11/2011, at 9:08 AM, Brian Wilson wrote:
>
>> Hi all,
>>
>> Quick question: what happens if you do not specify an initial
>> temperature condition, but just boundary conditions? Is there a
>> capability for GALE to initially compute a steady state temperature
>> solution?
>>
>> Thanks,
>> Brian
>> _______________________________________________
>> CIG-LONG mailing list
>> CIG-LONG at geodynamics.org
>> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>
> Associate Professor Patrice F. Rey
> EarthByte Research Group
> School of Geosciences
> The University of Sydney
> +61 2 9351 2067
> patrice.rey at mac.com
> patrice.rey at sydney.edu.au
>
>
>
>
>
>
>