[CIG-LONG] Gale help

[Walter Landry]

"谢祖军" <zjxie at asch.whigg.ac.cn> wrote:
> Dear professors,
> 
> I am Zujun Xie, a postgraduate student of institute of Geodesy and Geophysics,
> Chinese Academy of Sciences. I am interesting in Gale, a code for the
> long-term tectonics which could be available from the Gale website. I have
> tried to use it to simulate the upper mantle convection and crustal dynamics
> in Tibetan Plateau and its surroundings. I want to utilize the seismic
> tomography result to obtain the distribution of density which will be used as
> material input to my calculations. However when I want to do this, I find it
> is quite hard to give the material properties on the elements, because in the
> examples of Gale, most of the material properties are given to a body. For my
> work, I shall give the material properties to each of the element, because in
> tomography models, each element has its own density. I am wondering how to add
> the material properties on the elements.

The simplest way is to make a new shape for each element.  You should
write a small program that will take your tomography output and, for
each element, write out something like

    <struct name="ShapeN">
      <param name="Type">Box</param>
      <param name="startX">minX</param>
      <param name="endX">maxX</param>
      <param name="startY">minY</param>
      <param name="endY">10000</param>
      <param name="startZ">minZ</param>
      <param name="endZ">maxZ</param>
    </struct>

    <struct name="ViscosityN">
      <param name="Type">FrankKamenetskii</param>
      <param name="TemperatureField">TemperatureField</param>
      <param name="eta0">eta</param>
      <param name="theta">theta</param>
    </struct>

    <struct name="MaterialN">
      <param name="Type">RheologyMaterial</param>
      <param name="Shape">ShapeN</param>
      <param name="density">density</param>
      <param name="alpha">alpha</param>
      <param name="diffusivity">density</param>
      <list name="heatingElements">
        <struct>
          <param name="Q">Q</param>
          <param name="lambda">lambda</param>
        </struct>
      </list>
      <list name="Rheology">
        <param>ViscosityN</param>
        <param>storeViscosity</param>
        <param>storeStress</param>
      </list>
    </struct>

You use this to generate a huge input file.  I had another user do
this, so it should work.

Cheers,
Walter Landry
walter at geodynamics.org