gael/bulkMatlab/loop_1x1_sst.m


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


1	gforget	1.1
2			domaine_global_def; domaine;
3			mask=squeeze(tmask3D(:,:,1)); recl=jpijpj4;
4
5
6			rep_in='/net/ross/raid0/gforget/DATAbin/forcing/SST/';
7			fid_runoff=fopen([rep_in '../CORE/runoffannual1x1'],'r','b');
8			runoff=fread(fid_runoff,[jpi jpj],'float32').*mask;
9			fclose(fid_runoff);
10
11			%SST and ICE data:
12			fid_sst=fopen('/net/ross/raid0/gforget/DATAbin/ICECONC/HADLEY/HadISST1_SST_9201_monthly','r','b');
13			field_sst=fread(fid_sst,jpijpj12,'float32'); fclose(fid_sst); field_sst=reshape(field_sst,jpi,jpj,12);
14			fid_ice=fopen('/net/ross/raid0/gforget/DATAbin/ICECONC/HADLEY/HadISST1_ICE_9201_monthly','r','b');
15			field_ice=fread(fid_ice,jpijpj12,'float32'); fclose(fid_ice); field_ice=reshape(field_ice,jpi,jpj,12);
16			%for interannual sst:
17			rep_sst='/net/ross/raid2/king/data_1x1_92-03/obs/'; pref_sst='SST_monthly_r2_'; suff_sst='';
18
19			cen2kel=273.15;
20			rhoConstFresh=999.8;
21			stefanBoltzmann = 5.670e-8;
22			ocean_emissivity=5.50e-8 / 5.670e-8;
23			albedo=0.1;
24
25
26			rep_out=rep_in; suff_out=['.9206daily.reynolds.' num2str(nbtcur) '.bin']; list_tcur=[1:365]; list_year=1992+nbtcur-1;
27			doInitFiles=1; doWriteMean=0;
28			rep_out='./';
29
30			choice_plot=0;
31
32			for ycur=list_year
33			for tcur=list_tcur
34			for ttcur=1:4
35
36			hcur=(tcur-1)*4+ttcur;
37
38			%for interannual sst:
39			[sst]=sst_load_field(rep_sst,suff_sst,pref_sst,ycur,hcur);
40
41			%part 3: time average
42			writeFlds=struct('sst',sst);
43			averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);
44			doInitFiles=0;
45
46			end%for ttcur=1:4
47			doWriteMean=1;
48			averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);
49			doWriteMean=0;
50			end%for ycur=1992:2001
51			end%for tcur=1:365
52
53			doInitFiles=-1; doWriteMean=0;
54			averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,writeFlds);
55
56
57