Annotation 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	edhill	1.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]);