Ignore:
Timestamp:
08/11/20 19:43:31 (5 years ago)
Author:
Mathieu Morlighem
Message:

CHG: moving inputs2 back to inputs

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp

    r25317 r25379  
    107107        return 2;
    108108}/*}}}*/
    109 void StressbalanceSIAAnalysis::UpdateElements(Elements* elements,Inputs2* inputs2,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/
     109void StressbalanceSIAAnalysis::UpdateElements(Elements* elements,Inputs* inputs,IoModel* iomodel,int analysis_counter,int analysis_type){/*{{{*/
    110110
    111111        /*Fetch data needed: */
    … …  
    127127                if(iomodel->my_elements[i]){
    128128                        Element* element=(Element*)elements->GetObjectByOffset(counter);
    129                         element->Update(inputs2,i,iomodel,analysis_counter,analysis_type,P1Enum);
     129                        element->Update(inputs,i,iomodel,analysis_counter,analysis_type,P1Enum);
    130130                        /*Need to know the type of approximation for this element*/
    131131                        if(iomodel->Data("md.flowequation.element_equation")){
    132                                 inputs2->SetInput(ApproximationEnum,counter,IoCodeToEnumElementEquation(reCast<int>(iomodel->Data("md.flowequation.element_equation")[i])));
     132                                inputs->SetInput(ApproximationEnum,counter,IoCodeToEnumElementEquation(reCast<int>(iomodel->Data("md.flowequation.element_equation")[i])));
    133133                        }
    134134                        counter++;
    … …  
    142142        switch(frictionlaw){
    143143                case 1:
    144                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.coefficient",FrictionCoefficientEnum);
    145                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.p",FrictionPEnum);
    146                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.q",FrictionQEnum);
     144                        iomodel->FetchDataToInput(inputs,elements,"md.friction.coefficient",FrictionCoefficientEnum);
     145                        iomodel->FetchDataToInput(inputs,elements,"md.friction.p",FrictionPEnum);
     146                        iomodel->FetchDataToInput(inputs,elements,"md.friction.q",FrictionQEnum);
    147147                        break;
    148148                case 2:
    149                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.C",FrictionCEnum);
    150                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.m",FrictionMEnum);
     149                        iomodel->FetchDataToInput(inputs,elements,"md.friction.C",FrictionCEnum);
     150                        iomodel->FetchDataToInput(inputs,elements,"md.friction.m",FrictionMEnum);
    151151                        break;
    152152                case 6:
    153                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.C",FrictionCEnum);
    154                         iomodel->FetchDataToInput(inputs2,elements,"md.friction.m",FrictionMEnum);
    155                         iomodel->FetchDataToInput(inputs2,elements,"md.initialization.pressure",PressureEnum);
    156                         iomodel->FetchDataToInput(inputs2,elements,"md.initialization.temperature",TemperatureEnum);
     153                        iomodel->FetchDataToInput(inputs,elements,"md.friction.C",FrictionCEnum);
     154                        iomodel->FetchDataToInput(inputs,elements,"md.friction.m",FrictionMEnum);
     155                        iomodel->FetchDataToInput(inputs,elements,"md.initialization.pressure",PressureEnum);
     156                        iomodel->FetchDataToInput(inputs,elements,"md.initialization.temperature",TemperatureEnum);
    157157                        break;
    158158                default:
    … …  
    160160        }
    161161
    162         iomodel->FetchDataToInput(inputs2,elements,"md.geometry.thickness",ThicknessEnum);
    163         iomodel->FetchDataToInput(inputs2,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum);
     162        iomodel->FetchDataToInput(inputs,elements,"md.geometry.thickness",ThicknessEnum);
     163        iomodel->FetchDataToInput(inputs,elements,"md.mask.ocean_levelset",MaskOceanLevelsetEnum);
    164164        if(ismovingfront){
    165165                if(iomodel->domaintype!=Domain2DhorizontalEnum)
    166                         iomodel->FetchDataToInput(inputs2,elements,"md.mesh.vertexonbase",MeshVertexonbaseEnum); // required for updating active nodes
     166                        iomodel->FetchDataToInput(inputs,elements,"md.mesh.vertexonbase",MeshVertexonbaseEnum); // required for updating active nodes
    167167        }
    168168
    … …  
    326326        IssmDouble  gravity    = element->FindParam(ConstantsGEnum);
    327327        IssmDouble  B,n;
    328         Input2* B_input         = element->GetInput2(MaterialsRheologyBbarEnum);_assert_(B_input);
    329         Input2* n_input         = element->GetInput2(MaterialsRheologyNEnum);   _assert_(n_input);
    330         Input2* slopex_input    = element->GetInput2(SurfaceSlopeXEnum);        _assert_(slopex_input);
    331         Input2* slopey_input    = element->GetInput2(SurfaceSlopeYEnum);        _assert_(slopey_input);
    332         Input2* thickness_input = element->GetInput2(ThicknessEnum);            _assert_(thickness_input);
    333         Input2* surface_input   = element->GetInput2(SurfaceEnum);              _assert_(surface_input);
    334         Input2* drag_input      = NULL;
     328        Input* B_input         = element->GetInput(MaterialsRheologyBbarEnum);_assert_(B_input);
     329        Input* n_input         = element->GetInput(MaterialsRheologyNEnum);   _assert_(n_input);
     330        Input* slopex_input    = element->GetInput(SurfaceSlopeXEnum);        _assert_(slopex_input);
     331        Input* slopey_input    = element->GetInput(SurfaceSlopeYEnum);        _assert_(slopey_input);
     332        Input* thickness_input = element->GetInput(ThicknessEnum);            _assert_(thickness_input);
     333        Input* surface_input   = element->GetInput(SurfaceEnum);              _assert_(surface_input);
     334        Input* drag_input      = NULL;
    335335        if(frictionlaw!=5 && frictionlaw!=1){
    336                 drag_input = element->GetInput2(FrictionCoefficientEnum);  _assert_(drag_input);
     336                drag_input = element->GetInput(FrictionCoefficientEnum);  _assert_(drag_input);
    337337        }
    338338
    … …  
    416416        IssmDouble  gravity    = element->FindParam(ConstantsGEnum);
    417417        IssmDouble B,n;
    418         Input2* B_input         = element->GetInput2(MaterialsRheologyBEnum);   _assert_(B_input);
    419         Input2* n_input         = element->GetInput2(MaterialsRheologyNEnum);   _assert_(n_input);
    420         Input2* surface_input   = element->GetInput2(SurfaceEnum);              _assert_(surface_input);
    421         Input2* slopex_input    = element->GetInput2(SurfaceSlopeXEnum);        _assert_(slopex_input);
    422         Input2* slopey_input    = element->GetInput2(SurfaceSlopeYEnum);        _assert_(slopey_input);
    423         Input2* thickness_input = element->GetInput2(ThicknessEnum);            _assert_(thickness_input);
    424         Input2* drag_input      = NULL;
     418        Input* B_input         = element->GetInput(MaterialsRheologyBEnum);   _assert_(B_input);
     419        Input* n_input         = element->GetInput(MaterialsRheologyNEnum);   _assert_(n_input);
     420        Input* surface_input   = element->GetInput(SurfaceEnum);              _assert_(surface_input);
     421        Input* slopex_input    = element->GetInput(SurfaceSlopeXEnum);        _assert_(slopex_input);
     422        Input* slopey_input    = element->GetInput(SurfaceSlopeYEnum);        _assert_(slopey_input);
     423        Input* thickness_input = element->GetInput(ThicknessEnum);            _assert_(thickness_input);
     424        Input* drag_input      = NULL;
    425425        Friction* friction     = NULL;
    426426        if(frictionlaw!=5 && frictionlaw!=1){
    427                 drag_input = element->GetInput2(FrictionCoefficientEnum);  _assert_(drag_input);
     427                drag_input = element->GetInput(FrictionCoefficientEnum);  _assert_(drag_input);
    428428        }
    429429        else if(frictionlaw==5){
    … …  
    535535        /*Get dof list and inputs */
    536536        element->GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
    537         Input2* vx_input=element->GetInput2(VxEnum); _assert_(vx_input);
    538         Input2* vy_input=element->GetInput2(VyEnum); _assert_(vy_input);
     537        Input* vx_input=element->GetInput(VxEnum); _assert_(vx_input);
     538        Input* vy_input=element->GetInput(VyEnum); _assert_(vy_input);
    539539
    540540        /*Ok, we have the velocities in inputs, fill in solution */
    … …  
    622622
    623623        /*Add vx and vy as inputs to the tria element: */
    624         element->AddInput2(VxEnum,vx,P1Enum);
    625         element->AddInput2(VyEnum,vy,P1Enum);
    626         element->AddInput2(VelEnum,vel,P1Enum);
    627         element->AddInput2(PressureEnum,pressure,P1Enum);
     624        element->AddInput(VxEnum,vx,P1Enum);
     625        element->AddInput(VyEnum,vy,P1Enum);
     626        element->AddInput(VelEnum,vel,P1Enum);
     627        element->AddInput(PressureEnum,pressure,P1Enum);
    628628
    629629        /*Free ressources:*/
Note: See TracChangeset for help on using the changeset viewer.