parameter(m=288,n=120,mnth=12*6,ny=6) !m:lon; n:lat; nt:time; ny:year,2004-2009
! m: Longitudes:  288 bins centered on 179.375W to 179.375E (1.25 degree steps)
! n: Latitudes:  120 bins centered on -59.5S to 59.5N (1.00 degree steps)
    real OZONE(M,N),OZONE2(M,N),O3MX(m,n,mnth),O3T(m,n,mnth)
    character chmn(12)*3 
	character chmd*3 !FOR READ FILES
    CHARACTER CHYR*2 !FOR WRITE FILES
    data chmn/'jan','feb','mar','apr','may','jun','jul','aug','sep','oct','nov','dec'/
  
  OPEN(10,file='d:\Tropozone\TROPOZONE2005-2009.dat', form='unformatted', &
        access='direct',recl=M*N)
  IREC = 0
  OZONE=9999.;OZONE2=9999.
  DO IYEAR=05,09
    WRITE(CHYR(1:2),'(i2.2)')IYEAR
  do imn=1,12
    WRITE(chmd(1:3),'(A3)')chmn(imn)
    OPEN(1000,ERR=3,file='d:\PPB\L3_tropo_ozone_vmr_'//chmd//''//chyr//'',status='old')
    PRINT *,'d:\PPB\L3_tropo_ozone_vmr_'//chmd//''//chyr//''
	DO I=1,3;READ (1000, *);ENDDO
	  DO J=1,N
	  DO II=1,11
	  READ(1000,200,err=3) (OZONE(I,J),I=(II-1)*25+1,(II-1)*25+25)
	  ENDDO
	  READ(1000,201,err=3) (OZONE(I,J),I=276,288)
	  ENDDO
3	CLOSE(1000)	
    OPEN(1100,ERR=4,file='d:\TCO\L3_tropo_ozone_column_'//chmd//''//chyr//'',status='old')
    PRINT *,'d:\TCO\L3_tropo_ozone_column_'//chmd//''//chyr//''
	DO I=1,3;READ (1100, *);ENDDO
	  DO J=1,N
	  DO II=1,11
	  READ(1100,200,err=4) (OZONE2(I,J),I=(II-1)*25+1,(II-1)*25+25)
	  ENDDO
	  READ(1100,201,err=4) (OZONE2(I,J),I=276,288)
	  ENDDO
4	CLOSE(1100)	
!   PRINT *,'./PPB/L3_tropo_ozone_vmr_'//chmd//''//chyr//'IN Jlat,IIraw=',J,II
!    do il=1,288
!	do jl=1,120
!	if(ozone(il,jl).ge.600.) then
!	 print *,'.\PPB\L3_tropo_ozone_vmr_'//chmd//''//chyr//'',ozone(il,jl),il,jl
!    endif
!	enddo
!	enddo
	IREC=IREC+1
	WRITE(10,REC=IREC)((OZONE(i,j), i = 1,m),j=1,n)
	IREC=IREC+1
	WRITE(10,REC=IREC)((OZONE2(i,j), i = 1,m),j=1,n)
  enddo
  ENDDO

200	FORMAT(1X,25(F3.0))
201	FORMAT(1X,13(F3.0))
202 FORMAT(288(I3,1X))

	END