<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:΢ź
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Hi all,<br><br>Thank you for your reply, which is very inspirational and helpful. <br><br>I now realized that, in CitcomS, viscosity structure of a slab can be prescribed by a given temperature perterbation, in which viscosity is set to be strongly temperature-dependent. As an example, in the case of the Mid-Ocean Ridge, one can use half-space 
cooling temperture profile to implement its unique viscosity structure 
(see attachment, from Liu's <a href="http://www.sciencedirect.com/science/article/pii/S0012821X11005474" target="_blank">EPSL paper</a>).<br><br>But it's not flexble enough. I think being able to <i>assign arbitrarily complex geometries to temperature, composition, and viscosity material</i> with Netcdf/GMT grd files as mentioned by Dr. Becker is an attractive way, although I have no idea about how so far. <br><br>Perhaps I will need more time digging into the code like viscosity_structures.c as well the solution of Netcdf/GMT files input. <br><br>Possible numerical diffusion when it comes to convection of temperature and composition field is also what I concern.<br><br>Thank you for your time!<br><br>Best,<br>Qingwen<br><br>                                         </div></body>
</html>