<div class="gmail_quote">Hi Chris,<div><br></div><div>My group has implemented angular momentum in 2D and 1D simulations by - as Sean mentioned - adding a mass scalar which represents angular momentum.  If you just use velz in cylindrical or spherical coordinates, you end up creating some ambiguity about how total and kinetic energy should be calculated.  We incorporate the centrifugal force term in the gravity_accelOneRow .  We also include an alpha-viscosity treatment of angular momentum transport.  We edited the diffuse_species module and a few other spots to do this.  We describe our implementation and some tests of this we did in some of our recent papers: <a href="http://arxiv.org/PS_cache/arxiv/pdf/1108/1108.1415v1.pdf" target="_blank">http://arxiv.org/PS_cache/arxiv/pdf/1108/1108.1415v1.pdf</a>  <a href="http://arxiv.org/PS_cache/arxiv/pdf/0910/0910.4989v1.pdf" target="_blank">http://arxiv.org/PS_cache/arxiv/pdf/0910/0910.4989v1.pdf</a> .  Let me know if you run into any trouble or have some more questions about it.</div>


<div><br></div><font color="#888888"><div>Chris Lindner<div><a href="mailto:lindner@astro.as.utexas.edu" target="_blank">lindner@astro.as.utexas.edu</a></div><div><div>NSF Graduate Fellow</div><div>Department of Astronomy</div>

<div>The University of Texas at Austin</div>
</div></div>
</font></div><br>