<div dir="ltr">Just to clarify, the tracers behave like Eric said, but they become "active" whenever an actual chemical buoyancy number is assigned to them, as per the common usage of these terms. Moreover, they can be assigned a different viscosity. <div>
<br></div><div>If you want a static lateral viscosity variation field, the way to go is indeed to modify the viscosity routines so that the pre-factor is based on some geometric formulae. Alternative, the ggrd extension can read in GMT/Netcdf grid files of prefactor fields, though this is not well documented. </div>
<div><br></div><div>Cheers</div><div><br></div><div style>T</div><div style><br></div></div><div class="gmail_extra"><br clear="all"><div>Thorsten W Becker - <a href="http://geodynamics.usc.edu/~becker" target="_blank">http://geodynamics.usc.edu/~becker</a></div>
<br><br><div class="gmail_quote">On Tue, Jul 2, 2013 at 7:15 PM, Scott Tarlow <span dir="ltr"><<a href="mailto:scott_tarlow@yahoo.com" target="_blank">scott_tarlow@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word"><div>Hey Eric,</div><div><br></div><div>Thanks for the quick response. I think you are right, in general it would be better to hardcode in a low viscosity region. If I were to create this region, I do have a few questions:</div>
<div><br></div><div>In the low_viscosity_channel_factor function it finds the indices for the minimum and maximum radius:</div><div><br></div><div><div style="margin:0px;font-size:11px;font-family:Menlo"><span style="color:#bb2ca2">for</span>(m=<span style="color:#272ad8">1</span>; m<=E->sphere.caps_per_proc; m++) {</div>
<div style="margin:0px;font-size:11px;font-family:Menlo;color:rgb(0,132,0)"><span style="color:#000000"> </span>/* find index of radius corresponding to lv_min_radius */</div><div style="margin:0px;font-size:11px;font-family:Menlo">
<span style="color:#bb2ca2">for</span>(e=<span style="color:#272ad8">1</span>; e<=E->lmesh.elz; e++) {</div><div style="margin:0px;font-size:11px;font-family:Menlo"> rad_mean = <span style="color:#272ad8">0.5</span> * (E->sx[m][<span style="color:#272ad8">3</span>][E->ien[m][e].node[<span style="color:#272ad8">1</span>]] +</div>
<div style="margin:0px;font-size:11px;font-family:Menlo"> E->sx[m][<span style="color:#272ad8">3</span>][E->ien[m][e].node[<span style="color:#272ad8">8</span>]]);</div><div style="margin:0px;font-size:11px;font-family:Menlo">
<span style="color:#bb2ca2">if</span>(rad_mean >= E->viscosity.lv_min_radius) <span style="color:#bb2ca2">break</span>;</div><div style="margin:0px;font-size:11px;font-family:Menlo"> }</div><div style="margin:0px;font-size:11px;font-family:Menlo">
nz_min[m] = e;</div></div><div><br></div><div>To modify this for the x and y directions can i just change <span style="font-family:Menlo;font-size:11px">lmesh.elz</span> to <span style="font-family:Menlo;font-size:11px">lmesh.elx along with changing the looping parameter and if statement parameters? Does sx[m][3] refer to the radial direction, meaning sx[m][2] would be y and sx[m][1] would be x? </span></div>
<div><span style="font-family:Menlo;font-size:11px"><br></span></div><div><font face="Menlo"><span style="font-size:11px">How I think the best way to implement it is by creating a new static void, finding the start and ending indices for the x,y, and z direction, and then looping over that volume and applying the viscosity reduction factor. Doing this, where would the best place be to call the function in the Viscosity_structures.c file? </span></font></div>
<div><font face="Menlo"><span style="font-size:11px"><br></span></font></div><div><font face="Menlo"><span style="font-size:11px">Thank you for your time,</span></font></div><div><font face="Menlo"><span style="font-size:11px"><br>
</span></font></div><div><font face="Menlo"><span style="font-size:11px">Scott Tarlow</span></font></div><div><br></div><br><div><div>On Jul 2, 2013, at 6:15 PM, Eric Heien <<a href="mailto:emheien@ucdavis.edu" target="_blank">emheien@ucdavis.edu</a>> wrote:</div>
<br><blockquote type="cite">Hi Scott,<br><br>Keep in mind the tracers in themselves do not impart any properties on the simulation (they are passive). Even if you stop them from moving, you have to couple their existence (or flavor) to whatever viscosity property you want. Since this change in viscosity is just placed on certain nodes, it might be simpler to define a specified viscosity on whatever nodes you want to be a weak zone. This way there isn't any confusion about how the tracers are interacting with the simulation and it's clear what the size and properties of the weak zone are.<br>
<br>I don't know the details of what you're doing so your current solution might work fine, I just want to make sure there's no confusion about it. In any case let us (or the mailing list) know if you have any more questions. Good luck,<br>
<br>-Eric<br><br>On Jul 2, 2013, at 12:22 PM, scott tarlow wrote:<br><br><blockquote type="cite">I'd be happy to.<br><br>One of the necessary requirements for the project I'm working on is that the overriding and subducting plate are decoupled. I think there are a few ways of doing this: One is to write some sort of extension to the Viscosity_structures.C file and change a specific volume of nodes to a much reduced viscosity, ie "weak" nodes. The other solution is using the tracers in conjunction with a Low Viscosity Channel. This worked great, and seems much easier to use with the exception that because the tracers move, the Low Viscosity Channel also moves. This causes the slab to move horizontally, I think. If I am able to stop the tracers from moving, I am essentially reproducing "weak" nodes without significant coding changes to the file. <br>
<br>Thank you for your time,<br><br>Scott Tarlow<br><br><br>From: Louise Kellogg <<a href="mailto:kellogg@ucdavis.edu" target="_blank">kellogg@ucdavis.edu</a>><br>To: scott tarlow <<a href="mailto:scott_tarlow@yahoo.com" target="_blank">scott_tarlow@yahoo.com</a>> <br>
Cc: Eric Heien <<a href="mailto:emheien@ucdavis.edu" target="_blank">emheien@ucdavis.edu</a>> <br>Sent: Tuesday, July 2, 2013 12:24 PM<br>Subject: Re: [CIG-MC] Removing Tracer Advection<br><br>Hi Scott,<br><br>I'm curious what scientific problem you're working on that would need the tracers to stay in one place. The tracers are passive and massless, so locking them in place will not influence the flow - but it also defeats the usual purpose of the tracers. <br>
<br>Can you tell us a bit more about what you're working on? <br><br>Best wishes,<br><br>Louise<br><br><br>On Jul 1, 2013, at 2:01 PM, scott tarlow <<a href="mailto:scott_tarlow@yahoo.com" target="_blank">scott_tarlow@yahoo.com</a>> wrote:<br>
<br><blockquote type="cite">Hello,<br><br>I am looking to modify the of the files : Tracer_setup.c and Regional_tracer_advection.c (in CitcomS) in order to force the tracers to stay at their original positions instead of moving with the flow field. <br>
<br>My original approach was to find in the code where where the nodal velocities are pulled and then assigned to a tracer and set them to 0 instead of Vx,Vy,Vz respectably. This did not work, so my next solution was deleting Regional_tracer_advection.c from the lib file. This hasn't worked either. <br>
<br>If Regional_tracer_advection.c is on my desktop instead of the lib folder, CitcomS wont read it correct? <br><br>I'm not really sure where to go from here. Perhaps someone has a suggestion ? <br><br><br>Thank you for your time,<br>
<br>Scott Tarlow<br>_______________________________________________<br>CIG-MC mailing list<br><a href="mailto:CIG-MC@geodynamics.org" target="_blank">CIG-MC@geodynamics.org</a><br><a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc</a><br>
</blockquote><br>*********************************************<br>Louise Kellogg<br>Professor, Department of Geology<br>One Shields Avenue, University of California, Davis, CA 95616<br>Phone: <a href="tel:%28530%29752-3690" value="+15307523690" target="_blank">(530)752-3690</a><br>
<a href="http://mygeologypage.ucdavis.edu/kellogg/" target="_blank">http://mygeologypage.ucdavis.edu/kellogg/</a><br><a href="http://www.keckcaves.org" target="_blank">http://www.keckcaves.org</a><br>*********************************************<br>
<br><br><br><br><br></blockquote><br></blockquote></div><br></div><br>_______________________________________________<br>
CIG-MC mailing list<br>
<a href="mailto:CIG-MC@geodynamics.org">CIG-MC@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc</a><br></blockquote></div><br></div>