Before reading the atomic data from /aalib/atomic. Before reading the plot parameters. &PLOTPAR APLOT = ./bravo/prd5.5_fig.par / Before reading the data files. Note we can plot many density profiles, but only one abundance file as specified by iabund. bravo: Before zeroing and reading in. icell vel (cm/s) rho (cgs) m_i (m_sun) ratio sum X Ni-56 1 3.1912E+07 3.1014E-08 1.3690E-03 76.3608 1.0000 0.0000 2 4.1374E+07 2.9585E-08 2.9090E-03 58.8975 1.0000 0.0000 3 4.9155E+07 2.6023E-08 4.6037E-03 49.5743 1.0000 0.0000 4 5.6938E+07 2.0870E-08 6.4694E-03 42.7979 1.0009 0.0027 5 6.3063E+07 2.2835E-08 8.5229E-03 38.6412 0.9999 0.0074 6 6.9360E+07 2.0056E-08 1.0781E-02 35.1330 1.0000 0.0000 7 7.7804E+07 1.3320E-08 1.3265E-02 31.3201 1.0003 0.0041 8 8.4049E+07 1.6388E-08 1.5997E-02 28.9929 0.9999 0.0065 9 9.1714E+07 1.2464E-08 1.9005E-02 26.5698 0.9993 0.0210 10 9.9954E+07 1.0734E-08 2.2317E-02 24.3795 0.9993 0.0211 11 1.0929E+08 8.7437E-09 2.5961E-02 22.2969 0.9995 0.0290 12 1.1952E+08 7.3405E-09 2.9968E-02 20.3884 0.9997 0.0357 13 1.3076E+08 6.1453E-09 3.4379E-02 18.6359 1.0002 0.0571 14 1.4487E+08 4.4394E-09 3.9231E-02 16.8208 0.9996 0.1042 15 1.6227E+08 3.1890E-09 4.4568E-02 15.0171 1.0009 0.1438 16 1.8369E+08 2.2461E-09 5.0442E-02 13.2660 1.0006 0.2373 17 2.1052E+08 1.5183E-09 5.6900E-02 11.5753 0.9997 0.3659 18 2.4118E+08 1.1131E-09 6.4004E-02 10.1038 0.9999 0.5715 19 2.7252E+08 9.2599E-10 7.1815E-02 8.9418 1.0000 0.7064 20 3.0360E+08 8.1655E-10 8.0404E-02 8.0264 0.9998 0.7783 21 3.3276E+08 7.8461E-10 8.9849E-02 7.3231 0.9997 0.8239 22 3.6038E+08 7.6770E-10 1.0023E-01 6.7618 1.0001 0.8417 23 3.8719E+08 7.4764E-10 1.1165E-01 6.2936 0.9997 0.8540 24 4.1303E+08 7.4442E-10 1.2420E-01 5.8999 0.9995 0.8598 25 4.3799E+08 7.4920E-10 1.3800E-01 5.5637 0.9999 0.8771 26 4.6084E+08 7.3369E-10 1.5180E-01 5.2878 1.0003 0.8864 27 4.8175E+08 7.2908E-10 1.6560E-01 5.0583 0.9996 0.8856 28 5.0163E+08 7.0457E-10 1.7940E-01 4.8578 1.0004 0.8953 29 5.2014E+08 7.0095E-10 1.9320E-01 4.6849 0.9998 0.8912 30 5.3735E+08 7.0384E-10 2.0700E-01 4.5349 1.0004 0.8951 31 5.5384E+08 6.8991E-10 2.2080E-01 4.3999 0.9999 0.8918 32 5.6946E+08 6.8730E-10 2.3460E-01 4.2792 1.0004 0.8887 33 5.8440E+08 6.8103E-10 2.4840E-01 4.1698 0.9997 0.8888 34 5.9924E+08 6.5155E-10 2.6220E-01 4.0665 0.9999 0.8886 35 6.1302E+08 6.6894E-10 2.7600E-01 3.9751 1.0000 0.8845 36 6.2635E+08 6.6160E-10 2.8980E-01 3.8905 0.9999 0.8671 37 6.3958E+08 6.3892E-10 3.0360E-01 3.8100 0.9997 0.8716 38 6.5160E+08 6.7601E-10 3.1740E-01 3.7398 0.9995 0.8544 39 6.6425E+08 6.1848E-10 3.3120E-01 3.6685 1.0000 0.8581 40 6.7605E+08 6.3907E-10 3.4500E-01 3.6045 0.9999 0.8466 41 6.8727E+08 6.4960E-10 3.5880E-01 3.5457 1.0003 0.8490 42 6.9900E+08 6.0095E-10 3.7260E-01 3.4862 1.0004 0.8411 43 7.0995E+08 6.2320E-10 3.8640E-01 3.4324 1.0002 0.8321 44 7.1997E+08 6.6122E-10 4.0020E-01 3.3846 1.0002 0.8296 45 7.3118E+08 5.7386E-10 4.1400E-01 3.3327 0.9998 0.8271 46 7.4069E+08 6.5754E-10 4.2780E-01 3.2899 0.9997 0.8175 47 7.5118E+08 5.8023E-10 4.4160E-01 3.2440 1.0002 0.8173 48 7.6076E+08 6.1859E-10 4.5540E-01 3.2031 1.0002 0.7948 49 7.7013E+08 6.1690E-10 4.6920E-01 3.1642 0.9998 0.8035 50 7.7959E+08 5.9627E-10 4.8300E-01 3.1258 1.0003 0.7847 51 7.8926E+08 5.6918E-10 4.9680E-01 3.0875 0.9995 0.8002 52 7.9757E+08 6.4741E-10 5.1060E-01 3.0553 1.0004 0.7883 53 8.0670E+08 5.7652E-10 5.2440E-01 3.0207 0.9996 0.7892 54 8.1530E+08 5.9874E-10 5.3820E-01 2.9889 1.0000 0.7846 55 8.2463E+08 5.3989E-10 5.5200E-01 2.9551 0.9997 0.7714 56 8.3248E+08 6.2845E-10 5.6580E-01 2.9272 0.9997 0.7616 57 8.4110E+08 5.6110E-10 5.7960E-01 2.8972 1.0001 0.7481 58 8.4917E+08 5.8756E-10 5.9340E-01 2.8697 1.0000 0.7414 59 8.5700E+08 5.9434E-10 6.0720E-01 2.8434 1.0005 0.7426 60 8.6572E+08 5.2347E-10 6.2100E-01 2.8148 1.0002 0.7082 61 8.7356E+08 5.7120E-10 6.3480E-01 2.7895 1.0000 0.7168 62 8.8108E+08 5.8512E-10 6.4860E-01 2.7657 0.9997 0.7129 63 8.8929E+08 5.2647E-10 6.6240E-01 2.7402 1.0004 0.7050 64 8.9720E+08 5.3662E-10 6.7620E-01 2.7160 1.0004 0.6963 65 9.0455E+08 5.6776E-10 6.9000E-01 2.6940 1.0000 0.6955 66 9.1218E+08 5.3794E-10 7.0380E-01 2.6714 1.0001 0.6968 67 9.2036E+08 4.9315E-10 7.1760E-01 2.6477 0.9997 0.6783 68 9.2910E+08 4.5315E-10 7.3140E-01 2.6228 1.0001 0.6596 69 9.3668E+08 5.1340E-10 7.4520E-01 2.6016 1.0002 0.6552 70 9.4423E+08 5.0718E-10 7.5900E-01 2.5808 1.0001 0.6505 71 9.5204E+08 4.8238E-10 7.7280E-01 2.5596 1.0001 0.6364 72 9.6001E+08 4.6493E-10 7.8660E-01 2.5383 1.0003 0.6427 73 9.6877E+08 4.1569E-10 8.0040E-01 2.5154 0.9998 0.6409 74 9.7635E+08 4.7237E-10 8.1420E-01 2.4959 0.9994 0.6115 75 9.8466E+08 4.2392E-10 8.2800E-01 2.4748 1.0004 0.6249 76 9.9227E+08 4.5549E-10 8.4180E-01 2.4558 1.0001 0.6274 77 1.0013E+09 3.7748E-10 8.5560E-01 2.4337 1.0004 0.6138 78 1.0095E+09 4.0859E-10 8.6940E-01 2.4139 1.0002 0.6183 79 1.0185E+09 3.6599E-10 8.8320E-01 2.3926 0.9996 0.6201 80 1.0268E+09 3.9017E-10 8.9700E-01 2.3732 0.9997 0.6003 81 1.0354E+09 3.7041E-10 9.1080E-01 2.3535 0.9999 0.6005 82 1.0448E+09 3.3305E-10 9.2460E-01 2.3323 1.0003 0.6006 83 1.0540E+09 3.3428E-10 9.3840E-01 2.3120 0.9995 0.5857 84 1.0637E+09 3.1142E-10 9.5220E-01 2.2909 1.0000 0.5822 85 1.0714E+09 3.8594E-10 9.6600E-01 2.2744 1.0001 0.5738 86 1.0815E+09 2.8938E-10 9.7980E-01 2.2532 1.0003 0.5774 87 1.0917E+09 2.8122E-10 9.9360E-01 2.2321 0.9999 0.5758 88 1.1010E+09 3.0297E-10 1.0074E+00 2.2133 0.9996 0.5817 89 1.1119E+09 2.5380E-10 1.0212E+00 2.1916 1.0006 0.5619 90 1.1223E+09 2.6095E-10 1.0350E+00 2.1713 0.9999 0.5674 91 1.1327E+09 2.5616E-10 1.0488E+00 2.1513 0.9995 0.5556 92 1.1446E+09 2.1951E-10 1.0626E+00 2.1290 0.9998 0.5558 93 1.1568E+09 2.0965E-10 1.0764E+00 2.1065 0.9998 0.5575 94 1.1686E+09 2.1231E-10 1.0902E+00 2.0853 1.0004 0.5460 95 1.1821E+09 1.8160E-10 1.1040E+00 2.0614 0.9994 0.5574 96 1.1954E+09 1.8020E-10 1.1178E+00 2.0385 0.9997 0.5538 97 1.2094E+09 1.6732E-10 1.1316E+00 2.0149 0.9996 0.5431 98 1.2240E+09 1.5670E-10 1.1454E+00 1.9909 1.0002 0.5434 99 1.2401E+09 1.3858E-10 1.1592E+00 1.9650 0.9995 0.5473 100 1.2572E+09 1.2703E-10 1.1730E+00 1.9383 0.9998 0.5356 101 1.2752E+09 1.1735E-10 1.1868E+00 1.9109 0.9994 0.5302 102 1.2943E+09 1.0743E-10 1.2006E+00 1.8827 1.0002 0.5281 103 1.3156E+09 9.3370E-11 1.2144E+00 1.8523 0.9995 0.5445 104 1.3380E+09 8.5885E-11 1.2282E+00 1.8212 1.0003 0.5195 105 1.3625E+09 7.5819E-11 1.2420E+00 1.7885 0.9994 0.5166 106 1.3894E+09 6.6499E-11 1.2558E+00 1.7539 1.0001 0.5155 107 1.4162E+09 5.8400E-11 1.2684E+00 1.7207 1.0003 0.5127 108 1.4440E+09 4.9292E-11 1.2798E+00 1.6876 1.0002 0.5015 109 1.4708E+09 4.4751E-11 1.2901E+00 1.6568 0.9999 0.5164 110 1.4997E+09 3.6377E-11 1.2996E+00 1.6249 0.9998 0.5090 111 1.5254E+09 3.5812E-11 1.3082E+00 1.5975 0.9998 0.4965 112 1.5531E+09 2.9242E-11 1.3160E+00 1.5690 1.0007 0.5022 113 1.5805E+09 2.5904E-11 1.3231E+00 1.5418 1.0007 0.4777 114 1.6086E+09 2.2189E-11 1.3296E+00 1.5149 1.0009 0.4990 115 1.6381E+09 1.8530E-11 1.3354E+00 1.4876 1.0012 0.4638 116 1.6678E+09 1.6150E-11 1.3408E+00 1.4611 1.0004 0.4953 117 1.6983E+09 1.3777E-11 1.3456E+00 1.4349 1.0000 0.4731 118 1.7292E+09 1.1926E-11 1.3500E+00 1.4092 0.9996 0.4589 119 1.7585E+09 1.1045E-11 1.3540E+00 1.3857 1.0004 0.4534 120 1.7911E+09 8.6992E-12 1.3577E+00 1.3605 0.9997 0.4502 121 1.8265E+09 7.0135E-12 1.3610E+00 1.3342 1.0004 0.4562 122 1.8623E+09 6.0655E-12 1.3640E+00 1.3085 1.0003 0.4578 123 1.8999E+09 5.0539E-12 1.3667E+00 1.2826 0.9999 0.4064 124 1.9410E+09 4.0151E-12 1.3692E+00 1.2554 0.9990 0.4269 125 1.9828E+09 3.4472E-12 1.3715E+00 1.2290 1.0003 0.4143 126 2.0333E+09 2.4706E-12 1.3735E+00 1.1985 1.0008 0.4225 127 2.0918E+09 1.8440E-12 1.3754E+00 1.1649 1.0000 0.3902 128 2.1661E+09 1.2315E-12 1.3771E+00 1.1250 1.0002 0.3728 129 2.2792E+09 6.7628E-13 1.3786E+00 1.0692 1.0003 0.3221 130 2.8338E+09 9.2396E-14 1.3800E+00 0.8599 1.0006 0.3146 Getting the exponential fitted vefold and rhoc. Fitted vefold and rhoc for 1 day. vefold,rhoc,xmass_gram,xke_erg 2.826134790116200E+08 7.501981739819784E-09 2.744958000000000E+33 1.315445006873390E+51 Checking the Ni/Co-56, IPE mass, and v_ipe. Ni/Co56 mass,IPE mass,g-factor,v_ipe,vel(k) 0.8567140493515171 0.9953949814950958 0.8606774850971438 1.092909661442347E+9 1.101E+9 Si,S,IPE,Ni-56 at vel(k) 0.076 0.0526 0.672109752 0.543 1.101E+9 Si,S,IPE at vel(k+1) 0.0824 0.0569 0.6492619230000001 0.523 1.1119E+9 i,iatom,abund(i,8),awt(8) 1 30 0.49 16.0 After the model read-in. i,ifile,iatom,icell,iradio 1 1 30 130 4 rhoc,vefold 7.501981739819784E-9 2826.1347901162007 After logging density. Before density plot. Small tick spacing must be <= big spacing Before sm_limits. velmin,velmax,rhomin,rhomax 0.0E+0 26000.0 -13.1549015 -7.39794 Before sm_grelocate. Before connecting the curves. tmp1,tmp2 -8.124824 -12.120265 After connecting the curves. In sim_table In sim_table: x,y 8000.0 -8.522879 8000.0 -8.780113 The density plot is completed. Before abundance plot for file 1 ./bravo/prd5.5_mod.dat At the moment I allow only one model abundances to plotted at at time. Small tick spacing must be <= big spacing Before sm_limits. velmin_abund,velmax_abund,abundmin,abundmax 0.0E+0 26000.0 -5.0 0.0E+0 Before sm_grelocate. Before connecting the curves. After connecting the curves. The abundance plot is completed.