| 28 | |
| 29 | |
| 30 | {{{ |
| 31 | CALL CreateOutflow(Particle%OutflowObj) |
| 32 | call addtracer(particle%outflowobj%iTracer,'Outflow_Tracer') |
| 33 | Particle%OutflowObj%duration=3e9/timescale |
| 34 | Particle%OutflowObj%radius=REAL(radiusw,qPREC)*sink_dx |
| 35 | Particle%OutflowObj%thickness=REAL(radiusw,qPREC)*sink_dx |
| 36 | Particle%OutflowObj%open_angle=Pi/2 |
| 37 | Particle%OutflowObj%theta=0 |
| 38 | Particle%OutflowObj%phi=0 |
| 39 | Particle%OutflowObj%density=outflowdensity/rscale |
| 40 | Particle%OutflowObj%temperature=30000/tempscale |
| 41 | Particle%OutflowObj%speed=0/velscale |
| 42 | Particle%OutflowObj%fade=0d0 |
| 43 | Particle%OutflowObj%velocity(1:nDim)=Particle%Q(imom(1:nDim)) |
| 44 | Particle%OutflowObj%position= Particle%xloc |
| 45 | Particle%OutflowObj%begin=-1d0 |
| 46 | CALL UpdateOutflow(Particle%OutflowObj) |
| 47 | }}} |