Changeset 5184


Ignore:
Timestamp:
08/12/10 11:06:19 (15 years ago)
Author:
seroussi
Message:

added the possibility to spc hutter

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateConstraintsDiagnosticHutter.cpp

    r4236 r5184  
    3131
    3232        /*Fetch data: */
     33        IoModelFetchData(&iomodel->spcvelocity,NULL,NULL,iomodel_handle,"spcvelocity");
    3334        IoModelFetchData(&iomodel->gridonhutter,NULL,NULL,iomodel_handle,"gridonhutter");
    3435
     
    4849                                count++;
    4950                        }
     51                        else{
     52                                if ((int)iomodel->spcvelocity[6*i+0]){
     53                                        constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,1,*(iomodel->spcvelocity+6*i+3)/iomodel->yts,DiagnosticHutterAnalysisEnum)); //add count'th spc, on node i+1, setting dof 1 to vx.
     54                                        count++;
     55                                }
     56
     57                                if ((int)iomodel->spcvelocity[6*i+1]){
     58                                        constraints->AddObject(new Spc(iomodel->constraintcounter+count+1,iomodel->nodecounter+i+1,2,*(iomodel->spcvelocity+6*i+4)/iomodel->yts,DiagnosticHutterAnalysisEnum)); //add count'th spc, on node i+1, setting dof 2 to vy
     59                                        count++;
     60                                }
     61                        }
    5062                } //if((my_vertices[i]))
    5163        }
     
    5365        /*Free data: */
    5466        xfree((void**)&iomodel->gridonhutter);
     67        xfree((void**)&iomodel->spcvelocity);
    5568
    5669        cleanup_and_return:
Note: See TracChangeset for help on using the changeset viewer.