C PROGRAM YAG C CALCULATES YEARLY AVERAGES OF GAMMA implicit real*8(A-H,O-Z) DIMENSION YEAR(30),D0(30),E0(30),D(30,30),ED(30,30),P(30) # ,GAMM(30,20),EGAMM(30,20),DY(30),GA(20),WGA(30) # ,GASO(6,12),DMF(6),DML(6),AL(6,12) C ************************************************************************ NY=24 MY=5 OPEN(UNIT=8,FILE='mart',STATUS='OLD') OPEN(UNIT=9,FILE='DATASETS',STATUS='OLD') OPEN(UNIT=10,FILE='GBC',STATUS='OLD') OPEN(UNIT=11,FILE='EGBC',STATUS='OLD') OPEN(UNIT=12,FILE='GA',STATUS='OLD') OPEN(UNIT=13,FILE='EGA',STATUS='OLD') OPEN(UNIT=17,FILE='AGAM',STATUS='UNKNOWN') OPEN(UNIT=18,FILE='MADA',STATUS='UNKNOWN') DO M=1,MY READ(8,DMF(M),DML(M),(AL(M,K),K=1,10) ENDDO GA0=0. WGA0=0. DO K=1,9 GA(K)=0. WGA(K)=0. ENDDO M=1 Y=1996 DO J=1,NY READ(9,*) NJ,NS,YM,DY(J) READ(10,*)YEAR(J),D0(J) READ(11,*)NSE,NPA,E0(J) READ(12,*)YEAR(J),(D(J,K),K=1,9) READ(13,*)NSE,NMO,(ED(J,K),K=1,9) EDAT=YEAR(J)+DY(J) BDAT=YEAR(J)-DY(J) DE=DML(M)-BDAT IF(DE.LT.0.OR.BDAT.GT.Y) THEN GASO(M,1)=GA/WGA0 GA0=0. WGA0=0. DO K=1,9 GASO(M,K+1)=GA(K)/WGA(K) GA(K)=0. WGA(K)=0. ENDDO M=M+1 Y=Y+1 IF(M.GT.5) GO TO 100 ENDIF DE=DML(M)-BDAT DB=EDAT-DMF(M) IF(DB.GT.0) THEN H=1 IF(DB.LT.DY(J)) H=DB/DY(J) IF(DE.LT.DY(J)) H=DE/DY(J) WL=H/E0(J)**2 GA0=GA0+D0(J)*WL WGAO=WGA0+WL DO K=1,18 WL=H/ED(J,K)**2 GA(K)=GA(K)+D(J,K)*WL WGA(K)=WGA(K)+WL ENDDO ENDIF 100 CONTINUE GX=0 AX=0 DO M=1,5 IF(GASO(M,4).GT.GX) GX=GASO(M,4) IF(AL(M,4).GT.GX) GX=AL(M,4) ENDDO FAG=GX/AX DO M=1,5 WRITE(17,60) 1995+M,(GASO(M,K),K=1,10) WRITE(18,60) 1995+M,(AL(M,K)*FAG,K=1,NZO) ENDDO 60 FORMAT(I6,10F7.3) STOP END