anchor.diff:                 iarg = iarg + 4
               case('-o')
                  call get_command_argument(iarg+1, outfile)
  -               outfile_set = .true.
                  iarg = iarg + 2
               case('-rmean')
                  rmean = .true.
  
positive.code:                call usage
           end select
        enddo
        if (.not.outfile_set) outfile = trim(infile) // '.nc'
        call get_command(command)
     end subroutine get_args
  
  end program rad3toncf
  !-------------------------------------------------------------------------------
  
negative.code: !===============================================================================
  program rad3voltoncf
  !===============================================================================
  ! Convert a set of rad3 files to NetCDF format.  Supply a file containing lines
  !     file, y_offset, x_offset
  