Parent Directory
| 
Revision Log
| 
Revision Graph
Checking in work done with Dimitri on high-resolution cube gridding and parallel
communications.
o code is in a contrib experiment for now so we can continue collaborating
on it. However most code is general and will be moved into main branch once
it is fully hardened.
o There are README files in the contrib root and in the subdirectories that
explain the contents
| 1 | % read3d(filename,nx,ny,nz) returns an array size(nx,ny,nz) read |
| 2 | % from a f77_unformatted the binary file (double precision) named |
| 3 | % by the string filename. |
| 4 | % |
| 5 | % eg. temp=read3d('T.bin',160,10,16); |
| 6 | % |
| 7 | function [arr] = read3d(file,nx,ny,nz) |
| 8 | fid=fopen(file,'r','b'); |
| 9 | q=fread(fid,1,'uint32'); % read the initial 4 bytes used in f77_unformatted |
| 10 | arr=fread(fid,nx*ny*nz,'float64'); |
| 11 | arr=reshape(arr,nx,ny,nz); |
| 12 | st=fclose(fid); |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |