dv_rho.cc

Go to the documentation of this file.

 
#include "dv_rho.h"
#include "domain.h"
 
 
dv_rho::dv_rho(domain    *line,
               const      string s,
               const bool Lt,
               const bool Lo) {
 
    domn          = line;
    var_name      = s;
    L_transported = Lt;
    L_output      = Lo;
    d             = vector<double>(domn->ngrd, domn->pram->rho0);
 
}
 
 
void dv_rho::merge2cells(const int    imrg,
                         const double m1,
                         const double m2,
                         const bool   LconstVolume) {
 
    domn->domc->setGasStateAtPt(imrg);
    d.at(imrg) = domn->gas->density();
    d.erase(d.begin() + imrg+1);
 
}
 
 
void dv_rho::setVar(const int ipt){
 
    d.resize(domn->ngrd, domn->pram->rho0);
    if(ipt == -1)
        for(int i=0; i<domn->ngrd; i++) {
            domn->domc->setGasStateAtPt(i);
            d.at(i) = domn->gas->density();
        }
    else {
        domn->domc->setGasStateAtPt(ipt);
        d.at(ipt) = domn->gas->density();
    }
}