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

Index: ../trunk-jpl/src/m/classes/taoinversion.m
===================================================================
--- ../trunk-jpl/src/m/classes/taoinversion.m   (revision 18667)
+++ ../trunk-jpl/src/m/classes/taoinversion.m   (revision 18668)
@@ -25,6 +25,7 @@
                vz_obs                      = NaN
                vel_obs                     = NaN
                thickness_obs               = NaN
+               surface_obs               = NaN
        end
        methods
                function obj = taoinversion(varargin) % {{{
@@ -123,6 +124,7 @@
                        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,6 +159,7 @@
                        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);