hermite now in critical section to avoid segfault, that should be fixed in future

utils/TextFunctions.h

@@ -344,7 +344,7 @@ static void parseStrA(const std::string& data, int pos1, int pos2, int& ret, int
         try {
             checkBounds(ret, min, max, strict);
         }
-        catch (Exp& e) {
+        catch (const std::exception& e) {
             const std::string errstr = pos2 ? (ToStr(pos1) + ":" + ToStr(pos2)) : (ToStr(pos1) + ":-");
             throw Exp() << "bad substr format in line [" + data + "] at position" << errstr << " " << part << ". " << e.what();
         }
@@ -370,7 +370,7 @@ static void parseStrA(const std::string& data, int pos1, int pos2, double& ret,
         try {
             checkBounds(ret, (double*)min, (double*)max, strict);
         }
-        catch (Exp& e) {
+        catch (const std::exception& e) {
             const std::string errstr = pos2 ? (ToStr(pos1) + ":" + ToStr(pos2)) : (ToStr(pos1) + ":-");
             throw Exp() << "bad substr format in line [" + data + "] at position" << errstr << " " << part << ". " << e.what();
         }

utils/ephinterp.cpp

@@ -87,16 +87,19 @@ int cinterp_pvac(
         if (j < 3) {
             if (fabs(ephi.outv[j + 3][0]) > 1.e-10) {
                 /* Use velocities if they exists (hermite) */
+#pragma omp critical (pos)
                 hermite(pindex, 2, ephi.jds, &ephi.outv[j][0], &ephi.outv[j + 3][0], ephi.nv, 10, itime, &outvec[j], &vz, ircode);
             }
             else {
                 /* Otherwise, ignore them (fvel = -1) or create them (fvel = 1) */
+#pragma omp critical (pos)
                 hermite(pindex, -fvel, ephi.jds, &ephi.outv[j][0], &z, ephi.nv, 10, itime, &outvec[j], &outvec[j + 3], ircode);
             }
         }
             /* Velocity */
             /* Interpolate if they exists */
         else if (fabs(ephi.outv[j][0]) > 1.e-10) {
+#pragma omp critical (pos)
             hermite(pindex, 1, ephi.jds, &ephi.outv[j][0], &z, ephi.nv, 10, itime, &outvec[j], &vz, ircode);
         }