The following is an example of power law analysis 

'REM'    SAMPLE FILE FOR POWER LAW ANALYSIS 
'SER'    0.01   /SYSTEMATIC ERROR
'CND'    300     3.00      0.001  / MAXSTP, CHISQ-MIN, VARIATION
'BIN'    01111111111111111111111111111110
'MDL'    POWL(1,2)
'PAR'    1          1         0.1         10
'PAR'    2         3.0       2.0         4.0
'OUT'    './HXS_RESULTS.DATA'
'FILE'   './HXS_SPECTRUM.DATA'


The following is an example of broken power law + Gauss analysis 

'REM'      SAMPLE FILE FOR COUTOUR MAP (XAXIS INDEX1,YAXIS INDEX2)
'TRACE'    1                     / 1 OFF   0    ON
'SER'      0.01   /SYSTEMATIC ERROR
'CND'      300     3.00      0.001  / MAXSTP, CHISQ-MIN, VARIATION
'BIN'      01111111111111111111111111111110
'MDL'      BRPW(1,2,3,4) + GAUS (5,6,7)
'PAR'      1          1           0.1          10
'XPAR'     2          3.0        4.0          20
'YPAR'     3          4.0        5.0          20
'PAR'      4           70         50         100 
'PAR'      5          3.0        0.3          30
'PAR'      6          511         0            0
'PAR'      7           2          0.2          20
'XRMK'     'INDEX1  ( <Eb)'
'YRMK'     'INDEX2  ( >Eb)'
'TITLE'    'CONTOUR MAP INDEX 1/2'
'CHI2'     5.0    / MAXIMUM ALLOWED RED_CHI2 VALUE.
'LEVEL'    99.0     95.0      90.0     68.0    3.0   /MAX:(10)
'OUT'      './HXS_RESULTS.DATA'
'FILE'     './HXS_SPECTRUM.DATA'


 Here we give the explanation. 

'REM'      Comment statement
'TRACE'    Switching of program trace,   1:  OFF, 0: ON
             (DEFAULT: 1:OFF)
'SER'      Set systematic error. We recommend systematic error of 
           0.01 to 0.02.
'CMD'      Set convergence condition code  
           We recommend the following values:
           300 :   max step of simplex
           3.0 :   minimum chi-square value for judgement of 
                   convergence           
           0.001 : relative standard deviation of chi-square value 
                   in simplex for judgement of convergence 
'BIN'      Set the channel numbers used for spectral fitting by 32
           characters. We do not use both the 0 and 31th channel data
           for spectral fitting, because their channel widths are not 
           extractly determined.
           For example,
           01111111111111111111111111223330
           Both the 0 and 31th channels are not use for spectral fitting. 
           Each of 1 - 25th channel data is used for spectral fitting.
           sum of 26 and 27th channel data is used for spectral fitting. 
           sum of 28, 29  and 30th channel data is used for spectral fitting. 
'MDL'      Set spectral fitting models.  
           Available models and their name and parameters
           POWL(i,j) : power law (E**-A)
             i: flux at 50 keV; j: spectral index
           BRPW(i,j,k,l) : broken power law
             i: flux at 50 keV; j: index1 (<Eb); k: index2 (>Eb);
             l: Eb (break energy:keV) 
           BRP2(i,j,k,l,m,n) : broken power law (two breaks)
             i: flux at 50 keV; j: index1 (<Eb1); 
             k: index2 (Eb2>E>Eb1); l: index3 (>Eb2);
             m: Eb1 (keV); n: Eb2 (keV)
           THB1(i,j) : thermal (Crannel: ApJ 1978)
             i: EM(10^45/cm^3); j: kT (keV)
           GAUS(i,j,k) : Gaussian
             i: intensity; j: center energy; 
             k: width (keV) in 1 sigma
             For example, POWL(1,2) + THB1(3,4)
               combination of power law and thermal
'PAR'      Set the initial values of free parameters for fitting program.
           Three initial values of center, lower and upper are set
           in that order.
'OUT'      Input the output filename of fitting result
           './HXS_RESULTS.DATA' (default)
'FILE'     Input the filename of observed data for spectral fitting 
           './HXS_SPECTRUM.DATA' (default)
           This statement always should be bottom in modelfile. 

The following options are for chi-square map.

Warning: In the case of no free parameter, the fitting routine can not 
         work. "not converged" appears in the output file. 

'XPAR'     Set values of X-axis parameters for chi-square map.
           Three values are set for spectral parameters
           in the following order:
             lower       upper      number of division
           Calculation of chi-squared value is carried out between 
           lower and upper values.
'YPAR'     Set values of Y-axis parameters for chi-square map.
           Three values are set for spectral parameters
           in the following order:
             lower       upper      number of division
           Calculation of chi-squared value is carried out between 
           lower and upper values.
'XRMK'     Input the title of X-axis for chi-square map.
'YRMK'     Input the title of Y-axis for chi-square map.
'TITLE'    Input the title for chi-square map.
'CHI2'     Set the reduced chi-squared value at which the simplex
           method is switched to the nonlinear least square method.
'LEVEL'    Set confidence level (%) for chi-square map. 
           The maximum setnumber is 10.