Index: /issm/trunk/test/Validation/ThermalTests/Simpleconduction/runme.m
===================================================================
--- /issm/trunk/test/Validation/ThermalTests/Simpleconduction/runme.m	(revision 512)
+++ /issm/trunk/test/Validation/ThermalTests/Simpleconduction/runme.m	(revision 513)
@@ -12,15 +12,12 @@
 md=geography(md,'Shelf.exp','');
 md=parameterize(md,'Square.par');
-md=extrude(md,11,1);
+md=extrude(md,11,2);
 md=setelementstype(md,'Pattyn','all');
 thermalboundarycondition;
 
 %analytical results
-%d2T/dz2=0 T(bed)=10 T(surface)=0  => T=10/thickness*z-10*surface/thickness
+%d2T/dz2=0 T(bed)=10 T(surface)=0  => T=0*(z-bed)/thickness+10*(surface-z)/thickness
 %each layer of the 3d mesh must have a constant value
-Tlayer=10:-1:0;
-for layer=1:md.numlayers
-	md.observed_temperature((layer-1)*md.numberofgrids2d+1:layer*md.numberofgrids2d)=Tlayer(layer);
-end
+md.observed_temperature=10*(md.surface-md.z)./md.thickness;
 
 %modeled  results