[Dart-dev] DART/branches Revision: 11147

[dart at ucar.edu]

thoar at ucar.edu
2017-02-23 11:32:23 -0700 (Thu, 23 Feb 2017)
212
Removing a horribly broken routine that should never have made it out of the development branch: plot_coverage.m.

Last of the observation-space matlab tools to be pruned of the third-party netCDF dependencies.




Deleted: DART/branches/rma_trunk/diagnostics/matlab/plot_coverage.m
===================================================================
--- DART/branches/rma_trunk/diagnostics/matlab/plot_coverage.m	2017-02-23 18:09:28 UTC (rev 11146)
+++ DART/branches/rma_trunk/diagnostics/matlab/plot_coverage.m	2017-02-23 18:32:23 UTC (rev 11147)
@@ -1,257 +0,0 @@
-function obs = plot_coverage(fname)
-%% plot_coverage examines the spatial and temporal coverage for an observation type.
-%  obs_seq_coverage must create the netCDF file read by this routine.
-% 
-% fname = 'obs_coverage.nc';
-%
-% plot_coverage(fname)
-
-%% DART software - Copyright UCAR. This open source software is provided
-% by UCAR, "as is", without charge, subject to all terms of use at
-% http://www.image.ucar.edu/DAReS/DART/DART_download
-%
-% DART $Id$
-
-if (exist(fname,'file') ~= 2)
-   error('%s does not exist.',fname)
-end
-
-locs    = nc_varget(fname,'location');
-
-obs.fname         = fname;
-obs.ObsTypeString = nc_read_att(fname,nc_global,'obs_of_interest');
-obs.stations      = nc_varget(fname,'stations');
-obs.ntimes        = nc_varget(fname,'ntimes');
-obs.times         = nc_varget(fname,'time');
-obs.timeunits     = nc_read_att(fname,'time','units');
-
-obs.lons       = locs(:,1);
-obs.lats       = locs(:,2);
-obs.vert       = locs(:,3);
-obs.region     = [ min(obs.lons) max(obs.lons) min(obs.lats) max(obs.lats) ];
-timebase       = sscanf(obs.timeunits,'%*s%*s%d%*c%d%*c%d'); % YYYY MM DD
-obs.timeorigin = datenum(timebase(1),timebase(2),timebase(3));
-obs.times      = obs.times + obs.timeorigin;
-
-allt     = obs.times(:);
-obs.tmin = datestr(min(allt));
-obs.tmax = datestr(max(allt));
-clear allt;
-
-inds = find(obs.ntimes == 8);
-fprintf('There are %d stations with precisely 8 observation times.\n',length(inds))
-
-%% Now pack the data in the same fashion as the cell array of column labels.
-
-obs.lonindex  = 1;
-obs.latindex  = 2;
-obs.zindex    = 3;
-obs.ntindex   = 4;
-obs.Iindex    = 5;
-
-obs.colnames{obs.lonindex}  = 'longitude';
-obs.colnames{obs.latindex}  = 'latitude';
-obs.colnames{obs.zindex}    = 'vertical';
-obs.colnames{obs.ntindex}   = 'num times';
-obs.colnames{obs.Iindex}    = 'stationID';
-
-global obsmat
-obsmat = zeros(length(obs.lons),5);
-obsmat(:,obs.lonindex ) = obs.lons;
-obsmat(:,obs.latindex ) = obs.lats;
-obsmat(:,obs.zindex   ) = obs.vert;
-obsmat(:,obs.ntindex  ) = obs.ntimes;
-obsmat(:,obs.Iindex   ) = obs.stations;
-
-global timemat
-timemat = obs.times;
-
-%% create the linked plots
-
-figure1 = figure(1); clf(figure1); orient tall; wysiwyg
-
-%----------------------------------------------------------------------
-% Top half is a 2D scatterplot
-%----------------------------------------------------------------------
-
-axes0 = axes('Parent',figure1,'OuterPosition',[0 0.4 1 0.57],'FontSize',16);
-view(axes0,[0 90]);
-grid(axes0,'on');
-hold(axes0,'all');
-
-xstring = sprintf('obsmat(:,%d)',obs.lonindex);
-ystring = sprintf('obsmat(:,%d)',obs.latindex);
-
-scatter(obsmat(:,obs.lonindex), obsmat(:,obs.latindex), 18, obsmat(:,obs.ntindex), 'filled', ...
-             'Parent',axes0,'DisplayName','observation locations', ...
-             'XDataSource',xstring, ...
-             'YDataSource',ystring);
-% myworldmap(obs);
-
-xlabel(obs.colnames{obs.lonindex});
-ylabel(obs.colnames{obs.latindex});
-
-hc = colorbar;
-set(get(hc,'YLabel'),'String','# of obs times','Interpreter','none')
-