Program interpolation
       character*28 pname    
       dimension tempa0(144,73),tempa1(144,72)
       dimension temp0(144,73),temp1(144,72)
       dimension tempa(145,73),temp(145,73)
       pname='/tc2/yangbo/lukehenew_Dec06/'
       open(333,file='caforecast.dat',form='unformatted')
       open(334,file='caforecastn.dat',form='unformatted')

       do lead=1,6
       read(333) ((tempa0(i,j),i=1,144),j=1,73)
       read(333) ((temp0(i,j),i=1,144),j=1,73)

           do i=1,144
           do j=1,73
            tempa(i,74-j)=tempa0(i,j)
            temp(i,74-j)=temp0(i,j)
           enddo
           enddo


         do j=1,73
         tempa(145,j)=tempa(1,j)
         temp(145,j)=temp(1,j)
         enddo
        
         do i=1,144
         do j=1,72
         tempa1(i,j)=(tempa(i,j)+tempa(i,j+1)+tempa(i+1,j)
     +  +tempa(i+1,j+1))/4.
         temp1(i,j)=(temp(i,j)+temp(i,j+1)+temp(i+1,j)
     +  +temp(i+1,j+1))/4.
         enddo
         enddo
 
      
       write(334) ((tempa1(i,j),i=1,144),j=1,72)
       write(334) ((temp1(i,j),i=1,144),j=1,72)
         
        enddo
         stop
         end