Changeset 17356
- Timestamp:
- 02/26/14 11:57:45 (11 years ago)
- Location:
- issm/trunk-jpl/src/c
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
r17355 r17356 191 191 } 192 192 } 193 193 194 /*Clean up and return*/ 194 195 xDelete<IssmDouble>(xyz_list); … … 293 294 return pe; 294 295 }/*}}}*/ 295 296 296 void HydrologyDCEfficientAnalysis::GetSolutionFromInputs(Vector<IssmDouble>* solution,Element* element){/*{{{*/ 297 297 element->GetSolutionFromInputsOneDof(solution,EplHeadEnum); 298 298 }/*}}}*/ 299 300 299 void HydrologyDCEfficientAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/ 301 300 … … 336 335 if(meshtype!=Mesh2DhorizontalEnum){basalelement->DeleteMaterials(); delete basalelement;}; 337 336 } /*}}}*/ 338 339 337 void HydrologyDCEfficientAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/ 340 338 /*Default, do nothing*/ … … 351 349 return rho_freshwater*g*epl_porosity*(water_compressibility+(epl_compressibility/epl_porosity)); 352 350 }/*}}}*/ 353 354 351 IssmDouble HydrologyDCEfficientAnalysis::SedimentStoring(Element* element){/*{{{*/ 355 352 IssmDouble rho_freshwater = element->GetMaterialParameter(MaterialsRhoFreshwaterEnum); … … 361 358 return rho_freshwater*g*sediment_porosity*sediment_thickness*(water_compressibility+(sediment_compressibility/sediment_porosity)); 362 359 }/*}}}*/ 363 364 360 IssmDouble HydrologyDCEfficientAnalysis::GetHydrologyDCInefficientHmax(Element* element, Gauss* gauss){/*{{{*/ 365 361 int hmax_flag; … … 440 436 return transfer; 441 437 }/*}}}*/ 442 443 438 IssmDouble HydrologyDCEfficientAnalysis::GetHydrologyPVectorTransfer(Element* element, Gauss* gauss){/*{{{*/ 444 439 … … 486 481 return transfer; 487 482 }/*}}}*/ 488 489 483 void HydrologyDCEfficientAnalysis::GetB(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/ 490 484 /*Compute B matrix. B=[B1 B2 B3] where Bi is of size 3*NDOF2. -
issm/trunk-jpl/src/c/solutionsequences/solutionsequence_hydro_nonlinear.cpp
r17353 r17356 177 177 ug_epl->Copy(ug_epl_sub_iter); 178 178 179 /* {{{ *//*Retriev the EPL head slopes and compute EPL Thickness*/179 /* {{{ *//*Retrieve the EPL head slopes and compute EPL Thickness*/ 180 180 if(VerboseSolution()) _printf0_("computing EPL Head slope...\n"); 181 181 femmodel->SetCurrentConfiguration(L2ProjectionEPLAnalysisEnum);
Note:
See TracChangeset
for help on using the changeset viewer.