1 |
gforget |
1.1 |
|
2 |
|
|
%select kBudget: |
3 |
|
|
if ~isempty(setDiagsParams); |
4 |
|
|
kBudget=setDiagsParams{1}; |
5 |
|
|
else; |
6 |
|
|
kBudget=1; |
7 |
|
|
end; |
8 |
|
|
|
9 |
gforget |
1.9 |
doMoreBudgetOutput=0; |
10 |
gforget |
1.16 |
%doMoreBudgetOutput=1; |
11 |
gforget |
1.9 |
|
12 |
gforget |
1.2 |
%override default file name: |
13 |
|
|
%--------------------------- |
14 |
|
|
tmp1=setDiags; |
15 |
|
|
if kBudget>1; |
16 |
|
|
tmp1=sprintf('D%02i',kBudget); |
17 |
gforget |
1.1 |
end; |
18 |
gforget |
1.2 |
fileMat=['diags_set_' tmp1]; |
19 |
gforget |
1.1 |
|
20 |
|
|
if userStep==1;%diags to be computed |
21 |
|
|
listDiags=['glo_vol_ocn glo_vol_tot glo_vol_ice glo_bp']; |
22 |
|
|
listDiags=[listDiags ' north_vol_ocn north_vol_tot north_vol_ice north_bp']; |
23 |
|
|
listDiags=[listDiags ' south_vol_ocn south_vol_tot south_vol_ice south_bp']; |
24 |
|
|
listDiags=[listDiags ' glo_heat_ocn glo_heat_tot glo_heat_ice']; |
25 |
|
|
listDiags=[listDiags ' north_heat_ocn north_heat_tot north_heat_ice']; |
26 |
|
|
listDiags=[listDiags ' south_heat_ocn south_heat_tot south_heat_ice']; |
27 |
|
|
listDiags=[listDiags ' glo_salt_ocn glo_salt_tot glo_salt_ice']; |
28 |
|
|
listDiags=[listDiags ' north_salt_ocn north_salt_tot north_salt_ice']; |
29 |
|
|
listDiags=[listDiags ' south_salt_ocn south_salt_tot south_salt_ice']; |
30 |
gforget |
1.2 |
|
31 |
gforget |
1.1 |
elseif userStep==2;%input files and variables |
32 |
gforget |
1.15 |
dirSnap=fullfile(dirModel,'diags',filesep,'BUDG',filesep); |
33 |
|
|
if ~isdir(dirSnap); dirSnap=fullfile(dirModel,'diags',filesep); end; |
34 |
|
|
tmp1=fullfile(dirSnap,'budg2d_snap_set2*meta'); |
35 |
gforget |
1.11 |
test3d=isempty(dir(tmp1)); |
36 |
|
|
% |
37 |
gforget |
1.16 |
listFlds={ 'ETAN','SIheff','SIhsnow','THETA ','SALT ','PHIBOT','geothFlux'}; |
38 |
gforget |
1.1 |
listFlds={listFlds{:},'SIatmFW ','oceFWflx','SItflux','TFLUX','SFLUX','oceSPflx','SRELAX'}; |
39 |
|
|
listFlds={listFlds{:},'oceQnet ','SIatmQnt','SIaaflux','SIsnPrcp','SIacSubl'}; |
40 |
|
|
listFlds={listFlds{:},'TRELAX','WTHMASS','WSLTMASS','oceSflux','oceQsw','oceSPtnd'}; |
41 |
gforget |
1.11 |
if kBudget>1|test3d; |
42 |
gforget |
1.1 |
listFlds={listFlds{:},'ADVr_TH','DFrE_TH','DFrI_TH','ADVr_SLT','DFrE_SLT','DFrI_SLT','WVELMASS'}; |
43 |
|
|
end; |
44 |
gforget |
1.4 |
listFlds={listFlds{:},'SDIAG1','SDIAG2','SDIAG3'}; |
45 |
gforget |
1.1 |
listFlds={listFlds{:},'UVELMASS','VVELMASS','AB_gT','AB_gS'}; |
46 |
|
|
listFlds={listFlds{:},'ADVx_TH ','ADVy_TH ','DFxE_TH ','DFyE_TH '}; |
47 |
|
|
listFlds={listFlds{:},'ADVx_SLT','ADVy_SLT','DFxE_SLT','DFyE_SLT'}; |
48 |
|
|
listFlds={listFlds{:},'ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF'}; |
49 |
|
|
listFlds={listFlds{:},'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW'}; |
50 |
|
|
listFldsNames=deblank(listFlds); |
51 |
|
|
% |
52 |
|
|
listFiles={'rate_budg2d_snap_set1','budg2d_hflux_set1','budg2d_zflux_set1','budg2d_zflux_set2'}; |
53 |
gforget |
1.11 |
if test3d; |
54 |
|
|
listFiles={listFiles{:},'rate_budg3d_snap_set1','budg3d_hflux_set1','budg3d_zflux_set1'}; |
55 |
|
|
elseif kBudget==1; |
56 |
gforget |
1.16 |
listFiles={listFiles{:},'rate_budg2d_snap_set2','budg2d_hflux_set2','geothermalFlux'}; |
57 |
gforget |
1.1 |
else; |
58 |
|
|
tmp1=sprintf('rate_budg2d_snap_set3_%02i',kBudget); |
59 |
|
|
tmp2=sprintf('budg2d_zflux_set3_%02i',kBudget); |
60 |
|
|
tmp3=sprintf('budg2d_hflux_set3_%02i',kBudget); |
61 |
gforget |
1.16 |
tmp4=sprintf('geothermalFlux_%02i',kBudget); |
62 |
|
|
listFiles={listFiles{:},tmp1,tmp2,tmp3,tmp4}; |
63 |
gforget |
1.1 |
end; |
64 |
gforget |
1.15 |
listSubdirs={[dirMat 'BUDG/' ],[dirMat '../BUDG/' ],dirSnap}; |
65 |
gforget |
1.2 |
|
66 |
gforget |
1.1 |
elseif userStep==3;%computational part; |
67 |
gforget |
1.2 |
|
68 |
|
|
%preliminary tests |
69 |
|
|
test1=isempty(dir([dirModel 'diags/BUDG/budg2d_snap_set1*'])); |
70 |
gforget |
1.5 |
test2=isempty(dir([dirMat 'BUDG/rate_budg2d_snap_set1*']))&... |
71 |
|
|
isempty(dir([dirMat '../BUDG/rate_budg2d_snap_set1*'])); |
72 |
gforget |
1.2 |
if (strcmp(setDiags,'D')&test1&test2); |
73 |
|
|
fprintf('\n abort : global and regional budgets, due to missing \n'); |
74 |
|
|
fprintf(['\n ' dirModel 'diags/BUDG/budg2d_snap_set1* \n']); |
75 |
|
|
return; |
76 |
|
|
end; |
77 |
|
|
|
78 |
|
|
if (strcmp(setDiags,'D')&test2); |
79 |
|
|
fprintf('\n abort : global and regional budgets, due to missing \n'); |
80 |
|
|
fprintf(['\n ' dirModel 'diags/BUDG/rate_budg2d_snap_set1* \n']); |
81 |
|
|
return; |
82 |
|
|
end; |
83 |
gforget |
1.1 |
|
84 |
|
|
%override default file name: |
85 |
|
|
%--------------------------- |
86 |
|
|
tmp1=setDiags; |
87 |
|
|
if kBudget>1; |
88 |
|
|
tmp1=sprintf('D%02i',kBudget); |
89 |
|
|
end; |
90 |
|
|
fileMat=['diags_set_' tmp1 '_' num2str(tt) '.mat']; |
91 |
|
|
|
92 |
|
|
%fill in optional fields: |
93 |
|
|
%------------------------ |
94 |
gforget |
1.11 |
if isempty(who('TRELAX')); TRELAX=0; end; |
95 |
|
|
if isempty(who('SRELAX')); SRELAX=0; end; |
96 |
|
|
if isempty(who('AB_gT')); AB_gT=0; end; |
97 |
|
|
if isempty(who('AB_gS')); AB_gS=0; end; |
98 |
|
|
if isempty(who('oceSPtnd')); oceSPtnd=0; end; |
99 |
|
|
if isempty(who('oceSPflx')); oceSPflx=0; end; |
100 |
|
|
if isempty(who('PHIBOT')); PHIBOT=0; end; |
101 |
gforget |
1.16 |
if isempty(who('geothFlux')); geothFlux=0; end; |
102 |
gforget |
1.4 |
|
103 |
|
|
%aliases from development phase (applies to 2012 core runs) |
104 |
|
|
%--------------------------------------------------------- |
105 |
|
|
if ~isempty(who('SDIAG1')); SRELAX=SDIAG1; end; |
106 |
|
|
if ~isempty(who('SDIAG2')); SIatmFW=SDIAG2; end; |
107 |
|
|
if ~isempty(who('SDIAG3')); SItflux=SDIAG3; end; |
108 |
|
|
|
109 |
gforget |
1.17 |
%=======indexing and sign convention====== |
110 |
|
|
|
111 |
|
|
%- MITgcm: fluxes are >0 downward, k=1 start at free surface |
112 |
|
|
%- here, similarly: >0 downward, k=1 free surface k=2 sea floor |
113 |
|
|
if ~test3d; |
114 |
|
|
budgO.specs.top='free surface'; |
115 |
|
|
if kBudget>1; budgO.specs.top=['interface no. ' num2str(kBudget)]; end; |
116 |
|
|
budgO.specs.bottom='sea floor'; |
117 |
|
|
else; |
118 |
|
|
budgO.specs.top='interface k'; |
119 |
|
|
budgO.specs.bottom='interface k+1'; |
120 |
|
|
end; |
121 |
|
|
budgI.specs.top='ocn-ice to atm interface'; |
122 |
|
|
budgI.specs.bottom='free surface'; |
123 |
gforget |
1.1 |
|
124 |
gforget |
1.17 |
%here we output tendencies and fluxes in kg/s |
125 |
|
|
budgMo=budgO; budgMi=budgI; |
126 |
|
|
budgMo.specs.units='kg/s';%ocean only |
127 |
|
|
budgMi.specs.units='kg/s';%ice only |
128 |
|
|
%here we output tendencies and fluxes in Watts |
129 |
|
|
budgHo=budgO; budgHi=budgI; |
130 |
|
|
budgHo.specs.units='W';%ocean only |
131 |
|
|
budgHi.specs.units='W';%ice only |
132 |
|
|
%here we output tendencies and fluxes in g/s |
133 |
|
|
budgSo=budgO; budgSi=budgI; |
134 |
|
|
budgSo.specs.units='g/s';%ocean only |
135 |
|
|
budgSi.specs.units='g/s';%ice only |
136 |
|
|
|
137 |
gforget |
1.1 |
%=======MASS========= |
138 |
gforget |
1.17 |
|
139 |
gforget |
1.19 |
[budgMo,budgMi,budgMoi]=calc_budget_mass(kBudget); |
140 |
gforget |
1.1 |
|
141 |
|
|
%bottom pressure for comparison: |
142 |
|
|
bp=myparms.rhoconst/9.81*PHIBOT; |
143 |
|
|
|
144 |
|
|
%compute global integrals: |
145 |
|
|
%------------------------- |
146 |
|
|
msk=mygrid.mskC(:,:,kBudget); |
147 |
gforget |
1.19 |
glo_vol_tot=calc_budget_mean_mask(budgMoi,msk); |
148 |
|
|
glo_vol_ocn=calc_budget_mean_mask(budgMo,msk); |
149 |
|
|
glo_vol_ice=calc_budget_mean_mask(budgMi,msk); |
150 |
gforget |
1.1 |
glo_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC); |
151 |
|
|
|
152 |
|
|
%compute northern hemisphere integrals: |
153 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0); |
154 |
gforget |
1.19 |
north_vol_tot=calc_budget_mean_mask(budgMoi,msk); |
155 |
|
|
north_vol_ocn=calc_budget_mean_mask(budgMo,msk); |
156 |
|
|
north_vol_ice=calc_budget_mean_mask(budgMi,msk); |
157 |
gforget |
1.1 |
north_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC); |
158 |
|
|
|
159 |
|
|
%and southern hemisphere integrals: |
160 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0); |
161 |
gforget |
1.19 |
south_vol_tot=calc_budget_mean_mask(budgMoi,msk); |
162 |
|
|
south_vol_ocn=calc_budget_mean_mask(budgMo,msk); |
163 |
|
|
south_vol_ice=calc_budget_mean_mask(budgMi,msk); |
164 |
gforget |
1.1 |
south_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC); |
165 |
|
|
|
166 |
|
|
%=======HEAT======= |
167 |
gforget |
1.11 |
|
168 |
gforget |
1.19 |
[budgHo,budgHi,budgHoi]=calc_budget_heat(kBudget); |
169 |
gforget |
1.17 |
|
170 |
gforget |
1.1 |
%compute global integrals: |
171 |
|
|
%------------------------- |
172 |
|
|
msk=mygrid.mskC(:,:,kBudget); |
173 |
gforget |
1.19 |
glo_heat_tot=calc_budget_mean_mask(budgHoi,msk); |
174 |
|
|
glo_heat_ocn=calc_budget_mean_mask(budgHo,msk); |
175 |
|
|
glo_heat_ice=calc_budget_mean_mask(budgHi,msk); |
176 |
gforget |
1.1 |
|
177 |
|
|
%compute northern hemisphere integrals: |
178 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0); |
179 |
gforget |
1.19 |
north_heat_tot=calc_budget_mean_mask(budgHoi,msk); |
180 |
|
|
north_heat_ocn=calc_budget_mean_mask(budgHo,msk); |
181 |
|
|
north_heat_ice=calc_budget_mean_mask(budgHi,msk); |
182 |
gforget |
1.1 |
|
183 |
|
|
%and southern hemisphere integrals: |
184 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0); |
185 |
gforget |
1.19 |
south_heat_tot=calc_budget_mean_mask(budgHoi,msk); |
186 |
|
|
south_heat_ocn=calc_budget_mean_mask(budgHo,msk); |
187 |
|
|
south_heat_ice=calc_budget_mean_mask(budgHi,msk); |
188 |
gforget |
1.11 |
|
189 |
gforget |
1.1 |
%=======SALT======= |
190 |
gforget |
1.17 |
|
191 |
gforget |
1.19 |
[budgSo,budgSi,budgSoi]=calc_budget_salt(kBudget); |
192 |
gforget |
1.17 |
|
193 |
gforget |
1.1 |
%compute global integrals: |
194 |
|
|
%------------------------- |
195 |
|
|
msk=mygrid.mskC(:,:,kBudget); |
196 |
gforget |
1.19 |
glo_salt_tot=calc_budget_mean_mask(budgSoi,msk); |
197 |
|
|
glo_salt_ocn=calc_budget_mean_mask(budgSo,msk); |
198 |
|
|
glo_salt_ice=calc_budget_mean_mask(budgSi,msk); |
199 |
gforget |
1.1 |
|
200 |
|
|
%compute northern hemisphere integrals: |
201 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0); |
202 |
gforget |
1.19 |
north_salt_tot=calc_budget_mean_mask(budgSoi,msk); |
203 |
|
|
north_salt_ocn=calc_budget_mean_mask(budgSo,msk); |
204 |
|
|
north_salt_ice=calc_budget_mean_mask(budgSi,msk); |
205 |
gforget |
1.1 |
|
206 |
|
|
%and southern hemisphere integrals: |
207 |
|
|
msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0); |
208 |
gforget |
1.19 |
south_salt_tot=calc_budget_mean_mask(budgSoi,msk); |
209 |
|
|
south_salt_ocn=calc_budget_mean_mask(budgSo,msk); |
210 |
|
|
south_salt_ice=calc_budget_mean_mask(budgSi,msk); |
211 |
gforget |
1.2 |
|
212 |
gforget |
1.9 |
if doMoreBudgetOutput; |
213 |
|
|
%list of budgets to output |
214 |
|
|
listbudg={'budgMo','budgHo','budgSo'}; |
215 |
|
|
if kBudget==1; listbudg={listbudg{:},'budgMi','budgHi','budgSi'}; end; |
216 |
|
|
%the actual output |
217 |
|
|
for iibudg=1:length(listbudg); |
218 |
|
|
%set directory name |
219 |
|
|
dirbudg=dirMat; |
220 |
|
|
if ~isempty(strfind(dirMat,['diags_set_' setDiags '/'])) |
221 |
|
|
dirbudg=fullfile(dirMat,'..',filesep); |
222 |
|
|
end; |
223 |
|
|
sufbudg=''; |
224 |
|
|
if kBudget>1; sufbudg=num2str(kBudget); end; |
225 |
|
|
dirbudg=fullfile(dirbudg,['diags_set_' listbudg{iibudg} sufbudg],filesep); |
226 |
|
|
% |
227 |
|
|
if ~isdir(dirbudg); mkdir(dirbudg); end; |
228 |
|
|
%set file name |
229 |
|
|
filebudg=[listbudg{iibudg} '_' num2str(tt) '.mat']; |
230 |
|
|
%output to file |
231 |
|
|
eval(['tmpbudg=' listbudg{iibudg} ';']); |
232 |
|
|
save([dirbudg filebudg],'-struct','tmpbudg'); |
233 |
|
|
end; |
234 |
|
|
end; |
235 |
gforget |
1.13 |
|
236 |
gforget |
1.2 |
%===================== COMPUTATIONAL SEQUENCE ENDS =========================% |
237 |
|
|
%===================== PLOTTING SEQUENCE BEGINS =========================% |
238 |
|
|
|
239 |
|
|
elseif userStep==-1;%plotting |
240 |
|
|
|
241 |
|
|
if isempty(setDiagsParams); |
242 |
|
|
choicePlot={'all'}; |
243 |
|
|
elseif isnumeric(setDiagsParams{1})&length(setDiagsParams)==1; |
244 |
|
|
choicePlot={'all'}; |
245 |
|
|
elseif isnumeric(setDiagsParams{1}); |
246 |
|
|
choicePlot={setDiagsParams{2:end}}; |
247 |
|
|
else; |
248 |
|
|
choicePlot=setDiagsParams; |
249 |
|
|
end; |
250 |
|
|
|
251 |
gforget |
1.3 |
tt=[1:length(alldiag.listTimes)]; |
252 |
|
|
TT=alldiag.listTimes(tt); |
253 |
|
|
nt=length(TT); |
254 |
|
|
|
255 |
gforget |
1.2 |
if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'mass'))); |
256 |
|
|
|
257 |
|
|
%1.1) ocean+seaice mass budgets |
258 |
|
|
%------------------------------ |
259 |
|
|
figureL; |
260 |
|
|
%global volume budget: |
261 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_tot,'kg/m^2','Global Mean Mass (incl. ice)'); |
262 |
gforget |
1.2 |
%add bp: |
263 |
|
|
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp); |
264 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
265 |
|
|
%northern hemisphere budget: |
266 |
heimbach |
1.6 |
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_tot,'kg/m^2','Northern Mean Mass (incl. ice)'); |
267 |
gforget |
1.2 |
%add bp: |
268 |
|
|
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp); |
269 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
270 |
|
|
%southern hemisphere budget: |
271 |
heimbach |
1.6 |
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_tot,'kg/m^2','Southern Mean Mass (incl. ice)'); |
272 |
gforget |
1.2 |
%add bp: |
273 |
|
|
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp); |
274 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
275 |
|
|
%add to tex file |
276 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
277 |
heimbach |
1.6 |
myCaption={myCaption{:},'mass budget (ocean+ice) in kg/m$^2$.'}; |
278 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
279 |
|
|
|
280 |
|
|
%1.2) ice mass budgets |
281 |
|
|
%--------------------- |
282 |
|
|
figureL; |
283 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_ice,'kg/m^2','Global Mean Mass (only ice)'); |
284 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp); |
285 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
286 |
heimbach |
1.6 |
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_ice,'kg/m^2','Northern Mean Mass (only ice)'); |
287 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp); |
288 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
289 |
heimbach |
1.6 |
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_ice,'kg/m^2','Southern Mean Mass (only ice)'); |
290 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp); |
291 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
292 |
|
|
%add to tex file |
293 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
294 |
heimbach |
1.6 |
myCaption={myCaption{:},'mass budget (ice only) in kg/m$^2$.'}; |
295 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
296 |
|
|
|
297 |
|
|
end; |
298 |
|
|
|
299 |
|
|
if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'mass'))); |
300 |
|
|
|
301 |
|
|
%1.3) ocean mass budgets |
302 |
|
|
%----------------------- |
303 |
|
|
figureL; |
304 |
|
|
%global volume budget: |
305 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_ocn,'kg/m^2','Global Mean Mass (only ocean)'); |
306 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp); |
307 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
308 |
heimbach |
1.6 |
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_ocn,'kg/m^2','Northern Mean Mass (only ocean)'); |
309 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp); |
310 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
311 |
heimbach |
1.6 |
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_ocn,'kg/m^2','Southern Mean Mass (only ocean)'); |
312 |
gforget |
1.2 |
dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp); |
313 |
|
|
plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal'); |
314 |
|
|
%add to tex file |
315 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
316 |
heimbach |
1.6 |
myCaption={myCaption{:},'mass budget (ocean only) in kg/m$^2$.'}; |
317 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
318 |
|
|
|
319 |
|
|
end; |
320 |
|
|
|
321 |
|
|
if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'heat'))); |
322 |
|
|
|
323 |
|
|
%2.1) ocean+seaice heat budgets |
324 |
|
|
%------------------------------ |
325 |
|
|
figureL; |
326 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_heat_tot,'J/m^2','Global Mean Ocean Heat (incl. ice)'); |
327 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_tot,'J/m^2','Northern Mean Ocean Heat (incl. ice)'); |
328 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_tot,'J/m^2','Southern Mean Ocean Heat (incl. ice)'); |
329 |
gforget |
1.2 |
%add to tex file |
330 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
331 |
heimbach |
1.6 |
myCaption={myCaption{:},'heat budget (ocean+ice) in J/m$^2$.'}; |
332 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
333 |
|
|
|
334 |
|
|
%2.2) ice heat budgets |
335 |
|
|
%--------------------- |
336 |
|
|
figureL; |
337 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_heat_ice,'J/m^2','Global Mean Ocean Heat (only ice)'); |
338 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_ice,'J/m^2','Northern Mean Ocean Heat (only ice)'); |
339 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_ice,'J/m^2','Southern Mean Ocean Heat (only ice)'); |
340 |
gforget |
1.2 |
%add to tex file |
341 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
342 |
heimbach |
1.6 |
myCaption={myCaption{:},'heat budget (ice only) in J/m$^2$.'}; |
343 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
344 |
|
|
|
345 |
|
|
end; |
346 |
|
|
|
347 |
|
|
if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'heat'))); |
348 |
|
|
|
349 |
|
|
%2.3) ocean heat budgets |
350 |
|
|
%----------------------- |
351 |
|
|
figureL; |
352 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_heat_ocn,'J/m^2','Global Mean Ocean Heat (only ocean)'); |
353 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_ocn,'J/m^2','Northern Mean Ocean Heat (only ocean)'); |
354 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_ocn,'J/m^2','Southern Mean Ocean Heat (only ocean)'); |
355 |
gforget |
1.2 |
%add to tex file |
356 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
357 |
heimbach |
1.6 |
myCaption={myCaption{:},'heat budget (ocean only) in J/m$^2$.'}; |
358 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
359 |
|
|
|
360 |
|
|
end; |
361 |
|
|
|
362 |
|
|
if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'salt'))); |
363 |
|
|
|
364 |
|
|
%3.1) ocean+seaice salt budgets |
365 |
|
|
%------------------------------ |
366 |
|
|
figureL; |
367 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_salt_tot,'g/m^2','Global Mean Ocean Salt (incl. ice)'); |
368 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_tot,'g/m^2','Northern Mean Ocean Salt (incl. ice)'); |
369 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_tot,'g/m^2','Southern Mean Ocean Salt (incl. ice)'); |
370 |
gforget |
1.2 |
%add to tex file |
371 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
372 |
heimbach |
1.6 |
myCaption={myCaption{:},'salt budget (ocean+ice) in g/m$^2$.'}; |
373 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
374 |
|
|
|
375 |
|
|
%2.2) ice salt budgets |
376 |
|
|
%--------------------- |
377 |
|
|
figureL; |
378 |
heimbach |
1.6 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_salt_ice,'g/m^2','Global Mean Ocean Salt (only ice)'); |
379 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_ice,'g/m^2','Northern Mean Ocean Salt (only ice)'); |
380 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_ice,'g/m^2','Southern Mean Ocean Salt (only ice)'); |
381 |
gforget |
1.2 |
%add to tex file |
382 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
383 |
heimbach |
1.6 |
myCaption={myCaption{:},'salt budget (ice only) in g/m$^2$.'}; |
384 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
385 |
|
|
|
386 |
|
|
end; |
387 |
|
|
|
388 |
|
|
|
389 |
|
|
if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'salt'))); |
390 |
|
|
|
391 |
|
|
%3.3) ocean salt budgets |
392 |
|
|
%----------------------- |
393 |
|
|
figureL; |
394 |
gforget |
1.8 |
subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_salt_ocn,'g/m^2','Global Mean Ocean Salt (only ocean)'); |
395 |
|
|
subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_ocn,'g/m^2','Northern Mean Ocean Salt (only ocean)'); |
396 |
|
|
subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_ocn,'g/m^2','Southern Mean Ocean Salt (only ocean)'); |
397 |
gforget |
1.2 |
%add to tex file |
398 |
|
|
myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '}; |
399 |
heimbach |
1.6 |
myCaption={myCaption{:},'salt budget (ocean only) in g/m$^2$.'}; |
400 |
gforget |
1.2 |
if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end; |
401 |
|
|
|
402 |
|
|
end; |
403 |
|
|
|
404 |
gforget |
1.1 |
end; |