gael/bulkMatlab/coeff_MonthlyAtlasInterp.m

function [mcur,mcurW]=coeff_MonthlyAtlasInterp(dcur);
%=> sst=(1-mcurW)*field_sst(:,:,mcur(1))+mcurW*field_sst(:,:,mcur(2));

if dcur<1|dcur>365; 
fprintf('error in coeff_MonthlyAtlasInterp: dcur must be in [1 365] \n'); 
mcur=[]; mcurW=[]; return; end;

tmp1=greg_0h(jul_0h([1994 1 0 0 0 0])+dcur);
tmp2=jul_0h(tmp1)-jul_0h([1994 tmp1(2)+0.5 0 0 0 0]);
if tmp2>=0; tmp3=[tmp2 jul_0h([1994 tmp1(2)+1.5 0 0 0 0])-jul_0h(tmp1)];
mcur=greg_0h(jul_0h([1994 tmp1(2)+1.5 0 0 0 0])); mcur=[tmp1(2) mcur(2)];
else; tmp3=[-tmp2 jul_0h(tmp1)-jul_0h([1994 tmp1(2)-0.5 0 0 0 0])];
mcur=greg_0h(jul_0h([1994 tmp1(2)-0.5 0 0 0 0])); mcur=[tmp1(2) mcur(2)];
end;
mcurW=tmp3(1)/sum(tmp3);

1	function [mcur,mcurW]=coeff_MonthlyAtlasInterp(dcur);
2	%=> sst=(1-mcurW)field_sst(:,:,mcur(1))+mcurWfield_sst(:,:,mcur(2));
3
4	if dcur<1\|dcur>365;
5	fprintf('error in coeff_MonthlyAtlasInterp: dcur must be in [1 365] \n');
6	mcur=[]; mcurW=[]; return; end;
7
8	tmp1=greg_0h(jul_0h([1994 1 0 0 0 0])+dcur);
9	tmp2=jul_0h(tmp1)-jul_0h([1994 tmp1(2)+0.5 0 0 0 0]);
10	if tmp2>=0; tmp3=[tmp2 jul_0h([1994 tmp1(2)+1.5 0 0 0 0])-jul_0h(tmp1)];
11	mcur=greg_0h(jul_0h([1994 tmp1(2)+1.5 0 0 0 0])); mcur=[tmp1(2) mcur(2)];
12	else; tmp3=[-tmp2 jul_0h(tmp1)-jul_0h([1994 tmp1(2)-0.5 0 0 0 0])];
13	mcur=greg_0h(jul_0h([1994 tmp1(2)-0.5 0 0 0 0])); mcur=[tmp1(2) mcur(2)];
14	end;
15	mcurW=tmp3(1)/sum(tmp3);
16