function [field_cur]=gpcp_load_fields(list_years);
%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/forcing/GPCP/';

jpig=144; jpjg=72;
lon2D_g=[1.25:2.5:358.75]'*ones(1,jpjg); lat2D_g=ones(jpig,1)*[-88.75:2.5:88.75];
mask_g=ones(jpig,jpjg); recl=jpig*jpjg*4;

[x1,y1] = meshgrid([lon2D_g(end,1)-360  lon2D_g(:,1)' lon2D_g(1:2,1)'+360],lat2D_g(1,:));
[x2,y2] = meshgrid(lon2D_t(:,1),lat2D_t(1,:));


nb_cur=zeros(1,12); field_cur=zeros(jpi,jpj,12);
for ycur=list_years 
fid=fopen([rep_in 'gpcp_v2_psg.' num2str(ycur)],'r','b');
status=fseek(fid,576,'bof'); %skip header
for mcur=1:12;
field_in=flipdim(fread(fid,[jpig jpjg],'float32'),2);
%interpolate:
tmp1=[field_in(end,:);field_in;field_in(1,:);field_in(2,:)]';
field_out=interp2(x1,y1,tmp1,x2,y2)';
%units:
field_out=field_out.*mask/86400;
%sum:
field_cur(:,:,mcur)=field_cur(:,:,mcur)+field_out; nb_cur(mcur)=nb_cur(mcur)+1;
end;
fclose(fid);
end%for ycur=list_years

%sum2average:
for tcur=1:12; field_cur(:,:,tcur)=field_cur(:,:,tcur)/nb_cur(tcur); end;