Index: /issm/trunk/test/Par/79North.par
===================================================================
--- /issm/trunk/test/Par/79North.par	(revision 9424)
+++ /issm/trunk/test/Par/79North.par	(revision 9425)
@@ -11,6 +11,6 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
 
Index: /issm/trunk/test/Par/Pig.par
===================================================================
--- /issm/trunk/test/Par/Pig.par	(revision 9424)
+++ /issm/trunk/test/Par/Pig.par	(revision 9425)
@@ -15,8 +15,8 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
+md.temperature=md.temperature;
 
 %Friction
Index: /issm/trunk/test/Par/RoundSheetEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 9424)
+++ /issm/trunk/test/Par/RoundSheetEISMINT.par	(revision 9425)
@@ -15,6 +15,5 @@
 st=1.67*10^-2/1000; %k/m;
 radius=sqrt((md.x).^2+(md.y).^2);
-md.observed_temperature=(tmin+st*radius);
-md.temperature=md.observed_temperature;
+md.temperature=(tmin+st*radius);
 md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/RoundSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetShelf.par	(revision 9424)
+++ /issm/trunk/test/Par/RoundSheetShelf.par	(revision 9425)
@@ -40,8 +40,7 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
 
 %Surface mass balance and basal melting
Index: /issm/trunk/test/Par/RoundSheetStaticEISMINT.par
===================================================================
--- /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 9424)
+++ /issm/trunk/test/Par/RoundSheetStaticEISMINT.par	(revision 9425)
@@ -21,6 +21,5 @@
 tmin=238.15; %K
 st=1.67*10^-2/1000; %k/m;
-md.observed_temperature=(tmin+st*radius);
-md.temperature=md.observed_temperature;
+md.temperature=(tmin+st*radius);
 md.geothermalflux=4.2*10^-2*ones(md.numberofnodes,1);
 
Index: /issm/trunk/test/Par/SquareEISMINT.par
===================================================================
--- /issm/trunk/test/Par/SquareEISMINT.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareEISMINT.par	(revision 9425)
@@ -19,8 +19,7 @@
 
 disp('      creating temperature');
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
 
 disp('      creating flow law paramter');
-%md.B=paterson(md.observed_temperature);
 md.rheology_B=1.7687*10^8*ones(md.numberofnodes,1);
 md.rheology_n=3*ones(md.numberofelements,1);
Index: /issm/trunk/test/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 9425)
@@ -19,8 +19,7 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
 
 %Friction
Index: /issm/trunk/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetShelf.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareSheetShelf.par	(revision 9425)
@@ -22,8 +22,7 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
 
 %Accumulation and melting
Index: /issm/trunk/test/Par/SquareShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareShelf.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareShelf.par	(revision 9425)
@@ -19,8 +19,7 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
 
 %Friction
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 9425)
@@ -19,8 +19,7 @@
 
 %Materials
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
-md.rheology_B=paterson(md.observed_temperature);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
-md.temperature=md.observed_temperature;
 
 %Surface mass balance and basal melting
Index: /issm/trunk/test/Par/SquareThermal.par
===================================================================
--- /issm/trunk/test/Par/SquareThermal.par	(revision 9424)
+++ /issm/trunk/test/Par/SquareThermal.par	(revision 9425)
@@ -25,8 +25,8 @@
 
 disp('      creating temperatures');
-md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
+md.temperature=(273-20)*ones(md.numberofnodes,1);
 
 disp('      creating flow law paramter');
-md.rheology_B=paterson(md.observed_temperature);
+md.rheology_B=paterson(md.temperature);
 md.rheology_n=3*ones(md.numberofelements,1);
 
@@ -41,6 +41,5 @@
 
 disp('      boundary conditions for thermal model');
-md.temperature=md.observed_temperature;
-md.spctemperature(:)=md.observed_temperature;
+md.spctemperature(:)=md.temperature;
 md.geothermalflux=zeros(md.numberofnodes,1); 
 pos=find(md.elementonicesheet);md.geothermalflux(md.elements(pos,:))=1*10^-3; %1 mW/m^2