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

[scott tarlow]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Citcom.c
Type: text/x-csrc
Size: 8997 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0003.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Citcom_init.c
Type: text/x-csrc
Size: 3295 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0004.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: global_defs.h
Type: text/x-chdr
Size: 22008 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0001.h>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Viscosity_structures.c
Type: text/x-csrc
Size: 49114 bytes
Desc: not available
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0005.c>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: vis.txt
URL: <http://geodynamics.org/pipermail/cig-mc/attachments/20130807/b27da68d/attachment-0001.txt>