Index: /issm/trunk/src/m/classes/@bamggeom/bamggeom.m
===================================================================
--- /issm/trunk/src/m/classes/@bamggeom/bamggeom.m	(revision 4390)
+++ /issm/trunk/src/m/classes/@bamggeom/bamggeom.m	(revision 4391)
@@ -1,4 +1,4 @@
 function bg = bamggeom(varargin)
-%ICEFRONT - constructor for bamggeom object
+%BAMGGEOM - constructor for bamggeom object
 %
 %   Usage:
Index: /issm/trunk/src/m/classes/@bamgmesh/bamgmesh.m
===================================================================
--- /issm/trunk/src/m/classes/@bamgmesh/bamgmesh.m	(revision 4390)
+++ /issm/trunk/src/m/classes/@bamgmesh/bamgmesh.m	(revision 4391)
@@ -1,4 +1,4 @@
 function bm = bamgmesh(varargin)
-%ICEFRONT - constructor for bamgmesh object
+%BAMGMESH - constructor for bamgmesh object
 %
 %   Usage:
Index: /issm/trunk/src/m/classes/@bamgmesh/subsasgn.m
===================================================================
--- /issm/trunk/src/m/classes/@bamgmesh/subsasgn.m	(revision 4390)
+++ /issm/trunk/src/m/classes/@bamgmesh/subsasgn.m	(revision 4391)
@@ -1,3 +1,3 @@
-function icefront = subsasgn(icefront,index,val)
+function bamgmesh = subsasgn(bamgmesh,index,val)
 %SUBSASGN - handles indexed assignments to objects
 %
@@ -6,7 +6,7 @@
 %
 %   Usage:
-%      icefront = subsasgn(icefront,index,val)
+%      bamgmesh = subsasgn(bamgmesh,index,val)
 %
 %   See also SUBSREF
 
-icefront=builtin('subsasgn',icefront,index,val);
+bamgmesh=builtin('subsasgn',bamgmesh,index,val);
Index: /issm/trunk/src/m/classes/@bamgmesh/subsref.m
===================================================================
--- /issm/trunk/src/m/classes/@bamgmesh/subsref.m	(revision 4390)
+++ /issm/trunk/src/m/classes/@bamgmesh/subsref.m	(revision 4391)
@@ -1,10 +1,10 @@
-function icefront = subsref(icefront,index)
+function bamgmesh = subsref(bamgmesh,index)
 %SUBSREF - handles indexed references to objects
 %
 %   the first argument is the object and the second one a structure array
 %   Usage:
-%      icefront = subsref(icefront,index)
+%      bamgmesh = subsref(bamgmesh,index)
 % 
 %   See also SUBSASGN
 
-icefront=builtin('subsref',icefront,index);
+bamgmesh=builtin('subsref',bamgmesh,index);
Index: /issm/trunk/src/m/utils/Nightly/testsgetfields.m
===================================================================
--- /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4390)
+++ /issm/trunk/src/m/utils/Nightly/testsgetfields.m	(revision 4391)
@@ -26,8 +26,13 @@
 		vel=zeros(md.numberofgrids,1);
 		vel(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
-		j=zeros(md.numberofgrids,1);
-		j(md.results.DiagnosticSolution.J.index)=md.results.DiagnosticSolution.J.value;
+		j=md.results.DiagnosticSolution.J;
 		parameter=zeros(md.numberofgrids,1);
-		parameter(md.results.DiagnosticSolution.parameter.index)=md.results.DiagnosticSolution.parameter.value;
+		if strcmpi(md.control_type,'rheology_B'),
+			parameter(md.results.DiagnosticSolution.RheologyB.index)=md.results.DiagnosticSolution.RheologyB.value;
+		elseif strcmpi(md.control_type,'drag_coefficient'),
+			parameter(md.results.DiagnosticSolution.DragCoefficient.index)=md.results.DiagnosticSolution.DragCoefficient.value;
+		else
+			error('not implemented yet');
+		end
 		field_values={vel,j,parameter};
 	elseif qmu,