C==============================================================
C
c    calculation of longitudinal harmonics    
c
c    in NCAR/NCEP Geopotential height data
C
c==============================================================

c------------------------------ 3-hourly data for Nov-Apr      
      PARAMETER (IDAY=1448)  
      PARAMETER (ILON=64)
      PARAMETER (ILAT=36)
      parameter (kgit=23)
      parameter (nvar=1)

      REAL  av_wQBO(ILON,ILAT,kgit,nvar),av_eQBO(ILON,ILAT,kgit,nvar),
     $      aus1983(ILON,ILAT,kgit,nvar),aus1985(ILON,ILAT,kgit,nvar),
     $      aus1993(ILON,ILAT,kgit,nvar),aus1995(ILON,ILAT,kgit,nvar),
     $      aus1997(ILON,ILAT,kgit,nvar),aus1999(ILON,ILAT,kgit,nvar),
     $      aus2002(ILON,ILAT,kgit,nvar),aus2004(ILON,ILAT,kgit,nvar),
     $      aus2008(ILON,ILAT,kgit,nvar),aus2013(ILON,ILAT,kgit,nvar)
      REAL  ave_QBO(ILON,ILAT,kgit,nvar),
     $      aus1987(ILON,ILAT,kgit,nvar),aus1989(ILON,ILAT,kgit,nvar),
     $      aus1996(ILON,ILAT,kgit,nvar),aus1998(ILON,ILAT,kgit,nvar),
     $      aus2000(ILON,ILAT,kgit,nvar),aus2003(ILON,ILAT,kgit,nvar),
     $      aus2005(ILON,ILAT,kgit,nvar),aus2007(ILON,ILAT,kgit,nvar),
     $      aus2010(ILON,ILAT,kgit,nvar),aus2012(ILON,ILAT,kgit,nvar)
C....OUTPUTFILES
          OPEN (300, FILE='PHI_ave_QBO.dx',
     $      form='unformatted',
     $      access='direct', status='new',recl=4*ILON*ILAT*kgit*nvar)
          OPEN (301, FILE='PHI_av_wQBO.dx',
     $      form='unformatted',
     $      access='direct', status='new',recl=4*ILON*ILAT*kgit*nvar)
          OPEN (302, FILE='PHI_av_eQBO.dx',
     $      form='unformatted',
     $      access='direct', status='new',recl=4*ILON*ILAT*kgit*nvar)
c....INPUTFILE_____wQBO      
      OPEN (101,
     $FILE='e:/Merra2/wQBO/1983/merra2_H_1982_1983_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (102,
     $FILE='e:/Merra2/wQBO/1985/merra2_H_1984_1985_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (103,
     $FILE='e:/Merra2/wQBO/1993/merra2_H_1992_1993_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (104,
     $FILE='e:/Merra2/wQBO/1995/merra2_H_1994_1995_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (105,
     $FILE='e:/Merra2/wQBO/1997/merra2_H_1996_1997_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (106,
     $FILE='e:/Merra2/wQBO/1999/merra2_H_1998_1999_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (107,
     $FILE='e:/Merra2/wQBO/2002/merra2_H_2001_2002_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (108,
     $FILE='e:/Merra2/wQBO/2004/merra2_H_2003_2004_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (109,
     $FILE='e:/Merra2/wQBO/2008/merra2_H_2007_2008_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (110,
     $FILE='e:/Merra2/wQBO/2013/merra2_H_2012_2013_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
c....INPUTFILE_____eQBO 
      OPEN (201,
     $FILE='e:/Merra2/eQBO/1987/merra2_H_1986_1987_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (202,
     $FILE='e:/Merra2/eQBO/1989/merra2_H_1988_1989_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (203,
     $FILE='e:/Merra2/eQBO/1996/merra2_H_1995_1996_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (204,
     $FILE='e:/Merra2/eQBO/1998/merra2_H_1997_1998_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (205,
     $FILE='e:/Merra2/eQBO/2000/merra2_H_1999_2000_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (206,
     $FILE='e:/Merra2/eQBO/2003/merra2_H_2002_2003_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (207,
     $FILE='e:/Merra2/eQBO/2005/merra2_H_2004_2005_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (208,
     $FILE='e:/Merra2/eQBO/2007/merra2_H_2006_2007_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
      OPEN (209,
     $FILE='e:/Merra2/eQBO/2010/merra2_H_2009_2010_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
             OPEN (210,
     $FILE='e:/Merra2/eQBO/2012/merra2_H_2011_2012_x64y36z23.bin',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*kgit*nvar)
                                                         
c                
      DO  l=1,IDAY
      read(101,rec=l) aus1983
      read(102,rec=l) aus1985
      read(103,rec=l) aus1993
      read(104,rec=l) aus1995
      read(105,rec=l) aus1997
      read(106,rec=l) aus1999
      read(107,rec=l) aus2002
      read(108,rec=l) aus2004
      read(109,rec=l) aus2008
      read(110,rec=l) aus2013
c
      read(201,rec=l) aus1987
      read(202,rec=l) aus1989
      read(203,rec=l) aus1996
      read(204,rec=l) aus1998
      read(205,rec=l) aus2000
      read(206,rec=l) aus2003
      read(207,rec=l) aus2005
      read(208,rec=l) aus2007
      read(209,rec=l) aus2010
      read(210,rec=l) aus2012      
c
      do mm=1,nvar
      do  k=1,kgit
        do j=1,ILAT
          do i=1,ILON
      av_wQBO(i,j,k,mm)=(aus1983(i,j,k,mm)+aus1985(i,j,k,mm)+
     $                   aus1993(i,j,k,mm)+aus1995(i,j,k,mm)+
     $                   aus1997(i,j,k,mm)+aus1999(i,j,k,mm)+
     $                   aus2002(i,j,k,mm)+aus2004(i,j,k,mm)+
     $                   aus2008(i,j,k,mm)+aus2013(i,j,k,mm))/10.
      av_eQBO(i,j,k,mm)=(aus1987(i,j,k,mm)+aus1989(i,j,k,mm)+
     $                   aus1996(i,j,k,mm)+aus1998(i,j,k,mm)+
     $                   aus2000(i,j,k,mm)+aus2003(i,j,k,mm)+
     $                   aus2005(i,j,k,mm)+aus2007(i,j,k,mm)+
     $                   aus2010(i,j,k,mm)+aus2012(i,j,k,mm))/10.
      ave_QBO(i,j,k,mm)=(av_wQBO(i,j,k,mm)+
     $                   av_eQBO(i,j,k,mm))/2.
          end do
        end do
      end do
      end do
      write(300,rec=l) ave_QBO 
      write(301,rec=l) av_wQBO   
      write(302,rec=l) av_eQBO      
      END DO
       STOP
       END