Changeset 176

Timestamp:

05/01/09 14:54:07 (16 years ago)

Author:

Mathieu Morlighem

Message:

Added accumulation and melting fields

File:

• issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp (modified) (3 diffs)

issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp

@@ -228 +228 @@
                 ModelFetchData((void**)&model->B,NULL,NULL,model_handle,"B","Matrix","Mat");
                 ModelFetchData((void**)&model->n,NULL,NULL,model_handle,"n","Matrix","Mat");
+                ModelFetchData((void**)&model->accumulation,NULL,NULL,model_handle,"accumulation","Matrix","Mat");
+                ModelFetchData((void**)&model->melting,NULL,NULL,model_handle,"melting","Matrix","Mat");
                 
                 for (i=0;i<model->numberofelements;i++){
@@ -269 +271 @@
                         tria_p=model->p[i];
                         tria_q=model->q[i];
+
+                        /*meling and accumulation*/
+                        tria_melting[0]=*(model->melting+        ((int)*(model->elements+elements_width*i+0)-1));
+                        tria_melting[1]=*(model->melting+        ((int)*(model->elements+elements_width*i+1)-1));
+                        tria_melting[2]=*(model->melting+        ((int)*(model->elements+elements_width*i+2)-1));
+
+                        tria_accumulation[0]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+0)-1));
+                        tria_accumulation[1]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+1)-1));
+                        tria_accumulation[2]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+2)-1));
 
                         /*element on iceshelf?:*/
@@ -368 +379 @@
                                 penta_b[j]=*(model->bed+    ((int)*(model->elements+elements_width*i+j)-1)); 
                                 penta_k[j]=*(model->drag+        ((int)*(model->elements+elements_width*i+j)-1)); 
+                                tria_melting[j]=*(model->melting+        ((int)*(model->elements+elements_width*i+j)-1));
+                                tria_accumulation[j]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+j)-1));
                         }