The .spe files contain the spectra. Below: a cutout of one of my yorick codes that extracts stellar parameters from the filenames... // For abundances: factab=1.e4; // Teff teff = 0.; i1= (strfind("_t",mm))(2)+1; ndum=sread(strpart(mm,i1:i1+3),teff); // logg logg = 0.; i1= (strfind("_g",mm))(2)+2; ndum=sread(strpart(mm,i1:i1+3),logg); logg/=100.; // now standard logg(cgs gravity) // mass mass = 0.; i1= (strfind("_m",mm))(2)+1; ndum=sread(strpart(mm,i1:i1+3),mass); mass/=100.; // now in Msun // Oxygen abundance thisO = 0.; i1= (strfind("_o",mm))(2)+1; ndum=sread(strpart(mm,i1:i1+5),thisO); // Oxygen (12-log)*1e4 thisO /=factab ; // now standard 12+log(O/H) // Carbon abundance thisC = 0.; i1= (strfind("_c",mm))(2)+1; ndum=sread(strpart(mm,i1:i1+5),thisC); // Carbon (12-log)*1e4 thisC /=factab ; // now standard 12+log(C/H) // Iron abundance thisFe = 0.; i1= (strfind("_fe",mm))(2)+1; n=sread(strpart(mm,i1:i1+5),thisFe); thisFe /= factab ; // now standard 12+log(Fe/H) // Microturbulent velocity xi = 0.; i1= (strfind("_xi",mm))(2)+1; n=sread(strpart(mm,i1:i1+1),xi); xi /= 10.; // now in km/s // Derived quantities // C/O ratio // Actual expression : (C/H)/(O/H) = 10.^(thisC-12.) / 10.^(thisO-12.) CovO=10.^(thisC-thisO); // simplified exact expression // C excess : 12+log( (C-O)/H ) = 12+log( (C/H)-(O/H) ) // = log( 10.^thisC - 10.^thisO ) after calculus if (thisC > thisO ) { CmO= log10( 10.^(thisC) - 10.^(thisO) ); } else { CmO=-1.e-10; // default value when C