imaxrecord = INT(tmp(1))
print*,'read ',dataFname,imaxrecord
itotalrecord = itotalrecord + imaxrecord
do irec=2,imaxrecord+1