integer MAX_R8_BUFLEN
parameter( MAX_R8_BUFLEN=HEADER_SIZE+256*256 )
real*8 r8buf(MAX_R8_BUFLEN)
if (HEADER_SIZE+Ni*Nj.gt.MAX_R8_BUFLEN)
count=HEADER_SIZE+MAX_R8_BUFLEN
count=MAX_R8_BUFLEN