238 lines
13 KiB
Fortran
238 lines
13 KiB
Fortran
FUNCTION VERNER(FR,ITR)
|
|
C =======================
|
|
C
|
|
C Photoionization cross-sections for ground states
|
|
C of atoms and ions. Analytical fits from:
|
|
C Verner D.A. et al. 1996, ApJ 465
|
|
C Verner & Yakovlev 1995, A&AS 109, 125
|
|
C
|
|
C 10-July-1996: Version for H to Si, S, Ar, Ca and Fe.
|
|
C No test on threshold energy as given by Verner et al.
|
|
C to avoid inconsistencies with limits FR0 as read by Tlusty
|
|
C
|
|
C
|
|
INCLUDE 'IMPLIC.FOR'
|
|
INCLUDE 'BASICS.FOR'
|
|
INCLUDE 'ATOMIC.FOR'
|
|
C
|
|
PARAMETER (HHEV=H/1.6022D-12,T18=1.D-18,MVER=105)
|
|
DIMENSION S0(MVER),E0(MVER),EMX(MVER)
|
|
DIMENSION Y0(MVER),Y1(MVER),YW(MVER),YA(MVER),PV(MVER)
|
|
DIMENSION S95(MVER),E95(MVER),Y95(MVER),YW95(MVER),P95(MVER)
|
|
DIMENSION IV0(14)
|
|
C
|
|
DATA IV0/0,1,3,6,10,15,21,28,36,45,55,66,78,91/
|
|
DATA S0/5.475D4,9.492D2,1.369D4,6.245D1,3.201D2,6.083D3,2.932D5,
|
|
+ 2.678D2,5.458D2,3.422D3,5.466D0,1.859D4,5.393D1,2.846D2,
|
|
+ 2.190D3,5.027D2,8.709D0,1.539D4,1.068D2,2.344D2,1.521D3,
|
|
+ 8.235D2,1.944D0,9.375D-1,1.690D4,8.376D1,1.519D2,1.117D3,
|
|
+ 1.745D3,5.967D1,6.753D2,8.659D-1,1.642D4,6.864D1,1.329D2,
|
|
+ 8.554D2,3.803D3,8.013D1,1.541D2,3.165D0,4.690D-1,1.157D4,
|
|
+ 6.930D1,1.039D2,6.759D2,4.287D3,1.583D3,5.708D0,1.685D3,
|
|
+ 2.430D0,9.854D-1,1.198D4,5.631D1,6.695D1,5.475D2,1.601D0,
|
|
+ 1.04D3,1.885D3,2.33D0,2.346D1,7.101D1,1.609D0,7.215D3,
|
|
+ 4.729D1,3.995D1,4.525D2,1.372D8,3.278D0,5.377D2,1.394D3,
|
|
+ 1.728D0,2.185D0,3.104D0,6.344D-2,9.008D2,4.427D1,6.14D1,
|
|
+ 3.802D2,7.195D0,6.948D-2,4.915D0,1.513D1,2.925D2,1.962D-2,
|
|
+ 1.889D2,2.35D-1,4.982D-1,1.774D4,3.388D1,4.036D1,3.239D2,
|
|
+ 2.506D1,4.14D0,5.79D-4,6.083D0,8.863D-1,7.293D1,6.68D-2,
|
|
+ 3.477D-1,1.465D-1,1.95D-1,1.992D4,2.539D1,4.754D1,2.793D2/
|
|
DATA E0/0.4298,13.61,1.72,3.107,20.06,3.871,9.539,1.181,17.6,
|
|
+ 6.879,0.5213,2.869,1.041,33.36,10.75,2.144,0.4058,4.614,
|
|
+ 3.506,46.24,15.48,4.034,0.06128,0.242,5.494,4.471,69.43,
|
|
+ 21.08,1.24,1.386,0.1723,0.2044,2.854,7.824,87.09,27.54,
|
|
+ 12.97,1.763,2.542,0.7744,0.7286,4.008,2.563,113.1,34.85,
|
|
+ 4.87,12.47,0.7753,5.566,1.248,1.499,4.888,10.03,158.6,
|
|
+ 43.04,6.139,8.203,10.69,0.669,5.408,48.46,2.096,153.5,
|
|
+ 13.91,226.8,52.11,11.97,8.139,10.86,29.12,0.9762,1.711,
|
|
+ 3.57,0.4884,34.82,14.52,204.2,62.03,13.81,0.2048,10.27,
|
|
+ 3.13,24.14,0.3483,2.636,0.4866,1.842,8.044,23.55,273.8,
|
|
+ 72.81,23.17,2.556,0.1659,1.288,0.7761,63.05,0.3277,0.7655,
|
|
+ 0.3343,0.8787,12.05,35.6,275.2,84.47/
|
|
DATA EMX/5.D4,5.D4,5.D4,64.39,5.D4,5.D4,11.93,12.99,5.D4,5.D4,
|
|
+ 194.0,209.8,227.4,5.D4,5.D4,291.0,307.6,328.9,352.2,
|
|
+ 5.D4,5.D4,404.8,423.6,447.3,475.3,504.3,5.D4,5.D4,
|
|
+ 538.0,558.1,584.0,614.4,649.1,683.7,5.D4,5.D4,694.0,
|
|
+ 712.2,739.2,770.9,809.1,850.2,890.5,5.D4,5.D4,870.1,
|
|
+ 883.1,913.1,948.0,987.3,1031.,1078.,1125.,5.D4,5.D4,
|
|
+ 38.14,1074.,1110.,1143.,1185.,1230.,1281.,1335.,1386.,
|
|
+ 5.D4,5.D4,54.9,65.69,1317.,1356.,1400.,1449.,1501.,
|
|
+ 1558.,1618.,1675.,5.D4,5.D4,80.4,89.97,102.6,1588.,
|
|
+ 1634.,1688.,1739.,1799.,1862.,1929.,1992.,5.D4,5.D4,
|
|
+ 106.,118.6,131.1,146.6,1887.,1946.,2001.,2058.,2125.,
|
|
+ 2194.,2268.,2336.,5.D4,5.D4/
|
|
DATA Y0/0.,0.4434,0.,0.,0.,0.,8.278D-4,0.,0.,0.,1.319D1,4.96D-3,
|
|
+ 0.,0.,0.,1.133,4.929D1,4.378D-3,0.,0.,0.,8.598D-1,4.28D2,
|
|
+ 1.877D2,6.415D-3,0.,0.,0.,8.698,2.131D1,3.839D-3,3.328D2,
|
|
+ 3.036D-2,0.,0.,0.,1.701D-4,1.715D1,1.641D1,9.531D1,
|
|
+ 1.506D2,2.071D-2,0.,0.,0.,4.236D-2,1.52,7.654D1,5.149,
|
|
+ 9.169D1,1.042D2,2.536D-2,0.,0.,0.,0.,3.375,3.725,1.383D2,
|
|
+ 2.204D1,9.603D-4,9.94D1,0.1667,0.,0.,0.,0.,0.,4.86,
|
|
+ 9.402D-1,1.276D2,1.007D2,5.452D1,5.348D2,5.444,0.,0.,0.,
|
|
+ 2.041D-1,9.149D1,0.,3.994D1,3.495,5.675D-2,3.552D1,
|
|
+ 5.704D2,1.719D2,2.538D-2,0.,0.,0.,1.672D-5,6.634,9.613D1,
|
|
+ 0.,2.009D2,1.115,1.149D-2,3.85D2,1.036D3,4.528D2,1.99D-2,
|
|
+ 0.,0.,0./
|
|
DATA Y1/0.,2.136,0.,0.,0.,0.,1.269D-2,0.,0.,0.,4.556,3.4D-2,0.,
|
|
+ 0.,0.,1.607,3.234,2.528D-2,0.,0.,0.,2.325,2.03D1,3.999,
|
|
+ 1.937D-2,0.,0.,0.,1.271D-1,1.503D-2,4.569D-1,4.285D1,
|
|
+ 5.554D-2,0.,0.,0.,1.345D-2,7.724D-1,5.124,9.781,2.574D-1,
|
|
+ 3.998D-2,0.,0.,0.,5.873,1.084D-1,2.023,6.687,3.702D-1,
|
|
+ 1.435,4.417D-2,0.,0.,0.,0.,4.01,0.2279,4.26,0.7577,
|
|
+ 6.378D-3,3.278,1.766D-2,0.,0.,0.,0.,0.,3.722,0.1135,
|
|
+ 3.979,1.729,2.078,3.997D-3,7.918D-2,0.,0.,0.,0.4753,
|
|
+ 0.6565,0.,4.803,0.2701,0.2768,8.223D-3,0.158,6.595,
|
|
+ 1.203D-2,0.,0.,0.,0.4207,0.1272,0.6442,0.,4.537,8.051D-2,
|
|
+ 0.6396,8.999D-2,0.2936,1.015,1.007D-2,0.,0.,0./
|
|
DATA YW/0.,2.039,0.,0.,0.,0.,3.655D-1,0.,0.,0.,1.887D1,3.503,0.,
|
|
+ 0.,0.,9.157D-2,2.093,5.922,0.,0.,0.,9.097D-2,1.043D1,
|
|
+ 1.85,7.904,0.,0.,0.,7.589D-2,1.934D-2,1.191D-1,3.143,
|
|
+ 2.836D1,0.,0.,0.,2.17D-3,5.103D-1,1.115,6.812,2.57D-1,
|
|
+ 2.411D1,0.,0.,0.,2.434D-1,6.558D-2,4.633D-1,8.29D-3,
|
|
+ 6.855D-1,1.656,2.811D1,0.,0.,0.,0.,2.328,8.579D-2,0.7365,
|
|
+ 0.9275,1.285D-2,1.895,0.9121,0.,0.,0.,0.2805,0.,2.604,
|
|
+ 4.326D-2,0.809,0.6325,1.422,0.6666,2.751,0.,0.,0.,0.3166,
|
|
+ 0.4615,0.,5.342,0.1,8.839,1.836,0.2773,0.6945,29.53,0.,
|
|
+ 0.,0.,0.2837,1.57,0.8626,0.,1.303,2.989D-3,3.28,1.476D-3,
|
|
+ 1.646,0.4489,2.392D1,0.,0.,0./
|
|
DATA YA/3.288D1,1.469D0,3.288D1,1.501D1,7.391D0,3.288D1,4.301D-1,
|
|
+ 5.645D0,1.719D1,3.288D1,8.618D0,1.783D0,1.767D1,2.163D1,
|
|
+ 3.288D1,6.216D1,1.261D2,1.737D0,1.436D1,2.183D1,3.288D1,
|
|
+ 8.033D1,8.163D2,2.788D2,1.714D0,3.297D1,2.627D1,3.288D1,
|
|
+ 3.784D0,3.175D1,3.852D2,4.931D2,1.792D0,3.210D1,2.535D1,
|
|
+ 3.288D1,2.587D0,1.667D1,5.742D1,1.099D2,1.400D2,1.848D0,
|
|
+ 7.547D1,2.657D1,3.288D1,5.798D0,3.935D0,6.725D1,6.409D2,
|
|
+ 1.066D2,1.350D2,1.788D0,3.628D1,3.352D1,3.288D1,6.148D3,
|
|
+ 8.259D0,3.613D0,1.205D2,2.913D1,3.945D1,2.473D2,3.886D-1,
|
|
+ 3.889D1,5.315D1,3.288D1,2.228D-1,4.341D7,9.779D0,2.895D0,
|
|
+ 9.184D1,9.350D1,6.060D1,5.085D2,1.823D0,3.826D1,2.778D1,
|
|
+ 3.288D1,1.621D3,5.675D2,1.990D6,1.674D1,6.973D0,1.856D1,
|
|
+ 1.338D2,7.216D2,2.568D2,1.653D0,3.432D1,3.567D1,3.288D1,
|
|
+ 2.057D1,1.337D1,1.474D2,1.356D6,1.541D2,1.558D2,4.132D1,
|
|
+ 3.733D2,1.404D3,7.461D2,1.582D0,3.307D1,2.848D1,3.288D1/
|
|
DATA PV/2.963,31.88,2.963,4.895,2.916,2.963,10.52,11.7,3.157,
|
|
+ 2.963,17.28,16.18,9.54,2.624,2.963,5.101,8.578,15.93,
|
|
+ 7.457,2.581,2.963,3.928,8.773,9.156,17.06,6.003,2.315,
|
|
+ 2.963,17.64,8.943,6.822,8.785,26.47,5.495,2.336,2.963,
|
|
+ 7.275,10.5,6.614,9.203,9.718,24.46,6.448,2.255,2.963,
|
|
+ 8.355,7.81,10.05,3.056,8.999,8.836,25.5,5.585,2.002,
|
|
+ 2.963,3.839,7.362,9.803,9.714,8.26,2.832,7.681,8.476,
|
|
+ 5.265,1.678,2.963,15.74,3.61,7.117,6.487,10.06,9.202,
|
|
+ 8.857,9.385,14.44,5.46,2.161,2.963,3.642,9.049,3.477,
|
|
+ 11.8,6.724,20.84,6.204,8.659,8.406,26.55,5.085,1.915,
|
|
+ 2.963,3.546,11.91,13.36,3.353,9.98,2.4,16.06,8.986,
|
|
+ 8.503,8.302,24.25,4.728,2.135,2.963/
|
|
DATA S95/5.475D4,4.47D3,1.369D4,1.564D2,3.201D2,6.083D3,1.306D2,
|
|
+ 9.796D1,5.458D2,3.422D3,9.698D1,1.037D2,8.605D1,2.846D2,
|
|
+ 2.19D3,7.421D1,6.649D1,8.067D1,8.111D1,2.344D2,1.521D3,
|
|
+ 4.748D1,5.002D1,5.235D1,7.046D1,7.304D1,1.519D2,1.117D3,
|
|
+ 3.237D1,3.584D1,3.939D1,4.123D1,5.735D1,6.029D1,1.329D2,
|
|
+ 8.554D2,2.295D1,4.798D1,3.144D1,5.302D1,3.680D1,4.668D1,
|
|
+ 4.890D1,1.039D2,6.759D2,1.664D1,2.783D1,2.943D1,3.027D1,
|
|
+ 3.097D1,3.232D1,5.011D1,3.719D1,6.695D1,5.475D2,2.486D2,
|
|
+ 1.889D1,3.466D1,2.551D1,3.753D1,2.654D1,2.745D1,3.613D1,
|
|
+ 3.850D1,3.995D1,4.525D2,2.023D2,2.049D2,1.877D1,2.212D1,
|
|
+ 2.350D1,3.278D1,3.407D1,2.721D1,3.699D1,3.290D1,6.140D1,
|
|
+ 3.802D2,1.735D2,1.842D2,1.990D2,2.494D1,2.572D1,2.786D1,
|
|
+ 2.653D1,2.989D1,2.857D1,3.223D1,3.272D1,4.036D1,3.239D2,
|
|
+ 1.532D2,1.832D2,2.197D2,2.326D2,3.100D1,2.398D1,2.474D1,
|
|
+ 3.135D1,2.641D1,3.343D1,2.832D1,3.880D1,4.754D1,2.793D2/
|
|
DATA E95/0.4298,5.996,1.72,27.4,20.06,3.871,40.93,47.59,17.6,
|
|
+ 6.879,61.55,59.84,65.92,33.36,10.75,86.55,91.13,83.7,
|
|
+ 84.12,46.24,15.48,127.,124.2,122.,107.,106.,69.43,21.08,
|
|
+ 177.4,169.,162.,159.3,137.7,135.4,87.09,27.54,239.,166.,
|
|
+ 205.5,160.3,192.5,173.7,171.1,113.1,34.85,314.4,245.2,
|
|
+ 238.7,236.,234.3,230.8,188.6,219.3,158.6,43.04,36.55,
|
|
+ 326.6,242.6,284.1,236.,280.6,277.7,246.4,240.6,226.8,
|
|
+ 52.11,49.37,49.4,360.7,334.5,325.9,277.7,274.2,308.2,
|
|
+ 267.1,285.8,204.2,62.03,64.45,61.54,60.16,348.2,343.2,
|
|
+ 326.8,339.6,319.3,330.5,311.7,313.7,273.8,72.81,78.08,
|
|
+ 71.54,66.52,64.82,346.,379.9,375.9,344.2,367.9,335.6,
|
|
+ 359.9,315.,275.2,84.47/
|
|
DATA Y95/3.288D1,2.199D0,3.288D1,3.382D1,7.391D0,3.288D1,1.212D2,
|
|
+ 1.166D3,1.719D1,3.288D1,7.354D1,7.915D1,1.906D2,2.163D1,
|
|
+ 3.288D1,5.498D1,9.609D1,7.471D1,7.459D1,2.183D1,3.288D1,
|
|
+ 1.380D2,9.100D1,9.428D1,5.342D1,5.547D1,2.627D1,3.288D1,
|
|
+ 3.812D2,1.894D2,1.104D2,1.141D2,5.486D1,5.682D1,2.535D1,
|
|
+ 3.288D1,1.257D3,5.000D1,1.230D2,5.000D1,7.933D1,5.876D1,
|
|
+ 6.137D1,2.657D1,3.288D1,2.042D5,6.075D1,6.525D1,6.734D1,
|
|
+ 6.907D1,7.167D1,5.000D1,8.181D1,3.352D1,3.288D1,3.222D2,
|
|
+ 2.527D2,5.000D1,6.977D1,5.000D1,7.237D1,7.512D1,4.968D1,
|
|
+ 5.198D1,5.315D1,3.288D1,1.079D4,4.112D4,1.401D2,7.147D1,
|
|
+ 6.173D1,5.000D1,5.000D1,5.108D1,5.000D1,5.384D1,2.778D1,
|
|
+ 3.288D1,1.131D4,2.404D3,5.082D2,3.260D1,3.376D1,5.000D1,
|
|
+ 3.478D1,5.000D1,3.660D1,5.000D1,4.295D1,3.567D1,3.288D1,
|
|
+ 5.765D6,3.537D2,1.169D2,1.188D2,1.979D1,5.000D1,5.000D1,
|
|
+ 2.260D1,5.000D1,2.487D1,5.000D1,3.034D1,2.848D1,3.288D1/
|
|
DATA YW95/55*0.,0.1465,10*0.,1.463D-2,2.223D-5,10*0.,2.337D-2,
|
|
+ 7.839D-3,2.016D-2,10*0.,2.774D-4,2.87D-4,8.658D-4,
|
|
+ 8.417D-4,10*0./
|
|
DATA P95/2.963,6.098,2.963,1.49,2.916,2.963,1.348,1.022,3.157,
|
|
+ 2.963,1.438,1.436,1.21,2.624,2.963,1.503,1.338,1.442,
|
|
+ 1.428,2.581,2.963,1.252,1.335,1.335,1.552,1.538,2.315,
|
|
+ 2.963,1.083,1.185,1.289,1.287,1.540,1.533,2.336,2.963,
|
|
+ 0.9638,1.65,1.263,1.65,1.377,1.511,1.501,2.255,2.963,
|
|
+ 0.845,1.42,1.424,1.422,1.416,1.411,1.65,1.396,2.002,
|
|
+ 2.963,3.57,1.12,1.65,1.414,1.65,1.405,1.397,1.579,
|
|
+ 1.575,1.678,2.963,2.96,2.995,1.238,1.413,1.464,1.65,
|
|
+ 1.65,1.539,1.65,1.56,2.161,2.963,2.762,2.92,3.11,1.732,
|
|
+ 1.73,1.65,1.721,1.65,1.712,1.65,1.676,1.915,2.963,2.639,
|
|
+ 3.133,3.529,3.547,2.094,1.65,1.65,2.02,1.65,1.982,1.65,
|
|
+ 1.925,2.135,2.963/
|
|
C
|
|
VERNER=0.
|
|
E=HHEV*FR
|
|
C
|
|
II=ILOW(ITR)
|
|
N1=NFIRST(IEL(II))
|
|
IF(II.NE.N1)
|
|
+ CALL QUIT('Verner fits only for ground states',ii,n1)
|
|
IAT=NUMAT(IATM(II))
|
|
IZZ=IZ(IEL(II))
|
|
IF(IAT.GT.14) GO TO 10
|
|
IVER=IV0(IAT)+IZZ
|
|
C
|
|
C 1996 Expression
|
|
C
|
|
IF(E.LT.EMX(IVER)) THEN
|
|
XX=E/E0(IVER)-Y0(IVER)
|
|
YY=SQRT(XX*XX+Y1(IVER)*Y1(IVER))
|
|
AA=(XX-UN)*(XX-UN)+YW(IVER)*YW(IVER)
|
|
BB=YY**(HALF*PV(IVER)-5.5)
|
|
CC=(UN+SQRT(YY/YA(IVER)))**PV(IVER)
|
|
FY=AA*BB/CC
|
|
VERNER=S0(IVER)*T18*FY
|
|
ELSE
|
|
C
|
|
C 1995 Expression for high energies
|
|
C (ionization of inner shell electron)
|
|
C
|
|
YY=E/E95(IVER)
|
|
XL=0.
|
|
IF((IAT-IZZ).GE.10) XL=UN
|
|
Q=HALF*P95(IVER)-5.5-XL
|
|
AA=(YY-UN)*(YY-UN)+YW95(IVER)*YW95(IVER)
|
|
BB=YY**Q
|
|
CC=(UN+SQRT(YY/Y95(IVER)))**P95(IVER)
|
|
FY=AA*BB/CC
|
|
VERNER=S95(IVER)*T18*FY
|
|
END IF
|
|
RETURN
|
|
C
|
|
C Heavier elements
|
|
C
|
|
10 IF(IAT.EQ.26) THEN
|
|
VERNER=VERN26(E,IZZ)
|
|
ELSE IF(IAT.EQ.16) THEN
|
|
VERNER=VERN16(E,IZZ)
|
|
ELSE IF(IAT.EQ.18) THEN
|
|
VERNER=VERN18(E,IZZ)
|
|
ELSE IF(IAT.EQ.20) THEN
|
|
VERNER=VERN20(E,IZZ)
|
|
ELSE
|
|
CALL QUIT('VERNER - No data for this element',iat,izz)
|
|
ENDIF
|
|
C
|
|
RETURN
|
|
END
|