[CIG-LONG] Gale 2.0

[Patrice Rey]

G'day Walter

I am testing Gal 2.0 with a simple model of lithospheric extension. The model involves visco(non-Newtonian)-plastic rheology with internal radiogenic heating.

Boundary conditions involves wrapTop and staticBottom, and constant velocity is imposed of both vertical walls.

Mechanically the model evolves nicely but the temperature does not evolve as expected.  A negative temperature gradient develops in the viscinity of the Moho (cf. attached figure).  This suggests a problem with diffusivity which in my model is the same for all materials.

I have checked that Gale's time step (ca. 5e10 sec) is much smaller that dt_diff = courantFactor * (cellsize^2) / (maximum_diffusivity) = 2.5 10^11 sec.

I have tried to impose a smaller time step (ten times smaller)  via  <param_name="maxTimeStepSize">5e9sec</param>, but the problem persists.



0
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PastedGraphic-1.tiff
Type: image/tiff
Size: 228638 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20111110/71f67982/attachment-0001.tiff 
-------------- next part --------------




Cheers

-------------- next part --------------
A non-text attachment was scrubbed...
Name: AA-MCC-V-110b.xml
Type: application/xml
Size: 39293 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20111110/71f67982/attachment-0001.xml 
-------------- next part --------------


Patrice