/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m
ViewVC logotype

Annotation of /MITgcm_contrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.23 - (hide annotations) (download)
Wed Feb 10 15:04:34 2016 UTC (9 years, 5 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint65t
Changes since 1.22: +2 -0 lines 
- add diags_inifin_D.m to compute initial and final budget states
- add call to diags_inifin_D in diags_set_D
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 gforget 1.20
144 gforget 1.1 %compute global integrals:
145     %-------------------------
146     msk=mygrid.mskC(:,:,kBudget);
147 gforget 1.20 tmp1=calc_mskmean_T(budgMoi,msk,'extensive');
148 gforget 1.21 glo_vol_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
149 gforget 1.20 tmp1=calc_mskmean_T(budgMo,msk,'extensive');
150 gforget 1.21 glo_vol_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
151 gforget 1.20 tmp1=calc_mskmean_T(budgMi,msk,'extensive');
152 gforget 1.21 glo_vol_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
153 gforget 1.20
154 gforget 1.1 glo_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC);
155    
156     %compute northern hemisphere integrals:
157     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0);
158 gforget 1.20 tmp1=calc_mskmean_T(budgMoi,msk,'extensive');
159 gforget 1.21 north_vol_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
160 gforget 1.20 tmp1=calc_mskmean_T(budgMo,msk,'extensive');
161 gforget 1.21 north_vol_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
162 gforget 1.20 tmp1=calc_mskmean_T(budgMi,msk,'extensive');
163 gforget 1.21 north_vol_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
164 gforget 1.20
165 gforget 1.1 north_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC);
166    
167     %and southern hemisphere integrals:
168     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0);
169 gforget 1.20 tmp1=calc_mskmean_T(budgMoi,msk,'extensive');
170 gforget 1.21 south_vol_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
171 gforget 1.20 tmp1=calc_mskmean_T(budgMo,msk,'extensive');
172 gforget 1.21 south_vol_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
173 gforget 1.20 tmp1=calc_mskmean_T(budgMi,msk,'extensive');
174 gforget 1.21 south_vol_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
175 gforget 1.20
176 gforget 1.1 south_bp=nansum(bp.*msk.*mygrid.RAC)/nansum(msk.*mygrid.RAC);
177    
178     %=======HEAT=======
179 gforget 1.11
180 gforget 1.19 [budgHo,budgHi,budgHoi]=calc_budget_heat(kBudget);
181 gforget 1.17
182 gforget 1.1 %compute global integrals:
183     %-------------------------
184     msk=mygrid.mskC(:,:,kBudget);
185 gforget 1.20 tmp1=calc_mskmean_T(budgHoi,msk,'extensive');
186 gforget 1.21 glo_heat_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
187 gforget 1.20 tmp1=calc_mskmean_T(budgHo,msk,'extensive');
188 gforget 1.21 glo_heat_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
189 gforget 1.20 tmp1=calc_mskmean_T(budgHi,msk,'extensive');
190 gforget 1.21 glo_heat_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
191 gforget 1.1
192     %compute northern hemisphere integrals:
193     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0);
194 gforget 1.20 tmp1=calc_mskmean_T(budgHoi,msk,'extensive');
195 gforget 1.21 north_heat_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
196 gforget 1.20 tmp1=calc_mskmean_T(budgHo,msk,'extensive');
197 gforget 1.21 north_heat_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
198 gforget 1.20 tmp1=calc_mskmean_T(budgHi,msk,'extensive');
199 gforget 1.21 north_heat_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
200 gforget 1.1
201     %and southern hemisphere integrals:
202     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0);
203 gforget 1.20 tmp1=calc_mskmean_T(budgHoi,msk,'extensive');
204 gforget 1.21 south_heat_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
205 gforget 1.20 tmp1=calc_mskmean_T(budgHo,msk,'extensive');
206 gforget 1.21 south_heat_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
207 gforget 1.20 tmp1=calc_mskmean_T(budgHi,msk,'extensive');
208 gforget 1.21 south_heat_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
209 gforget 1.11
210 gforget 1.1 %=======SALT=======
211 gforget 1.17
212 gforget 1.19 [budgSo,budgSi,budgSoi]=calc_budget_salt(kBudget);
213 gforget 1.17
214 gforget 1.1 %compute global integrals:
215     %-------------------------
216     msk=mygrid.mskC(:,:,kBudget);
217 gforget 1.20 tmp1=calc_mskmean_T(budgSoi,msk,'extensive');
218 gforget 1.21 glo_salt_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
219 gforget 1.20 tmp1=calc_mskmean_T(budgSo,msk,'extensive');
220 gforget 1.21 glo_salt_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
221 gforget 1.20 tmp1=calc_mskmean_T(budgSi,msk,'extensive');
222 gforget 1.21 glo_salt_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
223 gforget 1.20
224 gforget 1.1 %compute northern hemisphere integrals:
225     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC>0);
226 gforget 1.20 tmp1=calc_mskmean_T(budgSoi,msk,'extensive');
227 gforget 1.21 north_salt_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
228 gforget 1.20 tmp1=calc_mskmean_T(budgSo,msk,'extensive');
229 gforget 1.21 north_salt_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
230 gforget 1.20 tmp1=calc_mskmean_T(budgSi,msk,'extensive');
231 gforget 1.21 north_salt_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
232 gforget 1.1
233     %and southern hemisphere integrals:
234     msk=mygrid.mskC(:,:,kBudget).*(mygrid.YC<=0);
235 gforget 1.20 tmp1=calc_mskmean_T(budgSoi,msk,'extensive');
236 gforget 1.21 south_salt_tot=[tmp1.tend;tmp1.hconv;tmp1.zconv];
237 gforget 1.20 tmp1=calc_mskmean_T(budgSo,msk,'extensive');
238 gforget 1.21 south_salt_ocn=[tmp1.tend;tmp1.hconv;tmp1.zconv];
239 gforget 1.20 tmp1=calc_mskmean_T(budgSi,msk,'extensive');
240 gforget 1.21 south_salt_ice=[tmp1.tend;tmp1.hconv;tmp1.zconv];
241 gforget 1.2
242 gforget 1.9 if doMoreBudgetOutput;
243 gforget 1.23 %initial and final state:
244     if ii==1; diags_inifin_D(kBudget,dirSnap,dirMat); end;
245 gforget 1.9 %list of budgets to output
246     listbudg={'budgMo','budgHo','budgSo'};
247     if kBudget==1; listbudg={listbudg{:},'budgMi','budgHi','budgSi'}; end;
248     %the actual output
249     for iibudg=1:length(listbudg);
250     %set directory name
251     dirbudg=dirMat;
252     if ~isempty(strfind(dirMat,['diags_set_' setDiags '/']))
253     dirbudg=fullfile(dirMat,'..',filesep);
254     end;
255     sufbudg='';
256     if kBudget>1; sufbudg=num2str(kBudget); end;
257     dirbudg=fullfile(dirbudg,['diags_set_' listbudg{iibudg} sufbudg],filesep);
258     %
259     if ~isdir(dirbudg); mkdir(dirbudg); end;
260     %set file name
261     filebudg=[listbudg{iibudg} '_' num2str(tt) '.mat'];
262     %output to file
263 gforget 1.22 eval(['tmpbudg=' listbudg{iibudg} '.fluxes;']);
264     eval(['tmpbudg.tend=' listbudg{iibudg} '.tend;']);
265     listterms=fieldnames(tmpbudg);
266     for iiterm=1:length(listterms);
267     tmp1=getfield(tmpbudg,listterms{iiterm});
268     tmp1=convert2gcmfaces(tmp1);
269     %tmp2=prod(size(tmp1));
270     fid=fopen([dirbudg listterms{iiterm} '.bin'],'a+','b');
271     %status=fseek(fid,(tt-1)*recl2D,'bof');
272     fwrite(fid,tmp1,'float64');
273     fclose(fid);
274     end;
275 gforget 1.9 end;
276     end;
277 gforget 1.13
278 gforget 1.2 %===================== COMPUTATIONAL SEQUENCE ENDS =========================%
279     %===================== PLOTTING SEQUENCE BEGINS =========================%
280    
281     elseif userStep==-1;%plotting
282    
283     if isempty(setDiagsParams);
284     choicePlot={'all'};
285     elseif isnumeric(setDiagsParams{1})&length(setDiagsParams)==1;
286     choicePlot={'all'};
287     elseif isnumeric(setDiagsParams{1});
288     choicePlot={setDiagsParams{2:end}};
289     else;
290     choicePlot=setDiagsParams;
291     end;
292    
293 gforget 1.3 tt=[1:length(alldiag.listTimes)];
294     TT=alldiag.listTimes(tt);
295     nt=length(TT);
296    
297 gforget 1.2 if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'mass')));
298    
299     %1.1) ocean+seaice mass budgets
300     %------------------------------
301     figureL;
302     %global volume budget:
303 heimbach 1.6 subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_tot,'kg/m^2','Global Mean Mass (incl. ice)');
304 gforget 1.2 %add bp:
305     dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp);
306     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
307     %northern hemisphere budget:
308 heimbach 1.6 subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_tot,'kg/m^2','Northern Mean Mass (incl. ice)');
309 gforget 1.2 %add bp:
310     dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp);
311     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
312     %southern hemisphere budget:
313 heimbach 1.6 subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_tot,'kg/m^2','Southern Mean Mass (incl. ice)');
314 gforget 1.2 %add bp:
315     dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp);
316     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
317     %add to tex file
318     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
319 heimbach 1.6 myCaption={myCaption{:},'mass budget (ocean+ice) in kg/m$^2$.'};
320 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
321    
322     %1.2) ice mass budgets
323     %---------------------
324     figureL;
325 heimbach 1.6 subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_ice,'kg/m^2','Global Mean Mass (only ice)');
326 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp);
327     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
328 heimbach 1.6 subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_ice,'kg/m^2','Northern Mean Mass (only ice)');
329 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp);
330     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
331 heimbach 1.6 subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_ice,'kg/m^2','Southern Mean Mass (only ice)');
332 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp);
333     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
334     %add to tex file
335     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
336 heimbach 1.6 myCaption={myCaption{:},'mass budget (ice only) in kg/m$^2$.'};
337 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
338    
339     end;
340    
341     if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'mass')));
342    
343     %1.3) ocean mass budgets
344     %-----------------------
345     figureL;
346     %global volume budget:
347 heimbach 1.6 subplot(3,1,1); disp_budget_mean_mask(TT,alldiag.glo_vol_ocn,'kg/m^2','Global Mean Mass (only ocean)');
348 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.glo_bp);
349     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
350 heimbach 1.6 subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_vol_ocn,'kg/m^2','Northern Mean Mass (only ocean)');
351 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.north_bp);
352     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
353 heimbach 1.6 subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_vol_ocn,'kg/m^2','Southern Mean Mass (only ocean)');
354 gforget 1.2 dt=median(diff(TT))*86400; bp=dt*cumsum(alldiag.south_bp);
355     plot(TT,bp,'k'); aa=legend; bb=get(aa,'String'); bb={bb{:},'bp'}; legend(bb,'Orientation','horizontal');
356     %add to tex file
357     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
358 heimbach 1.6 myCaption={myCaption{:},'mass budget (ocean only) in kg/m$^2$.'};
359 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
360    
361     end;
362    
363     if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'heat')));
364    
365     %2.1) ocean+seaice heat budgets
366     %------------------------------
367     figureL;
368 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)');
369     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_tot,'J/m^2','Northern Mean Ocean Heat (incl. ice)');
370     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_tot,'J/m^2','Southern Mean Ocean Heat (incl. ice)');
371 gforget 1.2 %add to tex file
372     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
373 heimbach 1.6 myCaption={myCaption{:},'heat budget (ocean+ice) in J/m$^2$.'};
374 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
375    
376     %2.2) ice heat budgets
377     %---------------------
378     figureL;
379 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)');
380     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_ice,'J/m^2','Northern Mean Ocean Heat (only ice)');
381     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_ice,'J/m^2','Southern Mean Ocean Heat (only ice)');
382 gforget 1.2 %add to tex file
383     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
384 heimbach 1.6 myCaption={myCaption{:},'heat budget (ice only) in J/m$^2$.'};
385 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
386    
387     end;
388    
389     if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'heat')));
390    
391     %2.3) ocean heat budgets
392     %-----------------------
393     figureL;
394 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)');
395     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_heat_ocn,'J/m^2','Northern Mean Ocean Heat (only ocean)');
396     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_heat_ocn,'J/m^2','Southern Mean Ocean Heat (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{:},'heat budget (ocean only) in J/m$^2$.'};
400 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
401    
402     end;
403    
404     if (kBudget==1)&(sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'salt')));
405    
406     %3.1) ocean+seaice salt budgets
407     %------------------------------
408     figureL;
409 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)');
410     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_tot,'g/m^2','Northern Mean Ocean Salt (incl. ice)');
411     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_tot,'g/m^2','Southern Mean Ocean Salt (incl. ice)');
412 gforget 1.2 %add to tex file
413     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
414 heimbach 1.6 myCaption={myCaption{:},'salt budget (ocean+ice) in g/m$^2$.'};
415 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
416    
417     %2.2) ice salt budgets
418     %---------------------
419     figureL;
420 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)');
421     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_ice,'g/m^2','Northern Mean Ocean Salt (only ice)');
422     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_ice,'g/m^2','Southern Mean Ocean Salt (only ice)');
423 gforget 1.2 %add to tex file
424     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
425 heimbach 1.6 myCaption={myCaption{:},'salt budget (ice only) in g/m$^2$.'};
426 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
427    
428     end;
429    
430    
431     if (sum(strcmp(choicePlot,'all'))|sum(strcmp(choicePlot,'salt')));
432    
433     %3.3) ocean salt budgets
434     %-----------------------
435     figureL;
436 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)');
437     subplot(3,1,2); disp_budget_mean_mask(TT,alldiag.north_salt_ocn,'g/m^2','Northern Mean Ocean Salt (only ocean)');
438     subplot(3,1,3); disp_budget_mean_mask(TT,alldiag.south_salt_ocn,'g/m^2','Southern Mean Ocean Salt (only ocean)');
439 gforget 1.2 %add to tex file
440     myCaption={myYmeanTxt,' global (upper) north (mid) and south (lower), '};
441 heimbach 1.6 myCaption={myCaption{:},'salt budget (ocean only) in g/m$^2$.'};
442 gforget 1.2 if addToTex&multiTimes; write2tex(fileTex,2,myCaption,gcf); elseif ~multiTimes; close; end;
443    
444     end;
445    
446 gforget 1.1 end;
  ViewVC Help
Powered by ViewVC 1.1.22