SUBROUTINE COLH(ID,T,COL) C ========================= C C Hydrogen collision rates C C All standard expressions are taken from Mihalas, Heasley, and C Auer, NCAR-TN/STR-104 (1975) C C New expressions (also from Mihalas) for collisional ionization C for first 10 levels taken from Klaus Werner. C C New standard expressions from Giovanardi et al. (1987, AAS, 70, 29) C for collisional excitation (valid from 3000K to 500000K) C C Meaning of ICOL: C a) for ionization - .ge.0 - standard expression C < 0 - non-standard, user suplied formula C b) for 1 - 2 transition =0 - standard theoretical formula C = 1 - experimental fit (formula quted in C Mihalas et al.) C = 2 - formula by Crandall et al (procedure C CEH12) C c) for all other line transitions C .ge.0 - standard expression C < 0 - non-standard, user supplie formula C INCLUDE 'IMPLIC.FOR' INCLUDE 'BASICS.FOR' INCLUDE 'ATOMIC.FOR' INCLUDE 'MODELQ.FOR' PARAMETER (CC0 = 5.465D-11, * CEX1 = -30.20581, * CEX2 = 3.8608704, * CEX3 = 305.63574, * CI1 = 0.3, * CI2 = 0.435, * CA1 = 5.444416D7, * CA2 = -2.8185937D4, * CA3 = 19.987261, * CA4 = -5.8906298D-5, * CB1 = 1.3935312D3, * CB2 = -1.6805859D2, * CB3 = -2.539D3, * CC1 = 2.0684609D3, * CC2 = -3.341582D2, * CC3 = -7.6440625D3, * CD1 = 3.2174844D3, * CD2 = -5.5882422D2, * CD3 = -6.86325D3, * CE1 = 5.759125D3, * CE2 = 81.75, * CE3 = -1.5163D3, * CF1 = 1.461475D4, * CF2 = 393.4, * CF3 = -4.8284D3, * ALF0 = 1.8, * ALF1 = 0.4, * BET0 = 3.0, * BET1 = 1.2, * O148 = 0.148, * CHMI = 5.59D-15) PARAMETER (EXPIA1=-0.57721566,EXPIA2=0.99999193, * EXPIA3=-0.24991055,EXPIA4=0.05519968, * EXPIA5=-0.00976004,EXPIA6=0.00107857, * EXPIB1=0.2677734343,EXPIB2=8.6347608925, * EXPIB3=18.059016973,EXPIB4=8.5733287401, * EXPIC1=3.9584969228,EXPIC2=21.0996530827, * EXPIC3=25.6329561486,EXPIC4=9.5733223454) DIMENSION COL(MTRANS),A(6,10) DIMENSION CCOOL(4,14,15),CHOT(4,14,15),XTT(4) DATA ((A(I,J),J=1,10),I=1,6) / * -86.7633398, 2632.8369 , 7478.9556 ,-4202.8442 ,-47995.930 , * -120942.89 ,-202300.81 ,-261373.03 ,-266337.91 ,-192293.20 , * 100.919188 ,-2738.7485 ,-8495.4590 , 1937.3763 , 45825.371 , * 122209.39 , 211928.67 , 285044.75 , 309455.47 , 258802.22 , * -45.7813807, 1121.3976 , 3794.6826 , 340.35764 ,-16617.055 , * -47390.313 ,-84973.688 ,-117833.95 ,-133243.61 ,-120363.95 , * 10.1978559 ,-224.30670 ,-822.83636 ,-290.10489 , 2905.7393 , * 8944.6025 , 16556.992 , 23544.543 , 27419.742 , 26002.143 , * -1.11223557, 21.923729 , 86.619110 , 48.840523 ,-246.99014 , * -828.41028 ,-1581.2722 ,-2297.9321 ,-2738.1743 ,-2686.4087 , * .0474198818,-.83974838,-3.5534720 ,-2.6097214 , 8.1972208 , * 30.267115 , 59.521984 , 88.178680 , 107.05288 , 107.73775 / C DATA ((CCOOL(I, 1, K),I=1,4),K=1,15)/ 4*0., & 5.742D-01, 1.818D-05,-1.093D-10, 8.687D-16, & 1.934D-01,-4.698D-07, 8.352D-11,-5.576D-16, & 6.323D-03, 2.237D-06,-1.620D-11, 8.955D-17, & 2.035D-02, 6.076D-07,-2.175D-13,-2.495D-18, & 1.136D-02, 3.428D-07,-1.467D-13,-1.300D-18, & 6.999D-03, 2.126D-07,-9.963D-14,-7.672D-19, & 4.624D-03, 1.410D-07,-6.969D-14,-4.927D-19, & 3.217D-03, 9.836D-08,-5.031D-14,-3.361D-19, & 2.329D-03, 7.135D-08,-3.737D-14,-2.400D-19, & 1.741D-03, 5.342D-08,-2.845D-14,-1.775D-19, & 1.336D-03, 4.103D-08,-2.213D-14,-1.351D-19, & 1.048D-03, 3.220D-08,-1.754D-14,-1.053D-19, & 8.369D-04, 2.574D-08,-1.413D-14,-8.368D-20, & 6.791D-04, 2.090D-08,-1.154D-14,-6.763D-20/ DATA ((CCOOL(I, 2, K),I=1,4),K=1,15)/ 8*0., & 2.253D+01, 9.350D-04, 1.215D-08,-9.969D-14, & 7.816D-01, 5.414D-04,-1.827D-09, 5.140D-17, & 1.459D+00, 2.858D-04,-2.207D-09, 9.028D-15, & 7.172D-01, 1.440D-04,-1.139D-09, 4.755D-15, & 4.107D-01, 8.360D-05,-6.699D-10, 2.823D-15, & 2.591D-01, 5.319D-05,-4.293D-10, 1.819D-15, & 1.747D-01, 3.608D-05,-2.925D-10, 1.243D-15, & 1.237D-01, 2.567D-05,-2.087D-10, 8.891D-16, & 9.097D-02, 1.893D-05,-1.539D-10, 6.585D-16, & 6.896D-02, 1.438D-05,-1.174D-10, 5.017D-16, & 5.356D-02, 1.119D-05,-9.150D-11, 3.913D-16, & 4.247D-02, 8.887D-06,-7.272D-11, 3.112D-16, & 3.425D-02, 7.176D-06,-5.877D-11, 2.516D-16/ DATA ((CCOOL(I, 3, K),I=1,4),K=1,15)/ 12*0., & -1.290D+01, 2.059D-02, 5.461D-08,-9.082D-13, & 3.562D+02, 7.337D-03,-9.622D-08, 5.596D-13, & 5.744D+00, 3.570D-03,-3.259D-08, 1.452D-13, & 2.968D+00, 1.813D-03,-1.703D-08, 7.744D-14, & 1.756D+00, 1.065D-03,-1.016D-08, 4.667D-14, & 1.135D+00, 6.865D-04,-6.601D-09, 3.053D-14, & 7.802D-01, 4.713D-04,-4.558D-09, 2.116D-14, & 5.615D-01, 3.390D-04,-3.292D-09, 1.532D-14, & 4.189D-01, 2.528D-04,-2.461D-09, 1.148D-14, & 3.213D-01, 1.939D-04,-1.891D-09, 8.833D-15, & 2.523D-01, 1.522D-04,-1.487D-09, 6.953D-15, & 2.018D-01, 1.218D-04,-1.192D-09, 5.576D-15/ DATA ((CCOOL(I, 4, K),I=1,4),K=1,15)/ 16*0., & 4.139D+03, 4.645D-01,-7.097D-06, 4.388D-11, & 1.794D+03, 4.443D-02,-6.484D-07, 3.936D-12, & 1.536D+01, 2.042D-02,-2.065D-07, 9.734D-13, & 8.730D+00, 1.033D-02,-1.074D-07, 5.161D-13, & 5.434D+00, 6.084D-03,-6.423D-08, 3.116D-13, & 3.628D+00, 3.938D-03,-4.196D-08, 2.048D-13, & 2.554D+00, 2.718D-03,-2.914D-08, 1.428D-13, & 1.873D+00, 1.967D-03,-2.119D-08, 1.041D-13, & 1.418D+00, 1.476D-03,-1.594D-08, 7.843D-14, & 1.102D+00, 1.138D-03,-1.232D-08, 6.075D-14, & 8.744D-01, 8.987D-04,-9.746D-09, 4.809D-14/ DATA ((CCOOL(I, 5, K),I=1,4),K=1,15)/ 20*0., & -9.122D+02, 1.260D+00,-1.070D-05, 4.290D-11, & 3.959D+01, 2.108D-01,-2.162D-06, 1.020D-11, & 3.691D+01, 7.806D-02,-8.485D-07, 4.166D-12, & 2.352D+01, 3.911D-02,-4.365D-07, 2.179D-12, & 1.542D+01, 2.296D-02,-2.601D-07, 1.310D-12, & 1.062D+01, 1.487D-02,-1.699D-07, 8.608D-13, & 7.642D+00, 1.029D-02,-1.183D-07, 6.014D-13, & 5.695D+00, 7.464D-03,-8.621D-08, 4.394D-13, & 4.368D+00, 5.617D-03,-6.508D-08, 3.323D-13, & 3.430D+00, 4.348D-03,-5.051D-08, 2.583D-13/ DATA ((CCOOL(I, 6, K),I=1,4),K=1,15)/ 24*0., & -3.431D+03, 4.116D+00,-3.853D-05, 1.679D-10, & 4.397D+01, 6.434D-01,-7.008D-06, 3.431D-11, & 8.927D+01, 2.325D-01,-2.667D-06, 1.350D-11, & 6.153D+01, 1.152D-01,-1.354D-06, 6.957D-12, & 4.165D+01, 6.729D-02,-8.024D-07, 4.156D-12, & 2.923D+01, 4.349D-02,-5.232D-07, 2.724D-12, & 2.130D+01, 3.008D-02,-3.641D-07, 1.902D-12, & 1.603D+01, 2.185D-02,-2.656D-07, 1.391D-12, & 1.239D+01, 1.647D-02,-2.008D-07, 1.054D-12/ DATA ((CCOOL(I, 7, K),I=1,4),K=1,15)/ 28*0., & -9.280D+03, 1.116D+01,-1.122D-04, 5.167D-10, & 6.658D+01, 1.651D+00,-1.884D-05, 9.487D-11, & 2.172D+02, 5.833D-01,-6.977D-06, 3.615D-11, & 1.535D+02, 2.858D-01,-3.499D-06, 1.838D-11, & 1.049D+02, 1.660D-01,-2.060D-06, 1.090D-11, & 7.412D+01, 1.070D-01,-1.339D-06, 7.118D-12, & 5.428D+01, 7.389D-02,-9.304D-07, 4.963D-12, & 4.103D+01, 5.366D-02,-6.786D-07, 3.629D-12/ DATA ((CCOOL(I, 8, K),I=1,4),K=1,15)/ 32*0., & -2.069D+04, 2.637D+01,-2.802D-04, 1.342D-09, & 2.055D+02, 3.731D+00,-4.420D-05, 2.276D-10, & 5.123D+02, 1.292D+00,-1.599D-05, 8.442D-11, & 3.578D+02, 6.265D-01,-7.922D-06, 4.235D-11, & 2.438D+02, 3.616D-01,-4.633D-06, 2.494D-11, & 1.721D+02, 2.322D-01,-3.000D-06, 1.622D-11, & 1.260D+02, 1.601D-01,-2.081D-06, 1.129D-11/ DATA ((CCOOL(I, 9, K),I=1,4),K=1,15)/ 36*0., & -4.032D+04, 5.614D+01,-6.231D-04, 3.073D-09, & 6.989D+02, 7.655D+00,-9.352D-05, 4.903D-10, & 1.141D+03, 2.605D+00,-3.313D-05, 1.777D-10, & 7.755D+02, 1.250D+00,-1.624D-05, 8.808D-11, & 5.234D+02, 7.175D-01,-9.437D-06, 5.153D-11, & 3.677D+02, 4.590D-01,-6.087D-06, 3.338D-11/ DATA ((CCOOL(I,10, K),I=1,4),K=1,15)/ 40*0., & -7.097D+04, 1.101D+02,-1.266D-03, 6.390D-09, & 2.018D+03, 1.455D+01,-1.824D-04, 9.708D-10, & 2.383D+03, 4.875D+00,-6.348D-05, 3.449D-10, & 1.569D+03, 2.319D+00,-3.081D-05, 1.691D-10, & 1.046D+03, 1.323D+00,-1.779D-05, 9.830D-11/ DATA ((CCOOL(I,11, K),I=1,4),K=1,15)/ 44*0., & -1.150D+05, 2.020D+02,-2.392D-03, 1.231D-08, & 4.988D+03, 2.601D+01,-3.334D-04, 1.797D-09, & 4.675D+03, 8.595D+00,-1.142D-04, 6.273D-10, & 2.986D+03, 4.054D+00,-5.491D-05, 3.046D-10/ DATA ((CCOOL(I,12, K),I=1,4),K=1,15)/ 48*0., & -1.737D+05, 3.511D+02,-4.263D-03, 2.227D-08, & 1.094D+04, 4.419D+01,-5.774D-04, 3.146D-09, & 8.673D+03, 1.442D+01,-1.950D-04, 1.082D-09/ DATA ((CCOOL(I,13, K),I=1,4),K=1,15)/ 52*0., & -2.459D+05, 5.829D+02,-7.233D-03, 3.830D-08, & 2.191D+04, 7.194D+01,-9.561D-04, 5.259D-09/ DATA ((CCOOL(I,14, K),I=1,4),K=1,15)/ 56*0., & -3.273D+05, 9.312D+02,-1.178D-02, 6.306D-08/ DATA ((CHOT(I, 1, K),I=1,4),K=1,15)/ 4*0., & 5.856D-01, 1.551D-05,-9.669D-12, 5.716D-19, & 1.537D-01, 3.548D-06,-3.224D-12, 7.626D-19, & 2.400D-02, 1.419D-06,-2.008D-12, 1.356D-18, & 2.002D-02, 6.325D-07,-7.070D-13, 4.096D-19, & 1.123D-02, 3.549D-07,-3.998D-13, 2.331D-19, & 6.940D-03, 2.194D-07,-2.483D-13, 1.453D-19, & 4.593D-03, 1.453D-07,-1.648D-13, 9.667D-20, & 3.199D-03, 1.012D-07,-1.150D-13, 6.758D-20, & 2.318D-03, 7.334D-08,-8.349D-14, 4.910D-20, & 1.727D-03, 5.493D-08,-6.270D-14, 3.695D-20, & 1.326D-03, 4.218D-08,-4.821D-14, 2.844D-20, & 1.040D-03, 3.310D-08,-3.786D-14, 2.236D-20, & 8.305D-04, 2.645D-08,-3.028D-14, 1.790D-20, & 6.740D-04, 2.147D-08,-2.460D-14, 1.455D-20/ DATA ((CHOT(I, 2, K),I=1,4),K=1,15)/ 8*0., & 1.710D+01, 1.530D-03,-2.553D-09, 1.924D-15, & 8.237D+00, 3.554D-04,-7.566D-10, 6.420D-16, & 5.932D+00, 1.301D-04,-2.912D-10, 2.535D-16, & 2.987D+00, 6.419D-05,-1.444D-10, 1.260D-16, & 1.733D+00, 3.689D-05,-8.324D-11, 7.267D-17, & 1.102D+00, 2.334D-05,-5.273D-11, 4.605D-17, & 7.472D-01, 1.576D-05,-3.567D-11, 3.116D-17, & 5.312D-01, 1.118D-05,-2.532D-11, 2.212D-17, & 3.919D-01, 8.232D-06,-1.865D-11, 1.630D-17, & 2.977D-01, 6.245D-06,-1.416D-11, 1.237D-17, & 2.315D-01, 4.855D-06,-1.101D-11, 9.622D-18, & 1.838D-01, 3.851D-06,-8.734D-12, 7.635D-18, & 1.484D-01, 3.108D-06,-7.050D-12, 6.164D-18/ DATA ((CHOT(I, 3, K),I=1,4),K=1,15)/ 12*0., & 1.940D+02, 1.949D-02,-3.832D-08, 3.137D-14, & 4.729D+02, 1.927D-03,-4.171D-09, 3.628D-15, & 6.741D+01, 1.315D-03,-3.145D-09, 2.814D-15, & 3.444D+01, 6.477D-04,-1.560D-09, 1.399D-15, & 2.031D+01, 3.744D-04,-9.054D-10, 8.130D-16, & 1.311D+01, 2.388D-04,-5.789D-10, 5.203D-16, & 9.007D+00, 1.629D-04,-3.955D-10, 3.556D-16, & 6.484D+00, 1.166D-04,-2.835D-10, 2.550D-16, & 4.837D+00, 8.666D-05,-2.108D-10, 1.896D-16, & 3.711D+00, 6.631D-05,-1.614D-10, 1.452D-16, & 2.914D+00, 5.194D-05,-1.265D-10, 1.138D-16, & 2.332D+00, 4.150D-05,-1.011D-10, 9.100D-17/ DATA ((CHOT(I, 4, K),I=1,4),K=1,15)/ 16*0., & 7.204D+03, 1.627D-01,-5.181D-07, 5.605D-13, & 2.507D+03, 9.370D-03,-2.091D-08, 1.842D-14, & 3.823D+02, 6.480D-03,-1.600D-08, 1.452D-14, & 1.950D+02, 3.161D-03,-7.869D-09, 7.157D-15, & 1.154D+02, 1.823D-03,-4.561D-09, 4.154D-15, & 7.486D+01, 1.165D-03,-2.924D-09, 2.665D-15, & 5.178D+01, 7.977D-04,-2.006D-09, 1.830D-15, & 3.752D+01, 5.737D-04,-1.444D-09, 1.318D-15, & 2.816D+01, 4.283D-04,-1.080D-09, 9.858D-16, & 2.174D+01, 3.293D-04,-8.307D-10, 7.587D-16, & 1.717D+01, 2.592D-04,-6.544D-10, 5.978D-16/ DATA ((CHOT(I, 5, K),I=1,4),K=1,15)/ 20*0., & 2.166D+04, 4.690D-01,-1.122D-06, 1.008D-12, & 3.874D+03, 6.443D-02,-1.596D-07, 1.452D-13, & 1.465D+03, 2.207D-02,-5.556D-08, 5.082D-14, & 7.410D+02, 1.062D-02,-2.698D-08, 2.476D-14, & 4.374D+02, 6.096D-03,-1.556D-08, 1.430D-14, & 2.841D+02, 3.889D-03,-9.962D-09, 9.167D-15, & 1.969D+02, 2.663D-03,-6.838D-09, 6.297D-15, & 1.431D+02, 1.918D-03,-4.935D-09, 4.547D-15, & 1.078D+02, 1.436D-03,-3.698D-09, 3.409D-15, & 8.353D+01, 1.107D-03,-2.854D-09, 2.632D-15/ DATA ((CHOT(I, 6, K),I=1,4),K=1,15)/ 24*0., & 7.146D+04, 1.379D+00,-3.346D-06, 3.023D-12, & 1.187D+04, 1.794D-01,-4.501D-07, 4.118D-13, & 4.380D+03, 5.990D-02,-1.527D-07, 1.405D-13, & 2.192D+03, 2.846D-02,-7.324D-08, 6.759D-14, & 1.288D+03, 1.621D-02,-4.197D-08, 3.881D-14, & 8.351D+02, 1.031D-02,-2.678D-08, 2.480D-14, & 5.790D+02, 7.050D-03,-1.837D-08, 1.702D-14, & 4.213D+02, 5.079D-03,-1.326D-08, 1.229D-14, & 3.179D+02, 3.804D-03,-9.944D-09, 9.226D-15/ DATA ((CHOT(I, 7, K),I=1,4),K=1,15)/ 28*0., & 1.954D+05, 3.426D+00,-8.397D-06, 7.624D-12, & 3.057D+04, 4.266D-01,-1.080D-06, 9.917D-13, & 1.103D+04, 1.392D-01,-3.582D-07, 3.309D-13, & 5.458D+03, 6.530D-02,-1.696D-07, 1.572D-13, & 3.189D+03, 3.693D-02,-9.653D-08, 8.966D-14, & 2.062D+03, 2.338D-02,-6.136D-08, 5.707D-14, & 1.429D+03, 1.595D-02,-4.200D-08, 3.910D-14, & 1.039D+03, 1.148D-02,-3.029D-08, 2.282D-14/ DATA ((CHOT(I, 8, K),I=1,4),K=1,15)/ 32*0., & 4.651D+05, 7.527D+00,-1.859D-05, 1.694D-11, & 6.930D+04, 9.038D-01,-2.302D-06, 2.121D-12, & 2.450D+04, 2.891D-01,-7.487D-07, 6.939D-13, & 1.200D+04, 1.340D-01,-3.505D-07, 3.260D-13, & 6.970D+03, 7.523D-02,-1.981D-07, 1.846D-13, & 4.493D+03, 4.741D-02,-1.254D-07, 1.170D-13, & 3.106D+03, 3.226D-02,-8.559D-08, 7.997D-14/ DATA ((CHOT(I, 9, K),I=1,4),K=1,15)/ 36*0., & 9.956D+05, 1.506D+01,-3.741D-05, 3.418D-11, & 1.425D+05, 1.754D+00,-4.489D-06, 4.146D-12, & 4.949D+04, 5.510D-01,-1.435D-06, 1.333D-12, & 2.401D+04, 2.526D-01,-6.645D-07, 6.196D-13, & 1.386D+04, 1.408D-01,-3.729D-07, 3.485D-13, & 8.904D+03, 8.835D-02,-2.350D-07, 2.200D-13/ DATA ((CHOT(I,10, K),I=1,4),K=1,15)/ 40*0., & 1.961D+06, 2.798D+01,-6.982D-05, 6.394D-11, & 2.715D+05, 3.175D+00,-8.158D-06, 7.551D-12, & 9.279D+04, 9.821D-01,-2.567D-06, 2.390D-12, & 4.460D+04, 4.458D-01,-1.178D-06, 1.100D-12, & 2.561D+04, 2.468D-01,-6.566D-07, 6.150D-13/ DATA ((CHOT(I,11, K),I=1,4),K=1,15)/ 44*0., & 3.613D+06, 4.898D+01,-1.227D-04, 1.125D-10, & 4.861D+05, 5.434D+00,-1.401D-05, 1.299D-11, & 1.638D+05, 1.658D+00,-4.348D-06, 4.054D-12, & 7.810D+04, 7.456D-01,-1.976D-06, 1.850D-12/ DATA ((CHOT(I,12, K),I=1,4),K=1,15)/ 48*0., & 6.300D+06, 8.163D+01,-2.051D-04, 1.884D-10, & 8.271D+05, 8.881D+00,-2.296D-05, 2.131D-11, & 2.753D+05, 2.676D+00,-7.037D-06, 6.571D-12/ DATA ((CHOT(I,13, K),I=1,4),K=1,15)/ 52*0., & 1.049D+07, 1.305D+02,-3.288D-04, 3.025D-10, & 1.348D+06, 1.396D+01,-3.617D-05, 3.361D-11/ DATA ((CHOT(I,14, K),I=1,4),K=1,15)/ 56*0., & 1.680D+07, 2.016D+02,-5.089D-04, 4.687D-10/ C HKT=HK/T CT=CC0*SQRT(T) TK=HKT/H t0=t X=LOG10(T) X2=X*X X3=X*X2 X4=X2*X2 X5=X3*X2 XTT(1)=1. XTT(2)=T XTT(3)=T*T XTT(4)=T*T*T SQT=SQRT(T) N0HN=NFIRST(IELH) N1H=NLAST(IELH) NKH=NKA(IATH) N0Q=NQUANT(N1H)+1 N1Q=ICUP(IELH) DO 200 II=N0HN,N1H I=II-N0HN+1 IT=ITRA(II,NKH) IF(IT.EQ.0) GO TO 100 C C *************** Collisional ionization C c for high temperature, use XSTAR formulae C if(t0.gt.1.e6) then rno=16. izc=1 call irc(i,t0,izc,rno,cs) col(it)=cs go to 100 end if C IC=ICOL(IT) U0=FR0(IT)*HKT IF(IC.LT.0) GO TO 90 if(ifwop(ii).lt.0) go to 95 GAM=I*I*I IF(I.GT.10) GO TO 80 GAM=A(1,I)+A(2,I)*X+A(3,I)*X2+A(4,I)*X3 * +A(5,I)*X4+A(6,I)*X5 80 COL(IT)=CT*EXP(-U0)*GAM GO TO 100 C C non-standard (user supplied) formula C 90 CALL CSPEC(II,NKH,IC,OSC0(IT),CPAR(IT),U0,T,COL(IT)) go to 100 c c ionization from the merged state c 95 sum1=0. sum2=0. ehk=eh/tk n00q=nquant(n1h-1)+1 n11q=nlmx do img=n00q,n11q xi=img xii=xi*xi sum1=sum1+xii*xii*xi*wnhint(img,id) sum2=sum2+xii*wnhint(img,id)*exp(ehk/xii) end do col(it)=ct*sum1/sum2 go to 200 C C ***************** Collisional excitation C 100 CONTINUE I1=I+1 XI=I VI=XI*XI ALF=ALF0-ALF1/VI BET=BET0-BET1/XI NHL=N1H-N0HN+1 IF(N1Q.GT.0) NHL=N1Q N1HC=N1H IF(IFWOP(N1H).LT.0) THEN NHL=NLMX N1HC=N1H-1 END IF CSUM=0. IF(I1.GT.NHL) GO TO 200 CSCA=8.63D-6/2./VI/SQT DO 190 J=I1,NHL XJ=J VJ=XJ*XJ IC=0 JJ=J+N0HN-1 IF(JJ.GT.N1HC) GO TO 150 ICT=ITRA(II,JJ) IF(ICT.EQ.0) GO TO 190 IC=ICOL(ICT) U0=FR0(ICT)*HKT E=U0/EH/TK C1=OSC0(ICT) IF(IC.LT.0) THEN CALL CSPEC(II,JJ,IC,C1,CPAR(ICT),U0,T,COL(ICT)) ELSE IF(IC.EQ.0) THEN GO TO 160 ELSE IF(IC.EQ.1) THEN COL(ICT)=CT*EXP(-U0)*(CEX1+CEX2*X+CEX3/X/X) ELSE IF(IC.GE.2) THEN COL(ICT)=CEH12(T) END IF GO TO 190 C C collisional excitations from level I to higher, non-explicit C levels are lumped into the collisional ionization rate C (the so-called modified collision ionization rate) C 150 CONTINUE E=UN/VI-UN/VJ U0=EH*E*TK IF(J.LE.20) C1=OSH(I,J) IF(J.GT.20) THEN C1=OSH(I,20)*((400.-VI)/20.*XJ/(VJ-VI))**3 end if 160 CONTINUE C IF(ICOLHN.EQ.1.AND.J.LE.7) GO TO 250 IF(ICOLHN.EQ.2.AND.J.LE.15) GO TO 260 C C Old standard formula for the collisional excitation rate - used for C rates in explicit transitions as well as for evaluation of the C modified collisional rate C IF(ICOLHN.EQ.1.AND.J.LE.7) GO TO 250 IF(ICOLHN.EQ.2.AND.J.LE.15) GO TO 260 CS=4.*CT*C1/E/E EX=EXP(-U0) IF(U0.LE.UN) THEN E1=-LOG(U0)+EXPIA1+U0*(EXPIA2+U0*(EXPIA3+U0*(EXPIA4+ * U0*(EXPIA5+U0*EXPIA6)))) ELSE E1=EXP(-U0)*((EXPIB1+U0*(EXPIB2+U0*(EXPIB3+ * U0*(EXPIB4+U0))))/(EXPIC1+U0*(EXPIC2+ * U0*(EXPIC3+U0*(EXPIC4+U0)))))/U0 END IF E5=E1 DO IX=1,4 E5=(EX-U0*E5)/IX END DO CS=CS*U0*(E1+O148*U0*E5) IF(J-I.NE.1) CS=CS*(BET+TWO*(ALF-BET)/(XJ-XI)) GO TO 180 C End of the old standard formula (Mihalas et al 1975) C c Butler new calculations c 250 call butler(i,j,t,u0,cs,ierr) go to 180 c C Giovanardi et al. 1987, AAS, 70, 269 C Cool: T<=60000K ; Hot: T>60000K C 260 IF(T.GT.60000.) GO TO 270 CS=CCOOL(1,I,J) DO ICA=2,4 CS=CS+CCOOL(ICA,I,J)*XTT(ICA) END DO GO TO 280 270 CS=CHOT(1,I,J) DO ICA=2,4 CS=CS+CHOT(ICA,I,J)*XTT(ICA) END DO 280 CS=CSCA*CS*EXP(-U0) C 180 IF(JJ.GT.N1HC) THEN CSUM=CSUM+CS ELSE COL(ICT)=CS END IF 190 CONTINUE IF(IT.NE.0.AND.N1Q.GT.0) COL(IT)=COL(IT)+CSUM ITH=ITRA(II,N1H) IF(IFWOP(N1H).LT.0.AND.ITH.GT.0) COL(ITH)=CSUM 200 CONTINUE C C special standard formula for collisional ionization of H- C IF(IELHM.EQ.0) RETURN IT=ITRA(NFIRST(IELHM),N0HN) IF(IT.EQ.0) RETURN IC=ICOL(IT) IF(IC.GE.0) THEN COL(IT)=CHMI*T*SQRT(T) ELSE C C if desired, non-standard, user supplied, formula for H- C U0=ENION(NFIRST(IELHM))*TK CALL CSPEC(NFIRST(IELHM),N0HN,IC,OSC0(IT),CPAR(IT),U0,T,CS) COL(IT)=CS END IF RETURN END