Example of routines in nebula.ijs: This run constructs a nebula of radius 0.05-0.10 parsec around a 120,000K (blackbody) central star. =================================================================== load'nebula.ijs' 55 250 16 4 SETUP'' Enter stellar temperature: 120000 Enter stellar luminosity (solar units): 100 Enter distance from star (in parsecs): 0.05 Enter Hydrogen density: 800 ┌──────┬───┬────┬───┐ │120000│100│0.05│800│ └──────┴───┴────┴───┘ Enter heavy element scale factor: 0.8 NH, Ne, new Ne = 800 880 911.525 tau0 = 0 --------------------------------------------- T = 12372.8 Ne = 919.317 0.0109209 0.989079 0 0 0 0 0.00257356 0.412586 0.58484 0 0 0 0.000357264 0.0949457 0.817316 0.0853684 0.00201268 3.44689e_17 0.00023905 0.108348 0.761156 0.128986 0.00126697 4.11093e_6 0.00158642 0.147297 0.722899 0.125278 0.0029361 3.55969e_6 2.92234e_5 0.0236192 0.869887 0.105945 0.000518769 7.21542e_7 ------- Fractional Cooling by Ion: ---------- 0.0805714 0 0 0 0 0 0 0 0 0 0 0 0.000124434 0.00815267 0.0617447 0.0085034 0 0 1.23638e_5 0.0207026 0.012115 0.000652124 0 0 0.000190844 0.0730318 0.53739 0.058753 0 0 0 0.000405177 0.133123 0.00419074 0.000337223 0 RUN'' Step size (parsecs): 0.0025 Number of steps: 20 --------------------------------------------- R = 0.0025 T = 12453.1 Ne = 911.801 0.0130916 0.986908 0 0 0 0 0.0031727 0.483071 0.513756 0 0 0 0.000433623 0.102178 0.831335 0.0649535 0.00109985 1.90427e_17 0.000325962 0.131977 0.770974 0.0959748 0.00074654 2.2868e_6 0.00244037 0.185202 0.728096 0.0825602 0.00170003 1.9831e_6 3.55243e_5 0.0269594 0.896936 0.0757754 0.000292776 3.97662e_7 --------------------------------------------- R = 0.005 T = 12623.8 Ne = 902.121 0.0162853 0.983715 0 0 0 0 0.00397658 0.569928 0.426096 0 0 0 0.000542067 0.110722 0.841458 0.0467561 0.000521918 9.19788e_18 0.000468498 0.16541 0.767471 0.0662683 0.000381187 1.0902e_6 0.0039789 0.236427 0.711223 0.0475402 0.000829921 9.23544e_7 4.32762e_5 0.0306255 0.919866 0.0493229 0.000141749 1.87589e_7 --------------------------------------------- R = 0.0075 T = 12744.1 Ne = 890.772 0.0209402 0.97906 0 0 0 0 0.00499168 0.662617 0.332392 0 0 0 0.000687717 0.120777 0.845792 0.0325205 0.000223035 3.99669e_18 0.000703531 0.211275 0.744422 0.0434245 0.00017418 4.58444e_7 0.00669866 0.301056 0.667543 0.0243534 0.00034853 3.66913e_7 5.27462e_5 0.0343565 0.935685 0.0298447 6.0855e_5 7.78592e_8 --------------------------------------------- R = 0.01 T = 12816.3 Ne = 877.826 0.0279984 0.972002 0 0 0 0 0.00624844 0.750707 0.243045 0 0 0 0.000891142 0.133437 0.843771 0.0218158 8.52863e_5 1.55271e_18 0.00111577 0.275695 0.696383 0.0267369 6.99646e_5 1.67035e_7 0.011712 0.381904 0.59528 0.0109805 0.000123864 1.22314e_7 6.47146e_5 0.0380783 0.945089 0.0167452 2.29605e_5 2.81805e_8 --------------------------------------------- R = 0.0125 T = 12849.5 Ne = 862.786 0.0390924 0.960908 0 0 0 0 0.00782914 0.823856 0.168315 0 0 0 0.00118761 0.150569 0.834195 0.0140193 2.90388e_5 5.3732e_19 0.00188181 0.366125 0.616843 0.0151262 2.40657e_5 5.15015e_8 0.0212484 0.479299 0.495099 0.00431639 3.67023e_5 3.37666e_8 8.09088e_5 0.0418695 0.949296 0.00874575 7.61742e_6 8.91546e_9 --------------------------------------------- R = 0.015 T = 12814.8 Ne = 843.556 0.0572597 0.94274 0 0 0 0 0.00999771 0.877323 0.11268 0 0 0 0.00163649 0.17551 0.814354 0.008491 8.70257e_6 1.6392e_19 0.00338505 0.488359 0.500716 0.00753314 6.79038e_6 1.2916e_8 0.0396739 0.586329 0.372537 0.00145168 8.83028e_6 7.54483e_9 0.000105556 0.0461333 0.949477 0.00428235 2.23549e_6 2.48578e_9 --------------------------------------------- R = 0.0175 T = 12663.6 Ne = 815.253 0.0883914 0.911609 0 0 0 0 0.0134007 0.911936 0.0746629 0 0 0 0.00233462 0.214248 0.7787 0.00471558 2.24317e_6 4.32432e_20 0.00647896 0.637451 0.352991 0.00307821 1.46756e_6 2.47435e_9 0.0750918 0.682353 0.242151 0.000402329 1.64213e_6 1.30776e_9 0.000148362 0.0516306 0.946243 0.0019779 5.83661e_7 6.1724e_10 --------------------------------------------- R = 0.02 T = 12357.8 Ne = 768.274 0.143931 0.856069 0 0 0 0 0.0195327 0.930345 0.0501218 0 0 0 0.00341179 0.277281 0.717002 0.00230447 4.79489e_7 9.61983e_21 0.0131783 0.785342 0.200542 0.000936873 2.21473e_7 3.34088e_10 0.141262 0.731067 0.127584 8.67115e_5 2.19558e_7 1.65752e_10 0.000232754 0.0595931 0.939308 0.000865927 1.36491e_7 1.39039e_10 --------------------------------------------- R = 0.0225 T = 11944.6 Ne = 688.176 0.2412 0.7588 0 0 0 0 0.0317073 0.933234 0.0350586 0 0 0 0.00488451 0.378676 0.615493 0.000946856 8.32131e_8 1.81583e_21 0.0286676 0.884976 0.0861598 0.00019623 2.2053e_8 3.08433e_11 0.255071 0.69388 0.0510353 1.39271e_5 2.02129e_8 1.50881e_11 .000417911 0.0719583 0.927257 0.000366881 2.98387e_8 3.04613e_11 --------------------------------------------- R = 0.025 T = 11444 Ne = 564.514 0.392236 0.607764 0 0 0 0 0.0568793 0.91688 0.0262408 0 0 0 0.00622691 0.524343 0.469118 0.000312719 1.18438e_8 3.05053e_22 0.067699 0.904512 0.0277608 2.81318e_5 1.52343e_9 2.11856e_12 0.420916 0.564249 0.0148334 1.68165e_6 1.34214e_9 1.06971e_12 0.000878346 0.0928879 0.906076 0.000158006 6.76774e_9 7.46655e_12 --------------------------------------------- R = 0.0275 T = 10877.5 Ne = 411.429 0.577896 0.422104 0 0 0 0 0.107352 0.871379 0.0212691 0 0 0 0.00640461 0.691418 0.302098 7.97561e_5 1.41884e_9 4.82766e_23 0.167209 0.825793 0.00699496 2.90457e_6 8.20818e_11 1.25056e_13 0.609395 0.387463 0.00314194 1.63176e_7 7.35657e_11 6.94606e_14 0.00218208 0.130863 0.866882 7.29668e_5 1.79164e_9 2.3705e_12 --------------------------------------------- R = 0.03 T = 10320.8 Ne = 265.909 0.749847 0.250153 0 0 0 0 0.202444 0.779286 0.0182698 0 0 0 0.00517985 0.836629 0.158176 1.50185e_5 1.41123e_10 7.1437e_24 0.377344 0.621296 0.00135963 2.11032e_7 3.46244e_12 6.39824e_15 0.773101 0.226393 0.000505862 1.36593e_8 3.73384e_12 4.64137e_15 0.00625651 0.199663 0.794045 3.6306e_5 5.69818e_10 1.00815e_12 --------------------------------------------- R = 0.0325 T = 9813.58 Ne = 154.401 0.872778 0.127222 0 0 0 0 0.363201 0.621161 0.0156389 0 0 0 0.0034015 0.930447 0.0661496 2.03436e_6 1.12314e_11 9.42825e_25 0.658793 0.341023 0.000183914 9.65736e_9 1.00904e_13 2.48149e_16 0.885665 0.11427 6.47522e_5 1.01781e_9 1.81273e_13 3.25834e_16 0.0188979 0.309538 0.671546 1.85174e_5 2.03307e_10 5.30777e_13 --------------------------------------------- R = 0.035 T = 9404.67 Ne = 82.358 0.94179 0.0582102 0 0 0 0 0.569614 0.417984 0.0124018 0 0 0 0.00191931 0.974362 0.0237187 2.28153e_7 8.01597e_13 1.22101e_25 0.858124 0.141857 1.89969e_5 3.27651e_10 2.32158e_15 8.20553e_18 0.947794 0.052199 7.43966e_6 7.40245e_11 9.04901e_15 2.53472e_17 0.0476379 0.430471 0.521881 9.35603e_6 7.6131e_11 3.17916e_13 --------------------------------------------- R = 0.0375 T = 9228.58 Ne = 44.8474 0.972022 0.0279784 0 0 0 0 0.732414 0.258563 0.0090228 0 0 0 0.00106622 0.989417 0.00951715 3.4921e_8 8.23675e_14 2.26416e_26 0.938946 0.061051 2.60351e_6 1.77895e_11 8.79546e_17 4.40697e_19 0.974929 0.0250703 1.14232e_6 7.60209e_12 6.54692e_16 2.79504e_18 0.0780389 0.501834 0.420122 5.10731e_6 3.15936e_11 2.08416e_13 --------------------------------------------- R = 0.04 T = 9052.67 Ne = 28.3912 0.983567 0.0164329 0 0 0 0 0.819551 0.174361 0.00608773 0 0 0 0.000681282 0.994262 0.00505677 9.68203e_9 1.55409e_14 6.62654e_27 0.966617 0.0333824 6.76172e_7 2.47597e_12 8.59854e_18 5.14868e_20 0.985251 0.0147483 3.19235e_7 1.44194e_12 8.76195e_17 4.81986e_19 0.0987802 0.538249 0.362968 3.03002e_6 1.42437e_11 1.25886e_13 --------------------------------------------- R = 0.0425 T = 8970.31 Ne = 22.1864 0.98732 0.0126802 0 0 0 0 0.857688 0.138028 0.00428471 0 0 0 0.00053448 0.99593 0.00353575 4.65629e_9 5.48399e_15 2.95971e_27 0.975456 0.0245433 3.33139e_7 8.45229e_13 2.19184e_18 1.32762e_20 0.98861 0.0113898 1.65056e_7 5.46001e_13 2.48497e_17 1.46845e_19 0.115793 0.560972 0.323233 1.98754e_6 7.40316e_12 7.28557e_14 --------------------------------------------- R = 0.045 T = 8909.74 Ne = 18.1399 0.989501 0.0104986 0 0 0 0 0.885028 0.112024 0.00294765 0 0 0 0.000438228 0.996978 0.00258384 2.44769e_9 2.1075e_15 1.37751e_27 0.980634 0.0193658 1.87102e_7 3.42322e_13 6.57675e_19 3.8047e_21 0.990564 0.00943593 9.72379e_8 2.34745e_13 7.93303e_18 4.7569e_20 0.135941 0.580567 0.28349 1.27733e_6 3.72041e_12 3.86229e_14 --------------------------------------------- R = 0.0475 T = 8865.62 Ne = 15.1895 0.99101 0.00899011 0 0 0 0 0.905841 0.092114 0.00204521 0 0 0 0.000367723 0.997689 0.00194374 1.36722e_9 8.75989e_16 6.7785e_28 0.984111 0.015889 1.13307e_7 1.5403e_13 2.22364e_19 1.20659e_21 0.991916 0.00808352 6.14295e_8 1.10185e_13 2.80167e_18 1.66742e_20 0.157251 0.595285 0.247463 8.30699e_7 1.90097e_12 2.03546e_14 --------------------------------------------- R = 0.05 T = 8827.28 Ne = 12.8607 0.992186 0.00781389 0 0 0 0 0.922502 0.0760675 0.00143085 0 0 0 0.000311901 0.998202 0.00148567 7.906e_10 3.82881e_16 3.47028e_28 0.986701 0.0132986 7.1409e_8 7.34654e_14 8.07913e_20 4.08866e_22 0.992971 0.00702853 4.01854e_8 5.44091e_14 1.05426e_18 6.16945e_21 0.17923 0.605539 0.215231 5.46507e_7 9.87756e_13 1.07973e_14 ---- Threshold Optical Depths by Ion: ------- 366.705 0 0 0 0 0 27.4095 11.1626 0 0 0 0 0.00114976 0.0815937 0.0209629 0.000224716 1.85607e_6 0 0.0476134 0.0273324 0.00383439 0.000174566 4.24756e_7 0 0.0546728 0.14705 0.0620894 0.000819778 1.03078e_5 0 0.00366092 0.0236907 0.039425 0.000669921 9.64732e_7 0 ------- Fractional Cooling by Ion: ---------- 0.577162 0 0 0 0 0 0 0 0 0 0 0 .000196984 0.0140541 4.03762e_5 1.57174e_11 0 0 0.0864179 0.00526717 6.95696e_9 6.6555e_17 0 0 0.185167 0.00445826 6.63255e_8 1.21657e_13 0 0 0 0.0416001 0.0856357 1.30282e_8 3.14829e_12 0 ======================================================================= Example of gas cooling for 60 years: SETUP'' Enter stellar temperature: 100000 Enter stellar luminosity (solar units): 100 Enter distance from star (in parsecs): 0.1 Enter Hydrogen density: 1e3 ┌──────┬───┬───┬────┐ │100000│100│0.1│1000│ └──────┴───┴───┴────┘ Enter heavy element scale factor: 1.5 NH, Ne, new Ne = 1000 1100 1058.66 tau0 = 0 --------------------------------------------- T = 9423.57 Ne = 1067.52 0.0451111 0.954889 0 0 0 0 0.0282478 0.840943 0.130809 0 0 0 0.00496136 0.368398 0.623776 0.00286317 2.07472e_6 5.46025e_24 0.00725667 0.690626 0.29956 0.00255684 8.07072e_7 1.21346e_10 0.0388183 0.769139 0.190826 0.00121448 1.72443e_6 8.04325e_11 0.00124452 0.17186 0.824021 0.00287456 2.1894e_7 8.15465e_12 ------- Fractional Cooling by Ion: ---------- 0.0196782 0 0 0 0 0 0 0 0 0 0 0 0.00219067 0.0165133 0.0186245 7.33164e_5 0 0 0.000384629 0.185101 0.00699341 3.06728e_6 0 0 0.00541043 0.34692 0.17898 0.00127562 0 0 0 0.00697661 0.210809 6.55138e_5 3.53413e_7 0 EVOL'' Step size (years): 2 Number of steps: 30 --------------------------------------------- t = 2 T = 8848.98 Ne = 1035.43 0.072705 0.927295 0 0 0 0 0 0.0528261 0.834964 0.11221 0 0 0 0 0.0270126 0.561249 0.410768 0.000968938 1.10991e_6 1.65242e_24 0 0.0336347 0.722163 0.243264 0.000938397 1.65307e_7 4.31622e_11 0 0.0559692 0.78298 0.160556 0.000495108 1.85198e_7 2.09888e_11 0 0.00453695 0.253554 0.739823 0.00208588 1.11851e_7 2.68766e_12 0 --------------------------------------------- t = 4 T = 8393.02 Ne = 1003.48 0.100347 0.899653 0 0 0 0 0 0.077968 0.826212 0.0958199 0 0 0 0 0.0562331 0.676108 0.267342 0.000316797 5.82989e_7 4.83712e_25 0 0.0617079 0.741746 0.196206 0.000340764 3.3138e_8 1.49018e_11 0 0.0738031 0.791742 0.134253 0.000201954 1.93959e_8 5.40942e_12 0 0.00920216 0.327571 0.661728 0.00149884 5.60446e_8 8.56315e_13 0 --------------------------------------------- t = 6 T = 8024.2 Ne = 970.955 0.128654 0.871346 0 0 0 0 0 0.104203 0.814672 0.0811253 0 0 0 0 0.0906616 0.739133 0.170108 9.71809e_5 2.95717e_7 1.32734e_25 0 0.0918396 0.751494 0.156547 0.000119032 6.2361e_9 4.85913e_12 0 0.092692 0.796069 0.111159 8.04663e_5 1.86288e_9 1.32811e_12 0 0.015365 0.395585 0.587993 0.00105755 2.70593e_8 2.56232e_13 0 --------------------------------------------- t = 8 T = 7719.78 Ne = 937.966 0.157511 0.842489 0 0 0 0 0 0.131406 0.800497 0.0680967 0 0 0 0 0.128468 0.765677 0.105826 2.79719e_5 1.44839e_7 3.41423e_26 0 0.1237 0.752618 0.123641 3.99436e_5 1.10196e_9 1.4962e_12 0 0.11255 0.796249 0.0911695 3.12933e_5 1.63837e_10 3.07844e_13 0 0.0230544 0.457164 0.519049 0.000732685 1.25876e_8 7.1822e_14 0 --------------------------------------------- t = 10 T = 7465.12 Ne = 905.912 0.185675 0.814325 0 0 0 0 0 0.158358 0.784563 0.0570793 0 0 0 0 0.16674 0.767619 0.065633 7.95873e_6 7.053e_8 8.70233e_27 0 0.155663 0.746699 0.0976242 1.3444e_5 1.96167e_10 4.58261e_13 0 0.132475 0.792827 0.0746863 1.23237e_5 1.45606e_11 7.25693e_14 0 0.031893 0.510064 0.457537 0.000505983 5.82125e_9 2.2179e_14 0 --------------------------------------------- -----> Skip ahead to the end <--------------------------- t = 58 T = 5624.6 Ne = 352.487 0.681674 0.318326 0 0 0 0 0 0.665744 0.333802 0.000453386 0 0 0 0 0.76658 0.23342 1.53663e_7 2.54372e_16 1.99544e_16 5.75828e_43 0 0.727201 0.272598 0.000200132 1.35591e_17 0 2.95749e_27 0 0.554678 0.444975 0.000346758 2.61358e_15 0 1.46249e_30 0 0.378012 0.608355 0.0136339 1.91436e_8 0 0 0 --------------------------------------------- t = 60 T = 5604.75 Ne = 338.453 0.694348 0.305652 0 0 0 0 0 0.679036 0.320596 0.000368083 0 0 0 0 0.778834 0.221166 8.81265e_8 1.09181e_16 8.54783e_17 1.1583e_43 0 0.74038 0.259466 0.000154048 1.88595e_17 0 7.29076e_28 0 0.567469 0.432255 0.00027556 6.99502e_16 0 2.98067e_31 0 0.391618 0.596675 0.0117067 1.23332e_8 0 0 0 ------- Fractional Cooling by Ion: ---------- 1.28031e_5 0 0 0 0 0 0 0 0 0 0 0 0.00168821 0.00590374 0.000313294 3.1853e_7 0 0 0.000135961 0.15168 0.036975 1.03199e_8 0 0 0.00306196 0.0814929 0.245516 0.0048223 0 0 0 0.0247115 0.443681 3.35613e_6 1.59197e_6 0 NB. plot logs of relative abundances of the ions of oxygen (4) NB. and nitrogen (3) as function of time. load'plot' 'pensize 3'plot tim; 10^.1e_4+ |: 4{"2 >zz 'pensize 3'plot tim; 10^.1e_4+ |: 3{"2 >zz