b2xyz.f

Go to the documentation of this file.

c> \ingroup core 
c> \if MANPAGE     
c> \page b2xyz
c> \endif      
c> 
c> Part of the NADCON5 \ref core , converts `*.b` grid to `xyz`
c>     
c> Program to convert standard "*.b" grid
c> formatted data to a binary xyz (lon, lat, value)
c> list, which can then be used by GMT for
c> various things (like running the GMT
c> routine "xyz2grd", to get a "*.grd" 
c> file, which is useful for plotting, etc)
c>     
c> ### Program arguments
c> Arguments are newline terminated and read from standard input
c>     
c> When run from the command line, the program prints a prompt for each argument    
c>     
c> They are enumerated here
c> \param infile Input File Name
c>      
c> ### Program Inputs:
c> 
c> -  Input File defined by `infile`
c>
c> ### Program Outputs:
c> 
c> - `temp.xyz`
c>
      program b2xyz
      implicit real*8(a-h,o-z)
      character*200 fnamein,fnameout,prefix
      real*4 data(100000)
      integer*4 idata(100000)
      real*4 xlo,xla,val
      equivalence(data(1),idata(1))

c ------------------------------------------------
c - User-supplied file name
c ------------------------------------------------
      read(5,'(a)')fnamein
      open(1,file=fnamein,status='old',form='unformatted')
      read(1)glamn,glomn,dla,dlo,nla,nlo,ikind
c     ll = len(trim(fnamein))
c     prefix = fnamein(1:ll-2)
c     fnameout = prefix//'.xyz'
      fnameout = 'temp.xyz'
      open(2,file=fnameout,status='new',form='binary')

c - Read, compute location, and write
      do 1 i=1,nla
        if(ikind.eq.1)read(1)( data(j),j=1,nlo)
        if(ikind.eq.0)read(1)(idata(j),j=1,nlo)
        xla = glamn + (i-1)*dla
        do 2 j=1,nlo
          xlo = glomn + (j-1)*dlo
          if(ikind.eq.1)val = data(j)
          if(ikind.eq.0)val = dble(idata(j))
          write(2)xlo,xla,val
    2   continue
    1 continue
      end