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