FUNCTION SBFOH(FR,T) C ==================== C C cross-section times partition function for OH C C from Kurucz ATLAS9 C INCLUDE 'PARAMS.FOR' parameter (fihu=500.,fihui=1./fihu, * twhu=200.,twhui=1./twhu, * tenl=2.30258509299405E0) DIMENSION CROSSOH(15,130),PARTOH(41),CROSSOHT(15) DIMENSION C1(150),C2(150),C3(150),C4(150),C5(150) DIMENSION C6(150),C7(150),C8(150),C9(150),C10(150) DIMENSION C11(150),C12(150),C13(150) EQUIVALENCE (CROSSOH(1, 1),C1(1)),(CROSSOH(1,11),C2(1)) EQUIVALENCE (CROSSOH(1,21),C3(1)),(CROSSOH(1,31),C4(1)) EQUIVALENCE (CROSSOH(1,41),C5(1)),(CROSSOH(1,51),C6(1)) EQUIVALENCE (CROSSOH(1,61),C7(1)),(CROSSOH(1,71),C8(1)) EQUIVALENCE (CROSSOH(1,81),C9(1)),(CROSSOH(1,91),C10(1)) EQUIVALENCE (CROSSOH(1,101),C11(1)) EQUIVALENCE (CROSSOH(1,111),C12(1)) EQUIVALENCE (CROSSOH(1,121),C13(1)) C DATA C1/-30.855,-29.121,-27.976,-27.166,-26.566,-26.106,-25.742, 1-25.448,-25.207,-25.006,-24.836,-24.691,-24.566,-24.457,-24.363, 2 -30.494,-28.760,-27.615,-26.806,-26.206,-25.745,-25.381, 2-25.088,-24.846,-24.645,-24.475,-24.330,-24.205,-24.097,-24.002, 3 -30.157,-28.425,-27.280,-26.472,-25.872,-25.411,-25.048, 3-24.754,-24.513,-24.312,-24.142,-23.997,-23.872,-23.764,-23.669, 4 -29.848,-28.117,-26.974,-26.165,-25.566,-25.105,-24.742, 4-24.448,-24.207,-24.006,-23.836,-23.692,-23.567,-23.458,-23.364, 5 -29.567,-27.837,-26.693,-25.885,-25.286,-24.826,-24.462, 5-24.169,-23.928,-23.727,-23.557,-23.412,-23.287,-23.179,-23.084, 6 -29.307,-27.578,-26.436,-25.628,-25.029,-24.569,-24.205, 6-23.912,-23.671,-23.470,-23.300,-23.155,-23.031,-22.922,-22.828, 7 -29.068,-27.341,-26.199,-25.391,-24.792,-24.332,-23.969, 7-23.676,-23.435,-23.234,-23.064,-22.920,-22.795,-22.687,-22.592, 8 -28.820,-27.115,-25.978,-25.172,-24.574,-24.115,-23.752, 8-23.459,-23.218,-23.017,-22.848,-22.703,-22.579,-22.470,-22.376, 9 -28.540,-26.891,-25.768,-24.968,-24.372,-23.914,-23.552, 9-23.259,-23.019,-22.818,-22.649,-22.504,-22.380,-22.272,-22.177, A -28.275,-26.681,-25.574,-24.779,-24.186,-23.729,-23.368, A-23.076,-22.836,-22.636,-22.467,-22.322,-22.198,-22.090,-21.996/ DATA C2/-27.993,-26.470,-25.388,-24.602,-24.014,-23.560,-23.200, 1-22.909,-22.669,-22.470,-22.301,-22.157,-22.033,-21.925,-21.831, 2 -27.698,-26.252,-25.204,-24.433,-23.851,-23.401,-23.043, 2-22.754,-22.515,-22.316,-22.148,-22.005,-21.881,-21.773,-21.679, 3 -27.398,-26.026,-25.019,-24.267,-23.696,-23.251,-22.896, 3-22.609,-22.372,-22.174,-22.007,-21.864,-21.741,-21.634,-21.540, 4 -27.100,-25.791,-24.828,-24.102,-23.543,-23.106,-22.756, 4-22.472,-22.238,-22.041,-21.875,-21.733,-21.611,-21.504,-21.411, 5 -26.807,-25.549,-24.631,-23.933,-23.391,-22.964,-22.621, 5-22.341,-22.109,-21.915,-21.751,-21.610,-21.488,-21.383,-21.290, 6 -26.531,-25.310,-24.431,-23.761,-23.238,-22.823,-22.488, 6-22.214,-21.986,-21.795,-21.633,-21.494,-21.374,-21.269,-21.178, 7 -26.239,-25.066,-24.225,-23.585,-23.082,-22.681,-22.356, 7-22.089,-21.866,-21.679,-21.520,-21.383,-21.265,-21.162,-21.072, 8 -25.945,-24.824,-24.017,-23.405,-22.923,-22.538,-22.223, 8-21.964,-21.748,-21.565,-21.410,-21.276,-21.160,-21.059,-20.970, 9 -25.663,-24.587,-23.810,-23.222,-22.761,-22.391,-22.088, 9-21.838,-21.629,-21.452,-21.300,-21.170,-21.057,-20.958,-20.872, A -25.372,-24.350,-23.603,-23.038,-22.596,-22.241,-21.950, A-21.710,-21.508,-21.337,-21.190,-21.064,-20.954,-20.858,-20.774/ DATA C3/-25.076,-24.111,-23.396,-22.853,-22.429,-22.088,-21.809, 1-21.578,-21.384,-21.220,-21.078,-20.957,-20.851,-20.758,-20.676, 2 -24.779,-23.870,-23.189,-22.669,-22.261,-21.934,-21.667, 2-21.445,-21.259,-21.101,-20.965,-20.848,-20.746,-20.656,-20.578, 3 -24.486,-23.629,-22.983,-22.486,-22.095,-21.781,-21.524, 3-21.311,-21.132,-20.980,-20.850,-20.737,-20.639,-20.553,-20.478, 4 -24.183,-23.382,-22.774,-22.302,-21.928,-21.627,-21.381, 4-21.177,-21.005,-20.859,-20.734,-20.625,-20.531,-20.449,-20.376, 5 -23.867,-23.127,-22.561,-22.116,-21.761,-21.474,-21.238, 5-21.043,-20.878,-20.738,-20.617,-20.513,-20.423,-20.344,-20.274, 6 -23.538,-22.862,-22.340,-21.926,-21.592,-21.320,-21.096, 6-20.909,-20.751,-20.617,-20.502,-20.402,-20.315,-20.239,-20.172, 7 -23.234,-22.604,-22.120,-21.734,-21.422,-21.166,-20.953, 7-20.776,-20.625,-20.497,-20.387,-20.291,-20.208,-20.135,-20.071, 8 -22.934,-22.347,-21.898,-21.541,-21.250,-21.010,-20.811, 8-20.643,-20.500,-20.378,-20.273,-20.182,-20.102,-20.033,-19.971, 9 -22.637,-22.092,-21.676,-21.345,-21.075,-20.853,-20.666, 9-20.508,-20.374,-20.259,-20.159,-20.073,-19.997,-19.931,-19.872, A -22.337,-21.835,-21.452,-21.147,-20.899,-20.693,-20.520, A-20.373,-20.247,-20.139,-20.046,-19.964,-19.892,-19.830,-19.774/ DATA C4/-22.049,-21.584,-21.230,-20.950,-20.721,-20.531,-20.372, 1-20.236,-20.119,-20.019,-19.931,-19.855,-19.788,-19.729,-19.676, 2 -21.768,-21.337,-21.011,-20.754,-20.544,-20.370,-20.223, 2-20.098,-19.991,-19.898,-19.817,-19.746,-19.683,-19.628,-19.579, 3 -21.494,-21.096,-20.796,-20.559,-20.367,-20.208,-20.074, 3-19.960,-19.861,-19.776,-19.701,-19.636,-19.578,-19.527,-19.482, 4 -21.233,-20.861,-20.585,-20.368,-20.193,-20.048,-19.926, 4-19.821,-19.732,-19.654,-19.586,-19.526,-19.473,-19.426,-19.384, 5 -20.983,-20.635,-20.380,-20.181,-20.021,-19.889,-19.778, 5-19.683,-19.602,-19.531,-19.469,-19.415,-19.367,-19.324,-19.286, 6 -20.743,-20.418,-20.182,-19.999,-19.853,-19.733,-19.633, 6-19.547,-19.474,-19.410,-19.354,-19.305,-19.261,-19.223,-19.189, 7 -20.515,-20.210,-19.991,-19.824,-19.690,-19.581,-19.490, 7-19.413,-19.347,-19.290,-19.240,-19.196,-19.157,-19.122,-19.092, 8 -20.297,-20.011,-19.808,-19.654,-19.532,-19.434,-19.352, 8-19.282,-19.223,-19.172,-19.127,-19.088,-19.054,-19.023,-18.996, 9 -20.090,-19.822,-19.633,-19.491,-19.381,-19.291,-19.218, 9-19.156,-19.103,-19.057,-19.018,-18.983,-18.952,-18.925,-18.901, A -19.893,-19.642,-19.467,-19.337,-19.236,-19.155,-19.089, A-19.034,-18.987,-18.946,-18.912,-18.881,-18.854,-18.831,-18.810/ DATA C5/-19.705,-19.472,-19.309,-19.190,-19.098,-19.025,-18.966, 1-18.917,-18.876,-18.840,-18.810,-18.783,-18.760,-18.739,-18.721, 2 -19.527,-19.310,-19.161,-19.051,-18.968,-18.903,-18.851, 2-18.807,-18.771,-18.740,-18.713,-18.690,-18.670,-18.653,-18.637, 3 -19.357,-19.159,-19.022,-18.922,-18.847,-18.789,-18.743, 3-18.704,-18.673,-18.646,-18.623,-18.603,-18.586,-18.571,-18.558, 4 -19.195,-19.016,-18.892,-18.803,-18.736,-18.684,-18.643, 4-18.610,-18.583,-18.560,-18.540,-18.523,-18.509,-18.496,-18.485, 5 -19.042,-18.883,-18.772,-18.693,-18.634,-18.589,-18.553, 5-18.525,-18.501,-18.481,-18.465,-18.451,-18.438,-18.428,-18.419, 6 -18.894,-18.758,-18.662,-18.593,-18.542,-18.503,-18.473, 6-18.448,-18.428,-18.412,-18.398,-18.386,-18.376,-18.367,-18.359, 7 -18.752,-18.639,-18.559,-18.501,-18.458,-18.426,-18.400, 7-18.380,-18.363,-18.350,-18.338,-18.328,-18.320,-18.313,-18.306, 8 -18.611,-18.523,-18.460,-18.415,-18.381,-18.355,-18.334, 8-18.318,-18.304,-18.293,-18.284,-18.276,-18.269,-18.263,-18.258, 9 -18.471,-18.408,-18.362,-18.329,-18.304,-18.285,-18.269, 9-18.257,-18.247,-18.238,-18.231,-18.224,-18.219,-18.214,-18.210, A -18.330,-18.290,-18.261,-18.239,-18.223,-18.211,-18.201, A-18.192,-18.185,-18.179,-18.174,-18.169,-18.165,-18.162,-18.159/ DATA C6/-18.190,-18.168,-18.154,-18.143,-18.135,-18.129,-18.124, 1-18.120,-18.116,-18.112,-18.109,-18.106,-18.104,-18.102,-18.100, 2 -18.055,-18.047,-18.043,-18.042,-18.040,-18.039,-18.039, 2-18.038,-18.037,-18.036,-18.035,-18.034,-18.033,-18.033,-18.032, 3 -17.929,-17.931,-17.935,-17.939,-17.943,-17.946,-17.948, 3-17.950,-17.952,-17.953,-17.955,-17.956,-17.957,-17.958,-17.959, 4 -17.818,-17.826,-17.834,-17.842,-17.849,-17.855,-17.860, 4-17.865,-17.869,-17.872,-17.875,-17.878,-17.881,-17.883,-17.886, 5 -17.724,-17.736,-17.747,-17.758,-17.767,-17.775,-17.782, 5-17.788,-17.793,-17.798,-17.803,-17.807,-17.811,-17.815,-17.819, 6 -17.651,-17.665,-17.678,-17.690,-17.701,-17.710,-17.718, 6-17.725,-17.732,-17.738,-17.744,-17.749,-17.755,-17.760,-17.765, 7 -17.601,-17.615,-17.629,-17.642,-17.653,-17.663,-17.672, 7-17.680,-17.688,-17.695,-17.701,-17.708,-17.714,-17.720,-17.726, 8 -17.572,-17.587,-17.602,-17.614,-17.626,-17.636,-17.645, 8-17.654,-17.662,-17.670,-17.677,-17.684,-17.691,-17.698,-17.704, 9 -17.565,-17.581,-17.595,-17.607,-17.619,-17.629,-17.638, 9-17.647,-17.656,-17.664,-17.671,-17.679,-17.686,-17.693,-17.700, A -17.580,-17.594,-17.608,-17.620,-17.630,-17.640,-17.650, A-17.658,-17.667,-17.675,-17.682,-17.690,-17.697,-17.704,-17.711/ DATA C7/-17.613,-17.626,-17.639,-17.649,-17.659,-17.669,-17.677, 1-17.686,-17.694,-17.701,-17.709,-17.716,-17.723,-17.730,-17.737, 2 -17.663,-17.675,-17.685,-17.695,-17.703,-17.711,-17.719, 2-17.727,-17.734,-17.741,-17.748,-17.755,-17.761,-17.768,-17.774, 3 -17.728,-17.737,-17.745,-17.752,-17.759,-17.766,-17.772, 3-17.778,-17.785,-17.791,-17.797,-17.803,-17.808,-17.814,-17.820, 4 -17.803,-17.809,-17.814,-17.818,-17.823,-17.828,-17.832, 4-17.837,-17.842,-17.847,-17.852,-17.856,-17.861,-17.866,-17.871, 5 -17.884,-17.886,-17.888,-17.889,-17.891,-17.893,-17.896, 5-17.899,-17.902,-17.905,-17.908,-17.912,-17.915,-17.919,-17.922, 6 -17.966,-17.964,-17.961,-17.959,-17.958,-17.958,-17.958, 6-17.959,-17.960,-17.961,-17.963,-17.964,-17.966,-17.968,-17.970, 7 -18.040,-18.034,-18.028,-18.023,-18.019,-18.016,-18.013, 7-18.012,-18.010,-18.010,-18.009,-18.009,-18.009,-18.009,-18.010, 8 -18.096,-18.087,-18.078,-18.071,-18.065,-18.059,-18.055, 8-18.051,-18.047,-18.045,-18.042,-18.040,-18.039,-18.037,-18.036, 9 -18.125,-18.115,-18.105,-18.097,-18.089,-18.082,-18.076, 9-18.070,-18.065,-18.061,-18.057,-18.053,-18.051,-18.048,-18.046, A -18.120,-18.112,-18.103,-18.095,-18.087,-18.079,-18.072, A-18.066,-18.060,-18.055,-18.050,-18.046,-18.042,-18.039,-18.036/ DATA C8/-18.083,-18.078,-18.071,-18.064,-18.057,-18.050,-18.044, 1-18.037,-18.032,-18.026,-18.022,-18.017,-18.014,-18.010,-18.007, 2 -18.025,-18.022,-18.017,-18.012,-18.006,-18.000,-17.994, 2-17.989,-17.984,-17.979,-17.975,-17.971,-17.968,-17.965,-17.963, 3 -17.957,-17.955,-17.952,-17.948,-17.943,-17.938,-17.934, 3-17.929,-17.925,-17.922,-17.918,-17.916,-17.913,-17.911,-17.910, 4 -17.890,-17.889,-17.886,-17.882,-17.879,-17.875,-17.871, 4-17.867,-17.864,-17.862,-17.860,-17.858,-17.857,-17.856,-17.855, 5 -17.831,-17.829,-17.826,-17.822,-17.819,-17.815,-17.812, 5-17.810,-17.807,-17.806,-17.804,-17.803,-17.803,-17.803,-17.803, 6 -17.786,-17.782,-17.777,-17.773,-17.769,-17.766,-17.763, 6-17.761,-17.759,-17.758,-17.757,-17.757,-17.757,-17.758,-17.759, 7 -17.753,-17.747,-17.741,-17.735,-17.731,-17.727,-17.724, 7-17.722,-17.721,-17.720,-17.720,-17.720,-17.721,-17.722,-17.724, 8 -17.733,-17.724,-17.716,-17.709,-17.703,-17.699,-17.696, 8-17.694,-17.693,-17.692,-17.692,-17.693,-17.694,-17.695,-17.697, 9 -17.723,-17.711,-17.700,-17.691,-17.685,-17.680,-17.676, 9-17.674,-17.673,-17.672,-17.673,-17.673,-17.675,-17.676,-17.678, A -17.718,-17.702,-17.689,-17.679,-17.672,-17.667,-17.663, A-17.660,-17.659,-17.659,-17.659,-17.660,-17.661,-17.663,-17.665/ DATA C9/-17.713,-17.695,-17.681,-17.670,-17.662,-17.656,-17.653, 1-17.650,-17.649,-17.649,-17.649,-17.650,-17.651,-17.653,-17.655, 2 -17.705,-17.686,-17.671,-17.660,-17.652,-17.647,-17.643, 2-17.641,-17.640,-17.640,-17.640,-17.641,-17.643,-17.645,-17.647, 3 -17.690,-17.671,-17.657,-17.647,-17.640,-17.635,-17.632, 3-17.630,-17.630,-17.630,-17.631,-17.632,-17.634,-17.636,-17.639, 4 -17.667,-17.649,-17.637,-17.629,-17.623,-17.619,-17.618, 4-17.617,-17.617,-17.618,-17.619,-17.621,-17.623,-17.626,-17.628, 5 -17.635,-17.621,-17.611,-17.605,-17.601,-17.600,-17.599, 5-17.599,-17.601,-17.602,-17.604,-17.607,-17.609,-17.612,-17.615, 6 -17.596,-17.585,-17.579,-17.576,-17.575,-17.575,-17.576, 6-17.578,-17.580,-17.582,-17.585,-17.588,-17.591,-17.595,-17.598, 7 -17.550,-17.544,-17.542,-17.542,-17.544,-17.546,-17.548, 7-17.552,-17.555,-17.558,-17.562,-17.566,-17.570,-17.573,-17.577, 8 -17.501,-17.500,-17.501,-17.504,-17.508,-17.513,-17.517, 8-17.521,-17.526,-17.530,-17.535,-17.539,-17.544,-17.548,-17.553, 9 -17.449,-17.452,-17.457,-17.463,-17.470,-17.476,-17.482, 9-17.488,-17.493,-17.499,-17.504,-17.509,-17.514,-17.519,-17.524, A -17.396,-17.403,-17.412,-17.420,-17.429,-17.437,-17.444, A-17.451,-17.458,-17.464,-17.470,-17.476,-17.481,-17.487,-17.492/ DATAC10/-17.344,-17.355,-17.366,-17.377,-17.387,-17.396,-17.405, 1-17.413,-17.420,-17.427,-17.434,-17.440,-17.446,-17.452,-17.458, 2 -17.295,-17.307,-17.321,-17.333,-17.345,-17.355,-17.365, 2-17.373,-17.382,-17.389,-17.397,-17.404,-17.410,-17.417,-17.423, 3 -17.249,-17.264,-17.278,-17.292,-17.304,-17.316,-17.326, 3-17.335,-17.344,-17.352,-17.360,-17.368,-17.375,-17.382,-17.389, 4 -17.209,-17.225,-17.241,-17.255,-17.268,-17.280,-17.291, 4-17.301,-17.310,-17.319,-17.327,-17.335,-17.343,-17.350,-17.357, 5 -17.177,-17.194,-17.210,-17.225,-17.239,-17.251,-17.262, 5-17.272,-17.282,-17.291,-17.300,-17.308,-17.316,-17.324,-17.331, 6 -17.154,-17.172,-17.189,-17.204,-17.218,-17.230,-17.242, 6-17.252,-17.262,-17.272,-17.280,-17.289,-17.298,-17.306,-17.314, 7 -17.144,-17.162,-17.179,-17.194,-17.208,-17.220,-17.232, 7-17.242,-17.253,-17.262,-17.271,-17.280,-17.289,-17.297,-17.306, 8 -17.146,-17.164,-17.181,-17.196,-17.210,-17.222,-17.234, 8-17.245,-17.255,-17.265,-17.274,-17.283,-17.292,-17.301,-17.309, 9 -17.163,-17.180,-17.197,-17.212,-17.225,-17.237,-17.249, 9-17.260,-17.270,-17.280,-17.289,-17.298,-17.307,-17.316,-17.325, A -17.193,-17.211,-17.227,-17.241,-17.254,-17.266,-17.277, A-17.288,-17.298,-17.308,-17.317,-17.327,-17.336,-17.345,-17.353/ DATAC11/-17.239,-17.256,-17.271,-17.284,-17.297,-17.309,-17.320, 1-17.330,-17.340,-17.350,-17.359,-17.369,-17.378,-17.387,-17.395, 2 -17.299,-17.315,-17.329,-17.342,-17.354,-17.365,-17.376, 2-17.386,-17.396,-17.405,-17.415,-17.424,-17.433,-17.442,-17.451, 3 -17.373,-17.388,-17.402,-17.414,-17.425,-17.436,-17.446, 3-17.456,-17.466,-17.475,-17.484,-17.493,-17.502,-17.511,-17.520, 4 -17.462,-17.476,-17.489,-17.500,-17.511,-17.521,-17.531, 4-17.541,-17.550,-17.559,-17.569,-17.578,-17.587,-17.595,-17.604, 5 -17.567,-17.581,-17.592,-17.603,-17.613,-17.623,-17.632, 5-17.641,-17.651,-17.660,-17.669,-17.678,-17.686,-17.695,-17.704, 6 -17.689,-17.701,-17.712,-17.722,-17.732,-17.741,-17.750, 6-17.759,-17.768,-17.777,-17.786,-17.795,-17.803,-17.812,-17.821, 7 -17.829,-17.840,-17.851,-17.860,-17.869,-17.878,-17.887, 7-17.896,-17.904,-17.913,-17.922,-17.930,-17.939,-17.948,-17.956, 8 -17.988,-18.000,-18.010,-18.019,-18.028,-18.036,-18.045, 8-18.053,-18.062,-18.070,-18.079,-18.087,-18.096,-18.104,-18.112, 9 -18.171,-18.183,-18.192,-18.201,-18.210,-18.218,-18.227, 9-18.235,-18.243,-18.252,-18.260,-18.268,-18.277,-18.285,-18.293, A -18.381,-18.393,-18.403,-18.413,-18.422,-18.430,-18.438, A-18.447,-18.455,-18.463,-18.471,-18.479,-18.487,-18.495,-18.503/ DATAC12/-18.625,-18.638,-18.650,-18.660,-18.669,-18.678,-18.687, 1-18.695,-18.703,-18.711,-18.719,-18.726,-18.734,-18.742,-18.750, 2 -18.912,-18.929,-18.943,-18.955,-18.966,-18.975,-18.984, 2-18.993,-19.001,-19.008,-19.016,-19.023,-19.031,-19.038,-19.045, 3 -19.260,-19.283,-19.303,-19.320,-19.333,-19.345,-19.355, 3-19.364,-19.372,-19.380,-19.387,-19.394,-19.400,-19.407,-19.413, 4 -19.704,-19.740,-19.771,-19.796,-19.816,-19.832,-19.845, 4-19.855,-19.863,-19.870,-19.876,-19.882,-19.887,-19.892,-19.897, 5 -20.339,-20.386,-20.424,-20.454,-20.476,-20.492,-20.502, 5-20.509,-20.513,-20.516,-20.518,-20.520,-20.521,-20.523,-20.524, 6 -21.052,-21.075,-21.093,-21.105,-21.114,-21.120,-21.123, 6-21.125,-21.126,-21.127,-21.128,-21.130,-21.131,-21.133,-21.135, 7 -21.174,-21.203,-21.230,-21.255,-21.278,-21.299,-21.320, 7-21.339,-21.357,-21.375,-21.392,-21.408,-21.424,-21.439,-21.454, 8 -21.285,-21.317,-21.346,-21.372,-21.395,-21.416,-21.435, 8-21.452,-21.468,-21.483,-21.497,-21.511,-21.524,-21.536,-21.548, 9 -21.396,-21.429,-21.459,-21.486,-21.511,-21.532,-21.551, 9-21.569,-21.585,-21.600,-21.614,-21.627,-21.640,-21.652,-21.663, A -21.516,-21.549,-21.580,-21.609,-21.635,-21.658,-21.678, A-21.696,-21.713,-21.728,-21.742,-21.755,-21.767,-21.779,-21.790/ DATAC13/-21.651,-21.681,-21.711,-21.738,-21.763,-21.785,-21.804, 1-21.821,-21.837,-21.851,-21.864,-21.876,-21.887,-21.898,-21.908, 2 -21.810,-21.831,-21.853,-21.874,-21.893,-21.910,-21.925, 2-21.938,-21.950,-21.961,-21.971,-21.980,-21.989,-21.998,-22.006, 3 -22.009,-22.016,-22.026,-22.037,-22.048,-22.058,-22.066, 3-22.074,-22.081,-22.088,-22.094,-22.099,-22.105,-22.111,-22.117, 4 -22.353,-22.317,-22.296,-22.284,-22.276,-22.270,-22.266, 4-22.262,-22.260,-22.258,-22.257,-22.257,-22.257,-22.258,-22.259, 5 -22.705,-22.609,-22.552,-22.515,-22.488,-22.468,-22.451, 5-22.438,-22.427,-22.418,-22.410,-22.405,-22.400,-22.397,-22.395, 6 -22.889,-22.791,-22.731,-22.690,-22.659,-22.634,-22.612, 6-22.594,-22.579,-22.566,-22.555,-22.546,-22.539,-22.533,-22.528, 7 -23.211,-23.109,-23.041,-22.989,-22.945,-22.906,-22.872, 7-22.842,-22.816,-22.793,-22.774,-22.757,-22.743,-22.732,-22.722, 8 -25.312,-24.669,-24.250,-23.959,-23.746,-23.587,-23.463, 8-23.366,-23.288,-23.225,-23.173,-23.131,-23.095,-23.066,-23.041, 9 -25.394,-24.752,-24.333,-24.041,-23.829,-23.669,-23.546, 9-23.449,-23.371,-23.308,-23.256,-23.214,-23.178,-23.149,-23.124, A -25.430,-24.787,-24.369,-24.077,-23.865,-23.705,-23.582, A-23.484,-23.407,-23.344,-23.292,-23.249,-23.214,-23.185,-23.160/ DATA PARTOH/ 1 145.979, 178.033, 211.618, 247.053, 284.584, 324.398, 2 366.639, 411.425, 458.854, 509.012, 561.976, 617.823, 3 676.626, 738.448, 803.363, 871.437, 942.735, 1017.330, 4 1095.284, 1176.654, 1261.510, 1349.898, 1441.875, 1537.483, 5 1636.753, 1739.733, 1846.434, 1956.883, 2071.080, 2189.029, 6 2310.724, 2436.155, 2565.283, 2698.103, 2834.571, 2974.627, 7 3118.242, 3265.366, 3415.912, 3569.837, 3727.077/ DATA FREQ1/0./ C SBFOH=0. IF(FR.EQ.FREQ1) GO TO 30 FREQ1=FR WAVENO=FR/2.99792458E10 EVOLT=WAVENO/8065.479 N=int(EVOLT*10.-20.) EN=FLOAT(N)*.1+2. IF(N.LE.0) RETURN IF(N.GE.130) RETURN DO IT=1,15 CROSSOHT(IT)=(CROSSOH(IT,N)+(CROSSOH(IT,N+1)-CROSSOH(IT,N))* * (EVOLT-EN)*10.) END DO c c interpolate to obtain partition function c 30 IF(T.GE.9000.) RETURN IF(N.LE.0) RETURN IF(N.GE.130) RETURN IT=int((T-1000.)*twhui+1.) IF(IT.LT.1) IT=1 TN=FLOAT(IT)*twhu+800. PART=PARTOH(IT)+(PARTOH(IT+1)-PARTOH(IT))*(T-TN)*twhui c c interpolate to obtain cross-section c IT=int((T-2000.)*fihui+1.) IF(IT.LT.1) IT=1 TN=FLOAT(IT)*fihu+1500. SBFOH=EXP((CROSSOHT(IT)+(CROSSOHT(IT+1)-CROSSOHT(IT))* c * (T-TN)*fihui)*tenl)*PART * (T-TN)*fihui)*tenl) RETURN END