[Dart-dev] DART/branches Revision: 12984
dart at ucar.edu
dart at ucar.edu
Tue Feb 12 15:31:17 MST 2019
thoar at ucar.edu
2019-02-12 15:31:16 -0700 (Tue, 12 Feb 2019)
216
Handles obs_diag_output files with no posterior information.
2018b complained about the wysiwyg() ... issued warnings about
figure being too large and getting cut off, etc.
Removing wysiwyg resolved the problem.
Modified: DART/branches/recam/diagnostics/matlab/two_experiments_evolution.m
===================================================================
--- DART/branches/recam/diagnostics/matlab/two_experiments_evolution.m 2019-02-12 18:44:53 UTC (rev 12983)
+++ DART/branches/recam/diagnostics/matlab/two_experiments_evolution.m 2019-02-12 22:31:16 UTC (rev 12984)
@@ -108,7 +108,7 @@
%% set up all the stuff that is common.
-commondata = check_compatibility(files, obsnames, copy);
+commondata = check_compatibility(files, prpo, obsnames, copy);
figuredata = setfigure(NumExp);
%%--------------------------------------------------------------------
@@ -119,48 +119,48 @@
for ivar = 1:nvars
fprintf('Working on %s ...\n',obsnames{ivar})
-
+
%------------------------------------------------------------------------
% Plot each region in a separate figure window.
%------------------------------------------------------------------------
-
+
for iregion = 1:commondata.nregions
-
+
figure(iregion);
clf(iregion);
orient(figuredata.orientation);
- wysiwyg;
-
+
%---------------------------------------------------------------------
% 1) Get the data for each experiment
% 2) plot the data
% 3) annotate
%---------------------------------------------------------------------
-
+
for iexp = 1:NumExp
-
- plotobj{iexp} = getvals(files{iexp}, obsnames{ivar}, copy, prpo, iregion, p.Results.level);
+
+ plotobj{iexp} = getvals(files{iexp}, commondata.targets{ivar}, copy, iregion, p.Results.level);
plotobj{iexp}.title = titles{iexp};
plotobj{iexp}.nregions = commondata.nregions;
plotobj{iexp}.region_names = commondata.region_names;
-
+ plotobj{iexp}.phase = commondata.phase;
+
end
-
+
myplot(plotobj, figuredata);
-
+
BottomAnnotation(plotobj)
-
+
psfname = sprintf('%s_%s_region%d_ilev%d_evolution_%dexp', ...
obsnames{ivar}, plotobj{1}.copystring, iregion, p.Results.level, NumExp);
print(iregion,'-dpdf',psfname)
-
+
end % of loop around regions
-
+
if ( ivar ~= nvars )
disp('Pausing, hit any key to continue ...')
pause
end
-
+
end % of loop around variable
@@ -171,7 +171,7 @@
-function common = check_compatibility(filenames, varnames, copystring)
+function common = check_compatibility(filenames, prpo, varnames, copystring)
%% Trying to prevent the comparison of apples and oranges.
% make sure the diagnostics were generated the same way.
@@ -182,18 +182,26 @@
mystat = 0;
nexp = length(filenames);
commondata = cell(1,nexp);
-priornames = struct([]);
-postenames = struct([]);
+targets = struct([]);
for i = 1:length(varnames)
- priornames{i} = sprintf('%s_guess',varnames{i});
- postenames{i} = sprintf('%s_analy',varnames{i});
+ switch lower(prpo)
+ case {'guess','forecast','prior'}
+ targets{i} = sprintf('%s_guess',varnames{i});
+ commondata{i}.phase = 'prior';
More information about the Dart-dev
mailing list