Changes between Version 2 and Version 3 of AstroBearObjectTracers

Timestamp:

10/22/12 12:43:49 (12 years ago)

Author:

Erica Kaminski

Comment:

—

AstroBearObjectTracers

@@ -13 +13 @@
 This call passes in a non-zero integer to index the tracer in q, which is assigned in the routine AddTracer, by identifying the next empty slot in q, as well as the name of the tracer, here "ClumpTracer". This name is then recognized by a visualization program such as Visit as a new variable.
 
-Next, the tracer is 
+Next, the tracer is distributed throughout the clump by the line in clumps.f90:
 
+IF (Clump%iTracer .ne. 0) s(Clump%iTracer) = s(1)
+
+This line says, if there is an assigned slot in q for the tracer, than assign as that tracer's value, s(1), where s(1) is defined earlier in clumps.f90 as:
+
+s(1)=(perturb-q(1))*f !fade into clump
+
+This source term takes the value of the clump density (perturb) and subtracts from it the density of the background density (q(1)), where background here means the ambient object "behind" the clump object. The f is a fading factor that is used within a smoothing region between clump and ambient object. If one is fully within the clump (and outside of the smoothing region), then f=1. The value of the density is then assigned by rho=q(1)+s(1)=perturb. If, on the other hand, one is fully outside of the clump, then f=0, and rho=q(1)+s(1)=q(1). Within the smoothing region, some value between 0 and 1 is used to fade the density between these two objects. 
+
+
+
+
+