Changeset 21809

Timestamp:

07/18/17 16:53:01 (8 years ago)

Author:

schlegel

Message:

CHG: add an extrude function for misfit output def

Location:

issm/trunk-jpl/src/m/classes

Files:

• SMBforcing.m (modified) (1 diff)
• misfit.m (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/SMBforcing.m

@@ -68 +68 @@
 
                         WriteData(fid,prefix,'name','md.smb.model','data',1,'format','Integer');
-                        WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
-                        %WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','CompressedMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
+                        %WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
+                        WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','CompressedMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
                         
                         %process requested outputs

issm/trunk-jpl/src/m/classes/misfit.m

@@ -19 +19 @@
                 %misfit
                 name                                    = '';
-                definitionstring                = ''; %string that identifies this output definition uniquely, from 'Outputdefinition[1-10]'
+                definitionstring                = ''; %string that identifies this output definition uniquely, from 'Outputdefinition[1-100]'
                 model_string                    = ''; %string for field that is modeled
                 observation                             = NaN; %observed field that we compare the model against
@@ -31 +31 @@
         
         methods
+                function self = extrude(self,md) % {{{
+                        if ~isnan(self.weights)
+                                self.weights=project3d(md,'vector',self.weights,'type','node');
+                        end
+                        if ~isnan(self.observation)
+                                self.observation=project3d(md,'vector',self.observation,'type','node');
+                        end
+                end % }}}
                 function self = misfit(varargin) % {{{
                         if nargin==0,