[CIG-MC] How to impose arbitrary viscosity and density structure to CitcomS?

[Eh Tan]

Hi Choon Ki,

To input arbitrary buoyancy, you can use tic_method=-1, as described in
Cookbook 5 in the manual. However, it is currently impossible to input
arbitrary viscosity. There are two possible approaches: 1) using
temperature and composition dependent viscosity to approximate your
viscosity structure; 2) modify the code to impose your viscosity
structure. If you decide to modify the code, I can help to point out
where to modify and explain the code.

Cheers,

Eh Tan



Choon-Ki Lee wrote:
> Hi all
>
> I am a new CitcomS user who is interested in only Stokes flow calculation.
> Looking over the manual, I didn't precisely understand the input
> parameters and files.
> I would like to calculate the Stokes flow driven by a density
> structure, such as a tomography model, and an arbitrary
> laterally-heterogeneous viscosity structure (not temperature dependent
> viscosity).
> How can I impose the density and viscosity structure as input
> parameters to CitcomS?
>
> If someone give me a hint, I will read the manual more carefully.
> Thanks.
>
> Choon Ki Lee
> GEST/UMBC 
> ------------------------------------------------------------------------
>
> _______________________________________________
> CIG-MC mailing list
> CIG-MC at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc
>   

-- 
Eh Tan
Staff Scientist
Computational Infrastructure for Geodynamics
California Institute of Technology, 158-79
Pasadena, CA 91125
(626) 395-1693
http://www.geodynamics.org