Hi YT users,<br>
<br>
I would like to make a (spherical and/or cylindrical) profile of the velocity dispersion.<br>
In order to do so, I need to compute the mean velocity (BulkVelocity) 
for a given shell (or ring), and pass this value to newly defined 
velocity dispersion field, and then make a profile for the velocity 
dispersion field.<br>
Is there a particular way in yt to handle this?<br>
<br>
Junhwan<br>
<br>
P.S.: A few weeks ago I posted a similar question and Matt gave me a following hint.<br>
@derived_field(name = "deltaLocalCircularVelocity")<br>
def delta_LocalCircularVelocity(fi<div id=":5e">eld, data):<br>
  fn, x, cv = data.get_field_parameter("local_circular_velocity")<br>
  cv0 = na.interp(data[fn].ravel(), x, cv)<br>
  cv0 = cv0.reshape(data[fn].shape)<br>
  delta_cv = cv0 - data["LocalCircularVelocity"]<br>
  return delta_cv<br>
<br>
You'll need to do something like this on the data object from which<br>
you are getting your data:<br>
<br>
  sp.set_field_parameter("local_circular_velocity",<br>
      (fn, prof[fn], prof["LocalCircularVelocity"]))<br>
<br>
However, it does not work for me and it is also complicate to understand and modify it.<div class="yj6qo ajU"><div id=":9f" class="ajR" tabindex="0"><img class="ajT" src="images/cleardot.gif"></div></div></div><br><br clear="all">
<br>-- <br>--------------------------------------------------------------<br>Jun-Hwan Choi, Ph.D.<br>Department of Physics and Astronomy, University of Kentucky<br>Tel: (859) 897-6737        Fax: (859) 323-2846<br>Email: <a href="mailto:jhchoi@pa.uky.edu" target="_blank">jhchoi@pa.uky.edu</a>   URL: <a href="http://www.pa.uky.edu/%7Ejhchoi" target="_blank">http://www.pa.uky.edu/~jhchoi</a><br>
--------------------------------------------------------------<br><br><br>