subroutine begin * * initialize variables for central boundary condition * implicit double precision(a-h,o-z) common/shells/ sa(400),ra(400),ba(400),pa(400),ta(400), 1 ea(400),xca(400),fca(400),s(400),r(400),b(400), 2 p(400),t(400),xc(400),e(400),sk(400), 3 rk(400),bk(400),pk(400),tk(400) common/contrl/ds,g,sm,wc,it,nite,ja,jb,j,k,l common/temp/s1,r1,s2,r2,b2,p2,t2,e2,xc2,xo2,fca2,f2,q2,w2,c common/fin/ifin j=2 k=0 l=2 s1=-10. r1=-10. s2=s(1) r2=r(1) b2=b(1) p2=p(1) t2=t(1) ifin = 0 return end ************************************************************************ block data * * block data: initiallizes some physical constants, and also the * constants used by the conductive opacity regimes. the numbers * have been taken from itoh et al (1983; 1984). * Also added Itoh and company Bremsstrahlung neutrino coefficients. * See Itoh etal 1989, ApJ, 339, 354. + erratum 360, 741. * implicit double precision(a-h,l,o-z) common/c/az(3),z(3),taux,ts,ps,sm,rs,err,te,gs,rho,ol1x,rl1x,el1x 1,atgx,rtgx,tl,pl,xmass,acu,am,bk,an0,pi * * the following are the fitting coefficients quoted by itoh * "cdats" ==> carbon, solid state * "cdatl" ==> carbon, liquid metal state * "odats" ==> oxygen, solid state * "odatl" ==> oxygen, liquid metal state * and so on. * common /hydatl/ alhy(4),blhy(3),clhy(3),dlhy(4),elhy(3),flhy(3) common /hedatl/ alhe(4),blhe(3),clhe(3),dlhe(4),elhe(3),flhe(3) common/hedats/ sahe(5),sbhe(3),sehe(5),sfhe(3),sihe(5),sjhe(3), 1 sphe(5), sqhe(3), alphahe(4), betahe(4), sche, 2 sdhe, sghe, shhe, skhe, slhe, srhe, sshe common /cdatl/ alc(4), blc(3), clc(3), dlc(4), elc(3), flc(3) common/cdats/ sac(5), sbc(3), sec(5), sfc(3), sic(5), sjc(3), 1 spc(5), sqc(3), alphac(4), betac(4), scc, 2 sdc, sgc, shc, skc, slc, src, ssc common /odatl/ alo(4), blo(3), clo(3), dlo(4), elo(3), flo(3) common/odats/ sao(5), sbo(3), seo(5), sfo(3), sio(5), sjo(3), 1 spo(5), sqo(3), alphao(4), betao(4), sco, 2 sdo, sgo, sho, sko, slo, sro, sso * * the following are the Bremsstrahlung neutrino rate fitting * coefficients quoted by Itoh. Common'd blocks are: * Helium: soldath, phondath, liqdath, quandath * Carbon: soldatc, phondatc, liqdatc, quandatc * Oxygen: soldato, phondato, liqdato, quandato * Iron: soldatf, phondatf, liqdatf, quandatf * common/soldath/ snah(5), snbh(3), sneh(5), snfh(3), snih(5), 1 snjh(3), snph(5), snqh(3), alphasnh(4), 2 betasnh(4), snch, sndh, sngh, snhh, snkh, 3 snlh, snrh, snsh common/phondath/ pnah(5), pnbh(3), pnih(5), pnjh(3), 1 alphapnh(4), betapnh(4), pnch, pndh, pnkh, pnlh common/liqdath/ lnah(6), lnbh(4), lneh(6), lnfh(4), lnih(6), 1 lnjh(4), lnph(6), lnqh(4), alphalnh(4), 2 betalnh(4), lnch, lndh, lngh, lnhh, lnkh, 3 lnlh, lnrh, lnsh common/quandath/ qnah(6), qnbh(4), qneh(6), qnfh(4), 1 alphaqnh(5), qnch, qndh, qngh, qnhh common/soldatc/ snac(5), snbc(3), snec(5), snfc(3), snic(5), 1 snjc(3), snpc(5), snqc(3), alphasnc(4), 2 betasnc(4), sncc, sndc, sngc, snhc, snkc, 3 snlc, snrc, snsc common/phondatc/ pnac(5), pnbc(3), pnic(5), pnjc(3), 1 alphapnc(4), betapnc(4), pncc, pndc, pnkc, pnlc common/liqdatc/ lnac(6), lnbc(4), lnec(6), lnfc(4), lnic(6), 1 lnjc(4), lnpc(6), lnqc(4), alphalnc(4), 2 betalnc(4), lncc, lndc, lngc, lnhc, lnkc, 3 lnlc, lnrc, lnsc common/quandatc/ qnac(6), qnbc(4), qnec(6), qnfc(4), 1 alphaqnc(5), qncc, qndc, qngc, qnhc common/soldato/ snao(5), snbo(3), sneo(5), snfo(3), snio(5), 1 snjo(3), snpo(5), snqo(3), alphasno(4), 2 betasno(4), snco, sndo, sngo, snho, snko, 3 snlo, snro, snso common/phondato/ pnao(5), pnbo(3), pnio(5), pnjo(3), 1 alphapno(4), betapno(4), pnco, pndo, pnko, pnlo common/liqdato/ lnao(6), lnbo(4), lneo(6), lnfo(4), lnio(6), 1 lnjo(4), lnpo(6), lnqo(4), alphalno(4), 2 betalno(4), lnco, lndo, lngo, lnho, lnko, 3 lnlo, lnro, lnso common/quandato/ qnao(6), qnbo(4), qneo(6), qnfo(4), 1 alphaqno(5), qnco, qndo, qngo, qnho common/sfodatf/ snaf(5), snbf(3), snef(5), snff(3), snif(5), 1 snjf(3), snpf(5), snqf(3), alphasnf(4), 2 betasnf(4), sncf, sndf, sngf, snhf, snkf, 3 snlf, snrf, snsf common/phondatf/ pnaf(5), pnbf(3), pnif(5), pnjf(3), 1 alphapnf(4), betapnf(4), pncf, pndf, pnkf, pnlf common/liqdatf/ lnaf(6), lnbf(4), lnef(6), lnff(4), lnif(6), 1 lnjf(4), lnpf(6), lnqf(4), alphalnf(4), 2 betalnf(4), lncf, lndf, lngf, lnhf, lnkf, 3 lnlf, lnrf, lnsf common/quandatf/ qnaf(6), qnbf(4), qnef(6), qnff(4), 1 alphaqnf(5), qncf, qndf, qngf, qnhf data an0,bk,pi/6.022169e23,1.380622e-16,3.1415926535e0/ * * variables used in condlhy * data alhy / 0.34542, -.28157, 0.09184, -.03734 / data blhy / -.61919, 0.40004, -.16585 / data clhy / 0.35742, -.41151, 0.21552 / data dlhy / 0.21512, -.10843, -.00596, -.00950 / data elhy / -.36667, 0.14040, -.04588 / data flhy / 0.10493, -.09537, 0.04682 / * * variables used in condlhe * data alhe / 0.62199, -.16110, 0.15574, -.02893 / data blhe / -.65222, 0.48601, -.18266 / data clhe / 0.36580, -.52176, 0.26240 / data dlhe / 0.36090, 0.02576, 0.05061, 0.00015 / data elhe / -.40559, 0.15316, -.04058 / data flhe / 0.12140, -.11621, 0.04939 / * * variables used in condshe * NOTE: the values stored in a(1), e(1), i(1) and p(1) are 1/2 the * values listed in the paper. this is because they are used in the * formulae in the paper as a(1)/2, etc. (same for carbon and oxygen). * data sahe / -.01006, 0.05975, -.04290, -.00341, -.00338 / data sbhe / 0.11689, 0.01560, -.00713 / data sche, sdhe / -.16902, 1.4525 / data sehe / -.05144, 0.10844, -.04380, -.00663, -.00657 / data sfhe / 0.12937, 0.02888, -.00413 / data sghe, shhe / -.20406, 1.36576 / data sihe / 0.12429, -.04498, -.06418, -.00931, -.00582 / data sjhe / 0.11189, 0.00562, -.01176 / data skhe, slhe / -.22355, 1.70151 / data sphe / 0.09218, -.00644, -.06552, -.01190, -.00832 / data sqhe / 0.12575, 0.01598, -.00950 / data srhe, sshe / -.25532, 1.90219 / data alphahe / 1.5611, -20.6104, 247.0830, -1003.2600 / data betahe / 1.5874, -21.4772, 255.7040, -1028.9000 / * * variables used in condlc * data alc/ 0.9896, -0.1851, 0.1019, -0.0360/ data blc/ -0.8825, 0.6657, -0.3798/ data clc/ -0.0915, -1.5848, 1.1882/ data dlc/ 0.4406, -0.0161, -0.0093, -0.0028/ data elc/ -0.4821, 0.0826, -0.0557/ data flc/ -0.5193, -0.0830, 0.0147/ * * variables used in condsc * data sac / 0.03387, 0.01795, -0.04324, -0.00437, -0.00422/ data sbc / 0.04865, -0.00473, -0.00865/ data scc, sdc/-0.09378, 0.94180/ data sec / -0.00329, 0.08497, -0.06909, -0.00582, -0.00677/ data sfc / 0.12695, 0.00854, -0.01197/ data sgc, shc/-0.19961, 1.63031/ data sic / 0.177955,-0.08588, -0.07995, -0.00957, -0.00256/ data sjc / 0.03417, -0.03395, -0.02519/ data skc, slc/-0.10127, 2.10721/ data spc / 0.18606, -0.06547, -0.10440, -0.01172, -0.00448/ data sqc / 0.10304, -0.02601, -0.02922/ data src, ssc/-0.20326, 2.80479/ data alphac/0.4786, 16.4194, -139.1920, 184.7580/ data betac/ 0.6146, 11.7296, -89.8319, 32.0399/ * * variables used in condlo * data alo/ 1.0779, -0.1838, 0.1059, -.0290/ data blo/ -.9743, 0.6955, -.3966/ data clo/ -.1040, -1.7692, 1.3546 / data dlo/ 0.4486, -0.0160, -.0014, 0.0039 / data elo/ -.5193, 0.0822, -.0467 / data flo/ -.5403, -0.1022, 0.0416 / * * variables used in condso * data sao / 0.02312, 0.02087, -.03684, -.00366, -.00348 / data sbo / 0.04177, -.00712, -.00865 / data sco, sdo/ -.07178, 0.85073 / data seo / -.001185, 0.08552, -.07206, -.00554, -.00673 / data sfo / 0.12586, 0.00371, -.01360 / data sgo, sho/ -.19226, 1.66185 / data sio / .160905, -.08292, -.07377, -.00577, 0.00154/ data sjo / 0.03816, -.04206, -.02870 / data sko, slo/ -.06064, 2.14532 / data spo / .19600, -.07598, -.10917, -.00942, -.00142 / data sqo / 0.10403, -.03774, -.03483 / data sro, sso/ -.17568, 2.99352 / data alphao/ .4682, 18.0197, -171.0090, 313.8890 / data betao/ 0.6085, 12.9705, -115.3180, 136.3590 / * * variables used in bremsolh * data snah / -.01148, 0.01601, -.00433, 0.00015, -.00034 / data snbh / 0.01558, 0.00191, -.00055 / data snch, sndh / -.01694, 0.10649 / data sneh / -.01827, 0.02395, -.00448, -.00033, -.00088 / data snfh / 0.01730, 0.00402, -.00005 / data sngh, snhh / -.02222, 0.13969 / data snih / -.00323, 0.00440, -.00110, 0.00001, -.00007 / data snjh / 0.00294, 0.00059, -.00018 / data snkh, snlh / -.00337, 0.02116 / data snph / -.00469, 0.00610, -.00114, -.00010, -.00018 / data snqh / 0.00320, 0.00107, -.00008 / data snrh, snsh / -.00442, 0.02775 / data alphasnh / 1.6449, -23.2588, 272.1670, -1074.7000 / data betasnh / 1.6443, -23.2414, 272.0080, -1074.2500 / data pnah / -.00687, 0.00957, -.00204, -.00005, -.00003 / data pnbh / 0.00661, 0.00135, -.00035 / data pnch, pndh / -.00811, 0.05098 / data pnih / -.00169, 0.00231, -.00047, -.00003, 0.00000 / data pnjh / 0.00111, 0.00042, -.00010 / data pnkh, pnlh / -.00161, 0.01013 / data alphapnh / -.1394, 7.0680, -115.5940, 619.9170 / data betapnh / -.1394, 7.0664, -115.5800, 619.8790 / * * variables used in bremsolc * data snac / 0.01838, -.01066, -.00458, -.00177, -.00138 / data snbc / -.00244, -.00206, -.00037 / data sncc, sndc / -.01093, 0.12431 / data snec / 0.02360, -.01353, -.00619, -.00211, -.00176 / data snfc / 0.00456, -.00174, -.00031 / data sngc, snhc / -.02259, 0.20343 / data snic / 0.00053, -.00048, -.00022, 0.00019, -.00001 / data snjc / 0.00658, -.00180, 0.00036 / data snkc, snlc / -.00398, 0.02499 / data snpc / -.00024, 0.00063, -.00064, 0.00030, -.00006 / data snqc / 0.01013, -.00247, 0.00052 / data snrc, snsc / -.00650, 0.04087 / data alphasnc / 0.6252, 10.6819, -70.6879, -44.3349 / data betasnc / 0.6307, 10.4966, -68.7973, -50.0581 / data pnac / 0.01116, -.00589, -.00279, -.00073, -.00043 / data pnbc / -.00095, -.00059, 0.00002 / data pncc, pndc / -.00729, 0.06630 / data pnic / 0.00012, 0.00018, -.00028, 0.00012, -.00004 / data pnjc / 0.00339, -.00082, 0.00015 / data pnkc, pnlc / -.00212, 0.01332 / data alphapnc / 0.5481, -20.4731, 223.9220, -534.9400 / data betapnc / 0.5413, -20.2069, 220.7060, -524.1240 / * * variables used in bremsolo * data snao / 0.01616, -.00874, -.00413, -.00190, -.00139 / data snbo / -.00344, -.00261, -.00070 / data snco, sndo / -.00791, 0.13980 / data sneo / 0.02210, -.00883, -.00857, -.00257, -.00214 / data snfo / 0.00629, -.00210, -.00099 / data sngo, snho / -.02610, 0.26993 / data snio / 0.00100, -.00112, -.00003, 0.00014, 0.00001 / data snjo / 0.00745, -.00209, 0.00044 / data snko, snlo / -.00447, 0.02811 / data snpo / -.00056, 0.00110, -.00074, 0.00024, -.00004 / data snqo / 0.01286, -.00281, 0.00057 / data snro, snso / -.00861, 0.05414 / data alphasno / 0.4889, 16.1962, -138.4860, 185.7060 / data betasno / 0.5111, 15.4195, -130.1540, 159.6050 / data pnao / 0.00800, -.00191, -.00330, -.00075, -.00047 / data pnbo / 0.00088, -.00098, -.00036 / data pnco, pndo / -.00776, 0.08995 / data pnio / -.00009, 0.00055, -.00038, 0.00011, -.00003 / data pnjo / 0.00429, -.00088, 0.00014 / data pnko, pnlo / -.00287, 0.01803 / data alphapno / 0.3173, -14.4048, 186.9100, -476.8100 / data betapno / 0.3073, -13.7973, 176.9940, -438.7520 / * * variables used in bremsolf * data snaf / 0.02096, -.01768, -.00007, -.00241, -.00080 / data snbf / -.01705, -.00268, -.00141 / data sncf, sndf / 0.00818, 0.13629 / data snef / 0.04628, -.03290, -.00523, -.00539, -.00276 / data snff / -.02574, -.00630, -.00285 / data sngf, snhf / 0.00022, 0.31871 / data snif / 0.00489, -.00653, 0.00171, -.00024, 0.00017 / data snjf / 0.00869, -.00323, 0.00075 / data snkf, snlf / -.00439, 0.02766 / data snpf / 0.00732, -.00957, 0.00222, -.00022, 0.00025 / data snqf / 0.01867, -.00615, 0.00133 / data snrf, snsf / -.01023, 0.06442 / data alphasnf / 0.6798, 12.7527, -140.1800, 268.8290 / data betasnf / 0.7783, 10.2315, -124.2640, 241.3060 / data pnaf / 0.01292, -.00894, -.00097, -.00125, -.00045 / data pnbf / -.00739, -.00190, -.00062 / data pncf, pndf / 0.00167, 0.09950 / data pnif / 0.00240, -.00271, 0.00047, 0.00001, -.00001 / data pnjf / 0.00604, -.00197, 0.00044 / data pnkf, pnlf / -.00320, 0.02012 / data alphapnf / 0.2847, -13.0828, 192.5030, -543.1480 / data betapnf / 0.3221, -13.8640, 201.0700, -573.1160 / * * variables used in bremliqh * data lnah / 0.04518, -.03009, -.00564, -.00544, -.00290, 1 -.00224 / data lnbh / -.02148, -.00817, -.00300, -.00170 / data lnch, lndh / 0.00671, 0.28130 / data lneh / -.01003, 0.01790, -.00783, -.00021, 0.00024, 1 -.00014 / data lnfh / 0.00538, -.00175, -.00346, -.00031 / data lngh, lnhh / -.02199, 0.17300 / data lnih / 0.00096, -.00301, -.00073, 0.00182, 0.00037, 1 0.00116 / data lnjh / 0.01706, -.00753, 0.00066, -.00060 / data lnkh, lnlh / -.01021, 0.06417 / data lnph / -.00556, 0.00603, -.00149, 0.00047, 0.00040, 1 0.00028 / data lnqh / 0.00422, -.00009, -.00066, -.00003 / data lnrh, lnsh / -.00561, 0.03522 / data alphalnh / -.07980, 0.17057, 1.51980, -.61058 / data betalnh / -.05881, 0.00165, 1.82700, -.76993 / data qnah / -.02260, 0.02634, -.00424, 0.00055, 0.00017, 1 -.00042 / data qnbh / 0.00550, 0.00105, -.00179, 0.00028 / data qnch, qndh / -.02408, 0.28763 / data qneh / -.00539, 0.00617, -.00205, 0.00089, 0.00022, 1 0.00032 / data qnfh / 0.00615, -.00149, -.00055, -.00004 / data qngh, qnhh / -.00764, 0.53346 / data alphaqnh / 2.5344, -7.7188, -9.5605, -5.6433, 1.2671 / * * variables used in bremliqc * data lnac / 0.08973, -.05821, -.01089, -.01147, -.00656, 1 -.00519 / data lnbc / -.04969, -.01584, -.00504, -.00281 / data lncc, lndc / 0.00945, 0.34529 / data lnec / 0.03390, -.00944, -.01289, -.00589, -.00404, 1 -.00330 / data lnfc / -.02213, -.01136, -.00467, -.00131 / data lngc, lnhc / -.02342, 0.24819 / data lnic / 0.00383, -.00710, -.00028, 0.00232, 0.00044, 1 0.00158 / data lnjc / 0.02300, -.01078, 0.00118, -.00089 / data lnkc, lnlc / -.01259, 0.07917 / data lnpc / -.00384, 0.00356, -.00184, 0.00146, 0.00031, 1 0.00069 / data lnqc / 0.01052, -.00354, -.00014, -.00018 / data lnrc, lnsc / -.00829, 0.05211 / data alphalnc / -.05483, -.01946, 1.86310, -.78873 / data betalnc / -.06711, 0.06859, 1.74360, -.74498 / data qnac / -.01900, 0.02446, -.00419, -.00030, -.00047, 1 -.00099 / data qnbc / 0.00093, 0.00060, -.00176, 0.00018 / data qncc, qndc / -.02458, 0.25843 / data qnec / -.00690, 0.01046, -.00351, 0.00010, -.00003, 1 -.00026 / data qnfc / 0.00289, -.00082, -.00155, -.00008 / data qngc, qnhc / -.01194, 0.34437 / data alphaqnc / 0.9859, 2.2468, -10.6808, 11.4582, -4.0100 / * * variables used in bremliqo * data lnao / 0.10466, -.06740, -.01293, -.01352, -.00776, 1 -.00613 / data lnbo / -.05950, -.01837, -.00567, -.00310 / data lnco, lndo / 0.00952, 0.36029 / data lneo / 0.04652, -.01656, -.01489, -.00778, -.00520, 1 -.00418 / data lnfo / -.03076, -.01390, -.00522, -.00161 / data lngo, lnho / -.02513, 0.27480 / data lnio / 0.00476, -.00838, -.00011, 0.00244, 0.00046, 1 0.00168 / data lnjo / 0.02455, -.01167, 0.00132, -.00097 / data lnko, lnlo / -.01314, 0.08263 / data lnpo / -.00350, 0.00295, -.00184, 0.00166, 0.00032, 1 0.00082 / data lnqo / 0.01231, -.00445, 0.00002, -.00026 / data lnro, lnso / -.00921, 0.05786 / data alphalno / -.06597, 0.06048, 1.74860, -.74308 / data betalno / -.07356, 0.10865, 1.70150, -.73653 / data qnao / -.01792, 0.02395, -.00425, -.00055, -.00065, 1 -.00115 / data qnbo / -.00039, 0.00043, -.00177, 0.00015 / data qnco, qndo / -.02491, 0.25207 / data qneo / -.00866, 0.01344, -.00448, -.00002, -.00009, 1 -.00041 / data qnfo / 0.00291, -.00091, -.00200, -.00010 / data qngo, qnho / -.01541, 0.32800 / data alphaqno / 0.8978, 2.5443, -10.3549, 10.2665, -3.3537 / * * variables used in bremliqf * data lnaf / 0.17444, -.11076, -.02349, -.02283, -.01250, 1 -.00971 / data lnbf / -.10661, -.02860, -.00785, -.00385 / data lncf, lndf / 0.00766, 0.40991 / data lnef / 0.11580, -.05891, -.02531, -.01747, -.01021, 1 -.00778 / data lnff / -.07767, -.02516, -.00711, -.00260 / data lngf, lnhf / -.03076, 0.36908 / data lnif / 0.00946, -.01493, -.00125, 0.00262, 0.00055, 1 0.00209 / data lnjf / 0.03034, -.01519, 0.00204, -.00135 / data lnkf, lnlf / -.01494, 0.09395 / data lnpf / -.00006, -.00222, -.00104, 0.00225, 0.00037, 1 0.00140 / data lnqf / 0.02013, -.00883, 0.00090, -.00071 / data lnrf, lnsf / -.01249, 0.07850 / data alphalnf / -.07608, 0.11559, 1.75730, -.79677 / data betalnf / -.08034, 0.14368, 1.73140, -.79467 / data qnaf / -.01326, 0.02179, -.00484, -.00146, -.00135, 1 -.00175 / data qnbf / -.00567, -.00040, -.00174, 0.00012 / data qncf, qndf / -.02710, 0.23276 / data qnef / -.00430, 0.01443, -.00677, -.00142, -.00119, 1 -.00151 / data qnff / -.00396, -.00245, -.00278, -.00023 / data qngf, qnhf / -.02288, 0.26787 / data alphaqnf / 1.0379, 1.0195, -5.1242, 3.9369, -0.8701 / end ************************************************************************