[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