***** Main program Sfitting.f ******************************************************** * * Definition of variables * * NZ : index which controles (Z,A) values, 1…NZ…11. * NZ= 1 (1H), NZ= 2 (4He), NZ= 3 (12C), NZ= 4 (14N), NZ= 5 (16O), NZ=6(20Ne), * NZ= 7 (24Mg), NZ= 8 (28Si), NZ= 9 (32S), NZ=10 (40Ca), NZ=11 (56Fe). * Z : atomic number, 1… Z…26. * A : mass number, 1… A…56. * * Gam: ion coupling parameter Gamma, 0.1…Gam…180. * rho: mass density in g/cm3, 0…log10(rho)…12.8 * * S : Fitted value for the Coulomb logarithm (scattering integral) . * is defined by Eq.(2.8) in the paper by Itoh et al. arXiv:0708.2967[astro-ph]. * The fitting formulae are defined by Eqs.(4.1)-(4.14), and the fitting * coefficients are given in Tables 2-5 in the paper. * * inputs: rho, Gam, NZ * * output: S * ******************************************************************************************* implicit none real*8 Z,A,pi,alp,Gam,x,rhop,u,C1,RHO,vv,ww,zz,xGam real*8 Smin,Smin1,Smin2,Szero,Szero1,Szero2,Splu,Splu1,Splu2,S real*8 T2(100,100),T3(100,100),T4(100,100),T5(100,100) integer*4 NZ * Open a file for writing data (rho, Gam, Z, A, S) open (unit=10, file='Sfitting.csv') * inputs ********************************************************************************** * Specify a rho value within 0…log10(rho)…12.8 rho=1.d+8 * Specify a Gam value within 0.1…Gam…180 Gam=10.d0 * Specify a NZ value within 1…NZ…11 NZ = 3 ******************************************************************************************* if(NZ.eq.1)then Z=1.d0 A=1.d0 elseif(NZ.eq.2)then Z=2.d0 A=4.d0 elseif(NZ.eq.3)then Z=6.d0 A=12.d0 elseif(NZ.eq.4)then Z=7.d0 A=14.d0 elseif(NZ.eq.5)then Z=8.d0 A=16.d0 elseif(NZ.eq.6)then Z=10.d0 A=20.d0 elseif(NZ.eq.7)then Z=12.d0 A=24.d0 elseif(NZ.eq.8)then Z=14.d0 A=28.d0 elseif(NZ.eq.9)then Z=16.d0 A=32.d0 elseif(NZ.eq.10)then Z=20.d0 A=40.d0 elseif(NZ.eq.11)then Z=26.d0 A=56.d0 endif data T2(1,1),T2(1,2),T2(1,3),T2(1,4),T2(1,5),T2(1,6),T2(1,7), * T2(1,8),T2(1,9),T2(1,10),T2(1,11),T2(1,12),T2(1,13),T2(1,14), * T2(2,1),T2(2,2),T2(2,3),T2(2,4),T2(2,5),T2(2,6),T2(2,7), * T2(2,8),T2(2,9),T2(2,10),T2(2,11),T2(2,12),T2(2,13),T2(2,14), * T2(3,1),T2(3,2),T2(3,3),T2(3,4),T2(3,5),T2(3,6),T2(3,7), * T2(3,8),T2(3,9),T2(3,10),T2(3,11),T2(3,12),T2(3,13),T2(3,14), * T2(4,1),T2(4,2),T2(4,3),T2(4,4),T2(4,5),T2(4,6),T2(4,7), * T2(4,8),T2(4,9),T2(4,10),T2(4,11),T2(4,12),T2(4,13),T2(4,14), * T2(5,1),T2(5,2),T2(5,3),T2(5,4),T2(5,5),T2(5,6),T2(5,7), * T2(5,8),T2(5,9),T2(5,10),T2(5,11),T2(5,12),T2(5,13),T2(5,14), * T2(6,1),T2(6,2),T2(6,3),T2(6,4),T2(6,5),T2(6,6),T2(6,7), * T2(6,8),T2(6,9),T2(6,10),T2(6,11),T2(6,12),T2(6,13),T2(6,14), * T2(7,1),T2(7,2),T2(7,3),T2(7,4),T2(7,5),T2(7,6),T2(7,7), * T2(7,8),T2(7,9),T2(7,10),T2(7,11),T2(7,12),T2(7,13),T2(7,14), * T2(8,1),T2(8,2),T2(8,3),T2(8,4),T2(8,5),T2(8,6),T2(8,7), * T2(8,8),T2(8,9),T2(8,10),T2(8,11),T2(8,12),T2(8,13),T2(8,14), * T2(9,1),T2(9,2),T2(9,3),T2(9,4),T2(9,5),T2(9,6),T2(9,7), * T2(9,8),T2(9,9),T2(9,10),T2(9,11),T2(9,12),T2(9,13),T2(9,14), * T2(10,1),T2(10,2),T2(10,3),T2(10,4),T2(10,5),T2(10,6), * T2(10,7),T2(10,8),T2(10,9),T2(10,10),T2(10,11),T2(10,12), * T2(10,13),T2(10,14), * T2(11,1),T2(11,2),T2(11,3),T2(11,4),T2(11,5),T2(11,6), * T2(11,7),T2(11,8),T2(11,9),T2(11,10),T2(11,11),T2(11,12), * T2(11,13),T2(11,14)/ * 0.6496d0,0.0471d0,-0.0056d0,-0.0284d0,0.0054d0,0.0921d0, * 0.4531d0,0.0268d0,0.0012d0,0.0051d0,-0.0007d0,0.0018d0, * 0.0011d0,0.0621d0,0.7407d0,-0.0007d0,-0.0165d0,-0.0376d0, * 0.0114d0,0.1037d0,0.3959d0,0.2196d0,0.0006d0,0.044d0, * -0.0081d0,0.0155d0,0.0056d0,0.3641d0,0.8981d0,-0.0666d0, * -0.0071d0,-0.0558d0,0.0247d0,0.1068d0,0.404d0,0.4347d0, * 0.0084d0,0.0741d0,-0.0189d0,0.0228d0,0.0174d0,0.3604d0, * 0.9232d0,-0.0781d0,-0.0045d0,-0.0588d0,0.027d0,0.1064d0, * 0.4047d0,0.4753d0,0.0052d0,0.0796d0,-0.0224d0,0.0241d0, * 0.0184d0,0.367d0,0.9457d0,-0.0884d0,-0.0019d0,-0.0615d0, * 0.0291d0,0.1059d0,0.4053d0,0.5166d0,0.0006d0,0.0854d0, * -0.0263d0,0.0257d0,0.0189d0,0.3787d0,0.9848d0,-0.1065d0, * 0.0031d0,-0.0663d0,0.0326d0,0.1046d0,0.4063d0,0.593d0, * -0.0094d0,0.0962d0,-0.0338d0,0.0287d0,0.0194d0,0.4011d0, * 1.0181d0,-0.1221d0,0.0076d0,-0.0703d0,0.0356d0,0.1032d0, * 0.4069d0,0.652d0,-0.0176d0,0.1035d0,-0.0399d0,0.0305d0, * 0.02d0,0.4097d0,1.0471d0,-0.1357d0,0.0117d0,-0.0737d0, * 0.0382d0,0.1018d0,0.4074d0,0.6976d0,-0.0243d0,0.1082d0, * -0.0449d0,0.0316d0,0.0208d0,0.4097d0,1.0729d0,-0.1477d0, * 0.0155d0,-0.0767d0,0.0404d0,0.1004d0,0.4078d0,0.7358d0, * -0.0304d0,0.1116d0,-0.0494d0,0.0323d0,0.0214d0,0.4074d0, * 1.1171d0,-0.1684d0,0.0222d0,-0.0818d0,0.044d0,0.0977d0, * 0.4084d0,0.8004d0,-0.0423d0,0.1169d0,-0.0571d0,0.0333d0, * 0.0223d0,0.4036d0,1.169d0,-0.197d0,0.0297d0,-0.0869d0, * 0.0481d0,0.0946d0,0.4017d0,0.8856d0,-0.0612d0,0.1224d0, * -0.0679d0,0.0344d0,0.023d0,0.3986d0/ data T3(1,1),T3(1,2),T3(1,3),T3(1,4),T3(1,5),T3(1,6),T3(1,7), * T3(1,8),T3(1,9),T3(1,10),T3(1,11),T3(1,12),T3(1,13),T3(1,14), * T3(2,1),T3(2,2),T3(2,3),T3(2,4),T3(2,5),T3(2,6),T3(2,7), * T3(2,8),T3(2,9),T3(2,10),T3(2,11),T3(2,12),T3(2,13),T3(2,14), * T3(3,1),T3(3,2),T3(3,3),T3(3,4),T3(3,5),T3(3,6),T3(3,7), * T3(3,8),T3(3,9),T3(3,10),T3(3,11),T3(3,12),T3(3,13),T3(3,14), * T3(4,1),T3(4,2),T3(4,3),T3(4,4),T3(4,5),T3(4,6),T3(4,7), * T3(4,8),T3(4,9),T3(4,10),T3(4,11),T3(4,12),T3(4,13),T3(4,14), * T3(5,1),T3(5,2),T3(5,3),T3(5,4),T3(5,5),T3(5,6),T3(5,7), * T3(5,8),T3(5,9),T3(5,10),T3(5,11),T3(5,12),T3(5,13),T3(5,14), * T3(6,1),T3(6,2),T3(6,3),T3(6,4),T3(6,5),T3(6,6),T3(6,7), * T3(6,8),T3(6,9),T3(6,10),T3(6,11),T3(6,12),T3(6,13),T3(6,14), * T3(7,1),T3(7,2),T3(7,3),T3(7,4),T3(7,5),T3(7,6),T3(7,7), * T3(7,8),T3(7,9),T3(7,10),T3(7,11),T3(7,12),T3(7,13),T3(7,14), * T3(8,1),T3(8,2),T3(8,3),T3(8,4),T3(8,5),T3(8,6),T3(8,7), * T3(8,8),T3(8,9),T3(8,10),T3(8,11),T3(8,12),T3(8,13),T3(8,14), * T3(9,1),T3(9,2),T3(9,3),T3(9,4),T3(9,5),T3(9,6),T3(9,7), * T3(9,8),T3(9,9),T3(9,10),T3(9,11),T3(9,12),T3(9,13),T3(9,14), * T3(10,1),T3(10,2),T3(10,3),T3(10,4),T3(10,5),T3(10,6), * T3(10,7),T3(10,8),T3(10,9),T3(10,10),T3(10,11),T3(10,12), * T3(10,13),T3(10,14), * T3(11,1),T3(11,2),T3(11,3),T3(11,4),T3(11,5),T3(11,6), * T3(11,7),T3(11,8),T3(11,9),T3(11,10),T3(11,11),T3(11,12), * T3(11,13),T3(11,14)/ * 0.2781d0,0.0357d0,0.0224d0,-0.0072d0,0.0059d0,0.0303d0, * 0.3087d0,0.0225d0,0.0007d0,0.0045d0,-0.0006d0,0.0016d0, * 0.0007d0,0.0557d0,0.3281d0,0.0222d0,0.0249d0,-0.0134d0, * 0.0082d0,0.0323d0,0.2749d0,0.1881d0,-0.0008d0,0.0385d0, * -0.0071d0,0.0136d0,0.0041d0,0.32d0,0.4042d0,-0.0077d0, * 0.0396d0,-0.0245d0,0.0152d0,0.0269d0,0.279d0,0.2913d0, * 0.0011d0,0.0519d0,-0.0133d0,0.0164d0,0.0095d0,0.2571d0, * 0.417d0,-0.0131d0,0.0423d0,-0.0264d0,0.0163d0,0.0258d0, * 0.2794d0,0.3121d0,-0.0026d0,0.0551d0,-0.0156d0,0.0173d0, * 0.0093d0,0.2604d0,0.4286d0,-0.018d0,0.0448d0,-0.028d0, * 0.0174d0,0.0247d0,0.2797d0,0.3353d0,-0.0072d0,0.059d0, * -0.0182d0,0.0185d0,0.0087d0,0.2694d0,0.4489d0,-0.0266d0, * 0.0491d0,-0.0309d0,0.0191d0,0.0227d0,0.2802d0,0.3787d0, * -0.0166d0,0.0664d0,-0.0229d0,0.0208d0,0.0074d0,0.2873d0, * 0.4662d0,-0.0339d0,0.0528d0,-0.0333d0,0.0205d0,0.0209d0, * 0.2805d0,0.4085d0,-0.0233d0,0.0711d0,-0.0263d0,0.0221d0, * 0.0066d0,0.2923d0,0.4813d0,-0.0404d0,0.056d0,-0.0353d0, * 0.0217d0,0.0193d0,0.2808d0,0.4283d0,-0.0281d0,0.0736d0, * -0.0288d0,0.0228d0,0.0061d0,0.2897d0,0.4946d0,-0.046d0, * 0.0588d0,-0.037d0,0.0227d0,0.0179d0,0.281d0,0.4436d0, * -0.0321d0,0.0754d0,-0.0307d0,0.0231d0,0.0056d0,0.2856d0, * 0.5173d0,-0.0556d0,0.0635d0,-0.0397d0,0.0241d0,0.0154d0, * 0.2813d0,0.4687d0,-0.0395d0,0.0781d0,-0.0339d0,0.0236d0, * 0.0047d0,0.2799d0,0.5452d0,-0.0672d0,0.0684d0,-0.0426d0, * 0.0254d0,0.0127d0,0.2773d0,0.5053d0,-0.0505d0,0.0819d0, * -0.0383d0,0.0243d0,0.0033d0,0.2763d0/ data T4(1,1),T4(1,2),T4(1,3),T4(1,4),T4(1,5),T4(1,6),T4(1,7), * T4(1,8),T4(1,9),T4(1,10),T4(1,11),T4(1,12),T4(1,13),T4(1,14), * T4(2,1),T4(2,2),T4(2,3),T4(2,4),T4(2,5),T4(2,6),T4(2,7), * T4(2,8),T4(2,9),T4(2,10),T4(2,11),T4(2,12),T4(2,13),T4(2,14), * T4(3,1),T4(3,2),T4(3,3),T4(3,4),T4(3,5),T4(3,6),T4(3,7), * T4(3,8),T4(3,9),T4(3,10),T4(3,11),T4(3,12),T4(3,13),T4(3,14), * T4(4,1),T4(4,2),T4(4,3),T4(4,4),T4(4,5),T4(4,6),T4(4,7), * T4(4,8),T4(4,9),T4(4,10),T4(4,11),T4(4,12),T4(4,13),T4(4,14), * T4(5,1),T4(5,2),T4(5,3),T4(5,4),T4(5,5),T4(5,6),T4(5,7), * T4(5,8),T4(5,9),T4(5,10),T4(5,11),T4(5,12),T4(5,13),T4(5,14), * T4(6,1),T4(6,2),T4(6,3),T4(6,4),T4(6,5),T4(6,6),T4(6,7), * T4(6,8),T4(6,9),T4(6,10),T4(6,11),T4(6,12),T4(6,13),T4(6,14), * T4(7,1),T4(7,2),T4(7,3),T4(7,4),T4(7,5),T4(7,6),T4(7,7), * T4(7,8),T4(7,9),T4(7,10),T4(7,11),T4(7,12),T4(7,13),T4(7,14), * T4(8,1),T4(8,2),T4(8,3),T4(8,4),T4(8,5),T4(8,6),T4(8,7), * T4(8,8),T4(8,9),T4(8,10),T4(8,11),T4(8,12),T4(8,13),T4(8,14), * T4(9,1),T4(9,2),T4(9,3),T4(9,4),T4(9,5),T4(9,6),T4(9,7), * T4(9,8),T4(9,9),T4(9,10),T4(9,11),T4(9,12),T4(9,13),T4(9,14), * T4(10,1),T4(10,2),T4(10,3),T4(10,4),T4(10,5),T4(10,6), * T4(10,7),T4(10,8),T4(10,9),T4(10,10),T4(10,11),T4(10,12), * T4(10,13),T4(10,14), * T4(11,1),T4(11,2),T4(11,3),T4(11,4),T4(11,5),T4(11,6), * T4(11,7),T4(11,8),T4(11,9),T4(11,10),T4(11,11),T4(11,12), * T4(11,13),T4(11,14)/ * 0.1543d0,0.0202d0,0.0205d0,-0.0024d0,0.0057d0,0.0133d0, * 0.2293d0,0.0196d0,0.0004d0,0.0041d0,-0.0006d0,0.0015d0, * 0.0005d0,0.0507d0,0.1881d0,0.0137d0,0.0248d0,-0.0063d0, * 0.0074d0,0.0135d0,0.2069d0,0.1636d0,-0.0017d0,0.0339d0, * -0.0064d0,0.0121d0,0.0032d0,0.2834d0,0.238d0,-0.0068d0, * 0.0367d0,-0.014d0,0.0123d0,0.0084d0,0.2095d0,0.2027d0, * -0.0028d0,0.0375d0,-0.0099d0,0.0121d0,0.0051d0,0.1902d0, * 0.2466d0,-0.0105d0,0.0388d0,-0.0152d0,0.0131d0,0.0074d0, * 0.2097d0,0.2148d0,-0.0064d0,0.0397d0,-0.0115d0,0.0129d0, * 0.0044d0,0.1933d0,0.2544d0,-0.0138d0,0.0406d0,-0.0164d0, * 0.0137d0,0.0065d0,0.2099d0,0.2302d0,-0.0108d0,0.0428d0, * -0.0134d0,0.0139d0,0.0035d0,0.2021d0,0.2679d0,-0.0196d0, * 0.0438d0,-0.0183d0,0.0149d0,0.0049d0,0.2102d0,0.2597d0, * -0.0191d0,0.0487d0,-0.0168d0,0.0159d0,0.0017d0,0.2188d0, * 0.2794d0,-0.0244d0,0.0465d0,-0.0198d0,0.0158d0,0.0036d0, * 0.2104d0,0.2776d0,-0.0245d0,0.0521d0,-0.019d0,0.017d0, * 0.0005d0,0.2226d0,0.2893d0,-0.0287d0,0.0487d0,-0.0211d0, * 0.0165d0,0.0024d0,0.2105d0,0.2872d0,-0.0279d0,0.0536d0, * -0.0202d0,0.0174d0,-0.0002d0,0.2192d0,0.298d0,-0.0323d0, * 0.0506d0,-0.0222d0,0.017d0,0.0013d0,0.2106d0,0.2938d0, * -0.0305d0,0.0545d0,-0.0211d0,0.0175d0,-0.0008d0,0.2147d0, * 0.3126d0,-0.0384d0,0.0537d0,-0.0238d0,0.0178d0,-0.0005d0, * 0.2108d0,0.3049d0,-0.0352d0,0.0561d0,-0.0226d0,0.0177d0, * -0.0019d0,0.2089d0,0.3306d0,-0.0451d0,0.0569d0,-0.0254d0, * 0.0183d0,-0.0024d0,0.2082d0,0.3248d0,-0.0422d0,0.0589d0, * -0.0246d0,0.0181d0,-0.0035d0,0.2068d0/ data T5(1,1),T5(1,2),T5(1,3),T5(1,4),T5(1,5),T5(1,6),T5(1,7), * T5(1,8),T5(1,9),T5(1,10),T5(1,11),T5(1,12), * T5(2,1),T5(2,2),T5(2,3),T5(2,4),T5(2,5),T5(2,6),T5(2,7), * T5(2,8),T5(2,9),T5(2,10),T5(2,11),T5(2,12), * T5(3,1),T5(3,2),T5(3,3),T5(3,4),T5(3,5),T5(3,6),T5(3,7), * T5(3,8),T5(3,9),T5(3,10),T5(3,11),T5(3,12), * T5(4,1),T5(4,2),T5(4,3),T5(4,4),T5(4,5),T5(4,6),T5(4,7), * T5(4,8),T5(4,9),T5(4,10),T5(4,11),T5(4,12), * T5(5,1),T5(5,2),T5(5,3),T5(5,4),T5(5,5),T5(5,6),T5(5,7), * T5(5,8),T5(5,9),T5(5,10),T5(5,11),T5(5,12), * T5(6,1),T5(6,2),T5(6,3),T5(6,4),T5(6,5),T5(6,6),T5(6,7), * T5(6,8),T5(6,9),T5(6,10),T5(6,11),T5(6,12), * T5(7,1),T5(7,2),T5(7,3),T5(7,4),T5(7,5),T5(7,6),T5(7,7), * T5(7,8),T5(7,9),T5(7,10),T5(7,11),T5(7,12), * T5(8,1),T5(8,2),T5(8,3),T5(8,4),T5(8,5),T5(8,6),T5(8,7), * T5(8,8),T5(8,9),T5(8,10),T5(8,11),T5(8,12), * T5(9,1),T5(9,2),T5(9,3),T5(9,4),T5(9,5),T5(9,6),T5(9,7), * T5(9,8),T5(9,9),T5(9,10),T5(9,11),T5(9,12), * T5(10,1),T5(10,2),T5(10,3),T5(10,4),T5(10,5),T5(10,6), * T5(10,7),T5(10,8),T5(10,9),T5(10,10),T5(10,11),T5(10,12), * T5(11,1),T5(11,2),T5(11,3),T5(11,4),T5(11,5),T5(11,6), * T5(11,7),T5(11,8),T5(11,9),T5(11,10),T5(11,11),T5(11,12)/ * 0.4288d0,-0.5654d0,0.0769d0,0.0662d0,0.5283d0,-0.5546d0, * -0.0256d0,0.0559d0,0.6073d0,-0.5135d0,-0.1064d0,0.0143d0, * 0.1778d0,-0.5446d0,0.3337d0,0.0467d0,0.009d0,-0.7314d0, * 0.5152d0,0.2429d0,-4.1969d0,-2.82d0,4.9321d0,2.461d0, * 0.2634d0,-0.5073d0,0.2504d0,0.0083d0,0.3358d0,-0.4548d0, * 0.1868d0,-0.0421d0,0.474d0,-0.1732d0,0.0621d0,-0.3263d0, * 0.2661d0,-0.5144d0,0.2474d0,0.0155d0,0.3352d0,-0.4724d0, * 0.185d0,-0.0245d0,0.4662d0,-0.2319d0,0.0612d0,-0.2675d0, * 0.2619d0,-0.5252d0,0.2512d0,0.0264d0,0.3055d0,-0.509d0, * 0.2121d0,0.0124d0,0.3641d0,-0.3358d0,0.1539d0,-0.1632d0, * 0.2489d0,-0.5454d0,0.2638d0,0.0472d0,0.2147d0,-0.5917d0, * 0.3007d0,0.0973d0,-0.2351d0,-0.7498d0,0.7369d0,0.2631d0, * 0.2461d0,-0.5525d0,0.2674d0,0.0549d0,0.1748d0,-0.6101d0, * 0.3446d0,0.1193d0,-0.9551d0,-1.0234d0,1.4926d0,0.5803d0, * 0.2513d0,-0.5518d0,0.2632d0,0.0547d0,0.1905d0,-0.5803d0, * 0.3334d0,0.0915d0,-0.6419d0,-0.6405d0,1.2221d0,0.2096d0, * 0.2583d0,-0.5496d0,0.2568d0,0.0527d0,0.2221d0,-0.5438d0, * 0.3045d0,0.0556d0,-0.1657d0,-0.3169d0,0.7513d0,-0.1263d0, * 0.27d0,-0.548d0,0.2454d0,0.0513d0,0.2725d0,-0.5002d0, * 0.2539d0,0.0107d0,0.2634d0,-0.1219d0,0.2993d0,-0.3459d0, * 0.276d0,-0.5561d0,0.2389d0,0.0594d0,0.2749d0,-0.5122d0, * 0.2464d0,0.0208d0,0.2677d0,-0.1927d0,0.2586d0,-0.2894d0/ pi=4.d0*datan(1.d0) alp=1.d0/137.036d0 if (Gam.eq.0.1d0) then xGam=-1.d0 elseif (Gam.eq.180.d0) then xGam=1.d0 else xGam=0.61439d0*dlog10(Gam)-0.38561d0 endif vv=T5(NZ,1)+T5(NZ,2)*xGam+T5(NZ,3)*xGam**2+T5(NZ,4)*xGam**3 ww=T5(NZ,5)+T5(NZ,6)*xGam+T5(NZ,7)*xGam**2+T5(NZ,8)*xGam**3 zz=T5(NZ,9)+T5(NZ,10)*xGam+T5(NZ,11)*xGam**2+T5(NZ,12)*xGam**3 x=log10(rho) u=2.d0*pi*x/25.6d0 C1=(1.018d0*(Z/A)**(2.d0/3.d0)*(rho*1.d0/(10.d0**6))**(2.d0/3.d0)) * /(1.d0+1.018d0*(Z/A)**(2.d0/3.d0)*(rho*1.d0/(10.d0**6)) * **(2.d0/3.d0)) Smin1=T2(NZ,1)*dsin(u)+T2(NZ,2)*dsin(2.d0*u) * +T2(NZ,3)*dsin(3.d0*u)+T2(NZ,4)*dsin(4.d0*u) * +T2(NZ,5)*dsin(5.d0*u)+12.8d0/pi*T2(NZ,6)*u+T2(NZ,7) Smin2=T2(NZ,8)*dsin(u)+T2(NZ,9)*dsin(2.d0*u) * +T2(NZ,10)*dsin(3.d0*u)+T2(NZ,11)*dsin(4.d0*u) * +T2(NZ,12)*dsin(5.d0*u)+12.8d0/pi*T2(NZ,13)*u+T2(NZ,14) Szero1=T3(NZ,1)*dsin(u)+T3(NZ,2)*dsin(2.d0*u) * +T3(NZ,3)*dsin(3.d0*u)+T3(NZ,4)*dsin(4.d0*u) * +T3(NZ,5)*dsin(5.d0*u)+12.8d0/pi*T3(NZ,6)*u+T3(NZ,7) Szero2=T3(NZ,8)*dsin(u)+T3(NZ,9)*dsin(2.d0*u) * +T3(NZ,10)*dsin(3.d0*u)+T3(NZ,11)*dsin(4.d0*u) * +T3(NZ,12)*dsin(5.d0*u)+12.8d0/pi*T3(NZ,13)*u+T3(NZ,14) Splu1=T4(NZ,1)*dsin(u)+T4(NZ,2)*dsin(2.d0*u) * +T4(NZ,3)*dsin(3.d0*u)+T4(NZ,4)*dsin(4.d0*u) * +T4(NZ,5)*dsin(5.d0*u)+12.8d0/pi*T4(NZ,6)*u+T4(NZ,7) Splu2=T4(NZ,8)*dsin(u)+T4(NZ,9)*dsin(2.d0*u) * +T4(NZ,10)*dsin(3.d0*u)+T4(NZ,11)*dsin(4.d0*u) * +T4(NZ,12)*dsin(5.d0*u)+12.8d0/pi*T4(NZ,13)*u+T4(NZ,14) Smin=vv*Smin1+(1.d0-vv)*Smin2 Szero=ww*Szero1+(1.d0-ww)*Szero2 Splu=zz*Splu1+(1.d0-zz)*Splu2 S=Smin-C1*Splu+pi*Z*alp*dsqrt(C1)*(Szero-Splu) * output: Fitted value S for the Coulomb logarithm *********************************** write(10, 1000) rho, Gam, Z, A, S ******************************************************************************************* 1000 format(d12.5, 4(',', f8.4)) end