function [OK,params]=read_params (fname);
% MATLAB HELP:
% 
% function [OK,params]=read_params (fname);
%
% reads paramaters from the file
%
% inputs:
%  fname, the file name
% outputs:
%  OK, this is one if things went well
%  params, a structure containing the paramaters for the kernel
%  calculation that you want to do
%
% see codes.ps or codes.pdf for details
%
% aaron birch, may 21, 2002

fid=fopen(fname,'rt');
if (fid==-1)
     OK=-1;
     return;
end;
fgetl(fid);
params.ONED=fscanf(fid,'one-dimensional kernel = %d \n');
fgetl(fid);
params.model_name=fscanf(fid,'MODEL = %s \n');
fgetl(fid);
params.NMODES=fscanf(fid,'NMODES = %d \n');
fgetl(fid);
params.SKIP=fscanf(fid,'SKIP = %d \n');
fgetl(fid);
params.FILTER_TYPE=fscanf(fid,'FILTER_TYPE = %d \n');
fgetl(fid);
params.FILTER_PARAM1=fscanf(fid,'FILTER_PARAM1 = %f \n');
fgetl(fid);
params.FILTER_PARAM2=fscanf(fid,'FILTER_PARAM2 = %f \n');
fgetl(fid);
params.FILTER_PARAM3=fscanf(fid,'FILTER_PARAM3 = %f \n');
fgetl(fid);
params.FILTER_PARAM4=fscanf(fid,'FILTER_PARAM4 = %f \n');
fgetl(fid);
params.EXCLUDE_F=fscanf(fid,'EXCLUDE_F = %d \n');
fgetl(fid);
params.T_MIN=fscanf(fid,'T_MIN = %f \n');
fgetl(fid);
params.T_MAX=fscanf(fid,'T_MAX = %f \n');
fgetl(fid);
params.T_STEPS=fscanf(fid,'T_STEPS = %d \n');
fgetl(fid);
params.delta=fscanf(fid,'delta = %f \n');
fgetl(fid);
params.CUT=fscanf(fid,'CUT = %d \n');
fgetl(fid);
params.NPHI=fscanf(fid,'NPHI = %d \n');
fgetl(fid);
params.NR=fscanf(fid,'NR = %d \n');
fgetl(fid);
params.MINR=fscanf(fid,'MINR = %f \n');
fgetl(fid);
params.VERBOSE=fscanf(fid,'VERBOSE = %f \n');
fgetl(fid);
params.PLOTS=fscanf(fid,'PLOTS = %f \n');
fgetl(fid);
params.CALC_F_ONLY=fscanf(fid,'CALC_F_ONLY = %f \n');



fclose(fid);
OK=1;
return