Nlte and p1 model rebase

run.sh

@@ -0,0 +1,25 @@
+#!/bin/bash
+
+cd test
+
+rm -R output-plasma
+rm -R output-plasma-refine2
+
+# cp Restart/lowP/restart_output-plasma.LTE.h5 restart_output-plasma.sol.h5
+cp Restart/highP/restart_output-plasma-refine2.LTE.h5 restart_output-plasma-refine2.sol.h5
+# cp Restart/highP/restart_output-plasma-refine2.NLTE.h5 restart_output-plasma-refine2.sol.h5
+
+
+
+EXE=../src/tps
+RUNFILE="inputs/plasma.ini"
+# RUNFILE="inputs/input.4iters.cyl.ini"
+
+# RUNFILE="inputs/perfectGas.air.ini"
+
+NPROC=10
+
+mpirun -np $NPROC $EXE --runFile $RUNFILE
+# mpirun -np $NPROC libtool --mode=execute gdb --args $EXE --runFile $RUNFILE
+
+

src/BCintegrator.cpp

@@ -360,6 +360,8 @@ void BCintegrator::AssembleFaceVector(const FiniteElement &el1, const FiniteElem
   const IntegrationRule *ir = &intRules->Get(Tr.GetGeometryType(), intorder);
   const int numActiveSpecies = mixture->GetNumActiveSpecies();
   const int nvel = mixture->GetNumVels();
+  const bool twoTemperature_ = mixture->IsTwoTemperature();
+  const int iTe = num_equation - 1;
 
   for (int i = 0; i < ir->GetNPoints(); i++) {
     const IntegrationPoint &ip = ir->IntPoint(i);
@@ -390,6 +392,18 @@ void BCintegrator::AssembleFaceVector(const FiniteElement &el1, const FiniteElem
         funval1(eq) = sum;
       }
 
+      if (twoTemperature_ && eq == iTe) {
+        funval1(eq) = max(sum,0.0); 
+      }
+
+      if (twoTemperature_ && eq == iTe) {
+        funval1(eq) = max(sum,0.0); 
+      }
+
+      if (twoTemperature_ && eq == iTe) {
+        funval1(eq) = max(sum,0.0); 
+      }
+
       // interpolation gradients
       for (int d = 0; d < dim; d++) {
         sum = 0.;