Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 24020)
+++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 24021)
@@ -180,20 +180,16 @@
 			break;
 		case BasalforcingsIsmip6Enum:{
-			//iomodel->FetchDataToInput(elements,"md.basalforcings.tf",BasalforcingsIsmp6TfEnum);
 			iomodel->FetchDataToInput(elements,"md.basalforcings.basin_id",BasalforcingsIsmp6BasinIdEnum);
-			IssmDouble** array3d = NULL;
-			int* Ms = NULL;
-			int* Ns = NULL;
-			int K;
+			IssmDouble** array3d = NULL; int* Ms = NULL; int* Ns = NULL; int K;
 			iomodel->FetchData(&array3d,&Ms,&Ns,&K,"md.basalforcings.tf");
 			if(!array3d) _error_("md.basalforcings.tf not found in binary file");
 			for(int i=0;i<elements->Size();i++){
 				Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(i));
+				if(!element->IsOnBase()) continue;
 				for(int kk=0;kk<K;kk++){
 					element->DatasetInputAdd(BasalforcingsIsmp6TfEnum,array3d[kk],iomodel,Ms[kk],Ns[kk],1,BasalforcingsIsmp6TfEnum,7,kk);
 				}
 			}
-			xDelete<int>(Ms);
-			xDelete<int>(Ns);
+			xDelete<int>(Ms); xDelete<int>(Ns);
 			for(int i=0;i<K;i++) xDelete<IssmDouble>(array3d[i]);
 			xDelete<IssmDouble*>(array3d);
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 24020)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 24021)
@@ -364,5 +364,5 @@
 	if(!IsOnBase()){
 		//put zero
-		sigma_b->SetValue(id-1,0.0,INS_VAL);
+		//sigma_b->SetValue(id-1,0.0,INS_VAL);
 		return;
 	}
@@ -416,5 +416,5 @@
 
 	/*Add value to output*/
-	sigma_b->SetValue(id-1,value,INS_VAL);
+	//sigma_b->SetValue(id-1,value,INS_VAL);
 }
 /*}}}*/