Re: problem with kerrdisk model

Hi,

xspec does indeed seem to have a more liberal understanding (which I
don't know) of a parameter's hard limits.

However, isis' notion of hard limits (as really limiting a parameter's
reasonable range) has always made much more sense to me, such that I
decided to patch xspec's model.dat with the script attached below.

Cheers,

Manfred


P.S. I hope that you are using isis-1.*6*.2-x.
     Version 1.4.2 would be quite ancient...! ;-)


>     The kerrdisk model does not work in isis-1.4.2-4 and 1.4.2-7.
> 
> isis> model ("kerrdisk");
> Error: invalid range for parameter 2: value=6.4 is outside [0.1,0.1)
> Error: function definition failed
> isis> list_par ;
> kerrdisk
> 
>  Looks like a problem in lmodel.dat but it seems to work well with xspec12.


::::::::::::::
patch_xspec_model.dat
::::::::::::::
--- model.dat
+++ model.dat
_at_email.domain.hidden Rout     Rg     1e5   1e4    1e4      1e8      1e8        -1.

 kerrdisk       9  0.         1.e20            spin     add  0
-lineE   keV     6.4   -0.1    0.1     0.1       100.0     100.0
+lineE   keV     6.4   -0.1    0.1     10.       100.0     100.0
 Index1   " "    3.    -10.    -10.    10.       10.       -0.1
 Index2   " "    3.    -10.    -10.    10.       10.       -0.1
 r_br(g) " "    6.0   1.0     1.0     400.      400.      -0.1
_at_email.domain.hidden FeAbun   " "    1.0   0.      0.    10.    10.     -0.01
 FeKedge  KeV    7.11  7.      7.    10.    10.     -0.01
 PhoIndex " "    2.    -3.    -2.     9.        10.       0.01
-HighECut keV    95.   0.01    1.      50.       200.     -0.01
+HighECut keV    25.   0.01    1.      50.       200.     -0.01
 foldE    " "  100.  1.      1.      1.e6      1.e6     -10.
 acrit    " "    1.    0.0     0.0     1.0       1.0      -1.e-2
 *FAST    " "  0

::::::::::::::
patch_xspec_model.dat_apply
::::::::::::::
#!/bin/sh

# This script must be run in the directory  $HEADAS/.. ,
# where the subdirectories  Xspec  and  spectral  reside.

# It assumes that patch_xspec_model.dat is in the same directory:
patch=`dirname $0`/patch_xspec_model.dat
echo "patch: $patch"

echo

cd Xspec/src/manager
echo "in directory "`pwd`
patch < $patch
cd ../../..

echo

cd spectral/manager
echo "in directory "`pwd`
patch < $patch
cd ../..

----
You received this message because you are
subscribed to the isis-users list.
To unsubscribe, send a message to
isis-users-request_at_email.domain.hiddenwith the first line of the message as:
unsubscribe