Changeset 513

Timestamp:

05/20/09 08:08:53 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed test

File:

• issm/trunk/test/Validation/ThermalTests/Simpleconduction/runme.m (modified) (1 diff)

issm/trunk/test/Validation/ThermalTests/Simpleconduction/runme.m

@@ -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