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