matlab_class/ecco_v4/read_grace_fld.m

function [fld,fldsm]=read_grace_fld(fname);

%[fld]=read_grace_fld('GRACE_CSR_Error.asc');

fid=fopen(fname,'rt');
fgetl(fid);
fgetl(fid);

fld4columns=zeros(1e5,4);
ii=1;
while ~feof(fid);
fld4columns(ii,:)=str2num(fgetl(fid));
ii=ii+1;
end;
fclose(fid);

ii=find(fld4columns(:,2)>0&fld4columns(:,2)<360); 
jj=(fld4columns(ii,1)+89.5)*360+fld4columns(ii,2)+0.5;
fld=NaN*zeros(360,180); 
fld(jj)=fld4columns(ii,3);
fldsm=NaN*zeros(360,180);
fldsm(jj)=fld4columns(ii,4);


1	gforget	1.1	function [fld,fldsm]=read_grace_fld(fname);
2
3			%[fld]=read_grace_fld('GRACE_CSR_Error.asc');
4
5			fid=fopen(fname,'rt');
6			fgetl(fid);
7			fgetl(fid);
8
9			fld4columns=zeros(1e5,4);
10			ii=1;
11			while ~feof(fid);
12			fld4columns(ii,:)=str2num(fgetl(fid));
13			ii=ii+1;
14			end;
15			fclose(fid);
16
17			ii=find(fld4columns(:,2)>0&fld4columns(:,2)<360);
18			jj=(fld4columns(ii,1)+89.5)*360+fld4columns(ii,2)+0.5;
19			fld=NaN*zeros(360,180);
20			fld(jj)=fld4columns(ii,3);
21			fldsm=NaN*zeros(360,180);
22			fldsm(jj)=fld4columns(ii,4);
23
24