Changeset 10587

Timestamp:

11/10/11 15:35:07 (13 years ago)

Author:

schlegel

Message:

add spcthickness

File:

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

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

@@ -37 +37 @@
 
 %remesh
-md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.thickness],'splitcorner',1,'KeepVertices',0,'err',err);
+md2=bamg(mde,'hmin',hmin,'hmax',hmax,'field',[mde.inversion.vel_obs mde.geometry.surface],'splitcorner',1,'KeepVertices',0,'err',err);
 md2=setmask(md2,'','');
 
@@ -92 +92 @@
         end
 
-        %Read transient velocities looping through only the populated times
+        %Read transient velocities and thickness, looping through only the populated times
         spcx=[];
         spcy=[];
+        spct=[];
         steps=[];
         nsteps=length(md1.results.TransientSolution);
@@ -103 +104 @@
                         vx=PatchToVec(md1.results.TransientSolution(t).Vx);
                         vy=PatchToVec(md1.results.TransientSolution(t).Vy);
+                        thickness=PatchToVec(md1.results.TransientSolution(t).Thickness);
                         spcx=[spcx InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,vx,md2.mesh.x,md2.mesh.y)];
                         spcy=[spcy InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,vx,md2.mesh.x,md2.mesh.y)];
+                        spct=[spct InterpFromMeshToMesh2d(md1.mesh.elements,md1.mesh.x,md1.mesh.y,thickness,md2.mesh.x,md2.mesh.y)];
                         steps=[steps t*md1.timestepping.time_step];
                 end
@@ -115 +118 @@
                 md2.diagnostic.spcvy=md2.diagnostic.spcvy*ones(1,size(spcy,2));
                 md2.diagnostic.spcvz=md2.diagnostic.spcvz*ones(1,size(spcx,2));
+                md2.prognostic.spcthickness=md2.prognostic.spcthickness*ones(1,size(spct,2));
                 md2.diagnostic.spcvx(find(md2.mesh.vertexonboundary),:)=spcx(find(md2.mesh.vertexonboundary),:);
                 md2.diagnostic.spcvy(find(md2.mesh.vertexonboundary),:)=spcy(find(md2.mesh.vertexonboundary),:);
                 md2.diagnostic.spcvz(find(md2.mesh.vertexonboundary),:)=0;
+                md2.prognostic.spcthickness(find(md2.mesh.vertexonboundary),:)=spct(find(md2.mesh.vertexonboundary),:);
                 md2.diagnostic.spcvx=[md2.diagnostic.spcvx; steps];
                 md2.diagnostic.spcvy=[md2.diagnostic.spcvy; steps];
                 md2.diagnostic.spcvz=[md2.diagnostic.spcvz; steps];
+                md2.prognostic.spcthickness=[md2.prognostic.spcthickness; steps];
         end
 
@@ -143 +149 @@
                 md2.diagnostic.spcvy(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
                 md2.diagnostic.spcvz(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
+                md2.prognostic.spcthickness(unique(md2.diagnostic.icefront(:,1:2)),:)=nan;
         end