{
float tim[100];
float td[100];
float th[100];
float tmin[100];
float ar[100];
float nit[100];
float co2[100];
ifstream gas_test("FDC_gas_test.dat");
int n=0;
while (gas_test.good()){
gas_test>> td[n]>>th[n]>>tmin[n]>>ar[n]>>co2[n]>>nit[n];
cout<<" ar co2 nit = "<<ar[n]<<" "<<co2[n]<<" "<<nit[n]<<endl;
tim[n]=(tmin[n]+th[n]*60.+td[n]*60.*24.)/60.;
co2[n]/=100.;
ar[n]/=100.;
nit[n]/=100.;
n++;
}
n--;
cout<<" npoints="<<n<<endl;
TGraph *gar = new TGraph(n,tim,ar);
TF1 *eeexp = new TF1("eeexp","exp(-1.38429-0.0355847*x)",0.,920.);
//TH1F *axis = new TH1F("axis","",1,0.,1600.);
TH1F *axis = new TH1F("axis","",1,400.,1200.);
//TH1F *axis = new TH1F("axis","",1,1520.,1580.);
//TH1F *axis = new TH1F("axis","",1,1550.,1575.);
axis->SetMaximum(1.);
axis->SetMinimum(5.e-5);
axis->Draw();
gStyle->SetOptStat(0);
TAxis *xx = (TAxis*)axis->GetXaxis();
TAxis *yy = (TAxis*)axis->GetYaxis();
xx->SetTitle("time, hours");
yy->SetTitle("Ar,N fraction");
axis->Draw();
//eeexp->Draw("same");
gar->SetMarkerStyle(2);
gar->SetMarkerColor(1);
gar->Draw("pl");

TGraph *gnit = new TGraph(n,tim,nit);
gnit->SetMarkerStyle(2);
gnit->SetMarkerColor(2);
gnit->Draw("pl");
gar->SetName("far");
gnit->SetName("fnit");


}