Changeset 8947

Timestamp:

07/13/11 09:32:38 (14 years ago)

Author:

Eric.Larour

Message:

Pourave solution to a pourave problem

File:

• issm/trunk/src/c/objects/Materials/Matice.cpp (modified) (3 diffs)

issm/trunk/src/c/objects/Materials/Matice.cpp

@@ -574 +574 @@
         /*Intermediaries*/
         Element *element      = NULL;
+        Parameters* parameters= NULL;
+        int         dim;
 
         /*Recover element*/
@@ -591 +593 @@
                                         for (int i=0;i<3;i++) values[i]=vector[((Tria*)element)->nodes[i]->GetSidList()]; //use sid list, to index into serial oriented vector 
                                         this->inputs->AddInput(new TriaVertexInput(name,values));
+                                        /*Special case for rheology B in 2D: Pourave land for this solution{{{1*/
+                                        if(name==RheologyBEnum){
+                                                /*Are we in 2D?:*/
+                                                if(element->Enum()==TriaEnum){
+                                                        parameters=((Tria*)(element))->parameters;
+                                                }
+                                                else{
+                                                        parameters=((Penta*)(element))->parameters;
+                                                }
+                                                parameters->FindParam(&dim,DimEnum);
+                                                if(dim==2){
+                                                        /*Dupliacte rheology input: */
+                                                        this->inputs->AddInput(new TriaVertexInput(RheologyBbarEnum,values));
+                                                }
+                                        }
+                                        /*}}}*/
                                         return;
 
@@ -597 +615 @@
                 default: _error_("type %i (%s) not implemented yet",type,EnumToStringx(type));
         }
+
+
+
 }
 /*}}}*/