[CIG-LONG] CIG-LONG Digest, Vol 60, Issue 4

[Brian Wilson]

Sure, here it is. I've commented out the Stokes equations to turn them
off - I'm assuming that's the way to do it?

Brian

On 11/16/11, cig-long-request at geodynamics.org
<cig-long-request at geodynamics.org> wrote:
> Send CIG-LONG mailing list submissions to
> 	cig-long at geodynamics.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> 	http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
> or, via email, send a message with subject or body 'help' to
> 	cig-long-request at geodynamics.org
>
> You can reach the person managing the list at
> 	cig-long-owner at geodynamics.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of CIG-LONG digest..."
>
>
> Today's Topics:
>
>    1. Re: Behavior without initial temperature conditions?
>       (Brian Wilson)
>    2. Re: Behavior without initial temperature conditions?
>       (Walter Landry)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 15 Nov 2011 17:40:04 -0700
> From: Brian Wilson <bjwilson83 at gmail.com>
> Subject: Re: [CIG-LONG] Behavior without initial temperature
> 	conditions?
> To: cig-long <cig-long at geodynamics.org>
> Message-ID:
> 	<CAPRTR0FxGe9GxxSgXJOc6pqAJcixbeK-a=4e=_bBwJMt6ES4QQ at mail.gmail.com>
> Content-Type: text/plain; charset=ISO-8859-1
>
> 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
>>
>>
>>
>>
>>
>>
>>
>
>
> ------------------------------
>
> Message: 2
> Date: Tue, 15 Nov 2011 19:15:07 -0800 (PST)
> From: Walter Landry <walter at geodynamics.org>
> Subject: Re: [CIG-LONG] Behavior without initial temperature
> 	conditions?
> To: bjwilson83 at gmail.com
> Cc: cig-long at geodynamics.org
> Message-ID:
> 	<20111115.191507.1675203964983243693.walter at geodynamics.org>
> Content-Type: Text/Plain; charset=us-ascii
>
> Brian Wilson <bjwilson83 at gmail.com> wrote:
>> 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.
>
> Like I said, you should probably not depend on that.
>
>> However, when I tried removing the Stokes equations the code simply
>> crashed with a memory error,
>
> That should not happen.  You should be able to run the code without
> Stokes.  Can you send your input file?
>
> Cheers,
> Walter Landry
>
>
> ------------------------------
>
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>
>
> End of CIG-LONG Digest, Vol 60, Issue 4
> ***************************************
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: antarctica_scaled_3layer.xml
Type: text/xml
Size: 34400 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20111116/155f505a/attachment-0001.bin