<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div></div><div><span>Right, the output, printed in vis.txt is this:</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span><br>E->EVI[0][1][50111] = 100 <----- before viscosity change</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>E->EVI[0][1][50111] = 10 <------ after viscosity change.</span></div><div
style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>Like I said, it seems to be working, with the exception that it does not seem to be influencing the flow field calculations. </span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent;
font-style: normal;"><span>Thank you for your time,</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br><span></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><span>Scott Tarlow<br></span></div><div style="display: block;" class="yahoo_quoted"> <br> <br> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> On Wednesday, August 7, 2013 1:29 PM, Shijie Zhong <shijie.zhong@Colorado.EDU> wrote:<br> </font> </div> <div class="y_msg_container"><br clear="none">You can always print out EVI to see whether they are what you want to be, and I realize you may have
tried this.<br clear="none"><br clear="none"><br clear="none"><br clear="none">Shijie Zhong, Professor<br clear="none">Department of Physics<br clear="none">University of Colorado at Boulder<br clear="none">Boulder, CO 80309<br clear="none">Tel: 303-735-5095; Fax: 303-492-7935<br clear="none">Web: <a shape="rect" href="http://anquetil.colorado.edu/szhong" target="_blank">http://anquetil.colorado.edu/szhong</a><br clear="none">________________________________________<div class="yqt4641765633" id="yqtfd00613"><br clear="none">From: <a shape="rect" ymailto="mailto:cig-mc-bounces@geodynamics.org" href="mailto:cig-mc-bounces@geodynamics.org">cig-mc-bounces@geodynamics.org</a> [<a shape="rect" ymailto="mailto:cig-mc-bounces@geodynamics.org" href="mailto:cig-mc-bounces@geodynamics.org">cig-mc-bounces@geodynamics.org</a>] On Behalf Of scott tarlow [<a shape="rect" ymailto="mailto:scott_tarlow@yahoo.com"
href="mailto:scott_tarlow@yahoo.com">scott_tarlow@yahoo.com</a>]<br clear="none">Sent: Wednesday, August 07, 2013 12:17 PM<br clear="none">To: <a shape="rect" ymailto="mailto:cig-mc@geodynamics.org" href="mailto:cig-mc@geodynamics.org">cig-mc@geodynamics.org</a><br clear="none">Subject: [CIG-MC] Which viscosity array to modify? (CitcomS)<br clear="none"><br clear="none">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.<br clear="none"><br clear="none">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<br clear="none"><br clear="none">The command that actually reduces the original viscosity in
viscosity_structures.c is:<br clear="none"><br clear="none">E->EVI[lev][m][(i-1)*vpts + j] = E->EVI[lev][m][(i-1)*vpts + j] * E->weak_viscosity_zone.visc_reduction;<br clear="none"><br clear="none">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.<br clear="none"><br clear="none">Thank you for your time,<br clear="none"><br clear="none">Scott Tarlow<br clear="none"><br clear="none"></div><br><br></div> </div> </div> </div></div></body></html>