Index: ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp =================================================================== --- ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 16774) +++ ../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (revision 16775) @@ -1171,8 +1171,8 @@ for(i=0;iTransformSolutionCoord(FSvalues,cs_list); - element->TransformSolutionCoord(HOvalues,XYEnum); + element->TransformSolutionCoord(FSvalues,2*numnodes,cs_list); + element->TransformSolutionCoord(HOvalues,numnodes,XYEnum); /*Ok, we have vx and vy in values, fill in vx and vy arrays: */ for(i=0;iObjectEnum()==DoubleInputEnum) total_integrated_input=new DoubleInput(average_enum_type,0.0); else{ - printf("id %i \n",this->id); - printf("enum %i \n",average_enum_type); _error_("object " << EnumToStringx(original_input->ObjectEnum()) << " not supported yet"); } } @@ -3328,20 +3326,34 @@ return dt; }/*}}}*/ -/*FUNCTION Penta::TransformSolutionCoord{{{*/ +/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int transformenum){{{*/ void Penta::TransformSolutionCoord(IssmDouble* values,int transformenum){ ::TransformSolutionCoord(values,this->nodes,this->NumberofNodes(),transformenum); } /*}}}*/ -/*FUNCTION Penta::TransformSolutionCoord{{{*/ +/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int* transformenum_list){{{*/ void Penta::TransformSolutionCoord(IssmDouble* values,int* transformenum_list){ ::TransformSolutionCoord(values,this->nodes,this->NumberofNodes(),transformenum_list); } /*}}}*/ +/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){{{*/ +void Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){ + + ::TransformSolutionCoord(values,this->nodes,numnodes,transformenum); + +} +/*}}}*/ +/*FUNCTION Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){{{*/ +void Penta::TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){ + + ::TransformSolutionCoord(values,this->nodes,numnodes,transformenum_list); + +} +/*}}}*/ /*FUNCTION Penta::Update {{{*/ void Penta::Update(int index,IoModel* iomodel,int analysis_counter,int analysis_type,int finiteelement_type){ Index: ../trunk-jpl/src/c/classes/Elements/Penta.h =================================================================== --- ../trunk-jpl/src/c/classes/Elements/Penta.h (revision 16774) +++ ../trunk-jpl/src/c/classes/Elements/Penta.h (revision 16775) @@ -248,6 +248,8 @@ void SurfaceNormal(IssmDouble* surface_normal, IssmDouble xyz_list[3][3]); void TransformSolutionCoord(IssmDouble* values,int transformenum); void TransformSolutionCoord(IssmDouble* values,int* transformenum_list); + void TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum); /*Tiling only*/ + void TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list);/*Tiling only*/ #ifdef _HAVE_STRESSBALANCE_ ElementMatrix* CreateKMatrixCouplingSSAHO(void); Index: ../trunk-jpl/src/c/classes/Elements/Seg.h =================================================================== --- ../trunk-jpl/src/c/classes/Elements/Seg.h (revision 16774) +++ ../trunk-jpl/src/c/classes/Elements/Seg.h (revision 16775) @@ -153,6 +153,8 @@ IssmDouble TimeAdapt(){_error_("not implemented yet");}; void TransformSolutionCoord(IssmDouble* values,int transformenum){_error_("not implemented yet");}; void TransformSolutionCoord(IssmDouble* values,int* transformenum_list){_error_("not implemented yet");}; + void TransformSolutionCoord(IssmDouble* values,int numnodes,int transformenum){_error_("not implemented yet");}; + void TransformSolutionCoord(IssmDouble* values,int numnodes,int* transformenum_list){_error_("not implemented yet");}; void UpdateConstraintsExtrudeFromBase(){_error_("not implemented");}; void UpdateConstraintsExtrudeFromTop(){_error_("not implemented");}; void UpdateConstraintsL2ProjectionEPL(){_error_("not implemented");};