c       Program to determine cold frontal passages over South America 
c234567890123456789001234567890123456789012345678901234567890123456789012
      parameter (ndmax=10958,nyrs=30,numcf=12600)
      dimension ndmon(12),x(ndmax,17)
      data ndmon/31,28,31,30,31,30,31,31,30,31,30,31/
      character path*30
      data path/'/cpc/home/wd52vk/samer/fronts/'/
      open (11,file=path//'SA-cf-dates-swne-section-12mons1980-2009',
     *access='direct',form='unformatted',recl=28)
      open (51,file=path//'SA-fropa-dates-for-swne-section-y1980-2009',
     *form='formatted')
      do 900 nday=1,ndmax
      do 901 k=1,17
      x(nday,k)=0.
 901  continue
 900  continue
      nrec=0
      irec=0
 99   continue
      irec=irec+1 
      if(irec.gt.numcf)go to 999
      read(11,rec=irec)i,j,idate,ndate,sum,v,tot
      lpyr=1980
      nday=0
      do 1 iyr=1,nyrs
      iyear=1979+iyr
      do 2 mon=1,12
      maxd=ndmon(mon)
      if(iyear.eq.lpyr.and.mon.eq.2)maxd=maxd+1
      do 3 iday=1,maxd
      nday=nday+1
      kdate=iyear*10000+mon*100+iday
      if(idate.eq.kdate)go to 4
 3    continue
 2    continue
      if(iyear.eq.lpyr)lpyr=lpyr+4
 1    continue
 4    ii=13
      jj=7
      do 5 k=1,17
      if(i.eq.ii.and.j.eq.jj)x(nday,k)=1.
      ii=ii+1
      jj=jj+1
 5    continue
      go to 99
 999  continue
      lpyr=1980
      nday=0
      do 6 iyr=1,nyrs
      iyear=1979+iyr
      do 7 mon=1,12
      maxd=ndmon(mon)
      if(iyear.eq.lpyr.and.mon.eq.2)maxd=maxd+1
      do 8 iday=1,maxd
      nday=nday+1
      kdate=iyear*10000+mon*100+iday
      write(51,200)kdate,(x(nday,k),k=1,17)
 200  format(i10,17f3.0)
 8    continue
 7    continue
      if(iyear.eq.lpyr)lpyr=lpyr+4
 6    continue
      stop
      end