function [line0] = allco2lines(idgas,database,f1,f2,iso);

%%function [line0] = lineparameters(idgas,database,f1,f2,iso)
% input : idgas : gasID
%         database : HITRAN96, HITRAN2K etc??
%                    default /asl/data/hitran/h2k.oldiso/
%         f1,f2 : start and stop frequencies within which to search HITRAN
%         iso   : gas isotope

clf
path(path,'/home/sergio/SPECTRA');

%idgas = input('Enter gasID : ');
%f1    = input('Enter start freq for hitread : ');
%f2    = input('Enter stop freq for hitread : ');

hitran = database;
%fname  = ['/asl/data/hitran/h2k.oldiso/g' num2str(idgas) '.dat'];
fname  = [hitran 'g' num2str(idgas) '.dat'];
fprintf(1,'reading in HITRAN data %s \n',fname);
lineX = hitread(f1,f2,0,idgas,fname);

%iso   = input('Enter gas isotope : ');
%iL    = input('Enter gas lower quantum number  : ');
%iU    = input('Enter gas upper quantum number  : ');

ii=find(lineX.iso==iso); 

line0.stren   = lineX.stren(ii); 
line0.elower  = lineX.els(ii); 
line0.freq    = lineX.wnum(ii); 
line0.j_lower = lineX.ilsgq(ii); 
line0.j_upper = lineX.iusgq(ii); 
line0.p_shift = lineX.tsp(ii); 
line0.w       = lineX.abroad(ii); 
line0.w_s     = lineX.sbroad(ii); 
line0.w_temp  = lineX.abcoef(ii); 
