public double[] getInitConcns(){
        if(!this.psi0set && !this.sigmaSet)unpack();
        double[] conc = Conv.copy(this.initConcn);
        for(int i=0; i<this.numOfIons; i++)conc[i] *= 1e-3;
        return conc;
    }