domaine_global_def; domaine;
mask=squeeze(tmask3D(:,:,1)); recl=jpi*jpj*4;


rep_in='/net/ross/raid0/gforget/DATAbin/forcing/SST/';
fid_runoff=fopen([rep_in '../CORE/runoffannual1x1'],'r','b');
runoff=fread(fid_runoff,[jpi jpj],'float32').*mask;
fclose(fid_runoff);

%SST and ICE data:
fid_sst=fopen('/net/ross/raid0/gforget/DATAbin/ICECONC/HADLEY/HadISST1_SST_9201_monthly','r','b');
field_sst=fread(fid_sst,jpi*jpj*12,'float32'); fclose(fid_sst); field_sst=reshape(field_sst,jpi,jpj,12);
fid_ice=fopen('/net/ross/raid0/gforget/DATAbin/ICECONC/HADLEY/HadISST1_ICE_9201_monthly','r','b');
field_ice=fread(fid_ice,jpi*jpj*12,'float32'); fclose(fid_ice); field_ice=reshape(field_ice,jpi,jpj,12);
%for interannual sst:
rep_sst='/net/ross/raid2/king/data_1x1_92-03/obs/'; pref_sst='SST_monthly_r2_'; suff_sst='';

cen2kel=273.15;
rhoConstFresh=999.8;
stefanBoltzmann = 5.670e-8;
ocean_emissivity=5.50e-8 / 5.670e-8;
albedo=0.1;


rep_out=rep_in; suff_out=['.9206daily.reynolds.' num2str(nbtcur) '.bin']; list_tcur=[1:365]; list_year=1992+nbtcur-1;
doInitFiles=1; doWriteMean=0;
rep_out='./'; 

choice_plot=0; 

for ycur=list_year
for tcur=list_tcur
for ttcur=1:4

hcur=(tcur-1)*4+ttcur;

%for interannual sst:
[sst]=sst_load_field(rep_sst,suff_sst,pref_sst,ycur,hcur);

%part 3: time average
writeFlds=struct('sst',sst);
averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);
doInitFiles=0; 

end%for ttcur=1:4
doWriteMean=1;
averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);
doWriteMean=0;
end%for ycur=1992:2001
end%for tcur=1:365

doInitFiles=-1; doWriteMean=0; 
averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);