Changeset 15451


Ignore:
Timestamp:
07/07/13 09:39:55 (12 years ago)
Author:
Mathieu Morlighem
Message:

CHG: spc vz only for FS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp

    r15448 r15451  
    341341                                        count++;
    342342                                }
    343                                 if(!xIsNan<IssmDouble>(spcvx[v1]) && !xIsNan<IssmDouble>(spcvx[v2])){
     343                                if(!xIsNan<IssmDouble>(spcvy[v1]) && !xIsNan<IssmDouble>(spcvy[v2])){
    344344                                        constraints->AddObject(new SpcStatic(iomodel->constraintcounter+count+1,iomodel->nodecounter+iomodel->numberofvertices+i+1,
    345345                                                                        2,(spcvy[v1]+spcvy[v2])/yts,DiagnosticHorizAnalysisEnum));
    346346                                        count++;
    347347                                }
    348                                 if(!xIsNan<IssmDouble>(spcvx[v1]) && !xIsNan<IssmDouble>(spcvx[v2])){
    349                                         constraints->AddObject(new SpcStatic(iomodel->constraintcounter+count+1,iomodel->nodecounter+iomodel->numberofvertices+i+1,
    350                                                                         3,(spcvz[v1]+spcvz[v2])/yts,DiagnosticHorizAnalysisEnum));
    351                                         count++;
     348                                if (reCast<int,IssmDouble>(vertices_type[v1])==StokesApproximationEnum ||  (reCast<int,IssmDouble>(vertices_type[v1])==NoneApproximationEnum)){
     349                                        if(!xIsNan<IssmDouble>(spcvz[v1]) && !xIsNan<IssmDouble>(spcvz[v2])){
     350                                                constraints->AddObject(new SpcStatic(iomodel->constraintcounter+count+1,iomodel->nodecounter+iomodel->numberofvertices+i+1,
     351                                                                                3,(spcvz[v1]+spcvz[v2])/yts,DiagnosticHorizAnalysisEnum));
     352                                                count++;
     353                                        }
    352354                                }
    353355                        }
Note: See TracChangeset for help on using the changeset viewer.