[CIG-LONG] Thermally driven counterflow

[Walter Landry]

Walter Landry <walter at geodynamics.org> wrote:
> Karen Paczkowski <karen.paczkowski at yale.edu> wrote:
>> Hi,
>> 
>> Thank you for all the help.  I am unfortunately still getting a very
>> unexpected flow field.  The only velocities I have set are:
>> 
>> Vx = 0 and Vy = -Vslab on the left side
>> Vy=0 across the top
>> 
>> My input file is producing a flow field that flows in the opposite
>> direction as Vslab (up and to the right) and has a magnitude 2-3
>> orders of magnitude larger than Vslab.  This was not occurring prior
>> to adding the thermal field, so I am assuming it has to do with that,
>> but I'm really not sure what is causing this exactly.  I have attached
>> a picture of the velocity field and a copy of the input file.  It is
>> written for Gale 1.5.0.  I would really appreciate any advice.
> 
> The problem is that there is no maximum temperature in the
> HydrostaticTerm.  I am on vacation until September, so I will not be
> able to get to it until then.

Sorry for the long delay, but I finally implemented this.  You can get
the code from the repository.  I am attaching a new input file and a
picture of the velocity and temperature.  Please let me know if you
have any more questions.

Cheers,
Walter Landry
walter at geodynamics.org
-------------- next part --------------
A non-text attachment was scrubbed...
Name: VerticalSubduction_ThermalField.xml
Type: text/xml
Size: 29357 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20101014/a6ce5b22/attachment-0001.bin 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: vertical_slab.png
Type: image/png
Size: 34182 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20101014/a6ce5b22/attachment-0001.png