Changeset 21455

Timestamp:

12/26/16 18:37:55 (8 years ago)

Author:

Eric.Larour

Message:

CHG: more documentation and minor fixes.

File:

• issm/trunk-jpl/src/m/miscellaneous/diagnostics.m (modified) (2 diffs)

issm/trunk-jpl/src/m/miscellaneous/diagnostics.m

@@ -1 +1 @@
 function diagnostics(md,varargin) 
-%DIAGNOSTICS: output all sort of information as shape files.
+%DIAGNOSTICS - output diagnostics for use in qgis
+%
+%
+%   Usage:
+%      diagnostics(md,options)
+%   
+%      where options include: 
+%            'path': where are the diagnostics file output
+%            'mesh': 0 or 1 (output mesh?, default 0)
+%            'gl': 0 or 1 (output grounding line position?, default 0)
+%            'vel': defaults to md.initialization.vel (can suplly another field if need be)
+%                Suboptions include:
+%                 'velposting': resolution at which the tiff file will be printed.
+%            'hotspots' (0 or 1, default 0). spots where md.results.TransientSolution(i).Vel>threshold
+%                Suboptions include:
+%                 'threshold': vel value
+%                 'hotspotsi': index into results.TransientSolution
+%
+%           
+%
+%   Examples:
+%      diagnostics(md,'vel',md.initialization.vel,'gl',1,'hotspots',1,'hotspotsi',10,'threshold',500);
+
 
         %process options: 
@@ -31 +53 @@
                 vel=uint16(flipud(vel));
 
-                imwrite(vel,[path '/vel.tif'],'tiff');
+                imwrite(vel,[path '/veltemp.tif'],'tiff');
 
-                string=sprintf('!gdal_translate -a_srs EPSG:3031 -a_ullr %g %g %g %g %s/vel.tif %s/velg.tif',...
+                string=sprintf('!gdal_translate -a_srs EPSG:3031 -a_ullr %g %g %g %g %s/veltemp.tif %s/vel.tif',...
                 min(md.mesh.x),max(md.mesh.y),max(md.mesh.x),min(md.mesh.y),path,path);
                 eval(string);
-
+                delete([path '/veltemp.tif']);
         end