GIStemp STEP2_split_binary

The program split_binary.f

Here is the listing. The explanation will follow some months or weeks from now below, after the === bar.

   
 
      integer,allocatable,dimension(:) :: idata
      integer info(9)
      character title*80,name*36,fileo*7

      open(1,file='Ts.bin',form='unformatted')
      read(1,end=100) info,title

      do n=51,56
      write(fileo,'(a6,i1)') 'Ts.bin',n-50
      open(n,file=fileo,form='unformatted')
      write(n) info,title
      end do
      allocate (idata(info(4)))
   10 read(1,end=100) idata,Lat,Lon,ID,iht,name,m1,m2
!!    write(99,*) Lat,Lon,ID
!!    if(ID.eq.603550002) write(98,'(12I6)') idata
      mu=lat+899
      if(mu.lt.0) mu=0
      mu=56-mu/300
      write(mu) idata,Lat,Lon,ID,iht,name,m1,m2
      go to 10

  100 stop 0
      end


=========================================================

The analysis of this program will have to wait a fair while as I’m still finishing another step. It looks like it’s a well done library of utilities so I’ll likely leave it to the end.

Advertisements

About E.M.Smith

A technical managerial sort interested in things from Stonehenge to computer science. My present "hot buttons' are the mythology of Climate Change and ancient metrology; but things change...
This entry was posted in GISStemp Technical and Source Code and tagged , , , , , . Bookmark the permalink.