gael/bulkMatlab/sst_load_field.m

function [field_sst]=sst_load_field(rep_sst,suff_sst,pref_sst,ycur,hcur);
%loads the fields and does the interpolation to the 1x1 ECCO grid
domaine_global_def;

mask=squeeze(tmask3D(:,:,1));

%typically: rep_sst='/net/ross/raid2/king/data_1x1_92-03/obs/'; pref_sst='SST_monthly_r2_'; suff_sst='';

%interpolate from the monthly field:
%       coeffs etc: 
dcur=ceil((hcur)/4); [mcur,mcurW]=coeff_MonthlyAtlasInterp(dcur);
ycur2=[ycur ycur];
if ~isempty(find(mcur==1))&dcur>300; ycur2(find(mcur==1))=ycur+1; end;
if ~isempty(find(mcur==12))&dcur<100; ycur2(find(mcur==12))=ycur-1; end;
tmp1=find(ycur2==1991); ycur2(tmp1)=1992; mcur(tmp1)=1;
tmp1=find(ycur2==2007); ycur2(tmp1)=2006; mcur(tmp1)=12;
%       load the two fields:
fid=fopen([rep_sst pref_sst num2str(ycur2(1)) suff_sst],'r','b');
recl=jpi*jpj*4; position0=recl*(mcur(1)-1); status=fseek(fid,position0,'bof');
field1=fread(fid,[jpi jpj],'float32'); fclose(fid);
field1(field1==0)=NaN; field1(field1<-1.9)=NaN; %accomodate various mask cases
fid=fopen([rep_sst pref_sst num2str(ycur2(2)) suff_sst],'r','b');
recl=jpi*jpj*4; position0=recl*(mcur(2)-1); status=fseek(fid,position0,'bof');
field2=fread(fid,[jpi jpj],'float32'); fclose(fid);
field2(field2==0)=NaN; field2(field2<-1.9)=NaN; %accomodate various mask cases
%       weighted average -> time interpolation
field_sst=(1-mcurW)*field1+mcurW*field2;
field_sst=field_sst.*mask;


1	gforget	1.1	function [field_sst]=sst_load_field(rep_sst,suff_sst,pref_sst,ycur,hcur);
2			%loads the fields and does the interpolation to the 1x1 ECCO grid
3			domaine_global_def;
4
5			mask=squeeze(tmask3D(:,:,1));
6
7			%typically: rep_sst='/net/ross/raid2/king/data_1x1_92-03/obs/'; pref_sst='SST_monthly_r2_'; suff_sst='';
8
9			%interpolate from the monthly field:
10			% coeffs etc:
11			dcur=ceil((hcur)/4); [mcur,mcurW]=coeff_MonthlyAtlasInterp(dcur);
12			ycur2=[ycur ycur];
13			if ~isempty(find(mcur==1))&dcur>300; ycur2(find(mcur==1))=ycur+1; end;
14			if ~isempty(find(mcur==12))&dcur<100; ycur2(find(mcur==12))=ycur-1; end;
15			tmp1=find(ycur2==1991); ycur2(tmp1)=1992; mcur(tmp1)=1;
16			tmp1=find(ycur2==2007); ycur2(tmp1)=2006; mcur(tmp1)=12;
17			% load the two fields:
18			fid=fopen([rep_sst pref_sst num2str(ycur2(1)) suff_sst],'r','b');
19			recl=jpijpj4; position0=recl*(mcur(1)-1); status=fseek(fid,position0,'bof');
20			field1=fread(fid,[jpi jpj],'float32'); fclose(fid);
21			field1(field1==0)=NaN; field1(field1<-1.9)=NaN; %accomodate various mask cases
22			fid=fopen([rep_sst pref_sst num2str(ycur2(2)) suff_sst],'r','b');
23			recl=jpijpj4; position0=recl*(mcur(2)-1); status=fseek(fid,position0,'bof');
24			field2=fread(fid,[jpi jpj],'float32'); fclose(fid);
25			field2(field2==0)=NaN; field2(field2<-1.9)=NaN; %accomodate various mask cases
26			% weighted average -> time interpolation
27			field_sst=(1-mcurW)field1+mcurWfield2;
28			field_sst=field_sst.*mask;
29
30