!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!    edit parameters im & jm
!    to your models grid size
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      parameter (im=144,jm=89,num_mon=432)
      dimension ice(im,jm)
      dimension ihead(45)
      dimension rhead(19)
      data ihead/45*0/, rhead/19*0.0/

      open (30,file='interp_ice.gra',form="unformatted",
     +      access="direct",recl=im*jm*4)
      open (10,file='interp_ice.pp',form='unformatted')

      ihead(18)=jm
      ihead(19)=im

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!    edit rhead(14), rhead(15), rhead(16) and rhead(17)
!    rhead(14) - starting latitude
!    rhead(15) - latitude distance
!    rhead(16) - starting longitude
!    rhead(17) - longitude distance
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      rhead(14)=-90.00
      rhead(15)=  2.00
      rhead(16)= -1.25
      rhead(17)=  2.50
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

      do n=1,num_mon
!---------
        ly=1971+(n-1)/12
        mo=mod(n-1,12)+1
      ihead(1)=ly
      ihead(2)=mo
      print*,ly,mo
!---------
      read(30,rec=n)ice
!
      write(10)ihead,rhead
      write(10)ice
!---------
      enddo

      end