Changeset 3751 for issm/trunk/src/c/parallel/thermal_core.cpp
- Timestamp:
- 05/15/10 19:36:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/parallel/thermal_core.cpp
r3717 r3751 10 10 #include "../modules.h" 11 11 12 void thermal_core(DataSet* results,Model* model , ParameterInputs* inputs){12 void thermal_core(DataSet* results,Model* model){ 13 13 14 14 extern int my_rank; … … 65 65 66 66 if(verbose)_printf_("computing temperatures:\n"); 67 thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t, inputs,ThermalAnalysisEnum,NoneAnalysisEnum);68 inputs->Add("temperature",t_g[0],1,numberofnodes);69 inputs->Add("melting_offset",melting_offset);67 thermal_core_nonlinear(&t_g[0],&melting_offset,fem_t,ThermalAnalysisEnum,NoneAnalysisEnum); 68 model->UpdateInputsFromVector(t_g[0],TemperatureEnum,VertexEnum); 69 model->UpdateInputsFromConstant(melting_offset,MeltingOffsetEnum); 70 70 71 71 if(verbose)_printf_("computing melting:\n"); 72 diagnostic_core_linear(&m_g[0],fem_m, inputs,MeltingAnalysisEnum,NoneAnalysisEnum);72 diagnostic_core_linear(&m_g[0],fem_m,MeltingAnalysisEnum,NoneAnalysisEnum); 73 73 } 74 74 else{ … … 86 86 t_g[0]=SerialToVec(t_g_serial,numberofnodes); 87 87 xfree((void**)&t_g_serial); 88 fem_m->FindParam(&m_g_serial,NULL,NULL,Melting Enum);88 fem_m->FindParam(&m_g_serial,NULL,NULL,MeltingRateEnum); 89 89 m_g[0]=SerialToVec(m_g_serial,numberofnodes); 90 90 xfree((void**)&m_g_serial); … … 95 95 96 96 if(verbose)_printf_("computing temperatures:\n"); 97 inputs->Add("temperature",t_g[i],1,numberofnodes);98 thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t, inputs,ThermalAnalysisEnum,NoneAnalysisEnum);97 model->UpdateInputsFromVector(t_g[i],TemperatureEnum,VertexEnum); 98 thermal_core_nonlinear(&t_g[i+1],&melting_offset,fem_t,ThermalAnalysisEnum,NoneAnalysisEnum); 99 99 100 100 if(verbose)_printf_("computing melting:\n"); 101 inputs->Add("temperature",t_g[i+1],1,numberofnodes);102 inputs->Add("melting_offset",melting_offset);103 diagnostic_core_linear(&m_g[i+1],fem_m, inputs,MeltingAnalysisEnum,NoneAnalysisEnum);101 model->UpdateInputsFromVector(t_g[i+1],TemperatureEnum,VertexEnum); 102 model->UpdateInputsFromConstant(melting_offset,MeltingOffsetEnum); 103 diagnostic_core_linear(&m_g[i+1],fem_m,MeltingAnalysisEnum,NoneAnalysisEnum); 104 104 } 105 105 }
Note:
See TracChangeset
for help on using the changeset viewer.