SystemModule.cpp

Go to the documentation of this file.

// SystemModule.cpp: implementation of the SystemModule class.
//
 
#include "SystemModule.h"
#include "../Branch/BranchRecord.h"
#include "../Global/SimulationContext.h"
using namespace std;
 
 
// Construction/Destruction
 
SystemModule::SystemModule()
{
 
}
 
SystemModule::~SystemModule()
{
    this->clearSystemModiule();
}
void SystemModule::clearSystemModiule(){
    control_variable_array.clear();
    hydraulic_junction_array.clear();
    hydraulic_volume_array.clear();
    interactive_variable_array.clear();
    logical_trip_array.clear();
    variable_trip_array.clear();
    heat_structure_array.clear();
    system_map.clear();
}
 
//**********************************************************************************************
//readInfoFromFile: read the output of GUI from harddisk to build the components in the TH system
//return 0 if reading correct
//return -1 if file can not be read;  error message is stored in "Error Message.txt" in the directory
//of executable file
//**********************************************************************************************
typedef map<string, THElement *> sm_type;
int SystemModule::readInfoFromFile(){
    int n_comp, i;
    string comp_type, comp_name, comp_unit;
    string file_name;
    //
    file_name = g_sim.SScheduler.getSystemDirectory();
    file_name += "RELAP5_channels.txt";
    ifstream fin(file_name.c_str(), ios::in);
    if(!fin.is_open()){
        // A plant with no RELAP5 channel map cannot run meaningfully: reading on
        // from the closed stream would yield an empty plant model, so this is fatal.
        cout << "FATAL: File: " << file_name <<  " could not be found" << endl;
        ofstream fout("err messages.txt", ios::app);
        fout << r5pars_.timesno << "  FATAL: File: " << file_name <<  " could not be found\n";
        fout.close();
        std::exit(1);
    }
    //read hydraulic volumes
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<HydraulicVolume>();
        p->readInfoFromFile(fin);
        hydraulic_volume_array.push_back(std::move(p));
    }
    //read hydraulic junctions
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<HydraulicJunction>();
        p->readInfoFromFile(fin);
        hydraulic_junction_array.push_back(std::move(p));
    }
    //read variable trips
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<VariableTrip>();
        p->readInfoFromFile(fin);
        variable_trip_array.push_back(std::move(p));
    }
    //read logical trips
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<LogicalTrip>();
        p->readInfoFromFile(fin);
        logical_trip_array.push_back(std::move(p));
    }
    //read Interactive Variables (Channels)
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<InteractiveVariable>();
        p->readInfoFromFile(fin);
        p->setSequenceInTheArray(i);
        interactive_variable_array.push_back(std::move(p));
    }
    //read Heat Structure
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<HeatStructure>();
        p->readInfoFromFile(fin);
        heat_structure_array.push_back(std::move(p));
    }
    //read Control Variables
    fin >> comp_type >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto p = make_unique<ControlVariable>();
        p->readInfoFromFile(fin);
        control_variable_array.push_back(std::move(p));
    }
    fin.close();
    //build up system map====================================
    for(i = 0; i < hydraulic_volume_array.size(); i++)
        system_map.insert(sm_type::value_type(hydraulic_volume_array[i]->getName(), hydraulic_volume_array[i].get()));
    for(i = 0; i < hydraulic_junction_array.size(); i++)
        system_map.insert(sm_type::value_type(hydraulic_junction_array[i]->getName(), hydraulic_junction_array[i].get()));
    for(i = 0; i < variable_trip_array.size(); i++)
        system_map.insert(sm_type::value_type(variable_trip_array[i]->getName(), variable_trip_array[i].get()));
    for(i = 0; i < logical_trip_array.size(); i++)
        system_map.insert(sm_type::value_type(logical_trip_array[i]->getName(), logical_trip_array[i].get()));
    for(i = 0; i < interactive_variable_array.size(); i++)
        system_map.insert(sm_type::value_type(interactive_variable_array[i]->getName(), interactive_variable_array[i].get()));
    for(i = 0; i < heat_structure_array.size(); i++)
        system_map.insert(sm_type::value_type(heat_structure_array[i]->getName(), heat_structure_array[i].get()));
    for(i = 0; i < control_variable_array.size(); i++)
        system_map.insert(sm_type::value_type(control_variable_array[i]->getName(), control_variable_array[i].get()));
    return 0;
}
 
