Changeset 14102 for issm/trunk-jpl/test/NightlyRun/test1304.m

Timestamp:

12/05/12 14:35:11 (12 years ago)

Author:

jschierm

Message:

NEW: Working python tests 1301-1304 (plus cosmetic changes).

File:

• issm/trunk-jpl/test/NightlyRun/test1304.m (modified) (3 diffs)

issm/trunk-jpl/test/NightlyRun/test1304.m

@@ -10 +10 @@
 md=parameterize(md,'../Par/SquareThermal.par');
 md=extrude(md,11,1.);
-md=setflowequation(md,'Pattyn','all');
+md=setflowequation(md,'pattyn','all');
 
-pos2=find(md.mesh.elementonsurface); md.thermal.spctemperature(md.mesh.elements(pos2,4:6))=0;
+pos2=find(md.mesh.elementonsurface); md.thermal.spctemperature(md.mesh.elements(pos2,4:6))=0.;
 md.initialization.pressure=zeros(md.mesh.numberofvertices,1);
 md.basalforcings.geothermalflux(:)=0.1; %100mW/m^2
@@ -21 +21 @@
 md.initialization.temperature=-0.1/md.materials.thermalconductivity*(md.mesh.z-md.geometry.surface); %G=0.1 W/m2
 
-%modeled  results
+%modeled results
 md.cluster=generic('name',oshostname(),'np',2);
 md=solve(md,ThermalSolutionEnum());
@@ -27 +27 @@
 %plot results
 comp_temp=md.results.ThermalSolution.Temperature;
-relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100;
-relative(find(comp_temp==md.initialization.temperature))=0;
+relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100.;
+relative(find(comp_temp==md.initialization.temperature))=0.;
 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data',md.initialization.temperature,'view',3,...
         'title','Analytical temperature','view',3,'data',comp_temp-md.initialization.temperature,...
         'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
         'figposition','mathieu','FontSize#all',20)
-if printingflag, 
+if printingflag,
         set(gcf,'Color','w')
         printmodel('thermalgeothermalflux','png','margin','on','marginsize',25,'frame','off','resolution',0.7,'hardcopy','off');