Changeset 18727


Ignore:
Timestamp:
11/03/14 16:12:49 (10 years ago)
Author:
bdef
Message:

CHG: modification in mask handling

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r18726 r18727  
    18161816               
    18171817        HydrologyDCEfficientAnalysis* effanalysis =  new HydrologyDCEfficientAnalysis();
     1818        HydrologyDCInefficientAnalysis* inefanalysis =  new HydrologyDCInefficientAnalysis();
    18181819        /*Step 1: update mask, the mask might be extended by residual and/or using downstream sediment head*/
    18191820        mask=new Vector<IssmDouble>(this->nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
     
    18341835        InputUpdateFromVectorx(this,serial_mask,HydrologydcMaskEplactiveNodeEnum,NodeSIdEnum);
    18351836        xDelete<IssmDouble>(serial_mask);
    1836 
     1837        inefanalysis->ElementizeEplMask(this);
    18371838        /*Step 2: update node activity. If one element is connected to mask=1, all nodes are active*/
    18381839        active=new Vector<IssmDouble>(nodes->NumberOfNodes(HydrologyDCEfficientAnalysisEnum));
     
    18411842                effanalysis->HydrologyEPLGetActive(active,element);
    18421843        }
     1844
    18431845
    18441846        /*Assemble and serialize*/
     
    18541856                        if(serial_active[node->Sid()]==1.){
    18551857                                node->Activate();
    1856                                 counter++;
     1858                                if(!node->IsClone());counter++;
    18571859                        }
    18581860                        else{
     
    19001902                        if(serial_active[node->Sid()]==1.){
    19011903                                node->Activate();
    1902                                 counter++;
     1904                                if(!node->IsClone());counter++;
    19031905                        }
    19041906                        else{
Note: See TracChangeset for help on using the changeset viewer.