[CIG-LONG] Initial temperature conditions

Walter Landry walter at geodynamics.org
Tue Feb 21 10:16:43 PST 2012


Robert Gray <graywacke at gmail.com> wrote:
> Hello everyone,
> I'm trying to set up an initial temperature profile in Gale that increases
> steadily from 273 K to 1623 K in the top 120 km of the box, while
> everything below that is 1623 K.  To do this I've used the following
> condition:
> 
>  "value": "(y>=0.48e6) ? -1.125e-2*y + 7023.0 : 1623.0"
> 
> It works great for the first time step, but after that temperatures jump
> several orders of magnitude to unrealistic values.  That said, the
> incremental increase :
> 
> "value": "(-1350.0/0.6e6)*y + 1623.0"
> 
> works fine but does not produce the initial temperature profile I need.

Can you send the input file to the mailing list?

Cheers,
Walter Landry


More information about the CIG-LONG mailing list