[CIG-LONG] Initial temperature conditions
Robert Gray
graywacke at gmail.com
Tue Feb 21 09:14:35 PST 2012
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
More information about the CIG-LONG
mailing list