Mg XI He-like complex using hewind - isoporous fitting 14Oct08 atomdb says: 9.1687 9.2282 9.2312 9.3143 and I'm going to use (for now) Maurice's phiratio value of 131 same continuum as for non-porous model: n=2; norm=1.55e-3 and same range for fitting the lines: 9.08:9.40 Reproducing the non-porous model: OK Now, setting the rosseland switch to 1 and thawing h: ======================================================================== Model hewind<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 hewind q 0.0 frozen 2 1 hewind taustar 0.897503 +/- 0.981677 3 1 hewind u0 0.643834 +/- 3.25393E-02 4 1 hewind h 1.28901E-04 +/- 9.89678 5 1 hewind tau0star 0.0 frozen 6 1 hewind beta 1.00000 frozen 7 1 hewind betaSob 0.0 frozen 8 1 hewind G 0.722239 +/- 5.32577E-02 9 1 hewind numerica 0 frozen 10 1 hewind anisotro 0 frozen 11 1 hewind rosselan 1 frozen 12 1 hewind expansio 0 frozen 13 1 hewind thick 0 frozen 14 1 hewind Z 12 frozen 15 1 hewind shift "mA" 0.0 frozen 16 1 hewind velocity "km/s" 2250.00 frozen 17 1 hewind phiratio 131.000 frozen 18 1 hewind verbose 0 frozen 19 1 hewind norm 1.80923E-04 +/- 9.81199E-06 20 2 powerlaw PhoIndex 2.00000 frozen 21 2 powerlaw norm 1.55000E-03 frozen ________________________________________________________________________ C-statistic = 410.03 using 392 PHA bins and 387 degrees of freedom. XSPEC12>goodness 1000 nosim 39.80% of realizations are < best fit statistic 410.03 Best-fit h~0 (note a slightly (~0.3) better C_min value; we found this, too, for the non-porous fit; but didn't change the log file, since the model parameters were hardly any different at all). XSPEC12>error 1. 4 Parameter Confidence Range (1.0000000) ***Warning: Parameter pegged at hard limit: 0 Apparent non-monotonicity in statistic space detected. Current bracket values 0.0000087, 1.7081864 and delta stat 0.0000000, 1.0541892 but latest trial 1.6373960 gives 1.0708481 Suggest that you check this result using the steppar command. 4 0.000000 0.854098 (-0.000009,0.854089) This turns out to be wrong; see below: XSPEC12>steppar 4 0.0 1.5 30 C-Statistic Delta h C-Statistic 4 410.26740 -0.06137 0 0.00000 410.27930 -0.04947 1 0.05000 410.29178 -0.03699 2 0.10000 410.30467 -0.02410 3 0.15000 410.31797 -0.01080 4 0.20000 410.33170 0.00293 5 0.25000 410.34587 0.01710 6 0.30000 410.36049 0.03172 7 0.35000 410.37557 0.04680 8 0.40000 410.39113 0.06236 9 0.45000 410.40717 0.07840 10 0.50000 410.42372 0.09495 11 0.55000 410.44079 0.11202 12 0.60000 410.45838 0.12961 13 0.65000 410.47651 0.14774 14 0.70000 410.49520 0.16643 15 0.75000 410.51447 0.18570 16 0.80000 410.53432 0.20555 17 0.85000 410.55477 0.22600 18 0.90000 410.57584 0.24707 19 0.95000 410.59754 0.26877 20 1.00000 410.61989 0.29112 21 1.05000 410.64290 0.31413 22 1.10000 410.66660 0.33783 23 1.15000 410.69100 0.36223 24 1.20000 410.71611 0.38734 25 1.25000 410.74195 0.41318 26 1.30000 410.76854 0.43977 27 1.35000 410.79590 0.46713 28 1.40000 410.82404 0.49527 29 1.45000 410.85299 0.52422 30 1.50000 XSPEC12>steppar 4 1.5 3.0 30 C-Statistic Delta h C-Statistic 4 410.98883 0.66245 0 1.50000 410.88275 0.55637 1 1.55000 410.91335 0.58697 2 1.60000 410.94480 0.61842 3 1.65000 410.97712 0.65075 4 1.70000 411.01034 0.68396 5 1.75000 411.04446 0.71808 6 1.80000 411.07950 0.75312 7 1.85000 411.11548 0.78910 8 1.90000 411.15242 0.82604 9 1.95000 411.19033 0.86395 10 2.00000 411.22923 0.90285 11 2.05000 411.26914 0.94276 12 2.10000 411.31006 0.98368 13 2.15000 411.35203 1.02565 14 2.20000 411.39504 1.06866 15 2.25000 411.43912 1.11274 16 2.30000 411.48427 1.15789 17 2.35000 411.53052 1.20414 18 2.40000 411.57787 1.25149 19 2.45000 411.62634 1.29996 20 2.50000 411.67593 1.34955 21 2.55000 411.72666 1.40028 22 2.60000 411.77854 1.45216 23 2.65000 411.83157 1.50519 24 2.70000 411.88576 1.55938 25 2.75000 411.94113 1.61475 26 2.80000 411.99767 1.67129 27 2.85000 412.05539 1.72901 28 2.90000 412.11429 1.78791 29 2.95000 412.21373 1.88735 30 3.00000 This isn't even up to the 90% limit of Delta-C = 2.7