source: issm/oecreview/Archive/18296-19100/ISSM-18667-18668.diff

../trunk-jpl/src/m/classes/taoinversion.m

@@ -25 +25 @@
                 vz_obs                      = NaN
                 vel_obs                     = NaN
                 thickness_obs               = NaN
+                surface_obs               = NaN
         end
         methods
                 function obj = taoinversion(varargin) % {{{
@@ -123 +124 @@
                         fielddisplay(obj,'vy_obs','observed velocity y component [m/yr]');
                         fielddisplay(obj,'vel_obs','observed velocity magnitude [m/yr]');
                         fielddisplay(obj,'thickness_obs','observed thickness [m]');
+                        fielddisplay(obj,'surface_obs','observed surface elevation [m]');
                         disp('Available cost functions:');
                         disp('   101: SurfaceAbsVelMisfit');
                         disp('   102: SurfaceRelVelMisfit');
@@ -157 +159 @@
                         WriteData(fid,'object',obj,'class','inversion','fieldname','vy_obs','format','DoubleMat','mattype',1,'scale',1./yts);
                         WriteData(fid,'object',obj,'class','inversion','fieldname','vz_obs','format','DoubleMat','mattype',1,'scale',1./yts);
                         WriteData(fid,'object',obj,'class','inversion','fieldname','thickness_obs','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',obj,'class','inversion','fieldname','surface_obs','format','DoubleMat','mattype',1);
 
                         %process control parameters
                         num_control_parameters=numel(obj.control_parameters);