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