[CIG-MC] Which viscosity array to modify? (CitcomS)

[Shijie Zhong]

You can always print out EVI to see whether they are what you want to be, and I realize you may have tried this.



Shijie Zhong, Professor
Department of Physics
University of Colorado at Boulder
Boulder, CO 80309
Tel: 303-735-5095; Fax: 303-492-7935
Web: http://anquetil.colorado.edu/szhong
________________________________________
From: cig-mc-bounces at geodynamics.org [cig-mc-bounces at geodynamics.org] On Behalf Of scott tarlow [scott_tarlow at yahoo.com]
Sent: Wednesday, August 07, 2013 12:17 PM
To: cig-mc at geodynamics.org
Subject: [CIG-MC] Which viscosity array to modify? (CitcomS)

Hello, I am currently  designing a static geometry weak viscosity zone, but having trouble figuring out which viscosity array is the correct one to modify, as my current solution does not change the flow field at all.  Attached are the files I modified to achieve this weak viscosity zone.

The volume of the weak viscosity zone is defined in Citcom_init.c, and its implementation  is a static void apply_weak_viscosity_zone in viscosity_structures.c

The command that actually reduces the original viscosity in viscosity_structures.c is:

E->EVI[lev][m][(i-1)*vpts + j] = E->EVI[lev][m][(i-1)*vpts + j] * E->weak_viscosity_zone.visc_reduction;

and its output, in vis.txt (by using the command, citcoms ######.cfg > vis.txt), seems to be changing the viscosity as intended. But when I look at the resulting flow field, it seems as if the weak viscosity zone does not exist. I'm fairly certain I am finding the correct index's for the nodes I want to change, as the output seems reasonable. I was wondering if I needed to modify another array before or after EVI to have an effect on the viscosity structure.

Thank you for your time,

Scott Tarlow