| 101 | |
| 102 | === Cylindrical Symmetry (2.5D) === |
| 103 | This section gives some brief details on how you can take advantage of cylindrical coordinates if your problem has such symmetry. |
| 104 | |
| 105 | * '''''{{{iCylindrical}}}''''' Controls whether astrobear is to use cylindrical coordinates. The default is 0 which leaves everything in Cartesian, 1 is for no angular momentum, and 2 is for with angular momentum. |
| 106 | |
| 107 | Basically when iCylindrical is not 0, astrobear changes the differential operators and how each axis is treated. In other words, astrobear switches to cylindrical coordinates. The x-axis becomes the radial axis, the y-axis becomes the cylindrical z-axis, and the z-axis is nothing since this is only implemented for ndim = 2. |
| 108 | |
| 109 | Furthermore, if your problem has rotation then iCylindrical should be set to 2. This will change the ivz component to the angular momentum component. Otherwise, if iCylindrical=1 then there is only ivx (radial) and ivy (cylindrical z). |
| 110 | |
| 111 | One other important thing to note...if iCylindrical is not 0, then mydim (a global parameter in physics declarations) will be set to 3, otherwise it is set to ndim. |