*-- Author : Girish D. Patel 07/06/93 SUBROUTINE FDEADR ************************************************************************ * * * * ************************************************************************ PARAMETER (NX=48) PARAMETER (NY=36) DIMENSION CON(NX,NY),YMAX(3),YTOT(3),ZMAX(NY) CALL HUNPAK(40,CON,'HIST',1) CALL VZERO(ZMAX,NY) CALL VZERO(YMAX,3) CALL VZERO(YTOT,3) DO 10 K = 1 , 3 DO 10 J = (K-1)*12+1 , K*12 DO 10 I = 1 , NX IF(CON(I,J).GT.YMAX(K)) THEN YMAX(K) = CON(I,J) ENDIF IF(CON(I,J).GT.ZMAX(J)) THEN ZMAX(J) = CON(I,J) ENDIF YTOT(K) = YTOT(K) + CON(I,J) 10 CONTINUE DO 20 K = 1 , 3 YMEAN = YTOT(K)/NX/12. DO 20 J = (K-1)*12+1 , K*12 DO 20 I = 1 , NX * CON(I,J) = CON(I,J)/YMAX(K) CON(I,J) = CON(I,J)/ZMAX(J) IF(CON(I,J).LT.0.02) CALL SHD(27,0,FLOAT(I-1),FLOAT(J)) CALL SHDW(28,0,FLOAT(I-1),FLOAT(J),CON(I,J)) CALL SHS (29,0,CON(I,J)) 20 CONTINUE 100 RETURN END *