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

c------------------------------ Monthly data       
      PARAMETER (MoNTHS=744)  
      PARAMETER (ILON=64)
      PARAMETER (ILAT=36)
      REAL gh_La_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....INPUTFILE mean Jan 10 years 2010-2019
c    2010 (625)
c    2011 (637)
c    2012 (649)
c    2013 (661)
c    2014 (673)
c    2015 (685)
c    2016 (697)
c    2017 (709)
c    2018 (721)
c    2019 (733)
c...................................................
      OPEN (101,
     $FILE='jra_hgt_1958m01-12_2019m01-12_x64y36z1t744.bin',
     $      form='unformatted',
     $      access='direct', status='unknown',recl=4*ILON*ILAT*MONTHS)
c########################################################################
C....OUTPUTFILES
      OPEN (301, FILE='hgt_mean10_dec.dx',
     $      form='unformatted',
     $      access='direct', status='unknown',recl=4*ILON*ILAT)
c-----------------------------------------------------------------
      do j=1,ilat
      do i=1,ilon 
      gh_La_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,636)+aus(i,j,648)+aus(i,j,660)+
     $               aus(i,j,672)+aus(i,j,684)+aus(i,j,696)+
     $               aus(i,j,708)+aus(i,j,720)+aus(i,j,732)+
     $               aus(i,j,744))/10.

          end do
        end do
      write(301,rec=1)gh_La_Jan

       STOP
       END