c  written by Vern Kousky 10/29/13
c  Create a dataset for cloud brightness temps. 1.75 deg boxes 1980-2009
c23456789012345678901234567890123456789012345678901234567890123456789012
      parameter (imax=40,jmax=40,ndmax=365)
      dimension x(imax,jmax)
      character path*33
      data path/'/cpc/home/wd52vk/samer/satellite/'/
      open (11,file=path//'irwin-SA-daily-anoms-boxes-y1980' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (12,file=path//'irwin-SA-daily-anoms-boxes-y1981' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (13,file=path//'irwin-SA-daily-anoms-boxes-y1982' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (14,file=path//'irwin-SA-daily-anoms-boxes-y1983' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (15,file=path//'irwin-SA-daily-anoms-boxes-y1984' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (16,file=path//'irwin-SA-daily-anoms-boxes-y1985' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (17,file=path//'irwin-SA-daily-anoms-boxes-y1986' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (18,file=path//'irwin-SA-daily-anoms-boxes-y1987' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (19,file=path//'irwin-SA-daily-anoms-boxes-y1988' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (20,file=path//'irwin-SA-daily-anoms-boxes-y1989' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (21,file=path//'irwin-SA-daily-anoms-boxes-y1990' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (22,file=path//'irwin-SA-daily-anoms-boxes-y1991' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (23,file=path//'irwin-SA-daily-anoms-boxes-y1992' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (24,file=path//'irwin-SA-daily-anoms-boxes-y1993' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (25,file=path//'irwin-SA-daily-anoms-boxes-y1994' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (26,file=path//'irwin-SA-daily-anoms-boxes-y1995' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (27,file=path//'irwin-SA-daily-anoms-boxes-y1996' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (28,file=path//'irwin-SA-daily-anoms-boxes-y1997' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (29,file=path//'irwin-SA-daily-anoms-boxes-y1998' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (30,file=path//'irwin-SA-daily-anoms-boxes-y1999' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (31,file=path//'irwin-SA-daily-anoms-boxes-y2000' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (32,file=path//'irwin-SA-daily-anoms-boxes-y2001' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (33,file=path//'irwin-SA-daily-anoms-boxes-y2002' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (34,file=path//'irwin-SA-daily-anoms-boxes-y2003' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (35,file=path//'irwin-SA-daily-anoms-boxes-y2004' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (36,file=path//'irwin-SA-daily-anoms-boxes-y2005' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (37,file=path//'irwin-SA-daily-anoms-boxes-y2006' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (38,file=path//'irwin-SA-daily-anoms-boxes-y2007' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (39,file=path//'irwin-SA-daily-anoms-boxes-y2008' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (40,file=path//'irwin-SA-daily-anoms-boxes-y2009' 
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      open (51,file=path//'irwin-SA-daily-anoms-boxes-y1980-2009'
     *,access='direct',form='unformatted',recl=4*imax*jmax)
      lpyr=1980
      ifile=10
      irec=0
      do 99 iyr=1,30
      ifile=ifile+1
      iyear=1979+iyr
      nmax=ndmax
      if(iyear.eq.lpyr)nmax=nmax+1
      do 1 n=1,nmax
      read(ifile,rec=n)x
      irec=irec+1
      write(51,rec=irec)x
 1    continue
      write(6,200)iyear,irec
 200  format(2i10)
      if(iyear.eq.lpyr)lpyr=lpyr+4
 99   continue
      stop
      end