king/raw_files/gridread.m

fid=fopen('tpg_ref.av11-341f2','r');

junk = fread(fid,1,'long');
iz = fread(fid,1,'long')
jz = fread(fid,1,'long')
ylat = fread(fid,2,'float32')
xlon = fread(fid,2,'float32')
junk = fread(fid,2,'long');
izjz = iz*jz;

ssh = fread(fid,izjz,'float32'); ssh = reshape(ssh,iz,jz); ssh = ssh';
ssh(33,:)

missing=-9999.;
[xx,yy]=find(ssh==missing);
nn=length(xx);
for ii=1:nn
  ssh(xx(ii),yy(ii))=NaN;
end
ssh(:,181)=ssh(:,1); xlon(2)=xlon(2)+2;
yl=ylat(1):2:ylat(2); xl=xlon(1):2:xlon(2);
scon=-160:20:180; contour(xl,yl,ssh,scon)
1	edhill	1.1	fid=fopen('tpg_ref.av11-341f2','r');
2
3			junk = fread(fid,1,'long');
4			iz = fread(fid,1,'long')
5			jz = fread(fid,1,'long')
6			ylat = fread(fid,2,'float32')
7			xlon = fread(fid,2,'float32')
8			junk = fread(fid,2,'long');
9			izjz = iz*jz;
10
11			ssh = fread(fid,izjz,'float32'); ssh = reshape(ssh,iz,jz); ssh = ssh';
12			ssh(33,:)
13
14			missing=-9999.;
15			[xx,yy]=find(ssh==missing);
16			nn=length(xx);
17			for ii=1:nn
18			ssh(xx(ii),yy(ii))=NaN;
19			end
20			ssh(:,181)=ssh(:,1); xlon(2)=xlon(2)+2;
21			yl=ylat(1):2:ylat(2); xl=xlon(1):2:xlon(2);
22			scon=-160:20:180; contour(xl,yl,ssh,scon)