Contents of /MITgcm_contrib/timour_matlab/mscripts/rawio

      PARAMETER (
     &           sNx = 124,
     &           sNy = 124,
     &           OLx =   3,
     &           OLy =   3,
     &           nSx =   1,
     &           nSy =   1,
     &           nPx =   1,
     &           nPy =   1,
     &           Nx  = sNx*nSx*nPx,
     &           Ny  = sNy*nSy*nPy,
     &           Nz  =  30)


      _RS  taux       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

      OPEN(31,FILE='taux',FORM='unformatted')
      WRITE(31) taux
      CLOSE(31)


fid=fopen('taux','r','ieee-be')
fread(fid,1,'int')
phi=fread(fid,130*130,'float32');
th=reshape(phi,[130,130]);


fid=fopen('tauy','r','ieee-be')
fread(fid,1,'int')
phi=fread(fid,130*130,'float32');
th=reshape(phi,[130,130]);
1	PARAMETER (
2	& sNx = 124,
3	& sNy = 124,
4	& OLx = 3,
5	& OLy = 3,
6	& nSx = 1,
7	& nSy = 1,
8	& nPx = 1,
9	& nPy = 1,
10	& Nx = sNxnSxnPx,
11	& Ny = sNynSynPy,
12	& Nz = 30)
13
14
15	_RS taux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
16
17	OPEN(31,FILE='taux',FORM='unformatted')
18	WRITE(31) taux
19	CLOSE(31)
20
21
22
23	fid=fopen('taux','r','ieee-be')
24	fread(fid,1,'int')
25	phi=fread(fid,130*130,'float32');
26	th=reshape(phi,[130,130]);
27
28
29	fid=fopen('tauy','r','ieee-be')
30	fread(fid,1,'int')
31	phi=fread(fid,130*130,'float32');
32	th=reshape(phi,[130,130]);