[CIG-LONG] Initial temperature conditions

[Robert Gray]

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.

Any help on the matter is greatly appreciated.
Cheers,
Rob
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20120221/ca2dd970/attachment.htm