SUBROUTINE HISCTR(IH) INTEGER*4 I2N CALL HSTAT(IH,IXSIZ,IYSIZ,IXS,IX0,IYS,IY0) IYSIZ = IYSIZ - 1 IF(IYSIZ .LT. 0)IYSIZ = 0 IF(IXSIZ .EQ. 0)RETURN MINX = 2 MAXX = IXSIZ - 1 MINY = 1 MAXY = IYSIZ - 1 ORG = IY0 II = 0 DO 100 I = 0,IYSIZ II = II + 1 SUM = 0. SUMX = 0. SUMX2 = 0. XX = IX0 DO 50 J = MINX,MAXX XX = XX + IXS CALL HGET(IH,I2N,J+I*IXSIZ,1,IJ) SUM = SUM + I2N SUMX = SUMX + I2N*XX SUMX2 = SUMX2 + I2N*XX*XX 50 CONTINUE IF(SUM .LE. 0.)GO TO 100 CENTR = SUMX/SUM WRITE(6,1000)ORG,SUM,CENTR 100 ORG = ORG + IYS IF(IYSIZ .LT. 10)RETURN ORG = IX0 DO 200 I = 1,IXSIZ SUM = 0. SUMX = 0. SUMX2 = 0. XX = IY0 DO 150 J = MINY,MAXY XX = XX + IYS CALL HGET(IH,I2N,I + J*IXSIZ,1,IJ) SUM = SUM + I2N SUMX = SUMX + I2N*XX SUMX2 = SUMX2 + I2N*XX*XX 150 CONTINUE IF(SUM .EQ. 0.)GO TO 200 CENTR = SUMX/SUM WRITE(6,1001)ORG,SUM,CENTR 200 ORG = ORG + IXS 1000 FORMAT(T10'Y ='F12.0, T30'SUM ='F12.0, T50'CENTR.='1P,G12.5) 1001 FORMAT(T10'X ='F12.0, T30'SUM ='F12.0, T50'CENTR.='1P,G12.5) END