function [ustress_cur,vstress_cur]=quickcow_load_atlas();
%loads the fields and does the interpolation to the 1x1 ECCO grid
domaine_global_def;

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

rep_in='/net/ross/raid0/gforget/DATAbin/QSCATcow/';
list_names=strvcat('january','february','march','april','may','june','july','august','september','october','november','december');

for vcur=1:2;
if vcur==1; eval(['ncload ' rep_in 'wind_stress_zonal_monthly_maps.nc;']); ustress_cur=NaN*zeros(jpi,jpj,12);
else; eval(['ncload ' rep_in 'wind_stress_meridional_monthly_maps.nc;']); vstress_cur=NaN*zeros(jpi,jpj,12);
end;
for mcur=1:12
tmp1=zeros(jpi,jpj); tmp2=zeros(jpi,jpj);
eval(['tmp3=' deblank(list_names(mcur,:)) ''';']);
tmp4=1*(tmp3>-1000); tmp3(tmp3<-1000)=0;
for icur=1:4; for jcur=1:4;
tmp1(:,11:150)=tmp1(:,11:150)+tmp3(icur:4:end,jcur:4:end);
tmp2(:,11:150)=tmp2(:,11:150)+tmp4(icur:4:end,jcur:4:end);
end;end;
%tmp1(find(tmp2>0))=tmp1(find(tmp2>0))./tmp2(find(tmp2>0)); tmp1(find(tmp2==0))=NaN;
%tmp1(find(tmp2>8))=tmp1(find(tmp2>8))./tmp2(find(tmp2>8)); tmp1(find(tmp2<=8))=NaN;
tmp1(find(tmp2>15))=tmp1(find(tmp2>15))./tmp2(find(tmp2>15)); tmp1(find(tmp2<=15))=NaN;

if vcur==1; ustress_cur(:,:,mcur)=tmp1;
else; vstress_cur(:,:,mcur)=tmp1; end;
end%for mcur=1:12
end%for vcur=1:2;