home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol266 / sumry.for < prev    next >
Encoding:
Text File  |  1986-05-19  |  3.2 KB  |  99 lines
  1.         SUBROUTINE      SUMRY(NSET)
  2. C////////////////////////////////////////////////////////////////
  3. C/                                                              /
  4. C/      Program-id.     SUMRY                                   /
  5. C/      Date-written.   Jan. 16th 1984                          /
  6. C/      File-name.      SUMRY.FOR                               /
  7. C/      Remarks.        Subroutine SUMRY.FOR page 84.           /
  8. C/            Subroutine SUMRY is the basic output     /
  9. C/            routine of GASP II. It processes the    /
  10. C/            data collected in subroutine COLCT    /
  11. C/            TMST, and HISTO and then prints out     /
  12. C/            a data summary.                /
  13. C/                                                              /
  14. C////////////////////////////////////////////////////////////////
  15. C
  16. C    * ONE WORD INTEGERS EXCEPT NSET
  17. C       
  18.         INTEGER*4       NSET(6,1)
  19. C
  20.       COMMON /C1/ID,IM,INIT,JEVNT,JMNIT,MFA,MSTOP,MX,MXC,NCLCT,NHIST,
  21.      1           NOQ,NORPT,NOT,NPRMS,NRUN,NRUNS,NSTAT,OUT,SCALE,ISEED,
  22.      2           TNOW,TBEG,TFIN,MXX,NPRNT,NCRDR,NEP,VNQ(4)
  23. C
  24.       COMMON /C2/ATRIB(4),ENQ(4),INN(4),JCELS(5,22),KRANK(4),MAXNQ(4),
  25.      1           MFE(4),MLC(4),MLE(4),NCELS(5),NQ(4),PARAM(20,4),
  26.      2           QTIME(4),SSUMA(10,5),SUMA(10,5),NAME(6),NPROJ,MON,
  27.      3         NDAY,NYR,JCLR
  28. C
  29.     WRITE(NPRNT,21)
  30.    21    FORMAT(1H1,29X,'** GASP Summary Report ** '/)
  31.     WRITE(NPRNT,102) NPROJ,NAME,MON,NDAY,NYR,NRUN
  32.   102    FORMAT(20X,'Simulation Project No.',I4,2X,'on',2X,
  33.      1  6A2//,20X,'Date',I3,'/',I3,'/',I5,12X,'Run number',I5/)
  34.     IF (NPRMS) 147,147,146
  35.   146    DO  64  I=1,NPRMS
  36.    64    WRITE(NPRNT,107) I,(PARAM(I,J),J=1,4)
  37.   107    FORMAT(10X,' Parameter No.',I5,4F12.4)
  38.   147    IF (NCLCT) 5,60,66
  39.     5    WRITE(NPRNT,199)
  40.   199    FORMAT(///26X,'Error Exit, Type 98 Error.')
  41.     CALL    EXIT
  42.    66    WRITE(NPRNT,23)
  43.    23    FORMAT(//34X,'** Generated Data ** ',/17X,'Code',4X,'Mean',6X,
  44.      1    'STD.DEV.',5X,'Min.',7X,'Max.',5X,'OBS.'/)
  45. C
  46. C    --- Compute and print statistics gathered by CLCT
  47. C
  48.       DO 2 I=1,NCLCT
  49.     IF (SUMA(I,3)) 5,62,61
  50.    62    WRITE(NPRNT,63) I
  51.    63    FORMAT(17X,I3,10X,'No Values Recorded ')
  52.                     GO TO 2
  53.    61    XS = SUMA(I,1)
  54.     XSS = SUMA(I,2)
  55.     XN = SUMA(I,3)
  56.     AVG = XS / XN
  57.     STD = (((XN * XSS) - (XS * XS))/(XN * (XN - 1.0)))**0.5
  58.     N = XN
  59.     WRITE(NPRNT,24) I,AVG,STD,SUMA(I,4),SUMA(I,5),N
  60.    24    FORMAT(17X,I3,4F11.4,I7)
  61.     2 CONTINUE
  62.    60    IF (NSTAT) 5,67,4
  63.     4    WRITE(NPRNT,29)
  64.    29    FORMAT(/34X,'** Time Generated Data **'/,17X,'Code',4X,'Mean',
  65.      1  6X,'STD.DEV.',5X,'Min.',7X,'Max.',3X,'Total Time '/)
  66. C
  67. C    --- Compute and print statistics gathered by TMST
  68. C
  69.       DO 6 I=1,NSTAT
  70.     IF (SSUMA(I,1)) 5,71,72
  71.    71    WRITE(NPRNT,63) I
  72.                     GO TO 6
  73.    72    XT = SSUMA(I,1)
  74.     XS = SSUMA(I,2)
  75.     XSS = SSUMA(I,3)
  76.     AVG = XS / XT
  77.     STD = (XSS/XT - AVG*AVG) ** 0.5
  78.     WRITE(NPRNT,30) I,AVG,STD,SSUMA(I,4),SSUMA(I,5),XT
  79.    30    FORMAT(17X,I3,5F11.4)
  80.     6 CONTINUE
  81.    67    IF (NHIST) 5,75,9
  82.     9    WRITE(NPRNT,25)
  83.    25    FORMAT(/27X,'** Generated Frequency Distributions **',/17X,
  84.      1  'Code',20X,'Histograms')
  85. C
  86. C    --- Print histograms
  87. C
  88.       DO  12  I=1,NHIST
  89.     NCL = NCELS(I) + 2
  90.    12 WRITE(NPRNT,26) I,(JCELS(I,J),J=1,NCL)
  91.    26    FORMAT(/17X,I3,5X,11I4,/(25X,11I4))
  92. C
  93. C    --- Print files and file statistics
  94. C
  95.    75    DO  15  I=1,NOQ
  96.    15    CALL    PRNTQ(I,NSET)
  97.     RETURN
  98.     END
  99.