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