Initial commit.

1 files changed, 767 insertions, 0 deletions

diff --git a/sources/over1.cpp b/sources/over1.cpp
new file mode 100644
index 0000000..35b6ced
--- /dev/null
+++ b/sources/over1.cpp
@@ -0,0 +1,767 @@
+//-*- mode: c++; indent-tabs-mode: t; coding: utf-8; show-trailing-whitespace: t -*-
+
+// file over1.cpp
+
+#include "over1.hpp"
+
+namespace over1 {
+
+  // subroutine over1.
+  void over1(void)
+  {
+    long int i;
+    long int j;
+    long int ll1, ll8, ll80;
+    long int n1, n5, n6, n7, n9;
+    long int nIAdQQ;
+    long int nIDA;
+    long int nIFKC;
+    long int nIP;
+    long int nIsWent;
+    long int nIY;
+    long int nNFDBR;
+    long int nNFDHST;
+    long int nNFDPH;
+    long int nNFDPol;
+    long int nNGroup;
+    long int nNK;
+    long int nNMAuto;
+    long int nNTLin;
+    long int nNTurn;
+    long int nNumBCO;
+    long int nNumNam;
+    long int pLSTACS[] = { 20, 90, 100, 20, 30, 250, 350, 60 };
+    double d1, d2, d3, d13;
+    double nPU;
+    double nZNVRef;
+    void *pLU2 = nullptr;
+    void *pLU6 = nullptr;
+    std::string sBus;
+    std::string sDiskFileName;
+    std::string sText2("name  ");
+    std::string sText6("copy  ");
+    std::string sText1("tacs o");
+    std::string sText3("tacs h");
+    std::string sText4("tacs s");
+    std::string sText5("4     ");
+    std::string sText7("tacs  ");
+    //
+    auto start_again = [ & ](void)
+    {
+      long int j;
+      //
+      blkcom::pIAlter = blkcom::pLFiles[ 1 ];
+      if(blkcom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << "+Continue partially-completed data case." << std::endl;
+      if(blkcom::sTexCol[ 12 ] == sText5[ 0 ])
+	goto a2861;
+      blkcom::sABuff >> sDiskFileName;
+      for(j = 1; j <= 2; j++) {
+	n7 = sDiskFileName.find(',');
+	sDiskFileName.replace(0, n7, " ");
+      }
+      for(j = n7; j <= 80; j++)
+	if(sDiskFileName[ j - 1 ] != ' ')
+	  break;
+      std::cout << " name of old pl4 file  disk_file(j : 80) = " << sDiskFileName.substr(j, 80) << std::endl;
+      blkcom::pLFiles[ 76 ] = nullptr;
+      blkcom::pLFiles[ 76 ] = (void *) new std::ifstream;
+      if(blkcom::pLFiles[ 76 ]) {
+	reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ]) -> open(sDiskFileName.substr(j, 80).c_str(), std::ios::in);
+      }
+      std::cout << " 1st record.  numnam, numnvo, numbco, numbrn =" << nNumNam << dekplt::nNumNVO << nNumBCO << dekplt::nNumBrn << std::endl;
+      if(((blkcom::nLBus + blkcom::pSStat -> m_nSize7) >= nNumNam) && \
+	 (blkcom::pSStat -> m_nSize12 >= dekplt::nNumNVO) &&	      \
+	 (blkcom::nLBranch >= dekplt::nNumBrn))
+	goto a2859;
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " Temporary error stop in \"over1\".  Program dimensioning is inadequate.  numnam, numnvo, numbrn =" << std::setfill(' ') << std::setw(8) << blkcom::nNumName << blkcom::nNumNVO << blkcom::nNumBRN << std::endl;
+      if(blkcom::pLFiles[ 76 ]) {
+	reinterpret_cast<std::ifstream *>(blkcom::pLFiles) -> close();
+	delete blkcom::pLFiles[ 76 ];
+      }
+      stoptp();
+
+    a2859:
+      reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ]) -> seekg(0, std::ios::beg);
+      n8 = blkcom::nLBus;
+      if(blkcom::nLBus >= blkcom::nNumNam)
+	n8 = blkcom::nNumNam - 1;
+      std::cout << " New name logic. n8, numnam, lbus =" << n8 << blkcom::nNumNam << blkcom::nLBus << std::endl;
+      (*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) << datexx << tcloxx << blkcom::nNumNam << blkcom << nNumNVO << blkcom::nNumBCO << blkcom::nNumBRN;
+      for(j = 1; j <= n8; j++)
+	(*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> blkcom::sBus[ j - 1 ];
+      std::cout << " bus(1 : n8) =";
+      for(j = 1; j <= n8; j++)
+	std::cout << blkcom::sBus[ j - 1 ] << " ";
+      std::cout << std::endl;
+
+    a2859:
+      reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ]) -> seekg(0, std::ios::beg);
+      blkcom::nNum888 = blkcom::nNumNam - n8;
+      std::cout << " Ready for final, full read with  num888 =" << blkcom::nNum888 << std::endl;
+      (*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> datexx >> tcloxx >> blkcom::nNumNam >> blckom::nNumNVO >> blkcom::nNumBCO >> blkcom::nNumBRN;
+      for(j = 1; j <= n8; j++)
+	(*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> blkcom::sBus[ j - 1 ];
+      for(j = 1; j <= blkcom::nNum888; j++)
+	(*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> blkcom::pTexVec[ j - 1 ];
+      std::cout << " After all named are read.  Next, ibsout." << std::endl;
+      if(blkcom::nNumNVO > 0)
+	for(j = 1; j <= blkcom::nNumNVO)
+	  (*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> blkcom::pIBSOut[ j - 1 ];
+      std::cout << " After ibsout.  Next, ...." << std::endl;
+      if(blkcom::nNumBRN > 0)
+	for(j = 1; j <= blkcom::nNumBRN)
+	  (*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> blkcom::sMBus[ j - 1 ];
+      std::cout << " After kbus, mbus." << std::endl;
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << blkcom::sDate1 << blkcom::sTClock << blkcom::nNumNam << blkcom::nNumNVO << blkcom::nNumBCO << blkcom::nNumBRN;
+      for(j = 1; j <= n8; j++)
+	(*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << blkcom::sBus[ j - 1 ];
+      for(j = 1; j <= blkcom::nNum888; j++)
+	(*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << blkcom::sTexVec[ j - 1 ];
+      std::cout << " over1, lunit4 date and time =" << blkcom::pDate1[ 0 ] << blkcom::pDate1[ 1 ] << " " << blkcom::pTClock[ 0 ] << blkcom::pTClock[ 1 ] << std::endl;
+      if(blkcom::nNumNVO > 0)
+	for(j = 1; j <= blkcom::nNumNVO; j++)
+	  (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << blkcom::pIBSOut[ j - 1 ];
+      if(blkcom::nNumBRN > 0)
+	for(j = 1; j <= blkcom::nNumBRN; j++)
+	  (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << blkcom::pKBus[ j - 1 ] << blkcom::pMBus[ j - 1 ];
+      n18 = blkcom::nNumNVO + blkcom::nNumBRN + 1;
+      std::cout << " Enter loop over numbers.  n18 =" << n18 << std::endl;
+      j = 1;
+      FOREVER {
+	for(k = 1; k <= n18; k++) {
+	  (*reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ])) >> labcom::sR4[ k - 1 ];
+	  if(reinterpret_cast<std::ifstream *>(blkcom::pLFiles[ 76 ]) -> eof())
+	    goto a6539;
+	}
+	if((blkcom::nIprsUp == 7) || (blkcom::nIprsUp > 9))
+	  std::cout << " j, r4(1) =" << j << labcom::sR4[ 0 ];\
+	if(sR4[ 0 ] == -9999.0e0)
+	  break;
+	for(k = 1; k <= n18; k++)
+	  (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 3 ])) << labcom::sR4[ k - 1 ];
+	++j;
+      }
+
+    a6539:
+      std::cout << " Done transferring lun77 to lunit(4).  j =" << j << std::endl;
+      runtym(d1, d2);
+      n18 = pLocker[ 0 ];
+      n19 = pLocker[ 1 ];
+      pfatch();
+      tables();
+      std::cout << " n18, n19, locker(1), locker(2) =" << n18 << n19 << pLocker[ 0 ] << pLocker[ 1 ] << std::endl;
+      blkcom::pFLStat[ 0 ] = -d1;
+      blkcom::pFLStat[ 1 ] = -d2;
+      if((n18 == pLocker[ 0 ]) && (n19 == pLocker[ 1 ]))
+	goto a2863;
+      blkcom::nKill = 201;
+      blkcom::sStat.m_nTACST = 2861;
+    };
+    auto absolute_TACS_dimensions = [ & ](void)
+    {
+      if(blkcom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << "+Set absolute TACS list size limits." << std::endl;
+      // To do.
+    };
+    auto relative_TACS_dimensions = [ & ](void)
+    {
+      if(blkcom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pLFile6)) << "+Proportional allocation of total TACS storage." << std::endl;
+      cimage();
+      if(blkcom::nKolBeg > 0)
+	goto a7120;
+      expchk(ll1, ll80, ll8);
+      if(blkcom::nKill > 0)
+	goto a9200;
+      for(i = 1; i <= 10; i++)
+	blkcom::sABuff >> labcom::sVoltBC[ i - 1 ];
+      goto a7130;
+
+    a7120:
+      blkcom::nNFrFld = 10;
+      freeOne(sVoltBC);
+
+    a7130:
+      if(blkcom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) "+Relative list sizes." << labcom::sVoltBC[ 0 ] << labcom::sVoltBC[ 1 ] << labcom::sVoltBC[ 2 ] << std::endl;
+      d1 = 0.0e0;
+      for(i = 1; i <= 8; i++)
+	d1 += labcom::sVoltBC[ i - 1 ];
+      d1 = blckom::nLTACST * blkcom::pNByte[ 2 ] / d1;
+      pLSTACS[ 0 ] = (long int) (labcom::sVoltBC[ 0 ] * d1 / (4 * blkcom::pNByte[ 2 ] + 8 * blkcom::pNByte[ 3 ]));
+      pLSTACS[ 1 ] = (long int) (labcom::sVoltBC[ 1 ] * d1 / (2 * blkcom::pNByte[ 2 ] + blkcom::pNByte[ 3 ]));
+      pLSTACS[ 2 ] = (long int) (labcom::sVoltBC[ 2 ] * d1 / (2 * blkcom::pNByte[ 3 ]));
+      pLSTACS[ 3 ] = (long int) (labcom::sVoltBC[ 3 ] * d1 / (5 * blkcom::pNByte[ 2 ] + blkcom::pNByte[ 3 ]));
+      pLSTACS[ 4 ] = (long int) (labcom::sVoltBC[ 4 ] * d1 / (3 * blkcom::pNByte[ 3 ]));
+      pLSTACS[ 5 ] = (long int) (labcom::sVoltBC[ 5 ] * d1 / blkcom::pNByte[ 3 ]);
+      PLSTACS[ 6 ] = (long int) (labcom::sVoltBC[ 6 ] * d1 / blkcom::pNByte[ 2 ]);
+      pLSTACS[ 7 ] = (long int) (labcom::sVoltBC[ 7 ] * d1 / (6 * blkcom::pNByte[ 2 ] + 2 * blkcom::pNByte[ 3 ]));
+    };
+    auto read_input_card = [ & ](void)
+    {
+      cimage();
+      for(i = 1; i <= 14; i++)
+	blkcom::sABuff >> blkcom::sAUpper[ i - 1 ];
+      if(strcom::toLower(blkcom::sAUpper[ 0 ]) == sText1)
+	goto a2697;
+      if(strcom::toLower(blkcom::sAUpper[ 0 ]) == sText3)
+	goto a2697;
+      if(strcom::toLower(blkcom::sAUpper[ 0 ]) == sText4)
+	goto a2699;
+      if(strcom::toLower(blkcom::sAUpper[ 0 ]) == sText6)
+	goto a7722;
+      (*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << " Begin TACS. " << std::string(27, ' ') << std::endl;
+      blkcom::nNewTAC = 1;
+      blkcom::nNTCSEX = 1;
+      blkcom::nNIUNRS = 1;
+      std::cout << " Prepare to call ntacs1 from over1." << std::endl;
+      ntacs1();
+      std::cout << " Back from ntacs1, back in over1." << std::endl;
+      goto a4284;
+
+    a2697:
+      blkcom::nNTCSEX = 1;
+
+    a2699:
+      n1 = 1;
+      if(blkcom::nNOutPR != 0)
+	goto a2691;
+      if(blckom::nNTCSEX == 0)
+	goto a22699;
+      if(blckom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << "+TACS hyvrid setup.  TACS data cards follow." << std::endl;
+      // 1984 format(18x, i2)
+
+      char *pTemp[ 20 ];
+      blkcom::sABuff.read(pTemp, 18);
+      blkcom::sABuff >> blkcom::sLStat[ 51 ];
+      goto a2691;
+
+    a4281:
+      if(n1 == 0)
+	goto a4284;
+      move(blkcom::pLSTACS, &blkcom::pLStat[ 60 ], ll8);
+      blkcom::nKTab = 1;
+      tacs1();
+      if(blkcom::nKill > 0)
+	goto a9200;
+      if((blkcom::nNTCSEX + blkcom::nNSTACS) > 0)
+	goto a4276;
+      blkcom::nIndStp = 1;
+      blkcom::LimStp = blkcom::pKprChg[ 0 ];
+      runtym(d1, d2);
+      blkcom::pFLStat[ 0 ] += d1;
+      blkcom::pFLStat[ 1 ] += d2;
+      blkcom::pFLStat[ 6 ] -= d1;
+      blkcom::pFLStat[ 7 ] -= d2;
+      blkcom::nChain = 12;
+    };
+    //
+    // over1 main function body.
+    //
+    ll8 = LL8;
+    ll80 = LL80;
+    //
+    if(blkcom::nIprsUp >= 1)
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << "Begin module over1." << std::endl;
+    nPU = std::numeric_limits<double>::infinity();
+    movecopy::move0<long int>(blkcom::pKTRLSW, LL6);
+    movecopy::move0<long int>(blkcom::pIpnTV, LL11);
+    movecopy::move0<long int>(blkcom::pLoopSS, LL11);
+    blkcom::sMonCar.m_nUserSST = 0;
+    nNMAuto = 0;
+    umcom::pBUSUM[ 0 ] = sText2;
+    blkcom::nNRight = 0;
+    blkcom::nIofgnd = 0;
+    blkcom::nMolDat = 0;
+    blkcom::nIStep = 0;
+    blkcom::sMonCar.m_nKLoaEP = 0;
+    blkcom::nTEnerg = 1e+20;
+    blkcom::sBranch = sText2;
+    blkcom::sCopy = sText6;
+    blkcom::nIAlter = 0;
+    blkcom::nIsPrin = 0;
+    blkcom::nIsPlot = 0;
+    blkcom::nIndStp = 1;
+    blkcom::nNOutPR = 0;
+    blkcom::pLFiles[ 0 ] = reinterpret_cast<void *>(&std::cerr);
+    blkcom::pLFiles[ 4 ] = reinterpret_cast<void *>(&std::cin);
+    blkcom::pLFiles[ 5 ] = reinterpret_cast<void *>(&std::cout);
+    blkcom::pKFile6 = blkcom::pLFiles[ 5 ];
+    blkcom::pPeakND[ 0 ] = 0.0;
+    blkcom::nKBurro = 0;
+    // Assign default relative precision for 6 EMTP variable types.
+    // 1 --- alphanumeric                 2 --- complex
+    // 3 --- floating-point numeric       4 --- integer numeric
+    // 5 --- 3 of 'subr31'                6 --- 4 of 'subr31'
+    // setting all equal to unity means that all have equal length.
+    for(j = 1; j <= 6; j++)
+      blkcom::pNByte[ j - 1 ] = 1;
+    emtp::dimens(blkcom::sLStat,					\
+		 blkcom::nChain,					\
+		 blkcom::sBus1,						\
+		 blkcom::sBus2);
+    /*
+    blkcom::nLBus = blkcom::sLStat[ 0 ];
+    blkcom::nLBrnch = blkcom::sLStat[ 1 ];
+    blkcom::nLData = blkcom::sLStat[ 2 ];
+    blkcom::nLExct = blkcom::sLStat[ 3 ];
+    blkcom::nLYMat = blkcom::sLStat[ 4 ];
+    blkcom::nLSwtch = blkcom::sLStat[ 5 ];
+    blkcom::nLSize7 = blkcom::sLStat[ 6 ];
+    blkcom::nLPast = blkcom::sLStat[ 7 ];
+    blkcom::nLNonL = blkcom::sLStat[ 8 ];
+    blkcom::nLChar = blkcom::sLStat[ 9 ];
+    blkcom::nLSMOut = blkcom::sLStat[ 10 ];
+    blkcom::nLSize12 = blkcom::sLStat[ 11 ];
+    blkcom::nLFDep = blkcom::sLStat[ 12 ];
+    blkcom::nLWT = blkcom::sLStat[ 13 ];
+    blkcom::nLTails = blkcom::sLStat[ 14 ];
+    blkcom::nLIMass = blkcom::sLStat[ 15 ];
+    blkcom::nLSyn = blkcom::sLStat[ 16 ];
+    blkcom::nMaxPE = blkcom::sLStat[ 17 ];
+    blkcom::nLTCST = blkcom::sLStat[ 18 ];
+    blkcom::nLFSem = blkcom::sLStat[ 19 ];
+    blkcom::nLFD = blkcom::sLStat[ 20 ];
+    blkcom::nLHist = blkcom::sLStat[ 21 ];
+    blkcom::nLSize23 = blkcom::sLStat[ 22 ];
+    blkcom::nLComp = blkcom::sLStat[ 23 ];
+    blkcom::nLSPCum = blkcom::sLStat[ 24 ];
+    blkcom::nLSize26 = blkcom::sLStat[ 25 ];
+    blkcom::nLSize27 = blkcom::sLStat[ 26 ];
+    blkcom::nLSize28 = blkcom::sLStat[ 27 ];
+    */
+    //    memcpy((void *) &blkcom::sSStat,				\
+    //	   (void *) blkcom::sLStat.data(),				\
+    //sizeof(blkcom::tStat));
+    blkcom::pSStat = blkcom::sLStat.data()
+    // Assign 'n1' equal to the number of EMTP lists of variable
+    // dimensioning (of the solution overlays).
+    //n1 = 27;
+    //blkcom::nLTLabl = blkcom::sLStat[ n1 + 1 ];
+    comlock::pLocker[ 0 ] = blkcom::sBus1;
+    comlock::pLocker[ 1 ] = blkcom::sBus2;
+    d13 = blkcom::pSStat -> m_nTACST;
+    d13 /= 1600.;
+    for(j = 1; j <= 8; j++)
+      pLSTACS[ j - 1 ] *= d13;
+    n1 = -9999;
+    movecopy::copy<long int>(n1, blkcom::sLStat, LL60);
+    sysdep();
+    movecopy::move0<double>(blkcom::sFLStat, LL20);
+    date_time::runtym(d1, d2);
+    blkcom::sFLStat[ 0 ] -= d1;
+    blkcom::sFLStat[ 1 ] -= d2;
+    blkcom::nTenM6 = pow(blkcom::nTenM3, 2.0);
+    blkcom::nMaxZNO = 50;
+    blkcom::nEpsZNo = blkcom::nEpsilon;
+    blkcom::nEpWarn = blkcom::nTenM3;
+    blkcom::nEpStop = 0.1;
+    blkcom::pZNoLim[ 0 ] = 1.0;
+    blkcom::pZNoLim[ 1 ] = 1.5;
+    blkcom::sMonCar.m_nKBase = 0;
+    blkcom::sMonCar.m_nISW = 0;
+    blkcom::pKPrChg[ 0 ] = -7777;
+    blkcom::nJFLSOS = 0;
+    // rewind lunit1
+    // rewind lunit3
+    // rewind lunit4
+    blkcom::sMonCar.m_nKnt = 1;
+    blkcom::sMonCar.m_pMTape = blkcom::pLFiles[ 4 ];
+    blkcom::nIFDep = 0;
+    blkcom::nNTCSex = 0;
+    blkcom::nNSTACS = 0;
+    blkcom::nKanal = 0;
+
+    // Initialize KPARTB with nonzero value so later minus
+    // sign can be applied as flag of "renumber bypass" usage:
+    blkcom::nKPartB = 1000;
+
+    // SYSDEP will redefine if low-prec. Complex ("cable constants"):
+    nZNVRef = 0.0;
+    blkcom::nFMaxFS = 0.0;
+    blkcom::pBegMax[ 0 ] = 0.0;
+    blkcom::pBegMax[ 1 ] = 0.0;
+    blkcom::pBegMax[ 2 ] = std::numeric_limits<double>::infinity();
+    blkcom::nKTRef = 0;
+    blkcom::nNComp = 0;
+    blkcom::nNumSM = 0;
+    blkcom::nNumUM = 0;
+    blkcom::nNumOut = 0;
+    blkcom::nKill = 0;
+    blkcom::nNum99 = 0;
+    blkcom::nModOut = 0;
+    blkcom::pInEcho = nullptr;
+    sBus = blkcom::sBlank;
+    blkcom::nIAverg = 0;
+    movecopy::move0<int>(labcom::sISourc, blkcom::sSStat.m_nSwitch);
+    movecopy::move0<int>(labcom::sKodeBR, blkcom::sSStat.m_nBranch);
+    movecopy::move0<int>(labcom::sKodSem, blkcom::sSStat.m_nBranch);
+    movecopy::move0<int>(labcom::sLength, blkcom::sSStat.m_nBranch);
+    movecopy::move0<int>(labcom::sIndHst, blkcom::sSStat.m_nBranch);
+    movecopy::move0<double>(labcom::sBValue, blkcom::sSStat.m_nSize12);
+    nIsWent = 1;
+    blkcom::nOmega = 0.0;
+    blkcom::nDegMin = 0.0;
+    blkcom::nDegMax = 0.0;
+    blkcom::nSglFir = 2.0 * M_PI;
+    blkcom::nJST = 0;
+    blkcom::nJST1 = 1;
+    blkcom::nIFDep2 = 1;
+    blkcom::nIAdd = 0;
+    blkcom::nIFSem = 0;
+    //
+    com2::nIDM = 1;
+    com2::nIDQ = 1;
+    com2::nIDU = 1;
+    com2::nIDT = 1;
+    nIAdQQ = 1;
+    com2::nIQ = 0;
+    com2::nIDA = 0;
+    com2::nIDY = 0;
+    com2::nIFKC = 0;
+    nNK = 0;
+    nNGroup = 0;
+    nNTurn = 0;
+    nNFDBR = 0;
+    nNFDPH = 9;
+    nNFDHST = 5400;
+    // nNFDPOL = 25;
+    nNFDPol = 30;
+    nNTLin = 0;
+    // pLU2 = ...
+    pLU6 = blkcom::pLFiles[ 5 ];
+
+  a5223:
+    blkcom::nIFX = 0;
+    nIP = 2;
+    nIY = 2;
+    blkcom::nICat = 0;
+    if(blkcom::nNOutPR != 0)
+      goto a15;
+    if(blkcom::nKol132 == 132)
+      goto a6452;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " *********  Begin \"M40.\" EMTP solution.";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << "   Size /label/ =";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(7) << std::setfill(' ') << std::dec;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << blkcom::sSStat.m_nTLabl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << "  integer words." << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " list limits  1-10 :";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nBus;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nBranch;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nData;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nExct;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nYMat;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nSwitch;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nSize7;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nPast;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nNOnL;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nChar << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " list limits 11-20 :";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nSMOut;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nSize12;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nFDep;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nWT;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nTails;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nIMass;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nSyn;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nMaxPE;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nTACST;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(6) << std::dec << blkcom::sSStat.m_nFSem << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " list limits 21-end:";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nFD;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nHist;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize23;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nComp;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSPCum;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize26;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize27;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize28 << std::endl;
+    goto a15;
+
+  a6452:
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " Associated user documentation is the 864-page EMTP rule book dated  June, 1984.   Version M43.   Vardim time/date =";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << comlock::pLocker[ 0 ] << comlock::pLocker[ 1 ] << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " Independent list limits follow.    Total length of /label/  equals ";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(8) << std::dec << blkcom::sSStat.m_nTLabl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << "  integer words.   ";
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nBus;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nBranch;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nData;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nExct;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nYMat;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSwitch << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " " << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize7;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nPast;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nNOnL;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nChar;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSMOut;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize12;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nFDep;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nWT;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nTails;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nIMass;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSyn;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nMaxPE;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nTACST;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nFSem;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nFD;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nHist;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize23;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nComp;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSPCum;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize26;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize27;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setw(6) << std::setfill(' ') << std::dec << blkcom::sSStat.m_nSize28 << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " " << std::string(51, '-') << "+" << std::string(80, '-') << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " Descriptive interpretation of new-case input data  |" << std::flush;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " input data card images printed below, all 80 columns, character by character." << std::endl << std::flush;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " " << std::string(51, ' ') << "|" << std::flush;
+    for(i = 1; i <= 8; i++) {
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::string(9, ' ') << std::flush;
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::setfill(' ') << std::setw(1) << std::dec << i << std::flush;
+    }
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::endl;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " " << std::string(51, ' ') << "|" << std::flush;
+    for(i = 1; i <= 8; i++)
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::string(9, ' ') << "0" << std::flush;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << std::endl << std::flush;
+    (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << " " << std::string(51, '-') << "+" << std::string(80, '-') << std::endl;
+
+  a15:
+    FOREVER {
+      emtp::cimage();
+      if(blkcom::nKill > 0)
+	goto a9200;
+      blkcom::nNRight = -2;
+      n9 = blkcom::nKolBeg;
+      blkcom::nKolBeg = 1;
+      blkcom::sTexCol = blkcom::sABuff.str().substr(0, 80);
+      freedom::freeOne(d1);
+
+    a3247:
+      blkcom::nNRight = 0;
+      if(n9 == -std::numeric_limits<long int>::infinity())
+	blkcom::nKolBeg = n9;
+      if(blkcom::nNFrFld > 0)
+	goto a3280;
+      if(blkcom::nNOutPR == 0)
+	(*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << "+Blank termination-of-run card." << std::endl;
+      interp();
+      blkcom::nIVolt = 7777;
+
+    a3273:
+      blkcom::nKill = 9999;
+      blkcom::nChain = 31;
+      goto a9800;
+
+    a3280:
+      reques();
+      i = sStat.m_nMaxPE;
+      switch(i) {
+      case 0:
+
+      a2843:
+	if((blkcom::nNOutPR != 0) && (blkcom::nIprsUp > 0))
+	  blkcom::nNOutPR = 0;
+	blkcom::nXOpt = blkcom::nStatFR;
+	blkcom::nCOpt = blkcom::nStatFR;
+	blkcom::nKolBeg = n9;
+	if(blkcom::nKolBeg > 0)
+	  goto a4201;
+	utilities::expchk(ll1, ll80, ll8);
+	if(blkcom::nKill > 0)
+	  goto a9200;
+	blkcom::sABuff >> blkcom::nDeltaT >> blkcom::nTMax >> d1 >> d2 >> d3 >> blkcom::nTolMat >> blkcom::nT;
+	if(blkcom::nT == 0.0)
+	  blkcom::nT = 0.0e0;
+	goto a4202;
+
+      a4201:
+	blkcom::nNFrFld = 1;
+	blkcom::nNRight = 0;
+	freedom::freeOne(blkcom::nDeltaT);
+	freedom::freeOne(blkcom::nTMax);
+	freedom::freeOne(d1);
+	freedom::freeOne(d2);
+	freedom::freeOne(d3);
+	freedom::freeOne(blkcom::nTolMat);
+	freedom::freeOne(blkcom::nT);
+
+      a4202:
+	if(blkcom::nNOutPR == 0)
+	  (*reinterpret_cast<std::ostream *>(blkcom::pKFile6)) << "+Misc. data." << std::setfill(' ') << std::setw(12) << blkcom::nDeltaT << blkcom::nTMax << d1 << std::endl;
+	if(blkcom::nIofbnd != 33666)
+	  goto a4206;
+	blkcom::nChain = 41;
+	blkcom::nXOpt = d1;
+	goto a9800;
+	break;
+
+      case 11:
+	goto a3273;
+	break;
+
+      case 15:
+	start_again();
+	goto a9200;
+	break;
+
+      case 28:
+	goto a3247;
+	break;
+
+      case 32:
+	absolute_TACS_dimensions();
+	goto a15;
+	break;
+
+      case 33:
+	relative_TACS_dimensions();
+	goto a15;
+	break;
+
+      case 34:
+	goto a523;
+	break;
+
+      case 38:
+	goto a15;
+	break;
+
+      case 40:
+	goto a4308;
+	break;
+      }
+      if(blkcom::nKill > 0)
+	goto a9200;
+      if(blkcom::nChain != 1)
+	goto a9800;
+    }
+
+  a4276:
+    blkcom::nChain = 12;
+    goto a9800;
+
+  a4284:
+    blkcom::nChain = 2;
+    goto a9800;
+
+  a9200:
+    blkcom::nChain = 51;
+    blkcom::sLStat[ 17 ] = 1;
+
+  a9800:
+    blkcom::nLastOV = 1;
+    n5 = (size_t) &nIDA;
+    n6 = (size_t) &nIFKC;
+    if(blkcom::nIprsUp >= 1)
+      (*reinterpret_cast<std::ostream *>(blkcom::pLFiles[ 5 ])) << "Exit module over1." << std::endl;
+
+  a99999:
+    return;
+  }
+
+  // subroutine tacs1c.
+  void tacs1c(void)
+  {
+    // To do.
+  }
+
+  // subroutine swmodf.
+  void swmodf(void)
+  {
+    // To do.
+  }
+
+  // subroutine reques.
+  void reques(void)
+  {
+    // To do.
+  }
+
+  // subroutine sysdep.
+  void sysdep(void)
+  {
+    std::string sCol("");
+    std::string sColXXX("plt00000000000.pl4");
+    std::string sLettra("a");
+    std::string sLettrb("b");
+    std::string sLettrc("c");
+    std::string sBusNM1("      ");
+    std::string sBusNM2("........");
+    std::string sBusNM3("terra");
+    std::string sText1 = ",";
+    std::string sText2 = "$";
+    //
+    sCol = sColXXX;
+    // To do.
+    blkcom::sBlank = sBusNM1;
+    blkcom::sTrash = sBusNM2;
+    blkcom::sTerra = sBusNM3;
+    blkcom::cCSepar = sText1[ 0 ];
+    blkcom::sChCont = sText2;
+    blkcom::nTenM3 = 1.0e-3;
+    blkcom::nUnity = 1.0e0;
+    blkcom::nOneHaf = 0.5e0;
+    blkcom::pNByte[ 0 ] = 2;
+    blkcom::pNByte[ 1 ] = 2;
+    blkcom::pKTRLSW[ 5 ] = 0;
+    blkcom::pNByte[ 2 ] = 2;
+    blkcom::nEpsilon = 1.0e-8;
+    //blkcom::nTwoPI = 2.0 * M_PI;
+    // To do.
+  }
+
+  // subroutine midov1.
+  void midov1(void)
+  {
+    // To do.
+  }
+
+  // subroutine nmincr.
+  void nmincr(std::string &, int &)
+  {
+    // To do.
+  }
+
+  // subroutine tacs1.
+  void tacs1(void)
+  {
+    // To do.
+  }
+
+  // subroutine tacs1a.
+  void tacs1a(void)
+  {
+    // To do.
+  }
+
+  // subroutine tacs1b.
+  void tacs1b(void)
+  {
+    // To do.
+  }
+
+  // subroutine expchk.
+  void expchk(int &, int &, int &)
+  {
+    // To do.
+  }
+
+  // subroutine intchk.
+  void intchk(int &, int &, int &)
+  {
+    // To do.
+  }
+
+  // subroutine pfatch.
+  void pfatch(void)
+  {
+    // To do.
+  }
+
+}
+
+// end of file over1.cpp