C==============================================================
C
c    Lower boundary for La and El ENSO
C
c==============================================================

c------------------------------ Monthly data       
      PARAMETER (MoNTHS=674)  
      PARAMETER (ILON=64)
      PARAMETER (ILAT=36)
      REAL gh_La_Jan(ILON,ILAT),gh_El_Jan(ILON,ILAT)
c--------------------- input field     
      REAL  aus(ILON,ILAT,MONTHS)
c######################################
c      REAL  aus(ILON,ILAT,17,MONTHS)
c....INPUTFILE mean Jan 10 years
c    2004 (553)
c    2005 (565)
c    2006 (577)
c    2007 (589)
c    2008 (601)
c    2009 (613)
c    2010 (625)
c    2011 (637)
c    2012 (649)
c    2013 (661)
c    2014 (673)
c....INPUTFILE La
c    1962 ( 49)
c    1971 (157)
c    1974 (193)
c    1976 (217)
c    1989 (373)
c    1999 (493)
c    2000 (505)
c    2008 (601)
c    2011 (637)
c----------------
c....INPUTFILE El
c    1958 (  1)
c    1966 ( 97)
c    1973 (181)
c    1983 (301)
c    1987 (349)
c    1992 (409)
c    1993 (421)
c    1998 (481)
c    2010 (625)
c----------------
      OPEN (101,
     $FILE='air.mon.mean.nc64x36x01.bin',
     $      form='unformatted',
     $      access='direct', status='unknown',recl=4*ILON*ILAT*MONTHS)
c########################################################################
c      OPEN (101,
c     $FILE='air.mon.mean.nc64x36x01.bin',
c     $      form='unformatted',
c     $      access='direct', status='unknown',recl=4*ILON*ILAT*17*MONTHS)
C....OUTPUTFILES
      OPEN (301, FILE='t_mean10_mar_04_13.dx',
     $      form='unformatted',
     $      access='direct', status='unknown',recl=4*ILON*ILAT)
c      OPEN (302, FILE='GH_1000_El_Feb.dx',
c     $      form='unformatted',
c     $      access='direct', status='unknown',recl=4*ILON*ILAT)
c-----------------------------------------------------------------
      do j=1,ilat
      do i=1,ilon 
      gh_La_Jan(i,j)=0.
*      gh_El_Jan(i,j)=0.
      end do
      end do
      read(101,rec=1) aus
        do j=1,ILAT
          do i=1,ILON
      gh_La_Jan(i,j)=(aus(i,j,555)+aus(i,j,567)+aus(i,j,579)+
     $               aus(i,j,591)+aus(i,j,603)+aus(i,j,615)+
     $               aus(i,j,627)+aus(i,j,639)+aus(i,j,651)+
     $               aus(i,j,663))/10.
c      gh_La_Jan(i,j)=(aus(i,j,373)+aus(i,j,493)+
c     $               aus(i,j,505)+aus(i,j,601)+aus(i,j,637))/5.
c      gh_El_Jan(i,j)=(aus(i,j,  1)+aus(i,j, 97)+aus(i,j,181)+
c     $               aus(i,j,301)+aus(i,j,349)+aus(i,j,409)+
c     $               aus(i,j,421)+aus(i,j,481)+aus(i,j,625))/9.
          end do
        end do
      write(301,rec=1)gh_La_Jan
c      write(302,rec=1)gh_El_Jan
       STOP
       END