initlayer.cpp

#include "initlayer.h"

#include <vector>
#include <string>
#include <cstdint>
#include <fstream>
#include <iostream>

using std::ofstream;
using std::ios_base;
using std::vector;
using std::string;
using std::cout;
using std::endl;

extern uint32_t beginPoint;
extern uint32_t endPoint;

extern vector <double> initLayerTV;
extern vector <double> initLayerTF;
extern vector <double> initLayerCV;
extern vector <double> initLayerCF;
extern vector <double> initLayerTB_ACU;
extern vector <double> initLayerTB_EVAP;
extern vector <double> initLayerTFG;

//---------------Borders--------------//
extern double initLayerTV_0, initLayerTV_1;
extern double initLayerTF_0, initLayerTF_1;
extern double initLayerCV_0, initLayerCV_1;
extern double initLayerCF_0, initLayerCF_1;
extern double initLayerTB_ACU_0, initLayerTB_ACU_1;
extern double initLayerTB_EVAP_0, initLayerTB_EVAP_1;
extern double initLayerTFG_0, initLayerTFG_1;

extern uint32_t spaceParametrBP;
extern uint32_t spaceParametrTP;
extern uint32_t spaceParametrACU;
extern uint32_t spaceParametrEVAP;

extern double dt;
extern unsigned long long selectN;


void initialLayerTV(vector <vector <double> > &TV, unsigned int it)
{
    size_t i = 0, j = 0, k = 1;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        TV.at(i).at(beginPoint) = initLayerTV_0;
        TV.at(i).at(it-1) = initLayerTV_1;
    }

    // Initial values
    for(j = beginPoint + 1; j <= (it-2); ++j)
    {
        TV.at(0).at(j) = initLayerTV.at(j-k);
    }
}

void initialLayerTF(vector <vector <double> > &TF, unsigned int it)
{
    size_t i = 0, j = 0;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        TF.at(i).at(beginPoint) = initLayerTF_0;
        TF.at(i).at(it-1) = initLayerTF_1;
    }

    // Initial values
    for(j = beginPoint + 1; j <= (it-2); ++j)
    {
        TF.at(0).at(j) = initLayerTF.at(j-1);
    }
}

void initialLayerCV(vector <vector <double> > &CV, unsigned int it)
{
    size_t i = 0, j = 0;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        CV.at(i).at(beginPoint) = initLayerCV_0;
        CV.at(i).at(it-1) = initLayerCV_1;
    }

    // Initial values
    for(j = beginPoint + 1; j <= (it-2); ++j)
    {
        CV.at(0).at(j) = initLayerCV.at(j-1);
    }
}

void initialLayerCF(vector <vector <double> > &CF, unsigned int it)
{
    size_t i = 0, j = 0;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        CF.at(i).at(beginPoint) = initLayerCF_0;
        CF.at(i).at(it-1) = initLayerCF_1;
    }

    // Initial values
    for(j = beginPoint + 1; j <= (it-2); ++j)
    {
        CF.at(0).at(j) = initLayerCF.at(j-1);
    }
}

void initialLayerTB_ACU(vector <vector <double> > &TB_ACU)
{
    size_t i = 0, j = 0;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        TB_ACU.at(i).at(0) = initLayerTB_ACU_0;
        TB_ACU.at(i).at(spaceParametrACU-1) = initLayerTB_ACU_1;
    }

    // Initial values
    for(j = 1; j <= (spaceParametrACU-2); ++j)
    {
        TB_ACU.at(0).at(j) = initLayerTB_ACU.at(j-1);
    }
}

void initialLayerTB_EVAP(vector <vector <double> > &TB_EVAP)
{
    size_t i = 0, j = 0;

    // Borders
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        TB_EVAP.at(i).at(0) = initLayerTB_EVAP_0;
//        TB_EVAP.at(i).at(spaceParametrEVAP-1) = initLayerTB_EVAP_1;
    }

//    // Initial values
//    for(j = 1; j <= (spaceParametrEVAP-2); ++j)
//    {
//        TB_EVAP.at(0).at(j) = initLayerTB_EVAP.at(j-1);
//    }
}

void initialLayerTFG(vector <vector <double> > &TFG)
{
    size_t i = 0, j = 0;

    // Bordersuint32_t
    for(i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        TFG.at(i).at(0) = initLayerTFG_0;
        TFG.at(i).at(spaceParametrEVAP-1) = initLayerTFG_1;
    }

    // Initial values
    for(j = 1; j <= (spaceParametrEVAP-2); ++j)
    {
        TFG.at(0).at(j) = initLayerTFG.at(j-1);
    }
}

//?????----Why it here?????????
void toFileMM(vector <vector <double> > MMM, string nameModel)
{
    uint32_t sizeModel = 0;

    // In future - RegExp
    if(nameModel == "TP_ACU_TV")
    {
        sizeModel = spaceParametrTP + spaceParametrACU - 2;
    }
    else
    {
        sizeModel = spaceParametrTP - 1;
    }

    if( (nameModel == "TV_ACU") || (nameModel == "TB_ACU") )
    {
        sizeModel = spaceParametrACU - 1;
    }

    nameModel = "MM_" + nameModel + ".dat";

    ofstream foutMM(nameModel, ios_base::out | ios_base::trunc);

    for(size_t i = 0; i < static_cast <size_t>(selectN / dt); ++i)
    {
        for(uint32_t j = 1; j < (sizeModel); j++)
        {
            foutMM << MMM.at(i).at(j) << " ";
        }
        foutMM << endl;
    }
    foutMM.close();
}