Потоки и работа одновременно с несколькими файлами

double re, e, Latrad, Longrad, Rmer, Rnorm, Reqviv, x, y, z, lat, lon, acpK, acpT, range, limUp, limDown, limLeft, limRight, w1, w2, ticks, UTC,
    devX, devY;    // Ox,Oy,
    string time;
    
    ifstream  aip{"lat.txt"};
    ifstream  aip1{"lon.txt"};
    ifstream  aip2{"time.txt"};
    ifstream ai4{"acpK.txt"};
    ifstream ai6{"acpT.txt"};
    ifstream ai8{"ticks.txt"};
    ifstream ai10{"UTC.txt"};
ofstream f3;
        f3.open("ACP.xls");
while ( ai4 >> acpK && ai6 >> acpT && ai8 >> ticks && ai10 >> UTC)
    {
                        range = 800;
                        limUp = 2.5;
                        limDown = -2.5;
                        limRight = 45;
                        limLeft = -45;
                         w1 = (acpK/range)*((abs(limLeft) + limRight)/(limUp + abs(limDown)));
                         w2 = (acpT/range)*((abs(limLeft) + limRight)/(limUp + abs(limDown)));
                         
    ////////////////////////////////////////////////////////////////////////////////////
    devX = ticks * cos(w1);
    //Ox = x + devX;                   
                         
    f3.precision(12);  
    f3 << UTC << "\t" << w1 << "\t"<< w2 << "\t"<< ticks <<"\t" << "\t" << "\t" << "\t" << "\t" << "\t" << "\t" << devX << endl;
    }
    
       while ( aip >> lat && aip1 >> lon && aip2 >> time )
    {
    re = 6378137;                  // экваториальный радиус(большая полуось)
    e = 0.081819190842622;         // эксцентриситет
    Latrad = lat * (M_PI / 180);       //перевод широты (x) в радианы
    Longrad = lon * (M_PI / 180);     //перевод долготы (y) в радианы
    Rmer = (re*(1 - pow(e, 2))) / pow((1 - pow(e, 2)*pow(sin(Latrad), 2)), 3 / 2);    // радиус меридиана
    Rnorm = (re) / pow((1 - pow(e, 2)*pow(sin(Latrad), 2)), 1 / 2);                // нормальный радиус (радиус кривизны вертикала)
    Reqviv = sqrt(Rmer*Rnorm);     //локальный радиус Земли(геоида)
    x = (Rnorm + Reqviv)*cos(Latrad)*cos(Longrad);     //координата х
    y = (Rnorm + Reqviv)*cos(Latrad)*sin(Longrad);     //координата у
    z = ((1 - pow(e, 2))*Rnorm + Reqviv)*sin(Latrad);       //координата z  
    
    f3 << "\t" << "\t" << "\t"<< "\t" <<"\t"<< time << "\t" << x << "\t" << y << "\t" << z << "\t" << "\t" << "\t" << endl;
 
    }
    ////////////////////////////////////////////////////////////////////////////////////
    aip.close();
    aip1.close();
    aip2.close();
    ai4.close();
    ai6.close();
    ai8.close();
    ai10.close();
    f3.close();