Changeset 820

Timestamp:

06/05/09 15:00:52 (16 years ago)

Author:

seroussi

Message:

need a matpar in prognostic for transient solution

File:

• issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp (modified) (2 diffs)

issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp

@@ -339 +339 @@
                         elements->AddObject(penta);
         
-
                         #ifdef _PARALLEL_
                         /*Now that we are here, we can also start building the list of grids belonging to this node partition: we use 
@@ -396 +395 @@
         #endif
 
+        /*Add one constant material property to materials: */
+        matpar_mid=1; //put it at the end of the materials
+        matpar_g=model->g; 
+        matpar_rho_ice=model->rho_ice; 
+        matpar_rho_water=model->rho_water; 
+        matpar_thermalconductivity=model->thermalconductivity; 
+        matpar_heatcapacity=model->heatcapacity; 
+        matpar_latentheat=model->latentheat; 
+        matpar_beta=model->beta; 
+        matpar_meltingpoint=model->meltingpoint; 
+        matpar_mixed_layer_capacity=model->mixed_layer_capacity; 
+        matpar_thermal_exchange_velocity=model->thermal_exchange_velocity; 
+
+        /*Create matpar object using its constructor: */
+        matpar=new Matpar(matpar_mid,matpar_rho_ice,matpar_rho_water,matpar_heatcapacity,matpar_thermalconductivity,
+                        matpar_latentheat,matpar_beta,matpar_meltingpoint,matpar_mixed_layer_capacity,
+                        matpar_thermal_exchange_velocity,matpar_g);
+                
+        /*Add to materials datset: */
+        materials->AddObject(matpar);
+
         /*Partition penalties in 3d: */
         if(strcmp(model->meshtype,"3d")==0){