From 3e84d8a661e2f304c213e43442b1b03922780a3e Mon Sep 17 00:00:00 2001
From: Knut Morten Okstad <knut.morten.okstad@sintef.no>
Date: Thu, 16 Mar 2017 01:24:20 +0100
Subject: [PATCH] Added: An overloded writeGlvS method that writes a single
 scalar field. Equvialent to writeGlvV, but output as fringe plot instead of
 vectors.

---
 src/SIM/SIMoutput.C | 34 ++++++++++++++++++++++++++++++++++
 src/SIM/SIMoutput.h | 11 ++++++++++-
 2 files changed, 44 insertions(+), 1 deletion(-)

diff --git a/src/SIM/SIMoutput.C b/src/SIM/SIMoutput.C
index a16022f9..47be05ee 100644
--- a/src/SIM/SIMoutput.C
+++ b/src/SIM/SIMoutput.C
@@ -434,6 +434,40 @@ bool SIMoutput::writeGlvV (const Vector& vec, const char* fieldName,
 }
 
 
+bool SIMoutput::writeGlvS (const Vector& scl, const char* fieldName,
+                           int iStep, int& nBlock, int idBlock) const
+{
+  if (scl.empty())
+    return true;
+  else if (!myVtf)
+    return false;
+
+  Matrix field;
+  Vector lovec;
+  IntVec sID;
+
+  int geomID = myGeomID;
+  for (size_t i = 0; i < myModel.size(); i++)
+  {
+    if (myModel[i]->empty()) continue; // skip empty patches
+
+    if (msgLevel > 1)
+      IFEM::cout <<"Writing scalar field for patch "<< i+1 << std::endl;
+
+    myModel[i]->extractNodeVec(scl,lovec,1);
+    if (!myModel[i]->evalSolution(field,lovec,opt.nViz))
+      return false;
+
+    if (!myVtf->writeNres(field,++nBlock,++geomID))
+      return false;
+    else
+      sID.push_back(nBlock);
+  }
+
+  return myVtf->writeSblk(sID,fieldName,idBlock,iStep);
+}
+
+
 bool SIMoutput::writeGlvS (const Vector& psol, int iStep, int& nBlock,
                            double time, const char* pvecName,
                            int idBlock, int psolComps)
diff --git a/src/SIM/SIMoutput.h b/src/SIM/SIMoutput.h
index 1b3a2462..5dd90ff1 100644
--- a/src/SIM/SIMoutput.h
+++ b/src/SIM/SIMoutput.h
@@ -87,11 +87,20 @@ public:
   //! \param[in] fieldName Name identifying the vector field
   //! \param[in] iStep Load/time step identifier
   //! \param nBlock Running result block counter
-  //! \param[in] idBlock Starting value of result block numbering
+  //! \param[in] idBlock Result block ID number
   //! \param[in] ncmp Number of components in vector field
   bool writeGlvV(const Vector& vec, const char* fieldName,
                  int iStep, int& nBlock, int idBlock = 2, int ncmp = 0) const;
 
+  //! \brief Writes a scalar field for a given load/time step to the VTF-file.
+  //! \param[in] scl The scalar field to output (nodal values)
+  //! \param[in] fieldName Name identifying the scalar field
+  //! \param[in] iStep Load/time step identifier
+  //! \param nBlock Running result block counter
+  //! \param[in] idBlock Result block ID number
+  bool writeGlvS(const Vector& scl, const char* fieldName,
+                 int iStep, int& nBlock, int idBlock = 2) const;
+
   //! \brief Writes solution fields for a given load/time step to the VTF-file.
   //! \param[in] psol Primary solution vector
   //! \param[in] iStep Load/time step identifier