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]);