Changeset 9045

Timestamp:

07/18/11 10:51:43 (14 years ago)

Author:

seroussi

Message:

fixed spctemperature

Location:

issm/trunk/test/NightlyRun

Files:

• test1301.m (modified) (1 diff)
• test1302.m (modified) (1 diff)
• test1303.m (modified) (1 diff)
• test1304.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test1301.m

@@ -16 +16 @@
 md.observed_temperature=273.15*ones(md.numberofnodes,1);
 md.temperature=md.observed_temperature;
-md.spctemperature(:,2)=md.observed_temperature;
+pos=find(md.nodeonsurface);
+md.spctemperature(pos)=md.observed_temperature(pos);
 md.rheology_B=paterson(md.temperature);
 

issm/trunk/test/NightlyRun/test1302.m

@@ -12 +12 @@
 
 %Thermal boundary conditions
-pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;
-pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
+pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3))=10;
+pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
 md.vz=0.1*ones(md.numberofnodes,1);
 md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2);

issm/trunk/test/NightlyRun/test1303.m

@@ -11 +11 @@
 md=extrude(md,11,2);
 md=setelementstype(md,'Pattyn','all');
-pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3),1)=1; md.spctemperature(md.elements(pos1,1:3),2)=10;
-pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
+pos1=find(md.elementonbed);     md.spctemperature(md.elements(pos1,1:3))=10;
+pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
 md.pressure=zeros(md.numberofnodes,1);
 

issm/trunk/test/NightlyRun/test1304.m

@@ -12 +12 @@
 md=setelementstype(md,'Pattyn','all');
 
-pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6),1)=1; md.spctemperature(md.elements(pos2,4:6),2)=0;
+pos2=find(md.elementonsurface); md.spctemperature(md.elements(pos2,4:6))=0;
 md.pressure=zeros(md.numberofnodes,1);
 md.geothermalflux(:)=0.1; %100mW/m^2