mydata=OpenRead[ "/Users/your_account/test.dat",BinaryFormat->True] t=BinaryRead[mydata,"Real64"]; Print[t] 2.00067 ix=BinaryRead[mydata,"Integer32"]; 202 jx=BinaryRead[mydata,"Integer32"]; 202 x=BinaryReadList[mydata,"Real64", ix]; ListPlot[x] y=BinaryReadList[mydata,"Real64", jx]; ijx=ix*jx 40804 mymx=BinaryReadList[mydata,"Real64",ijx] (* Momentum x *); mymy=BinaryReadList[mydata,"Real64",ijx] (* Momentum y *); mymz=BinaryReadList[mydata,"Real64",ijx] (* Momentum z *); myen=BinaryReadList[mydata,"Real64",ijx] (* Energy *); myro=BinaryReadList[mydata,"Real64",ijx] (* Density *); mybx=BinaryReadList[mydata,"Real64",ijx] (* B x *); myby=BinaryReadList[mydata,"Real64",ijx] (* B y *); mybz=BinaryReadList[mydata,"Real64",ijx] (* B z *); myps=BinaryReadList[mydata,"Real64",ijx] (* Div B control *); myvx=BinaryReadList[mydata,"Real64",ijx] (* Velocity x *); myvy=BinaryReadList[mydata,"Real64",ijx] (* Velocity y *); myvz=BinaryReadList[mydata,"Real64",ijx] (* Velocity z *); mypr=BinaryReadList[mydata,"Real64",ijx] (* Pressure *); Close[mydata] ListContourPlot[ArrayReshape[mypr,{ix,jx}]] ListDensityPlot[ArrayReshape[mypr,{ix,jx}]] ListPlot3D[ArrayReshape[mypr,{ix,jx}]]