readrec[h_, npar_]:=
(Check[jter=Read[h, Number],
      WriteString[out,"read error before, at or after iteration
                  record ", jter, "\n"];
      Quit[1]];
If[jter==EndOfFile,
   WriteString[out,"unexpected end of simulation input file before,
               at or after iteration record ", jter, "\n"];
   Quit[1];
  ];
Check[wtlog=Read[h, Number],
      WriteString[out,"read error before, at or after iteration
                  record ", jter, "\n"];
      Quit[1]]; 
If[wtlog==EndOfFile,
   WriteString[out,"unexpected end of simulation input file before,
               at or after iteration record ", jter, "\n"];
   Quit[1];
  ];
Check[pri=Read[h, Number],
      WriteString[out,"read error before, at or after iteration
                  record ", jter, "\n"];
      Quit[1]]; 
If[pri==EndOfFile,
   WriteString[out,"unexpected end of simulation input file before,
               at or after iteration record ", jter, "\n"];
   Quit[1];
  ];
Check[pdat=Read[h, Number],
      WriteString[out,"read error before, at or after iteration
                  record ", jter, "\n"];
      Quit[1]]; 
If[pdat==EndOfFile,
   WriteString[out,"unexpected end of simulation input file before,
               at or after iteration record ", jter, "\n"];
   Quit[1];
  ];
Check[par=ReadList[ h, Number, npar ],
      WriteString[out,"read error in simulation input file iteration
                  record ", jter, "\n"];
      Quit[1]];
If[Length[par] < npar,
   WriteString[out,"unexpected end of simulation input file in",
               jter, "\n"];
   Quit[1];
  ];)