diff --git a/src/sme/sme_synth_faster.cpp b/src/sme/sme_synth_faster.cpp index 7f65054..3baf332 100644 --- a/src/sme/sme_synth_faster.cpp +++ b/src/sme/sme_synth_faster.cpp @@ -606,6 +606,11 @@ extern "C" char const *SME_DLL InputLineList(int n, void *arg[]) /* Read in line GWLG10 = 0.; ANSTEE[LINE] = 1; } + else + { + // GAMVW between 0 and 10 + ANSTEE[LINE] = 0; + } } else /* For hydrogen lines GAMQST & GAMVW have special meaning */ { @@ -866,7 +871,8 @@ extern "C" char const *SME_DLL InputModel(int n, void *arg[]) /* Read in model a arg_offset = 1; RADIUS = *(double *)arg[5]; } - else { + else + { printf("ERROR: MOTYPE not recognized, got: %s", motype); } @@ -1997,7 +2003,7 @@ void HMINOP_old(double *ahmin, int iH1, int iHmin) void HMINOP(double *ahmin, int iH1, int iHmin) { - //From Mathisen (1984), after Wishart (1979) and Broad & Reinhardt (1976) + // From Mathisen (1984), after Wishart (1979) and Broad & Reinhardt (1976) static double WBF[85] = {18.00, 19.60, 21.40, 23.60, 26.40, 29.80, 34.30, 40.40, 49.10, 62.60, 111.30, 112.10, 112.67, 112.95, 113.05, 113.10, 113.20, 113.23, 113.50, 114.40, 121.00, 139.00, 164.00, @@ -2048,8 +2054,8 @@ void HMINOP(double *ahmin, int iH1, int iHmin) {3.04, 3.65, 4.80, 5.86, 6.86, 7.79, 8.67, 9.50, 10.3, 13.2, 15.6}, // 30378 {6.79, 8.16, 10.7, 13.1, 15.3, 17.4, 19.4, 21.2, 23.0, 29.5, 35.0}, // 45567 {27.0, 32.4, 42.6, 51.9, 60.7, 68.9, 76.8, 84.2, 91.4, 117., 140.}, // 91134 - {42.3, 50.6, 66.4, 80.8, 94.5, 107., 120., 131., 142., 183., 219.}, //113918 - {75.1, 90.0, 118., 144., 168., 191., 212., 234., 253., 325., 388.}}; //151890 + {42.3, 50.6, 66.4, 80.8, 94.5, 107., 120., 131., 142., 183., 219.}, // 113918 + {75.1, 90.0, 118., 144., 168., 191., 212., 234., 253., 325., 388.}}; // 151890 double WFFLOG[22], FFLOG[11][22], FFTT[11], THETA[MOSIZE], FFTHETA[MOSIZE]; double WAVE[1], WAVELOG[1], XHMIN[MOSIZE], FFTLOG[1], H, HMINBF[1], HMINFF; @@ -3354,7 +3360,7 @@ double SI1OP_new(int J) /* Cross-section */ while (1) { - //Si II 3s2 3p 2P average + // Si II 3s2 3p 2P average ELIM = 65939.18e0; // 3s2 3p4d 3P @@ -3536,7 +3542,7 @@ double SI1OP_new(int J) /* Cross-section */ break; } - //Si II 3s2 3p 2P1/2 + // Si II 3s2 3p 2P1/2 ELIM = 65747.55e0; while (1) @@ -3657,7 +3663,7 @@ double SI1OP_new(int J) /* Cross-section */ break; } - //Si II 3s 3p2 4P1/2 + // Si II 3s 3p2 4P1/2 ELIM = 65747.5e0 + 42824.35e0; while (1) @@ -3730,9 +3736,9 @@ double SI1OP_new(int J) /* Cross-section */ double FE1OP(int J) { - /* - Cross-section time partition functions - This routine is based on R.L. Kurucz Atlas12 + /* + Cross-section time partition functions + This routine is based on R.L. Kurucz Atlas12 */ static double G[48] = {25., 35., 21., 15., 9., 35., 33., 21., 27., 49., 9., 21., 27., 9., 9., 25., 33., 15., 35., 3., 5., 11., 15., 13., @@ -3853,7 +3859,7 @@ double FE1OP_new(int J) BOLT = exp(-Ebin[i_en] / kT_eV); fe1op += ((GCROSS[i_wn + 1][i_en] - GCROSS[i_wn][i_en]) * FACTOR + GCROSS[i_wn][i_en]) * BOLT; } - return fe1op; ///PARTITION_FUNCTIONS[J][IXFE1]; + return fe1op; /// PARTITION_FUNCTIONS[J][IXFE1]; } double CHOP(int J) /* Cross-section for CH molecule */ @@ -4057,18 +4063,18 @@ double CHOP(int J) /* Cross-section for CH molecule */ -20.822, -20.782, -20.750, -20.724, -20.704, -20.687, -20.674, -20.663}, // 9.8 {-22.319, -21.937, -21.686, -21.510, -21.380, -21.282, -21.206, // 9.9 -21.147, -21.099, -21.061, -21.031, -21.006, -20.985, -20.968, -20.954}, // 9.9 - {-22.969, -22.561, -22.288, -22.092, -21.945, -21.832, -21.743, //10.0 - -21.672, -21.616, -21.570, -21.533, -21.503, -21.477, -21.457, -21.439}, //10.0 - {-24.001, -23.527, -23.199, -22.957, -22.772, -22.629, -22.516, //10.1 - -22.427, -22.355, -22.297, -22.250, -22.212, -22.180, -22.153, -22.131}, //10.1 - {-24.233, -23.774, -23.477, -23.273, -23.128, -23.022, -22.943, //10.2 - -22.883, -22.837, -22.802, -22.774, -22.752, -22.735, -22.721, -22.710}, //10.2 - {-24.550, -23.913, -23.521, -23.266, -23.094, -22.976, -22.893, //10.3 - -22.836, -22.796, -22.768, -22.750, -22.737, -22.730, -22.726, -22.725}, //10.3 - {-24.301, -23.665, -23.274, -23.019, -22.848, -22.730, -22.648, //10.4 - -22.591, -22.552, -22.525, -22.507, -22.495, -22.489, -22.485, -22.485}, //10.4 - {-24.519, -23.883, -23.491, -23.237, -23.065, -22.948, -22.866, //10.5 - -22.809, -22.770, -22.743, -22.724, -22.713, -22.706, -22.703, -22.702}}; //10.5 + {-22.969, -22.561, -22.288, -22.092, -21.945, -21.832, -21.743, // 10.0 + -21.672, -21.616, -21.570, -21.533, -21.503, -21.477, -21.457, -21.439}, // 10.0 + {-24.001, -23.527, -23.199, -22.957, -22.772, -22.629, -22.516, // 10.1 + -22.427, -22.355, -22.297, -22.250, -22.212, -22.180, -22.153, -22.131}, // 10.1 + {-24.233, -23.774, -23.477, -23.273, -23.128, -23.022, -22.943, // 10.2 + -22.883, -22.837, -22.802, -22.774, -22.752, -22.735, -22.721, -22.710}, // 10.2 + {-24.550, -23.913, -23.521, -23.266, -23.094, -22.976, -22.893, // 10.3 + -22.836, -22.796, -22.768, -22.750, -22.737, -22.730, -22.726, -22.725}, // 10.3 + {-24.301, -23.665, -23.274, -23.019, -22.848, -22.730, -22.648, // 10.4 + -22.591, -22.552, -22.525, -22.507, -22.495, -22.489, -22.485, -22.485}, // 10.4 + {-24.519, -23.883, -23.491, -23.237, -23.065, -22.948, -22.866, // 10.5 + -22.809, -22.770, -22.743, -22.724, -22.713, -22.706, -22.703, -22.702}}; // 10.5 double WAVENO, EVOLT, EN, TN, CROSSCHT[15], CHop; int N, IT; @@ -4359,108 +4365,108 @@ double OHOP(int J) -17.694, -17.693, -17.692, -17.692, -17.693, -17.694, -17.695, -17.697}, // 9.8 {-17.723, -17.711, -17.700, -17.691, -17.685, -17.680, -17.676, // 9.9 -17.674, -17.673, -17.672, -17.673, -17.673, -17.675, -17.676, -17.678}, // 9.9 - {-17.718, -17.702, -17.689, -17.679, -17.672, -17.667, -17.663, //10.0 - -17.660, -17.659, -17.659, -17.659, -17.660, -17.661, -17.663, -17.665}, //10.0 - {-17.713, -17.695, -17.681, -17.670, -17.662, -17.656, -17.653, //10.1 - -17.650, -17.649, -17.649, -17.649, -17.650, -17.651, -17.653, -17.655}, //10.1 - {-17.705, -17.686, -17.671, -17.660, -17.652, -17.647, -17.643, //10.2 - -17.641, -17.640, -17.640, -17.640, -17.641, -17.643, -17.645, -17.647}, //10.2 - {-17.690, -17.671, -17.657, -17.647, -17.640, -17.635, -17.632, //10.3 - -17.630, -17.630, -17.630, -17.631, -17.632, -17.634, -17.636, -17.639}, //10.3 - {-17.667, -17.649, -17.637, -17.629, -17.623, -17.619, -17.618, //10.4 - -17.617, -17.617, -17.618, -17.619, -17.621, -17.623, -17.626, -17.628}, //10.4 - {-17.635, -17.621, -17.611, -17.605, -17.601, -17.600, -17.599, //10.5 - -17.599, -17.601, -17.602, -17.604, -17.607, -17.609, -17.612, -17.615}, //10.5 - {-17.596, -17.585, -17.579, -17.576, -17.575, -17.575, -17.576, //10.6 - -17.578, -17.580, -17.582, -17.585, -17.588, -17.591, -17.595, -17.598}, //10.6 - {-17.550, -17.544, -17.542, -17.542, -17.544, -17.546, -17.548, //10.7 - -17.552, -17.555, -17.558, -17.562, -17.566, -17.570, -17.573, -17.577}, //10.7 - {-17.501, -17.500, -17.501, -17.504, -17.508, -17.513, -17.517, //10.8 - -17.521, -17.526, -17.530, -17.535, -17.539, -17.544, -17.548, -17.553}, //10.8 - {-17.449, -17.452, -17.457, -17.463, -17.470, -17.476, -17.482, //10.9 - -17.488, -17.493, -17.499, -17.504, -17.509, -17.514, -17.519, -17.524}, //10.9 - {-17.396, -17.403, -17.412, -17.420, -17.429, -17.437, -17.444, //11.0 - -17.451, -17.458, -17.464, -17.470, -17.476, -17.481, -17.487, -17.492}, //11.0 - {-17.344, -17.355, -17.366, -17.377, -17.387, -17.396, -17.405, //11.1 - -17.413, -17.420, -17.427, -17.434, -17.440, -17.446, -17.452, -17.458}, //11.1 - {-17.295, -17.307, -17.321, -17.333, -17.345, -17.355, -17.365, //11.2 - -17.373, -17.382, -17.389, -17.397, -17.404, -17.410, -17.417, -17.423}, //11.2 - {-17.249, -17.264, -17.278, -17.292, -17.304, -17.316, -17.326, //11.3 - -17.335, -17.344, -17.352, -17.360, -17.368, -17.375, -17.382, -17.389}, //11.3 - {-17.209, -17.225, -17.241, -17.255, -17.268, -17.280, -17.291, //11.4 - -17.301, -17.310, -17.319, -17.327, -17.335, -17.343, -17.350, -17.357}, //11.4 - {-17.177, -17.194, -17.210, -17.225, -17.239, -17.251, -17.262, //11.5 - -17.272, -17.282, -17.291, -17.300, -17.308, -17.316, -17.324, -17.331}, //11.5 - {-17.154, -17.172, -17.189, -17.204, -17.218, -17.230, -17.242, //11.6 - -17.252, -17.262, -17.272, -17.280, -17.289, -17.298, -17.306, -17.314}, //11.6 - {-17.144, -17.162, -17.179, -17.194, -17.208, -17.220, -17.232, //11.7 - -17.242, -17.253, -17.262, -17.271, -17.280, -17.289, -17.297, -17.306}, //11.7 - {-17.146, -17.164, -17.181, -17.196, -17.210, -17.222, -17.234, //11.8 - -17.245, -17.255, -17.265, -17.274, -17.283, -17.292, -17.301, -17.309}, //11.8 - {-17.163, -17.180, -17.197, -17.212, -17.225, -17.237, -17.249, //11.9 - -17.260, -17.270, -17.280, -17.289, -17.298, -17.307, -17.316, -17.325}, //11.9 - {-17.193, -17.211, -17.227, -17.241, -17.254, -17.266, -17.277, //12.0 - -17.288, -17.298, -17.308, -17.317, -17.327, -17.336, -17.345, -17.353}, //12.0 - {-17.239, -17.256, -17.271, -17.284, -17.297, -17.309, -17.320, //12.1 - -17.330, -17.340, -17.350, -17.359, -17.369, -17.378, -17.387, -17.395}, //12.1 - {-17.299, -17.315, -17.329, -17.342, -17.354, -17.365, -17.376, //12.2 - -17.386, -17.396, -17.405, -17.415, -17.424, -17.433, -17.442, -17.451}, //12.2 - {-17.373, -17.388, -17.402, -17.414, -17.425, -17.436, -17.446, //12.3 - -17.456, -17.466, -17.475, -17.484, -17.493, -17.502, -17.511, -17.520}, //12.3 - {-17.462, -17.476, -17.489, -17.500, -17.511, -17.521, -17.531, //12.4 - -17.541, -17.550, -17.559, -17.569, -17.578, -17.587, -17.595, -17.604}, //12.4 - {-17.567, -17.581, -17.592, -17.603, -17.613, -17.623, -17.632, //12.5 - -17.641, -17.651, -17.660, -17.669, -17.678, -17.686, -17.695, -17.704}, //12.5 - {-17.689, -17.701, -17.712, -17.722, -17.732, -17.741, -17.750, //12.6 - -17.759, -17.768, -17.777, -17.786, -17.795, -17.803, -17.812, -17.821}, //12.6 - {-17.829, -17.840, -17.851, -17.860, -17.869, -17.878, -17.887, //12.7 - -17.896, -17.904, -17.913, -17.922, -17.930, -17.939, -17.948, -17.956}, //12.7 - {-17.988, -18.000, -18.010, -18.019, -18.028, -18.036, -18.045, //12.8 - -18.053, -18.062, -18.070, -18.079, -18.087, -18.096, -18.104, -18.112}, //12.8 - {-18.171, -18.183, -18.192, -18.201, -18.210, -18.218, -18.227, //12.9 - -18.235, -18.243, -18.252, -18.260, -18.268, -18.277, -18.285, -18.293}, //12.9 - {-18.381, -18.393, -18.403, -18.413, -18.422, -18.430, -18.438, //13.0 - -18.447, -18.455, -18.463, -18.471, -18.479, -18.487, -18.495, -18.503}, //13.0 - {-18.625, -18.638, -18.650, -18.660, -18.669, -18.678, -18.687, //13.1 - -18.695, -18.703, -18.711, -18.719, -18.726, -18.734, -18.742, -18.750}, //13.1 - {-18.912, -18.929, -18.943, -18.955, -18.966, -18.975, -18.984, //13.2 - -18.993, -19.001, -19.008, -19.016, -19.023, -19.031, -19.038, -19.045}, //13.2 - {-19.260, -19.283, -19.303, -19.320, -19.333, -19.345, -19.355, //13.3 - -19.364, -19.372, -19.380, -19.387, -19.394, -19.400, -19.407, -19.413}, //13.3 - {-19.704, -19.740, -19.771, -19.796, -19.816, -19.832, -19.845, //13.4 - -19.855, -19.863, -19.870, -19.876, -19.882, -19.887, -19.892, -19.897}, //13.4 - {-20.339, -20.386, -20.424, -20.454, -20.476, -20.492, -20.502, //13.5 - -20.509, -20.513, -20.516, -20.518, -20.520, -20.521, -20.523, -20.524}, //13.5 - {-21.052, -21.075, -21.093, -21.105, -21.114, -21.120, -21.123, //13.6 - -21.125, -21.126, -21.127, -21.128, -21.130, -21.131, -21.133, -21.135}, //13.6 - {-21.174, -21.203, -21.230, -21.255, -21.278, -21.299, -21.320, //13.7 - -21.339, -21.357, -21.375, -21.392, -21.408, -21.424, -21.439, -21.454}, //13.7 - {-21.285, -21.317, -21.346, -21.372, -21.395, -21.416, -21.435, //13.8 - -21.452, -21.468, -21.483, -21.497, -21.511, -21.524, -21.536, -21.548}, //13.8 - {-21.396, -21.429, -21.459, -21.486, -21.511, -21.532, -21.551, //13.9 - -21.569, -21.585, -21.600, -21.614, -21.627, -21.640, -21.652, -21.663}, //13.9 - {-21.516, -21.549, -21.580, -21.609, -21.635, -21.658, -21.678, //14.0 - -21.696, -21.713, -21.728, -21.742, -21.755, -21.767, -21.779, -21.790}, //14.0 - {-21.651, -21.681, -21.711, -21.738, -21.763, -21.785, -21.804, //14.1 - -21.821, -21.837, -21.851, -21.864, -21.876, -21.887, -21.898, -21.908}, //14.1 - {-21.810, -21.831, -21.853, -21.874, -21.893, -21.910, -21.925, //14.2 - -21.938, -21.950, -21.961, -21.971, -21.980, -21.989, -21.998, -22.006}, //14.2 - {-22.009, -22.016, -22.026, -22.037, -22.048, -22.058, -22.066, //14.3 - -22.074, -22.081, -22.088, -22.094, -22.099, -22.105, -22.111, -22.117}, //14.3 - {-22.353, -22.317, -22.296, -22.284, -22.276, -22.270, -22.266, //14.4 - -22.262, -22.260, -22.258, -22.257, -22.257, -22.257, -22.258, -22.259}, //14.4 - {-22.705, -22.609, -22.552, -22.515, -22.488, -22.468, -22.451, //14.5 - 5 - 22.438, -22.427, -22.418, -22.410, -22.405, -22.400, -22.397, -22.395}, //14.5 - {-22.889, -22.791, -22.731, -22.690, -22.659, -22.634, -22.612, //14.6 - -22.594, -22.579, -22.566, -22.555, -22.546, -22.539, -22.533, -22.528}, //14.6 - {-23.211, -23.109, -23.041, -22.989, -22.945, -22.906, -22.872, //14.7 - -22.842, -22.816, -22.793, -22.774, -22.757, -22.743, -22.732, -22.722}, //14.7 - {-25.312, -24.669, -24.250, -23.959, -23.746, -23.587, -23.463, //14.8 - -23.366, -23.288, -23.225, -23.173, -23.131, -23.095, -23.066, -23.041}, //14.8 - {-25.394, -24.752, -24.333, -24.041, -23.829, -23.669, -23.546, //14.9 - -23.449, -23.371, -23.308, -23.256, -23.214, -23.178, -23.149, -23.124}, //14.9 - {-25.430, -24.787, -24.369, -24.077, -23.865, -23.705, -23.582, //15.0 - -23.484, -23.407, -23.344, -23.292, -23.249, -23.214, -23.185, -23.160}}; //15.0 + {-17.718, -17.702, -17.689, -17.679, -17.672, -17.667, -17.663, // 10.0 + -17.660, -17.659, -17.659, -17.659, -17.660, -17.661, -17.663, -17.665}, // 10.0 + {-17.713, -17.695, -17.681, -17.670, -17.662, -17.656, -17.653, // 10.1 + -17.650, -17.649, -17.649, -17.649, -17.650, -17.651, -17.653, -17.655}, // 10.1 + {-17.705, -17.686, -17.671, -17.660, -17.652, -17.647, -17.643, // 10.2 + -17.641, -17.640, -17.640, -17.640, -17.641, -17.643, -17.645, -17.647}, // 10.2 + {-17.690, -17.671, -17.657, -17.647, -17.640, -17.635, -17.632, // 10.3 + -17.630, -17.630, -17.630, -17.631, -17.632, -17.634, -17.636, -17.639}, // 10.3 + {-17.667, -17.649, -17.637, -17.629, -17.623, -17.619, -17.618, // 10.4 + -17.617, -17.617, -17.618, -17.619, -17.621, -17.623, -17.626, -17.628}, // 10.4 + {-17.635, -17.621, -17.611, -17.605, -17.601, -17.600, -17.599, // 10.5 + -17.599, -17.601, -17.602, -17.604, -17.607, -17.609, -17.612, -17.615}, // 10.5 + {-17.596, -17.585, -17.579, -17.576, -17.575, -17.575, -17.576, // 10.6 + -17.578, -17.580, -17.582, -17.585, -17.588, -17.591, -17.595, -17.598}, // 10.6 + {-17.550, -17.544, -17.542, -17.542, -17.544, -17.546, -17.548, // 10.7 + -17.552, -17.555, -17.558, -17.562, -17.566, -17.570, -17.573, -17.577}, // 10.7 + {-17.501, -17.500, -17.501, -17.504, -17.508, -17.513, -17.517, // 10.8 + -17.521, -17.526, -17.530, -17.535, -17.539, -17.544, -17.548, -17.553}, // 10.8 + {-17.449, -17.452, -17.457, -17.463, -17.470, -17.476, -17.482, // 10.9 + -17.488, -17.493, -17.499, -17.504, -17.509, -17.514, -17.519, -17.524}, // 10.9 + {-17.396, -17.403, -17.412, -17.420, -17.429, -17.437, -17.444, // 11.0 + -17.451, -17.458, -17.464, -17.470, -17.476, -17.481, -17.487, -17.492}, // 11.0 + {-17.344, -17.355, -17.366, -17.377, -17.387, -17.396, -17.405, // 11.1 + -17.413, -17.420, -17.427, -17.434, -17.440, -17.446, -17.452, -17.458}, // 11.1 + {-17.295, -17.307, -17.321, -17.333, -17.345, -17.355, -17.365, // 11.2 + -17.373, -17.382, -17.389, -17.397, -17.404, -17.410, -17.417, -17.423}, // 11.2 + {-17.249, -17.264, -17.278, -17.292, -17.304, -17.316, -17.326, // 11.3 + -17.335, -17.344, -17.352, -17.360, -17.368, -17.375, -17.382, -17.389}, // 11.3 + {-17.209, -17.225, -17.241, -17.255, -17.268, -17.280, -17.291, // 11.4 + -17.301, -17.310, -17.319, -17.327, -17.335, -17.343, -17.350, -17.357}, // 11.4 + {-17.177, -17.194, -17.210, -17.225, -17.239, -17.251, -17.262, // 11.5 + -17.272, -17.282, -17.291, -17.300, -17.308, -17.316, -17.324, -17.331}, // 11.5 + {-17.154, -17.172, -17.189, -17.204, -17.218, -17.230, -17.242, // 11.6 + -17.252, -17.262, -17.272, -17.280, -17.289, -17.298, -17.306, -17.314}, // 11.6 + {-17.144, -17.162, -17.179, -17.194, -17.208, -17.220, -17.232, // 11.7 + -17.242, -17.253, -17.262, -17.271, -17.280, -17.289, -17.297, -17.306}, // 11.7 + {-17.146, -17.164, -17.181, -17.196, -17.210, -17.222, -17.234, // 11.8 + -17.245, -17.255, -17.265, -17.274, -17.283, -17.292, -17.301, -17.309}, // 11.8 + {-17.163, -17.180, -17.197, -17.212, -17.225, -17.237, -17.249, // 11.9 + -17.260, -17.270, -17.280, -17.289, -17.298, -17.307, -17.316, -17.325}, // 11.9 + {-17.193, -17.211, -17.227, -17.241, -17.254, -17.266, -17.277, // 12.0 + -17.288, -17.298, -17.308, -17.317, -17.327, -17.336, -17.345, -17.353}, // 12.0 + {-17.239, -17.256, -17.271, -17.284, -17.297, -17.309, -17.320, // 12.1 + -17.330, -17.340, -17.350, -17.359, -17.369, -17.378, -17.387, -17.395}, // 12.1 + {-17.299, -17.315, -17.329, -17.342, -17.354, -17.365, -17.376, // 12.2 + -17.386, -17.396, -17.405, -17.415, -17.424, -17.433, -17.442, -17.451}, // 12.2 + {-17.373, -17.388, -17.402, -17.414, -17.425, -17.436, -17.446, // 12.3 + -17.456, -17.466, -17.475, -17.484, -17.493, -17.502, -17.511, -17.520}, // 12.3 + {-17.462, -17.476, -17.489, -17.500, -17.511, -17.521, -17.531, // 12.4 + -17.541, -17.550, -17.559, -17.569, -17.578, -17.587, -17.595, -17.604}, // 12.4 + {-17.567, -17.581, -17.592, -17.603, -17.613, -17.623, -17.632, // 12.5 + -17.641, -17.651, -17.660, -17.669, -17.678, -17.686, -17.695, -17.704}, // 12.5 + {-17.689, -17.701, -17.712, -17.722, -17.732, -17.741, -17.750, // 12.6 + -17.759, -17.768, -17.777, -17.786, -17.795, -17.803, -17.812, -17.821}, // 12.6 + {-17.829, -17.840, -17.851, -17.860, -17.869, -17.878, -17.887, // 12.7 + -17.896, -17.904, -17.913, -17.922, -17.930, -17.939, -17.948, -17.956}, // 12.7 + {-17.988, -18.000, -18.010, -18.019, -18.028, -18.036, -18.045, // 12.8 + -18.053, -18.062, -18.070, -18.079, -18.087, -18.096, -18.104, -18.112}, // 12.8 + {-18.171, -18.183, -18.192, -18.201, -18.210, -18.218, -18.227, // 12.9 + -18.235, -18.243, -18.252, -18.260, -18.268, -18.277, -18.285, -18.293}, // 12.9 + {-18.381, -18.393, -18.403, -18.413, -18.422, -18.430, -18.438, // 13.0 + -18.447, -18.455, -18.463, -18.471, -18.479, -18.487, -18.495, -18.503}, // 13.0 + {-18.625, -18.638, -18.650, -18.660, -18.669, -18.678, -18.687, // 13.1 + -18.695, -18.703, -18.711, -18.719, -18.726, -18.734, -18.742, -18.750}, // 13.1 + {-18.912, -18.929, -18.943, -18.955, -18.966, -18.975, -18.984, // 13.2 + -18.993, -19.001, -19.008, -19.016, -19.023, -19.031, -19.038, -19.045}, // 13.2 + {-19.260, -19.283, -19.303, -19.320, -19.333, -19.345, -19.355, // 13.3 + -19.364, -19.372, -19.380, -19.387, -19.394, -19.400, -19.407, -19.413}, // 13.3 + {-19.704, -19.740, -19.771, -19.796, -19.816, -19.832, -19.845, // 13.4 + -19.855, -19.863, -19.870, -19.876, -19.882, -19.887, -19.892, -19.897}, // 13.4 + {-20.339, -20.386, -20.424, -20.454, -20.476, -20.492, -20.502, // 13.5 + -20.509, -20.513, -20.516, -20.518, -20.520, -20.521, -20.523, -20.524}, // 13.5 + {-21.052, -21.075, -21.093, -21.105, -21.114, -21.120, -21.123, // 13.6 + -21.125, -21.126, -21.127, -21.128, -21.130, -21.131, -21.133, -21.135}, // 13.6 + {-21.174, -21.203, -21.230, -21.255, -21.278, -21.299, -21.320, // 13.7 + -21.339, -21.357, -21.375, -21.392, -21.408, -21.424, -21.439, -21.454}, // 13.7 + {-21.285, -21.317, -21.346, -21.372, -21.395, -21.416, -21.435, // 13.8 + -21.452, -21.468, -21.483, -21.497, -21.511, -21.524, -21.536, -21.548}, // 13.8 + {-21.396, -21.429, -21.459, -21.486, -21.511, -21.532, -21.551, // 13.9 + -21.569, -21.585, -21.600, -21.614, -21.627, -21.640, -21.652, -21.663}, // 13.9 + {-21.516, -21.549, -21.580, -21.609, -21.635, -21.658, -21.678, // 14.0 + -21.696, -21.713, -21.728, -21.742, -21.755, -21.767, -21.779, -21.790}, // 14.0 + {-21.651, -21.681, -21.711, -21.738, -21.763, -21.785, -21.804, // 14.1 + -21.821, -21.837, -21.851, -21.864, -21.876, -21.887, -21.898, -21.908}, // 14.1 + {-21.810, -21.831, -21.853, -21.874, -21.893, -21.910, -21.925, // 14.2 + -21.938, -21.950, -21.961, -21.971, -21.980, -21.989, -21.998, -22.006}, // 14.2 + {-22.009, -22.016, -22.026, -22.037, -22.048, -22.058, -22.066, // 14.3 + -22.074, -22.081, -22.088, -22.094, -22.099, -22.105, -22.111, -22.117}, // 14.3 + {-22.353, -22.317, -22.296, -22.284, -22.276, -22.270, -22.266, // 14.4 + -22.262, -22.260, -22.258, -22.257, -22.257, -22.257, -22.258, -22.259}, // 14.4 + {-22.705, -22.609, -22.552, -22.515, -22.488, -22.468, -22.451, // 14.5 + 5 - 22.438, -22.427, -22.418, -22.410, -22.405, -22.400, -22.397, -22.395}, // 14.5 + {-22.889, -22.791, -22.731, -22.690, -22.659, -22.634, -22.612, // 14.6 + -22.594, -22.579, -22.566, -22.555, -22.546, -22.539, -22.533, -22.528}, // 14.6 + {-23.211, -23.109, -23.041, -22.989, -22.945, -22.906, -22.872, // 14.7 + -22.842, -22.816, -22.793, -22.774, -22.757, -22.743, -22.732, -22.722}, // 14.7 + {-25.312, -24.669, -24.250, -23.959, -23.746, -23.587, -23.463, // 14.8 + -23.366, -23.288, -23.225, -23.173, -23.131, -23.095, -23.066, -23.041}, // 14.8 + {-25.394, -24.752, -24.333, -24.041, -23.829, -23.669, -23.546, // 14.9 + -23.449, -23.371, -23.308, -23.256, -23.214, -23.178, -23.149, -23.124}, // 14.9 + {-25.430, -24.787, -24.369, -24.077, -23.865, -23.705, -23.582, // 15.0 + -23.484, -23.407, -23.344, -23.292, -23.249, -23.214, -23.185, -23.160}}; // 15.0 double WAVENO, EVOLT, EN, TN, CROSSOHT[15], OHop; int N, IT; @@ -5790,7 +5796,7 @@ extern "C" char const *SME_DLL Transf(int n, void *arg[]) TABLE = (double *)arg[7]; /* Array for synthetic spectrum */ EPS1 = *(double *)arg[8]; /* Accuracy of the radiative transfer integration */ EPS2 = *(double *)arg[9]; /* Accuracy of the interpolation on wl grid */ - keep_lineop = *(short *)arg[10]; /* For several spectral segments there is no + keep_lineop = *(short *)arg[10]; /* For several spectral segments there is no point recomputing line opacities. This flag tells when recalculations are needed */ } @@ -6159,7 +6165,7 @@ int RKINTS_sph(double rhox[][2 * MOSIZE], int NMU, int NRHOXs[], double EPS1, do double *FCBLUE, double *FCRED, double *TABLE, int NWSIZE, int &NWL, double *WL, short long_continuum, int grazing[]) { - /* + /* THIS SUBROUTINE CALLS SUBROUTINE FCINTG TO INTEGRATE THE EMMERGING SPECIFIC INTENSITIES FOR CONTINUUM AT THE EDGES OF SPECTRAL INTERVAL (returned as "FC*") AND SUBROUTINE TBINTG FOR THE LINE @@ -7372,7 +7378,7 @@ void LINEOPAC(int LINE) pi*e^2 Line opacity is: ------ * gf * N_absorb * STIM m*c - + The Hydrogen line profiles are computed externally by Kurucz approximation (HLINOP) or by interpolation in Stehle's tables (HTABLE) and are area normalized! @@ -7680,11 +7686,11 @@ void OPMTRX(double WAVE, double *XK, double *XC, double *source_line, pi*e^2 Line opacity is: ------ * gf * N_absorb * STIM * f(wl-wl0) m*c - + where the line profile f(wl) is assumed to be nomalized so that: - + \integ f(wl-wl0) d wl = 1 - + This is true for Voigt, Hydrogen and (I hope) Fano profiles. 1 E.g., in case of Voigt profile f(wl-wl0)= -------------------- * H(a,v) @@ -7693,7 +7699,7 @@ void OPMTRX(double WAVE, double *XK, double *XC, double *source_line, where H(a,v) is the Voigt function with normalization: \integ H(a,v) d v = sqrt(pi) - + Two Hydrogen line profiles are computed externally by Kurucz approximation (HLINOP) or by interpolation in Stehle's tables (HTABLE) and are area normalized! @@ -7946,11 +7952,11 @@ void OPMTRX1(int LINE, double *XK) PER GRAMM OF STELLAR MATER (CM**2/GM) PER ANGSTROEM AT DEPTH #IM OF THE STANDARD MODEL DEPTH SCALE. WAVELENGTH IS THE CENTRAL WAVELENGTH OF LINE "LINE". - + For comments and constants description see OPMTRX above. Author: N.Piskunov - + C++ Version: January 15, 1999 LAST UPDATE: See OPMTRX above */