<div dir="ltr"><div><div><div><div><div><div><div><div>Hi yt-users!<br><br></div>So I want to calculate the turbulent energy in a stripped tail of gas.  This tail is moving with a bulk flow in the z-direction, which varies as a function of z, and is rotating, and v_theta will also likely change as a function of z (and while I am at it, might as well check v_r).  So I would like to do something like this:<br>
<br></div>def _Eturb(field,data):<br></div>    vrmean = value I calculate in the code<br></div>    vthetamean = value I calculate in the code<br></div>    vzmean = value I calculate in the code<br></div>    return(pow(data["cyl_RadialVelocity"]-vrmean,2) + pow(data["cyl_TangentialVelocity"]-vthetamean,2) + pow(data["z-velocity"]-vzmean,2))<br>
<br></div>add_field("Eturb",function=_Eturb,units="velocityunits")<br><br>totalEturb = np.zeros(100,float)<br>pf = load("id0/rps."+outstring[i]+"<div>.vtk",parameters={"TimeUnits":9.78262,"LengthUnits":26.,"DensityUnits":9.999e-30,"VelocityUnits":2.6e8})<br>
<br></div><div>alld = pf.h.all_data()<br></div><div>i = 0<br></div><div>while i < 100:<br></div><div>    disk = pf.h.disk([0.0,0.0,i*0.385/2.],[0.0,0.0,1.0],1.0,0.385)<br></div><div>    vrmean = disk.quantities["WeightedAverageQuantity"]("cyl_RadialVelocity","Ones")<br>
</div><div>    vthetamean = disk.quantities["WeightedAverageQuantity"]("cyl_TangentialVelocity","Ones")<br></div><div>    vzmean = disk.quantities["WeightedAverageQuantity"]("z-velocity","Ones")<br>
<br></div><div>    totalEturb[i] = disk.quantities["TotalQuantity"]["Eturb"]<br></div><div>    i = i+1<br><br></div>I am not sure if/how I can make the def statement read in vrmean, vthetamean, and vzmean?<br>
<br></div><div>Thanks!<br><br>Stephanie<br></div></div>