/*	ftread reads unformatted data that has been written by fortran
	or in fortran convention -- i.e. an integer with the number of
	bytes in the record, the record data, and another integer with
	the number of bytes in the record.  ftread does various error
	checks, writing a message to stderr and returning a negative 
	value if an error or warning occurs.  The call is identical to
	the standard i/o library routine fread.
*/

#include <stdio.h>

int ftread(ptr,size,nitems,stream)
char *ptr ;
unsigned size, nitems ;
FILE *stream ;

{
    int nbyte1, nbyte2, nitem1 ;
    int errno ;

    errno = 0 ;
    if ( fread(&nbyte1,sizeof(int),1,stream) != 1 )  {
	errno = -10 ;
	fprintf(stderr,"read error, file empty ? \n") ;
	}
    nitem1 = fread(ptr,size,nitems,stream) ;
    if ( nitem1 != nitems ) {
	errno = -20 ;
	fprintf(stderr,"read error, %d items expected, %d items read. \n",
		nitems,nitem1) ;
    }
    if ( fread(&nbyte2,sizeof(int),1,stream) != 1 )  {
	errno = -30 ;
	fprintf(stderr,"read error, file too short ? \n") ;
    }
    if ( nbyte1 != nbyte2 ) {
	errno = errno - 1 ;
	fprintf(stderr,
	"read warning, byte #s do not match, nbyte1 = %d, nbyte2 = %d \n",
	nbyte1,nbyte2) ;
    }
    if ( nbyte1 != size*nitems) {
	errno = errno - 2 ;
	fprintf(stderr,
      "read warning, byte # does not match item #, nbyte1 = %d, nitems = %d \n",
	nbyte1,nitems) ;
    }
    
    return(errno) ;

}