Changes between Version 6 and Version 7 of u/johannjc/scratchpad22
- Timestamp:
- 05/25/17 16:33:20 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/johannjc/scratchpad22
v6 v7 95 95 2D version of Tripathi without atmosphere 96 96 97 $\rho_c = 6\times 10^6 cm^{-3}$ 97 98 98 $ \rho_c = 6\times 10^8 cm^{-3}$99 $R_c = 1.5\times 10^{10} cm$ 99 100 100 $ R_c=1.5\times 10^{10} cm$101 $F = 2\times 10^{13} cm^{-2} s^{-1}$ 101 102 102 $F=2\times 10^{13} cm^{-2} s^{-1}$ 103 104 $T= 10^3 K$ 105 103 $T = 10^5 K$ 106 104 107 105 embedded in a lower density pressure balanced ambient with a density ratio 108 106 109 107 $\chi=100$ 108 109 For reference, calculations done in setting up grid: 110 111 {{{#!fortran 112 Clump%density=rho*chi/nscale 113 Clump%Temperature=Temp/TempScale/chi 114 Ambient%density=rho/nscale 115 Ambient%pressure=ambient%density*Temp/TempScale 116 Clump%position=(/0d0,0d0,0d0/) 117 Clump%radius=radius ! in CU already 118 }}} 110 119 111 120 With a cross section for photo-ionization of 6.3e-18 and our clump's column density of 9e18, we would expect the tau=1 surface to be 5% of the way in towards the clump center.