Changeset 10295

Timestamp:

10/25/11 13:26:58 (13 years ago)

Author:

schlegel

Message:

Added comments

File:

• issm/trunk/src/m/model/regionaltransient2d.m (modified) (5 diffs)

issm/trunk/src/m/model/regionaltransient2d.m

@@ -37 +37 @@
 
 %remesh
-%md2=bamg(md2,'hmax',hmax,'hmin',hmin,'field',[md2.inversion.vel_obs md2.geometry.thickness],'splitcorner',1,'KeepVertices',0);
 md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.thickness],'splitcorner',1,'KeepVertices',0,'err',err);
 md2=setmask(md2,'','');
@@ -51 +50 @@
                 fieldsize=size(field);
 
+                %copy field, interpolated to new mesh
                 if isobject(field), %recursive call
                         object_fields=fields(md1.(model_fields{i}));
@@ -92 +92 @@
         end
 
+        %Read transient velocities looping through only the populated times
         spcx=[];
         spcy=[];
@@ -109 +110 @@
         end
 
+        %As long as there are recorded time steps, spc the boundaries with velocities
         if nsteps > 0
                 md2.diagnostic.spcvx=md2.diagnostic.spcvx*ones(1,size(spcx,2));
@@ -121 +123 @@
         end
 
-        %Diagnostic
+        %Diagnostic.  Don't spc the icefront vertices.
         if ~isnan(md2.diagnostic.icefront)
                 md1s=modelextract(md1,area);