home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol266 / tmst.for < prev    next >
Encoding:
Text File  |  1986-05-19  |  1.6 KB  |  37 lines
  1.         SUBROUTINE      TMST(X,T,N,NSET)
  2. C////////////////////////////////////////////////////////////////
  3. C/                                                              /
  4. C/      Program-id.     TMST                                    /
  5. C/      Date-written.   Jan. 16th 1984                          /
  6. C/      File-name.      TMST.FOR                                /
  7. C/      Remarks.        Subroutine TMST.FOR page 76.            /
  8. C/                      This subroutine collects sample data    /
  9. C/                      on observations of a variable made over /
  10. C/                      a period of time.                       /
  11. C/                                                              /
  12. C////////////////////////////////////////////////////////////////
  13. C       
  14.         INTEGER*4       NSET(6,1)
  15. C
  16.       COMMON /C1/ID,IM,INIT,JEVNT,JMNIT,MFA,MSTOP,MX,MXC,NCLCT,NHIST,
  17.      1           NOQ,NORPT,NOT,NPRMS,NRUN,NRUNS,NSTAT,OUT,SCALE,ISEED,
  18.      2           TNOW,TBEG,TFIN,MXX,NPRNT,NCRDR,NEP,VNQ(4)
  19. C
  20.       COMMON /C2/ATRIB(4),ENQ(4),INN(4),JCELS(5,22),KRANK(4),MAXNQ(4),
  21.      1           MFE(4),MLC(4),MLE(4),NCELS(5),NQ(4),PARAM(20,4),
  22.      2           QTIME(4),SSUMA(10,5),SUMA(10,5),NAME(6),NPROJ,MON,
  23.      3         NDAY,NYR,JCLR
  24. C
  25. C
  26.         IF (N.GT.0) GO TO 20
  27.    10   CALL    ERROR(91,NSET)
  28.    20   IF (N .GT. NSTAT) GO TO 10
  29.         TT = T - SSUMA(N,1)
  30.         SSUMA(N,1) = SSUMA(N,1) + TT
  31.         SSUMA(N,2) = SSUMA(N,2) + X*TT
  32.         SSUMA(N,3) = SSUMA(N,3) + X*X*TT
  33.         SSUMA(N,4) = AMIN1(SSUMA(N,4),X)
  34.         SSUMA(N,5) = AMAX1(SSUMA(N,5),X)
  35.         RETURN
  36.         END
  37.