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

c------------------------------ 24-hourly data for 181 days      
      PARAMETER (IDAY=182)  
      PARAMETER (ILON=64)
      PARAMETER (ILAT=36)
      parameter (kgit=20)
      parameter (nvar=3)
c year        phi  uvt
c 1992         23   20(182)   
c 1993-2003    20   20
c 2004         23   23
c 2005         23   23
c 2006         23   23
c 2007         19   19
c 2008         23   23
c 2009         23   23
c 2010         28   28
c 2011         28   28
c 2012         28   28
c
      REAL  aus_ave(ILON,ILAT,kgit,nvar),aus1992(ILON,ILAT,20,nvar),
     $      aus1993(ILON,ILAT,20,nvar),aus1995(ILON,ILAT,20,nvar),
     $      aus1996(ILON,ILAT,20,nvar),aus1998(ILON,ILAT,20,nvar),
     $      aus1999(ILON,ILAT,20,nvar),aus2001(ILON,ILAT,20,nvar),
     $      aus2002(ILON,ILAT,20,nvar),aus2003(ILON,ILAT,20,nvar),
     $      aus2004(ILON,ILAT,23,nvar),aus2006(ILON,ILAT,23,nvar),
     $      aus2008(ILON,ILAT,23,nvar),aus2009(ILON,ILAT,23,nvar),
     $      aus2010(ILON,ILAT,28,nvar),aus2011(ILON,ILAT,28,nvar),
     $      aus2012(ILON,ILAT,28,nvar),aus2000(ILON,ILAT,20,nvar)
C....OUTPUTFILES
          OPEN (200, FILE='UVT_ave.dx',
     $      form='unformatted',
     $      access='direct', status='new',recl=4*ILON*ILAT*kgit*nvar)
c....INPUTFILE      
      OPEN (101,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1992/UVT91_92.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (102,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1993/UVT92_93.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (103,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1995/UVT_94_95.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (104,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1996/UVT_95_96.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (105,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1998/UVT97_98.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (106,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_1999/UVT_98_99.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (107,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2001/UVT_00_01.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (125,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2000/UVT_99_00.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (108,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2002/UVT01_02.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
      OPEN (109,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2003/UVT02_03.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*20*nvar)
c      OPEN (110,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2004/UVT03_04.dx',
c     $      form='unformatted',
c     $      access='direct', status='old',recl=4*ILON*ILAT*23*nvar)
c      OPEN (111,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2006/UVT_05_06.dx',
c     $      form='unformatted',
c     $      access='direct', status='old',recl=4*ILON*ILAT*23*nvar)
c      OPEN (112,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2007/UVT06_07.dx',
c     $      form='unformatted',
c     $      access='direct', status='old',recl=4*ILON*ILAT*19*nvar)
      OPEN (112,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2008/UVT07_08.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*23*nvar)  
      OPEN (113,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2009/UVT08_09.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*23*nvar)
      OPEN (114,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2010/UVT09_10.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*28*nvar)
      OPEN (115,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2011/UVT10_11.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*28*nvar)
      OPEN (116,FILE='c:/MUAM_O3lng_March/HP_ws/UKMO_2012/UVT_11_12.dx',                               UVT_11_12.dx',
     $      form='unformatted',
     $      access='direct', status='old',recl=4*ILON*ILAT*28*nvar)                      
c                
c
      DO  l=1,61
      read(101,rec=l+43) aus1992
      read(102,rec=l+83) aus1993
      read(103,rec=l+60) aus1995
      read(104,rec=l+79) aus1996
      read(105,rec=l+68) aus1998
      read(125,rec=l+68) aus2000
      read(106,rec=l+87) aus1999
      read(107,rec=l+59) aus2001
      read(108,rec=l+78) aus2002
      read(109,rec=l+45) aus2003
c      read(110,rec=l+35) aus2004
c      read(111,rec=l+57) aus2006
      read(112,rec=l+68) aus2008
      read(113,rec=l+53) aus2009
      read(114,rec=l+60) aus2010
      read(115,rec=l+61) aus2011
      read(116,rec=l+48) aus2012
      do mm=1,3
      do  k=1,kgit
        do j=1,ILAT
          do i=1,ILON
      aus_ave(i,j,k,mm)=(aus1992(i,j,k,mm)+aus1993(i,j,k,mm)+
     $                   aus1995(i,j,k,mm)+aus1996(i,j,k,mm)+
     $                   aus1998(i,j,k,mm)+aus1999(i,j,k,mm)+
     $                   aus2001(i,j,k,mm)+aus2002(i,j,k,mm)+
     $                   aus2003(i,j,k,mm)+aus2000(i,j,k,mm)+
     $                   aus2008(i,j,k,mm)+
     $                   aus2009(i,j,k,mm)+aus2010(i,j,k,mm)+
     $                   aus2011(i,j,k,mm)+aus2012(i,j,k,mm))/15.    
          end do
        end do
      end do
      end do
      write(200,rec=l) aus_ave    
      END DO
       STOP
       END