//**********************************************************************************************
//setConnectionToRELAP5: specify the number of TH elements (i.e., hraraulic volumes, ydraulic junctions,
//trips, interactive variables, control variables, and heat structures) and teh corresponding card numbers 
//which the user specified to interact with rest of ADS routines;  
//This reoutine must be called after the routine readInfoFromFile()
//return 0 if successful completion of the task
//return -1 if errors occur
//currently set the number of each type of information is 100; number of heat structure is 10
// edit systemParameters.for and fmodule.h for change 
//**********************************************************************************************
int SystemModule::setConnectionToRELAP5(){
    int i, max_channel, max_ht, max_mesh, n_return, n_mesh;
    n_return = 0;
    max_channel = 100;
    max_ht = 10;
    max_mesh = 100;
    //
    r5pars_.n_hydrov = hydraulic_volume_array.size();
    r5pars_.n_hydroj = hydraulic_junction_array.size();
    r5pars_.n_vartrip = variable_trip_array.size();
    r5pars_.n_logtrip = logical_trip_array.size();
    r5pars_.n_interact = interactive_variable_array.size();
    r5pars_.n_ctrl_var = control_variable_array.size();
    r5pars_.n_heatstr = heat_structure_array.size();
    if (r5pars_.n_hydrov > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_hydrov value " << r5pars_.n_hydrov 
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    if (r5pars_.n_hydroj > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_hydroj value " <<  r5pars_.n_hydroj
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    if (r5pars_.n_vartrip > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_vartrip value " <<  r5pars_.n_vartrip
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    if (r5pars_.n_logtrip > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_logtrip value " << r5pars_.n_logtrip 
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    if (r5pars_.n_interact > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_interact value " << r5pars_.n_interact 
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    if (r5pars_.n_ctrl_var > max_channel){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_ctrl_var value " << r5pars_.n_ctrl_var 
            << " exceeds the maximum allowance value " << max_channel << '\n';
        n_return = -1;
        fout.close();
    }
    // heat_str_cn/heat_str_obj/n_mesh are [10] in r5pars_, so heat structures
    // are checked against max_ht (10), not max_channel: more than 10 would
    // overwrite adjacent struct fields shared with Fortran.
    if (r5pars_.n_heatstr > max_ht){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: r5pars_.n_heatstr value " << r5pars_.n_heatstr
            << " exceeds the maximum allowance value " << max_ht << '\n';
        n_return = -1;
        fout.close();
    }
    if(n_return == -1){
        return n_return;
    }
 
    n_mesh = 0;
    for(i = 0; i < r5pars_.n_heatstr; i++){
        n_mesh += heat_structure_array[i]->getNumberOfMeshPoints();
    }
    if (n_mesh > max_mesh){
        ofstream fout("err messages.txt", ios_base::out | ios_base::app);
        fout << r5pars_.timesno << "  SystemModule::setConnectionToRELAP5: total mesh points" << n_mesh
            << " exceeds the maximum allowance value " << max_mesh << '\n';
        n_return = -1;
        fout.close();
    }
    //pass card numbers to systemParameter.for to retrieve RELAP5 data
    if (n_return != -1){
        for(i = 0; i < r5pars_.n_hydrov; i++){
            r5pars_.hydro_v_cn[i] = hydraulic_volume_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_hydroj; i++){
            r5pars_.hydro_j_cn[i] = hydraulic_junction_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_vartrip; i++){
            r5pars_.var_trip_cn[i] = variable_trip_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_logtrip; i++){
            r5pars_.log_trip_cn[i] = logical_trip_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_interact; i++){
            r5pars_.interact_cn[i] = interactive_variable_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_ctrl_var; i++){
            r5pars_.ctrl_var_cn[i] = control_variable_array[i]->getCardNumber();}
        for(i = 0; i < r5pars_.n_heatstr; i++){
            r5pars_.heat_str_cn[i] = heat_structure_array[i]->getCardNumber();
            r5pars_.n_mesh[i] = heat_structure_array[i]->getNumberOfMeshPoints();
        }
    }
    return n_return;
};
 
void SystemModule::updateParametersValues(){
    int i, j;
    double dbg_x1;
    decltype(hydraulic_volume_array)::iterator hv_itr;
    decltype(hydraulic_junction_array)::iterator hj_itr;
    decltype(interactive_variable_array)::iterator iv_itr;
    decltype(logical_trip_array)::iterator lt_itr;
    decltype(variable_trip_array)::iterator vt_itr;
    decltype(control_variable_array)::iterator ctrl_itr;
    decltype(heat_structure_array)::iterator ht_itr;
    //Update parameters' values of Hydraulic Volumes
    i = 0;
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        (*hv_itr)->setLiquidTemperature(r5pars_.liq_t[i]);
        (*hv_itr)->setVaporTemperature(r5pars_.vap_t[i]);
        (*hv_itr)->setSaturationTemperature(r5pars_.sat_t[i]);
        (*hv_itr)->setPressure(r5pars_.sys_pressure[i]);
        (*hv_itr)->setBoronDensity(r5pars_.boron_d[i]);
        (*hv_itr)->setEquilibriumQuality(r5pars_.quality[i]);
        (*hv_itr)->setFlowRegimeNumber(r5pars_.flow_reg[i]);
        i++;
    }
    //Update parameters' values of Hydraulic Junctions
    i = 0;
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        (*hj_itr)->setLiquidVelocity(r5pars_.liq_vel[i]);
        (*hj_itr)->setMassFlowRate(r5pars_.mass_fr[i]);
        (*hj_itr)->setVaporVelocity(r5pars_.vap_vel[i]);
        i++;
    }
    //Update parameters' values of Variable Trips
    i = 0;
    for(vt_itr = variable_trip_array.begin(); vt_itr != variable_trip_array.end(); vt_itr++){
        (*vt_itr)->setTripTime(r5pars_.vartrip_t[i]);
        i++;
    }
    //Update parameters' values of Logical Trips
    i = 0;
    for(lt_itr = logical_trip_array.begin(); lt_itr != logical_trip_array.end(); lt_itr++){
        (*lt_itr)->setTripTime(r5pars_.logtrip_t[i]);
        i++;
    }
    //Update parameters' values of Interactive Variables
    i = 0;
    for(iv_itr = interactive_variable_array.begin(); iv_itr != interactive_variable_array.end(); iv_itr++){
        (*iv_itr)->setDisplayValue(r5pars_.inter_val[i]);
        i++;
    }
    //Update parameters' values of Control Variables
    i = 0;
    for(ctrl_itr = control_variable_array.begin(); ctrl_itr != control_variable_array.end(); ctrl_itr++){
        (*ctrl_itr)->setControlValue(r5pars_.ctrl_var_val[i]);
        i++;
    }
    //Update parameters' values of Heat Structure
    j = 0; 
    for(ht_itr = heat_structure_array.begin(); ht_itr != heat_structure_array.end(); ht_itr++){
        for(i = 0; i < (*ht_itr)->getNumberOfMeshPoints(); i++){
            (*ht_itr)->setMeshPointTemperature(i,r5pars_.mesh_t[j]);
            dbg_x1 = r5pars_.mesh_t[j];
            j++;
        }
    }
}
 
void SystemModule::outputParametersLabelsToFile(ofstream &fout){
    int i;
    string name, str1, str2, str3, str4, str5, str6, str7;
    decltype(hydraulic_volume_array)::iterator hv_itr;
    decltype(hydraulic_junction_array)::iterator hj_itr;
    decltype(variable_trip_array)::iterator vt_itr;
    decltype(logical_trip_array)::iterator lt_itr;
    decltype(interactive_variable_array)::iterator iv_itr;
    decltype(control_variable_array)::iterator ctrl_itr;
    decltype(heat_structure_array)::iterator ht_itr;
    //Output Name
    //Hydraulic Volumes
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        name = (*hv_itr)->getName();
        for(i = 0; i < 7; i++){
            fout << name << '\t';
        }
    }
    //Hydraulic Junctions
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        name = (*hj_itr)->getName();
        for(i = 0; i < 3; i++){
            fout << name << '\t';
        }
    }
    //Variable Trips
    for(vt_itr = variable_trip_array.begin(); vt_itr != variable_trip_array.end(); vt_itr++){
        fout << (*vt_itr)->getName() << '\t';
    }
    //Logical Trips
    for(lt_itr = logical_trip_array.begin(); lt_itr != logical_trip_array.end(); lt_itr++){
        fout << (*lt_itr)->getName() << '\t';
    }
    //Interactive Variables
    for(iv_itr = interactive_variable_array.begin(); iv_itr != interactive_variable_array.end(); iv_itr++){
        fout << (*iv_itr)->getName() << '\t';
    }
    //Control Variables
    for(ctrl_itr = control_variable_array.begin(); ctrl_itr != control_variable_array.end(); ctrl_itr++){
        fout << (*ctrl_itr)->getName() << '\t';
    }
    //Heat Structure
    for(ht_itr = heat_structure_array.begin(); ht_itr != heat_structure_array.end(); ht_itr++){
        name = (*ht_itr)->getName();
        for(i = 0; i < (*ht_itr)->getNumberOfMeshPoints(); i++){
            fout << name << '\t';
        }
    }
    fout << '\n';
    /*
    c
c       liq_t(i) = Volume liquid temperature (K, oF)->
c       vap_t(i) = Volume vapor temperature (K, oF)->
c       sat_t(i) = Volume saturation temperature based on the partial pressure of steam (K, oF)->
c       sys_pressure(i): Volume pressure (Pa, lbf/in^2)
c       boron_d(i) = Spatial boron density, rb (kg/m3, lb/ft3)-> This is the volume liquid fraction (af) times the
c                      liquid density (rf) times the boron concentration (Cb)-> Boron concentration is used for
c                      hydrodynamic input, and boron density is used for minor edits and plots.
c       quality(i) = Volume equilibrium quality. This quality uses phasic enthalpies and mixture quality, with
c                      the mixture enthalpy calculated using the flow quality.
c       flow_reg(i) = Flow regime number;
 
    */
    //output parameters' Names and units;
    //Hydraulic Volumes
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        (*hv_itr)->setUnitOfLiquidTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfVaporTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfSaturationTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfPressure(r5pars_.si_unit);
        (*hv_itr)->setUnitOfBoronDensity(r5pars_.si_unit);
        str1 = "Liquid_T("; str1 += (*hv_itr)->getUnitOfLiquidTemperature(); str1 += ")";
        str2 = "Vapor_T("; str2 += (*hv_itr)->getUnitOfVaporTemperature(); str2 += ")";
        str3 = "Sat_T("; str3 += (*hv_itr)->getUnitOfSaturationTemperature(); str3 += ")";
        str4 = "Pressure("; str4 += (*hv_itr)->getUnitOfPressure(); str4 += ")";
        str5 = "Boron_density("; str5 += (*hv_itr)->getUnitOfBoronDensity(); str5 += ")";
        str6 = "Equilibrium";
        str7 = "Flow_regime";
        fout << str1 << '\t' << str2 << '\t' << str3 << '\t' << str4 << '\t' << str5 << '\t' << str6 << '\t' << str7 << '\t';
    }
    //Hydraulic Junctions
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        (*hj_itr)->setUnitOfLiquidVelocity(r5pars_.si_unit);
        (*hj_itr)->setUnitOfMassFlowRate(r5pars_.si_unit);
        (*hj_itr)->setUnitOfVaporVelocity(r5pars_.si_unit);
        str1 = "Liquid_vel("; str1 += (*hj_itr)->getUnitOfLiquidVelocity(); str1 += ")";
        str2 = "Mass_flow("; str2 += (*hj_itr)->getUnitOfMassFlowRate(); str2 += ")";
        str3 = "Vapor_vel("; str3 += (*hj_itr)->getUnitOfVaporVelocity(); str3 += ")";
        fout << str1 << '\t' << str2 << '\t' << str3 << '\t';
    }
    //Variable Trips
    for(vt_itr = variable_trip_array.begin(); vt_itr != variable_trip_array.end(); vt_itr++){
        fout << "Trip_time(s)\t";
    }
    //Logical Trips
    for(lt_itr = logical_trip_array.begin(); lt_itr != logical_trip_array.end(); lt_itr++){
        fout << "Trip_time(s)\t";
    }
    //Interactive Variables
    for(iv_itr = interactive_variable_array.begin(); iv_itr != interactive_variable_array.end(); iv_itr++){
        fout << "--\t";
    }
    //Control Variables
    for(ctrl_itr = control_variable_array.begin(); ctrl_itr != control_variable_array.end(); ctrl_itr++){
        fout << (*ctrl_itr)->getUnitOfControlValue() << '\t';
        i++;
    }
    //Heat Structure
    for(ht_itr = heat_structure_array.begin(); ht_itr != heat_structure_array.end(); ht_itr++){
        (*ht_itr)->setUnitOfMeshPointTemperature(r5pars_.si_unit);
        for(i = 0; i < (*ht_itr)->getNumberOfMeshPoints(); i++){
            str1 = "Mesh_T_";  
            str1 += to_string(i); 
            str1 += "("; 
            str1 += (*ht_itr)->getUnitOfMeshPointTemperature();
            str1 += ")";
            fout << str1 << '\t';
        }
    }
    fout << '\n';
}
//
void SystemModule::setUnits(){
    decltype(hydraulic_volume_array)::iterator hv_itr;
    decltype(hydraulic_junction_array)::iterator hj_itr;
    decltype(heat_structure_array)::iterator ht_itr;
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        (*hv_itr)->setUnitOfLiquidTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfVaporTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfSaturationTemperature(r5pars_.si_unit);
        (*hv_itr)->setUnitOfPressure(r5pars_.si_unit);
        (*hv_itr)->setUnitOfBoronDensity(r5pars_.si_unit);
    }
    //Hydraulic Junctions
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        (*hj_itr)->setUnitOfLiquidVelocity(r5pars_.si_unit);
        (*hj_itr)->setUnitOfMassFlowRate(r5pars_.si_unit);
        (*hj_itr)->setUnitOfVaporVelocity(r5pars_.si_unit);
    }
    //Heat Structure
    for(ht_itr = heat_structure_array.begin(); ht_itr != heat_structure_array.end(); ht_itr++){
        (*ht_itr)->setUnitOfMeshPointTemperature(r5pars_.si_unit);
    }
}
//  
void SystemModule::outputParametersValuesToFile(ofstream &fout){
    int i;
    decltype(hydraulic_volume_array)::iterator hv_itr;
    decltype(hydraulic_junction_array)::iterator hj_itr;
    decltype(variable_trip_array)::iterator vt_itr;
    decltype(logical_trip_array)::iterator lt_itr;
    decltype(interactive_variable_array)::iterator iv_itr;
    decltype(control_variable_array)::iterator ctrl_itr;
    decltype(heat_structure_array)::iterator ht_itr;
    //Hydraulic Volumes
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        fout << (*hv_itr)->getLiquidTemperature() << '\t'
            << (*hv_itr)->getVaporTemperature() << '\t'
            << (*hv_itr)->getSaturationTemperature() << '\t'
            << (*hv_itr)->getPressure() << '\t'
            << (*hv_itr)->getBoronDensity() << '\t'
            << (*hv_itr)->getEquilibriumQuality() << '\t'
            << (*hv_itr)->getFlowRegimeNumber() << '\t';
    }
    //Hydraulic Junctions
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        fout << (*hj_itr)->getLiquidVelocity() << '\t'
            << (*hj_itr)->getMassFlowRate() << '\t'
            << (*hj_itr)->getVaporVelocity() << '\t';
    }
    //Variable Trips
    for(vt_itr = variable_trip_array.begin(); vt_itr != variable_trip_array.end(); vt_itr++){
        fout << (*vt_itr)->getTripTime() << '\t';
    }
    //Logical Trips
    for(lt_itr = logical_trip_array.begin(); lt_itr != logical_trip_array.end(); lt_itr++){
        fout << (*lt_itr)->getTripTime() << '\t';
    }
    //Interactive Variables
    for(iv_itr = interactive_variable_array.begin(); iv_itr != interactive_variable_array.end(); iv_itr++){
        fout << (*iv_itr)->getDisplayValue() << '\t';
    }
    //Control Variables
    for(ctrl_itr = control_variable_array.begin(); ctrl_itr != control_variable_array.end(); ctrl_itr++){
        fout << (*ctrl_itr)->getControlValue() << '\t';
    }
    //Heat Structure
    for(ht_itr = heat_structure_array.begin(); ht_itr != heat_structure_array.end(); ht_itr++){
        for(i = 0; i < (*ht_itr)->getNumberOfMeshPoints(); i++){
            fout << (*ht_itr)->getMeshPointTemperature(i) << '\t';
        }
    }
    fout << '\n';
}
//typedef map<string, THElement *> sm_type;
THElement *SystemModule::getComponentFromSystemMap(string _str){
    THElement *ptr;
    ptr = 0;
    map<string, THElement *>::iterator itr;
    itr = system_map.find(_str);
    if (itr != system_map.end()){
        ptr = itr->second;
    }
    else{
        ofstream fout("err messages.txt", ios::app);
        cout << "SystemModule::getComponentFromSystemMap indicator " << _str << " could not be found in plant." << endl;
        fout << r5pars_.timesno << "  SystemModule::getComponentFromSystemMap indicator " << _str << " could not be found in plant.\n";
        fout.close();
        r5pars_.sim_err = true;
        g_sim.Branch_center.setTempESMessage("SystemModule::getComponentFromSystemMap err.");
        g_sim.Branch_center.setTempTypeOfConsequence(Consequence::AdsSimErr);
    }
    return ptr;
}
//
void SystemModule::initiationForSimulation(){
    this->readInfoFromFile();
    if (this->setConnectionToRELAP5() == -1) {
        cout << "FATAL: SystemModule::setConnectionToRELAP5 failed — channel counts exceed "
             << "the r5pars_ bridge capacity (see err messages.txt)." << endl;
        std::exit(1);
    }
    this->setUnits();
}
void SystemModule::initiationForSimulationClient(){
    if (this->setConnectionToRELAP5() == -1) {
        cout << "FATAL: SystemModule::setConnectionToRELAP5 failed — channel counts exceed "
             << "the r5pars_ bridge capacity (see err messages.txt)." << endl;
        std::exit(1);
    }
    this->setUnits();
}//
void SystemModule::executeActionsToInteractWithSystem(){
    decltype(interactive_variable_array)::iterator itr;
    for(itr = interactive_variable_array.begin(); itr != interactive_variable_array.end(); itr++){
        (*itr)->executeInteractiveAction();
    }
}
void SystemModule::oneThroughProcess(){
    this->updateParametersValues();
    this->executeActionsToInteractWithSystem();
}
void SystemModule::saveSystemModuleToFile(ofstream &fout){
    decltype(control_variable_array)::iterator cv_itr;//control_variable_array;
    decltype(heat_structure_array)::iterator hs_itr;// heat_structure_array;
    decltype(hydraulic_junction_array)::iterator hj_itr;// hydraulic_junction_array;
    decltype(hydraulic_volume_array)::iterator hv_itr;// hydraulic_volume_array;
    decltype(interactive_variable_array)::iterator iv_itr;// interactive_variable_array;
    decltype(logical_trip_array)::iterator lt_itr;// logical_trip_array;
    decltype(variable_trip_array)::iterator vt_itr;// variable_trip_array;
    //
    fout << "SystemModule\n" << hydraulic_volume_array.size() << '\n';
    for(hv_itr = hydraulic_volume_array.begin(); hv_itr != hydraulic_volume_array.end(); hv_itr++){
        (*hv_itr)->saveHydraulicVolumeToFile(fout);
    }
    fout << hydraulic_junction_array.size() << '\n';
    for(hj_itr = hydraulic_junction_array.begin(); hj_itr != hydraulic_junction_array.end(); hj_itr++){
        (*hj_itr)->saveHydraulicJunctionToFile(fout);
    }
    fout << variable_trip_array.size() << '\n';
    for(vt_itr = variable_trip_array.begin(); vt_itr != variable_trip_array.end(); vt_itr++){
        (*vt_itr)->saveVariableTripToFile(fout);;
    }
    fout << logical_trip_array.size() << '\n';
    for(lt_itr = logical_trip_array.begin(); lt_itr != logical_trip_array.end(); lt_itr++){
        (*lt_itr)->saveLogicalTripToFile(fout);
    }
    fout << interactive_variable_array.size() << '\n';
    for(iv_itr = interactive_variable_array.begin(); iv_itr != interactive_variable_array.end(); iv_itr++){
        (*iv_itr)->saveInteractiveVariableToFile(fout);
    }
    fout << heat_structure_array.size() << '\n';
    for(hs_itr = heat_structure_array.begin(); hs_itr != heat_structure_array.end(); hs_itr++){
        (*hs_itr)->saveHeatStructureToFile(fout);
    }
    fout << control_variable_array.size() << '\n';
    for(cv_itr = control_variable_array.begin(); cv_itr != control_variable_array.end(); cv_itr++){
        (*cv_itr)->saveControlVariableToFile(fout);
    }
}
void SystemModule::retrieveSystemModuleFromFile(ifstream &fin){
    int n_comp, i;
    string str;
    THElement *element_ptr;
    string file_name;
    //
 
    //read hydraulic volumes
    fin >> str >> n_comp;
    if(str != "SystemModule"){
        ofstream fout("err messages.txt", ios::app);
        cout << "SystemModule::retrieveSystemModuleFromFile File retrieval error in previous function" << endl;
        fout << r5pars_.timesno << "  SystemModule::retrieveSystemModuleFromFile File retrieval error in previous function\n";
        fout.close();
        r5pars_.sim_err = true;
        g_sim.Branch_center.setTempESMessage("SystemModule::retrieveSystemModuleFromFile Prior error");
        g_sim.Branch_center.setTempTypeOfConsequence(Consequence::NoKB);
        return;
    }
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<HydraulicVolume>();
        ptr->retrieveHydraulicVolumeFromFile(fin);
        hydraulic_volume_array.push_back(std::move(ptr));
    }
    //read hydraulic junctions
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<HydraulicJunction>();
        ptr->retrieveHydraulicJunctionFromFile(fin);
        hydraulic_junction_array.push_back(std::move(ptr));
    }
    //read variable trips
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<VariableTrip>();
        ptr->retrieveVariableTripFromFile(fin);
        variable_trip_array.push_back(std::move(ptr));
    }
    //read logical trips
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<LogicalTrip>();
        ptr->retrieveLogicalTripFromFile(fin);
        logical_trip_array.push_back(std::move(ptr));
    }
    //read Interactive Variables (Channels)
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<InteractiveVariable>();
        ptr->retrieveInteractiveVariableFromFile(fin);
        ptr->setSequenceInTheArray(i);
        interactive_variable_array.push_back(std::move(ptr));
    }
    //read Heat Structure
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<HeatStructure>();
        ptr->retrieveHeatStructureFromFile(fin);
        heat_structure_array.push_back(std::move(ptr));
    }
    //read Control Variables
    fin >> n_comp;
    for(i = 0; i < n_comp; i++){
        auto ptr = make_unique<ControlVariable>();
        ptr->retrieveControlVariableFromFile(fin);
        control_variable_array.push_back(std::move(ptr));
    }
    //build up system map====================================
    //hydraulic volumes
    for(i = 0; i < hydraulic_volume_array.size(); i++){
        system_map.insert(sm_type::value_type(hydraulic_volume_array[i]->getName(), hydraulic_volume_array[i].get()));
    }
    //hydraulic junctions
    for(i = 0; i < hydraulic_junction_array.size(); i++){
        system_map.insert(sm_type::value_type(hydraulic_junction_array[i]->getName(), hydraulic_junction_array[i].get()));
    }
    //variable trips
    for(i = 0; i < variable_trip_array.size(); i++){
        system_map.insert(sm_type::value_type(variable_trip_array[i]->getName(), variable_trip_array[i].get()));
    }
    //logical trips
    for(i = 0; i < logical_trip_array.size(); i++){
        system_map.insert(sm_type::value_type(logical_trip_array[i]->getName(), logical_trip_array[i].get()));
    }
    //Interactive Variables (Channels)
    for(i = 0; i < interactive_variable_array.size(); i++){
        system_map.insert(sm_type::value_type(interactive_variable_array[i]->getName(), interactive_variable_array[i].get()));
    }
    //Heat Structure
    for(i = 0; i < heat_structure_array.size(); i++){
        system_map.insert(sm_type::value_type(heat_structure_array[i]->getName(), heat_structure_array[i].get()));
    }
    //Control Variables
    for(i = 0; i < control_variable_array.size(); i++){
        system_map.insert(sm_type::value_type(control_variable_array[i]->getName(), control_variable_array[i].get()));
    }
}
void SystemModule::copySystemModule(SystemModule &_sm){
    int i;
//
    for(auto& src : _sm.control_variable_array){
        auto ptr = make_unique<ControlVariable>();
        ptr->copyControlVariable(*src);
        control_variable_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.heat_structure_array){
        auto ptr = make_unique<HeatStructure>();
        ptr->copyHeatStructure(*src);
        heat_structure_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.hydraulic_junction_array){
        auto ptr = make_unique<HydraulicJunction>();
        ptr->copyHydraulicJunction(*src);
        hydraulic_junction_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.hydraulic_volume_array){
        auto ptr = make_unique<HydraulicVolume>();
        ptr->copyHydraulicVolume(*src);
        hydraulic_volume_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.interactive_variable_array){
        auto ptr = make_unique<InteractiveVariable>();
        ptr->copyInteractiveVariable(*src);
        interactive_variable_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.logical_trip_array){
        auto ptr = make_unique<LogicalTrip>();
        ptr->copyLogicalTrip(*src);
        logical_trip_array.push_back(std::move(ptr));
    }
    for(auto& src : _sm.variable_trip_array){
        auto ptr = make_unique<VariableTrip>();
        ptr->copyVariableTrip(*src);
        variable_trip_array.push_back(std::move(ptr));
    }
    //build up system map====================================
    //hydraulic volumes
    for(i = 0; i < hydraulic_volume_array.size(); i++){
        system_map.insert(sm_type::value_type(hydraulic_volume_array[i]->getName(), hydraulic_volume_array[i].get()));
    }
    //hydraulic junctions
    for(i = 0; i < hydraulic_junction_array.size(); i++){
        system_map.insert(sm_type::value_type(hydraulic_junction_array[i]->getName(), hydraulic_junction_array[i].get()));
    }
    //variable trips
    for(i = 0; i < variable_trip_array.size(); i++){
        system_map.insert(sm_type::value_type(variable_trip_array[i]->getName(), variable_trip_array[i].get()));
    }
    //logical trips
    for(i = 0; i < logical_trip_array.size(); i++){
        system_map.insert(sm_type::value_type(logical_trip_array[i]->getName(), logical_trip_array[i].get()));
    }
    //Interactive Variables (Channels)
    for(i = 0; i < interactive_variable_array.size(); i++){
        system_map.insert(sm_type::value_type(interactive_variable_array[i]->getName(), interactive_variable_array[i].get()));
    }
    //Heat Structure
    for(i = 0; i < heat_structure_array.size(); i++){
        system_map.insert(sm_type::value_type(heat_structure_array[i]->getName(), heat_structure_array[i].get()));
    }
    //Control Variables
    for(i = 0; i < control_variable_array.size(); i++){
        system_map.insert(sm_type::value_type(control_variable_array[i]->getName(), control_variable_array[i].get()));
    }
};