Changeset 17555 for issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
- Timestamp:
- 03/26/14 20:18:48 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
r17553 r17555 194 194 iomodel->FetchDataToInput(elements,ThicknessEnum); 195 195 iomodel->FetchDataToInput(elements,SurfaceEnum); 196 iomodel->FetchDataToInput(elements,B edEnum);196 iomodel->FetchDataToInput(elements,BaseEnum); 197 197 iomodel->FetchDataToInput(elements,FrictionCoefficientEnum); 198 198 iomodel->FetchDataToInput(elements,FrictionPEnum); … … 809 809 if(numrifts){ 810 810 iomodel->FetchData(&riftinfo,&numriftsegments,NULL,RiftsRiftstructEnum); 811 iomodel->FetchData(5,RiftsRiftstructEnum,ThicknessEnum,B edEnum,SurfaceEnum,MaskGroundediceLevelsetEnum);811 iomodel->FetchData(5,RiftsRiftstructEnum,ThicknessEnum,BaseEnum,SurfaceEnum,MaskGroundediceLevelsetEnum); 812 812 for(i=0;i<numriftsegments;i++){ 813 813 if(iomodel->my_elements[reCast<int,IssmDouble>(*(riftinfo+RIFTINFOSIZE*i+2))-1]){ … … 816 816 } 817 817 } 818 iomodel->DeleteData(5,RiftsRiftstructEnum,ThicknessEnum,B edEnum,SurfaceEnum,MaskGroundediceLevelsetEnum);818 iomodel->DeleteData(5,RiftsRiftstructEnum,ThicknessEnum,BaseEnum,SurfaceEnum,MaskGroundediceLevelsetEnum); 819 819 xDelete<IssmDouble>(riftinfo); 820 820 } … … 1486 1486 /*Retrieve all inputs and parameters*/ 1487 1487 Input* thickness_input = element->GetInput(ThicknessEnum); _assert_(thickness_input); 1488 Input* b ed_input = element->GetInput(BedEnum); _assert_(bed_input);1488 Input* base_input = element->GetInput(BaseEnum); _assert_(base_input); 1489 1489 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoWaterEnum); 1490 1490 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); … … 1501 1501 gauss->GaussPoint(ig); 1502 1502 thickness_input->GetInputValue(&thickness,gauss); 1503 b ed_input->GetInputValue(&bed,gauss);1503 base_input->GetInputValue(&bed,gauss); 1504 1504 element->JacobianDeterminantSurface(&Jdet,xyz_list_front,gauss); 1505 1505 element->NodalFunctions(basis,gauss); … … 1938 1938 /*Retrieve all inputs and parameters*/ 1939 1939 Input* thickness_input = element->GetInput(ThicknessEnum); _assert_(thickness_input); 1940 Input* b ed_input = element->GetInput(BedEnum); _assert_(bed_input);1940 Input* base_input = element->GetInput(BaseEnum); _assert_(base_input); 1941 1941 IssmDouble rho_water = element->GetMaterialParameter(MaterialsRhoWaterEnum); 1942 1942 IssmDouble rho_ice = element->GetMaterialParameter(MaterialsRhoIceEnum); … … 1952 1952 gauss->GaussPoint(ig); 1953 1953 thickness_input->GetInputValue(&thickness,gauss); 1954 b ed_input->GetInputValue(&bed,gauss);1954 base_input->GetInputValue(&bed,gauss); 1955 1955 element->JacobianDeterminantSurface(&Jdet,xyz_list_front,gauss); 1956 1956 element->NodalFunctions(basis,gauss); … … 3592 3592 /*Retrieve all inputs and parameters*/ 3593 3593 element->GetVerticesCoordinatesBase(&xyz_list_base); 3594 Input* b ed_input=element->GetInput(BedEnum); _assert_(bed_input);3594 Input* base_input=element->GetInput(BaseEnum); _assert_(base_input); 3595 3595 IssmDouble rho_water=element->GetMaterialParameter(MaterialsRhoWaterEnum); 3596 3596 IssmDouble gravity =element->GetMaterialParameter(ConstantsGEnum); … … 3606 3606 element->NormalBase(&normal[0],xyz_list_base); 3607 3607 _assert_(normal[dim-1]<0.); 3608 b ed_input->GetInputValue(&bed, gauss);3608 base_input->GetInputValue(&bed, gauss); 3609 3609 water_pressure=gravity*rho_water*bed; 3610 3610
Note:
See TracChangeset
for help on using the changeset viewer.