program readata
       parameter (ny=64)
       parameter (ix=81, il=41)
c ny: 1955-2018

c MAM: 31 + 30 + 31 = 92 days
       real dt(ix,il,ny,92)
       real dh(ix,il,ny,92)

       real t(ix,il,ny)
       real c(ix,il),sdv(ix,il)

       real s(ix,il)


       open(11,file='td-PPH_1955-2018_MAM.92d.1x1.gr',
     1 form='unformatted',access='direct',recl=ix*il*64*92*4)

       open(12,file='hl-PPH_1955-2018_MAM.92d.1x1.gr',
     1 form='unformatted',access='direct',recl=ix*il*64*92*4)

       open(15,file='us.land.1x1.gr',
     1 form='unformatted',access='direct',recl=ix*il*4)

c       open(21,file='MAM.clm.sdv.1981-2010.eventday.gr',
c     1 form='unformatted',access='direct',recl=ix*il*4)


        read(11,rec=1) 
     1  ((((dt(i,j,iy,k),i=1,ix),j=1,il),k=1,92),iy=1,ny)
        read(12,rec=1) 
     1  ((((dh(i,j,iy,k),i=1,ix),j=1,il),k=1,92),iy=1,ny)

        read(15,rec=1) ((s(i,j),i=1,ix),j=1,il)

        do i=1,ix
        do j=1,il
        do k=1,92
        do iy=1,ny
         dt(i,j,iy,k)=dt(i,j,iy,k)*100.0
         dh(i,j,iy,k)=dh(i,j,iy,k)*100.0
        end do
        end do
        end do
        end do


        do iy=1,ny

         nn=0

        do k=1,92
        
        do i=1,ix
        do j=1,il
         if(dt(i,j,iy,k).ge.10.or.dh(i,j,iy,k).ge.30) then
          nn=nn+1
          goto 100
         end if
        end do
        end do
100      continue
        end do

        print *, iy+1954, nn
        end do

        
        stop
        end