astrobear - Blog
https://bluehound2.circ.rochester.edu/astrobear/blog
About blog postsen-USTrac 1.4.1Solution to streamlines issues.madamsWed, 03 Dec 2014 18:43:25 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/madams120320142
https://bluehound2.circ.rochester.edu/astrobear/blog/madams120320142<p>
Turns out we are calling the wrong position in the array. The expressions should be:
</p>
<p>
By_downx = array_decompose(projections,0)
</p>
<p>
Bz_downx = array_decompose(projections,1)
</p>
<p>
Bz_downy = array_decompose(projections,0)
</p>
<p>
Bx_downy = array_decompose(projections,1)
</p>
<p>
Byz_downx = {<By_downx>, <Bz_downx>}
</p>
<p>
Bzx_downy = {<Bz_downy>, <Bx_downy>}
</p>
<p>
This yields the following:
</p>
<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/madams120320142/B10S0_Streamlinesdownx_solved0000.png" style="padding:0; border:none"><img alt="down the barrel" crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/madams120320142/B10S0_Streamlinesdownx_solved0000.png" title="down the barrel" width="50%" /></a>
</p>
<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/madams120320142/B10S0_Streamlinesdowny_solvedbetter0000.png" style="padding:0; border:none"><img crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/madams120320142/B10S0_Streamlinesdowny_solvedbetter0000.png" width="50%" /></a>
</p>
<p>
Which makes much more sense.
</p>
visitstreamlinesCollidingFlowsPotential problems with projected streamlines?madamsWed, 03 Dec 2014 16:44:04 GMT
https://bluehound2.circ.rochester.edu/astrobear/blog/madams12032014
https://bluehound2.circ.rochester.edu/astrobear/blog/madams12032014<p>
So I am attempting to plot the streamlines for our colliding flows problem. Here is an example of the shear-0 case at 10.1 Myr (or frame 101) (also I did these under Erica' account on BH2 - hence the username haha).
</p>
<p>
The first image is down the barrel of the two flows (otherwise projecting the mass down the x-axis) thus the vertical axis is z and the horizontal is y.
</p>
<p>
The second image is a projection down the y-axis. Thus the vertical axis is x and the horizontal is z. This makes sense given that we've defined <a class="missing wiki">GxBounds</a> = 0d0,0d0,0d0,62.5d0,75d0,75d0. The two flows are colliding along x, so in the second image, they are coming in from top and bottom.
</p>
<p>
In both images I've plotted the column density maps for min = 60 and max = 1000. I did similarly for the min/max of the streamlines which are plotted on top of the column density maps. I also checked that they are scaled by magnitude. Now after talking with Erica we are not sure if these streamlines make any physical sense if we have defined a magnetic field along the flow axis (i.e. x). Ignore the visit axis labels as they are generic and don't define the dimensions of our problem.
</p>
<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/madams12032014/B10S0_streamlines_downx0000.png" style="padding:0; border:none"><img alt="down the barrel" crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/madams12032014/B10S0_streamlines_downx0000.png" title="down the barrel" width="50%" /></a>
</p>
<p>
<a href="https://bluehound2.circ.rochester.edu/astrobear/attachment/blog/madams12032014/B10S0_streamlines_downy0000.png" style="padding:0; border:none"><img alt="down y" crossorigin="anonymous" src="https://bluehound2.circ.rochester.edu/astrobear/raw-attachment/blog/madams12032014/B10S0_streamlines_downy0000.png" title="down y" width="50%" /></a>
</p>
<p>
In our problem.f90 we have defined the projections for streamlines like so:
</p>
<pre class="wiki"> !For 'projected' streamlines plot of the data down x:
CALL CreateProjection(projection)
Projection%dim=1
Projection%Field(1)%iD=By_field
Projection%Field(1)%component=BOTHCOMP
Projection%field(1)%name='By'
Projection%Field(2)%iD=Bz_field
Projection%Field(2)%component=BOTHCOMP
Projection%field(2)%name='Bz'
!For 'projected' streamlines plot of the data down y:
CALL CreateProjection(projection)
Projection%dim=2
Projection%Field(1)%iD=Bz_field
Projection%Field(1)%component=BOTHCOMP
Projection%field(1)%name='Bz'
Projection%Field(2)%iD=Bx_field
Projection%Field(2)%component=BOTHCOMP
Projection%field(2)%name='Bx'
</pre><p>
So in Visit I defined a few expressions to be able to plot the streamlines. For down the x-axis (which correspond to the mass1 CDMs):
</p>
<p>
By_downx = array_decompose(projections, 1)
</p>
<p>
Bz_downx = array_decompose(projections, 2)
</p>
<p>
which you can create the expression for the vector Byz_downx = {<By_downx, Bz_downx>} to plot the streamlines like I have above. The first component corresponds should correspond to the right axis if the horizontal component is truly y. Thus the second component will correspond to z if the vertical is truly z. So I think I have these lined up correctly? For projections down the y-axis (corresponding to mass2 CDMs):
</p>
<p>
Bz_downy = array_decompose(projections, 1)
</p>
<p>
Bx_downy = array_decompose(projections, 2)
</p>
<p>
you can create Bzx_downy = {<Bz_downy>, <Bx_downy>}. Clearly from the size of the box we know the horizonal component is z, and thus the first parameter in our vector should be Bz. Similarly for the second being x. However the streamlines don't seem right? Not sure what is going on.
</p>
visitstreamlinesCollidingFlows