Changeset 878

Timestamp:

06/10/09 10:21:39 (16 years ago)

Author:

seroussi

Message:

no update of B if no matice

File:

• issm/trunk/src/c/objects/Penta.cpp (modified) (1 diff)

issm/trunk/src/c/objects/Penta.cpp

@@ -1099 +1099 @@
         inputs->Recover("melting",&melting[0],1,dofs,6,(void**)nodes);
         inputs->Recover("accumulation_param",&accumulation[0],1,dofs,6,(void**)nodes);
-        
+
         //Update material if necessary
         if(inputs->Recover("temperature_average",&temperature_list[0],1,dofs,6,(void**)nodes)){
-                temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
-                B_average=Paterson(temperature_average);
-                matice->SetB(B_average);
+                if(matice){
+                        temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
+                        B_average=Paterson(temperature_average);
+                        matice->SetB(B_average);
+                }
         }
         
         if(inputs->Recover("B",&B_list[0],1,dofs,6,(void**)nodes)){
-                B_average=(B_list[0]+B_list[1]+B_list[2]+B_list[3]+B_list[4]+B_list[5])/6.0;
-                matice->SetB(B_average);
+                if(matice){
+                        B_average=(B_list[0]+B_list[1]+B_list[2]+B_list[3]+B_list[4]+B_list[5])/6.0;
+                        matice->SetB(B_average);
+                }
         }