Changeset 15986
- Timestamp:
- 08/27/13 16:12:21 (12 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 4 deleted
- 38 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r15982 r15986 5318 5318 /*update meltrate and watercolumn*/ 5319 5319 this->inputs->AddInput(new PentaInput(WatercolumnEnum, watercolumn, P1Enum)); 5320 this->inputs->AddInput(new PentaInput(Basal MeltrateEnum, meltrate, P1Enum));5320 this->inputs->AddInput(new PentaInput(BasalforcingsMeltingRateEnum, meltrate, P1Enum)); 5321 5321 5322 5322 /*Clean up and return*/ … … 11287 11287 this->inputs->AddInput(new PentaInput(SurfaceEnum,&s[0],P1Enum)); 11288 11288 this->inputs->AddInput(new PentaInput(BedEnum,&b[0],P1Enum)); 11289 this->inputs->AddInput(new BoolInput(MaskElementonfloatingiceEnum,floatingelement));11290 11289 11291 11290 /*Extrude inputs*/ 11292 11291 this->InputExtrude(SurfaceEnum,ElementEnum); 11293 11292 this->InputExtrude(BedEnum,ElementEnum); 11294 this->InputExtrude(MaskElementonfloatingiceEnum,ElementEnum);11295 11293 } 11296 11294 /*}}}*/ -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r15982 r15986 7402 7402 7403 7403 /*Update inputs*/ 7404 this->inputs->AddInput(new BoolInput(MaskElementonfloatingiceEnum,floatingelement));7405 7404 this->inputs->AddInput(new TriaInput(SurfaceEnum,&s[0],P1Enum)); 7406 7405 this->inputs->AddInput(new TriaInput(BedEnum,&b[0],P1Enum)); -
issm/trunk-jpl/src/c/classes/Loads/Riftfront.cpp
r15771 r15986 96 96 riftfront_friction=*(iomodel->Data(RiftsRiftstructEnum)+RIFTINFOSIZE*i+8); 97 97 riftfront_fractionincrement=*(iomodel->Data(RiftsRiftstructEnum)+RIFTINFOSIZE*i+10); 98 riftfront_shelf=reCast<bool,IssmDouble>(iomodel->Data(Mask VertexonfloatingiceEnum)[node1-1]);98 riftfront_shelf=reCast<bool,IssmDouble>(iomodel->Data(MaskGroundediceLevelsetEnum)[node1-1]<0.); 99 99 100 100 this->inputs->AddInput(new IntInput(RiftfrontTypeEnum,riftfront_type)); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateNodesBalancethickness.cpp
r15749 r15986 19 19 20 20 /*First fetch data: */ 21 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);21 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 22 22 if(stabilization!=3){ 23 23 CreateNodes(pnodes,iomodel,BalancethicknessAnalysisEnum,P1Enum); … … 26 26 CreateNodes(pnodes,iomodel,BalancethicknessAnalysisEnum,P1DGEnum); 27 27 } 28 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);28 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 29 29 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/UpdateElementsBalancethickness.cpp
r15941 r15986 36 36 iomodel->FetchDataToInput(elements,BedEnum); 37 37 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 38 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);39 38 iomodel->FetchDataToInput(elements,VxEnum); 40 39 iomodel->FetchDataToInput(elements,VyEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateNodesBedSlope.cpp
r15749 r15986 11 11 void CreateNodesBedSlope(Nodes** pnodes, IoModel* iomodel){ 12 12 13 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);13 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 14 14 CreateNodes(pnodes,iomodel,BedSlopeAnalysisEnum,P1Enum); 15 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);15 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 16 16 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateNodesEnthalpy.cpp
r15749 r15986 11 11 void CreateNodesEnthalpy(Nodes** pnodes, IoModel* iomodel){ 12 12 13 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);13 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 14 14 CreateNodes(pnodes,iomodel,EnthalpyAnalysisEnum,P1Enum); 15 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);15 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 16 16 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/UpdateElementsEnthalpy.cpp
r15979 r15986 35 35 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 36 36 iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum); 37 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);38 37 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 39 38 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceBase/CreateNodesFreeSurfaceBase.cpp
r15877 r15986 12 12 13 13 /*Create Nodes either DG or CG depending on stabilization*/ 14 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);14 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 15 15 CreateNodes(pnodes,iomodel,FreeSurfaceBaseAnalysisEnum,P1Enum); 16 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);16 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 17 17 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceBase/UpdateElementsFreeSurfaceBase.cpp
r15941 r15986 32 32 iomodel->FetchDataToInput(elements,SurfaceEnum); 33 33 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 34 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);35 34 iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum); 36 35 iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceTop/CreateNodesFreeSurfaceTop.cpp
r15877 r15986 12 12 13 13 /*Create Nodes either DG or CG depending on stabilization*/ 14 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);14 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 15 15 CreateNodes(pnodes,iomodel,FreeSurfaceTopAnalysisEnum,P1Enum); 16 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);16 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 17 17 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/FreeSurfaceTop/UpdateElementsFreeSurfaceTop.cpp
r15941 r15986 31 31 iomodel->FetchDataToInput(elements,SurfaceEnum); 32 32 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 33 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);34 33 iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum); 35 34 iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/CreateNodesHydrologyDCEfficient.cpp
r15749 r15986 21 21 if(!isefficientlayer) return; 22 22 23 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);23 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 24 24 CreateNodes(pnodes,iomodel,HydrologyDCEfficientAnalysisEnum,P1Enum); 25 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);25 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 26 26 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCEfficient/UpdateElementsHydrologyDCEfficient.cpp
r15941 r15986 36 36 iomodel->FetchDataToInput(elements,BedEnum); 37 37 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 38 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);39 38 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 40 39 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/CreateNodesHydrologyDCInefficient.cpp
r15749 r15986 18 18 if(hydrology_model!=HydrologydcEnum) return; 19 19 20 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);20 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 21 21 CreateNodes(pnodes,iomodel,HydrologyDCInefficientAnalysisEnum,P1Enum); 22 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);22 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 23 23 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyDCInefficient/UpdateElementsHydrologyDCInefficient.cpp
r15941 r15986 37 37 iomodel->FetchDataToInput(elements,BedEnum); 38 38 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 39 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);40 39 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 41 40 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/CreateNodesHydrologyShreve.cpp
r15749 r15986 18 18 if(hydrology_model!=HydrologyshreveEnum) return; 19 19 20 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);20 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 21 21 CreateNodes(pnodes,iomodel,HydrologyShreveAnalysisEnum,P1Enum); 22 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);22 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 23 23 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/HydrologyShreve/UpdateElementsHydrologyShreve.cpp
r15980 r15986 36 36 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 37 37 iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum); 38 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);39 38 iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum); 40 39 iomodel->FetchDataToInput(elements,WatercolumnEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/CreateNodesMasstransport.cpp
r15767 r15986 19 19 20 20 /*Create Nodes either DG or CG depending on stabilization*/ 21 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);21 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 22 22 if(stabilization!=3){ 23 23 CreateNodes(pnodes,iomodel,MasstransportAnalysisEnum,P1Enum); … … 26 26 CreateNodes(pnodes,iomodel,MasstransportAnalysisEnum,P1DGEnum); 27 27 } 28 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);28 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 29 29 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Masstransport/UpdateElementsMasstransport.cpp
r15952 r15986 46 46 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 47 47 iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum); 48 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);49 48 iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateEnum); 50 49 iomodel->FetchDataToInput(elements,BasalforcingsMeltingRateCorrectionEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateNodesMelting.cpp
r15749 r15986 11 11 void CreateNodesMelting(Nodes** pnodes, IoModel* iomodel){ 12 12 13 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);13 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 14 14 CreateNodes(pnodes,iomodel,MeltingAnalysisEnum,P1Enum); 15 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);15 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 16 16 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/UpdateElementsMelting.cpp
r15941 r15986 32 32 iomodel->FetchDataToInput(elements,FrictionQEnum); 33 33 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 34 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);35 34 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 36 35 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp
r15925 r15986 29 29 IssmDouble *nodeonFS = NULL; 30 30 IssmDouble *nodeonbed = NULL; 31 IssmDouble * nodeonicesheet= NULL;31 IssmDouble *groundedice_ls = NULL; 32 32 IssmDouble *vertices_type = NULL; 33 33 IssmDouble *surface = NULL; … … 96 96 iomodel->FetchData(&nodeonFS,NULL,NULL,FlowequationBorderFSEnum); 97 97 iomodel->FetchData(&nodeonbed,NULL,NULL,MeshVertexonbedEnum); 98 iomodel->FetchData(& nodeonicesheet,NULL,NULL,MaskVertexongroundediceEnum);98 iomodel->FetchData(&groundedice_ls,NULL,NULL,MaskGroundediceLevelsetEnum); 99 99 for(i=0;i<iomodel->numberofvertices;i++){ 100 100 if(iomodel->my_vertices[i]){ 101 if( reCast<int,IssmDouble>(nodeonbed[i]) && reCast<int,IssmDouble>(nodeonicesheet[i]) && reCast<int,IssmDouble>(nodeonFS[i])){101 if(nodeonbed[i]>0. && groundedice_ls[i]>0. && nodeonFS[i]>0.){ 102 102 if(vertices_type[i] == FSApproximationEnum){ 103 103 for(j=0;j<Nz;j++) spcvz[i*Nz+j] = 0.; … … 114 114 iomodel->DeleteData(nodeonFS,FlowequationBorderFSEnum); 115 115 iomodel->DeleteData(nodeonbed,MeshVertexonbedEnum); 116 iomodel->DeleteData( nodeonicesheet,MaskVertexongroundediceEnum);116 iomodel->DeleteData(groundedice_ls,MaskGroundediceLevelsetEnum); 117 117 118 118 /*Pressure spc*/ … … 125 125 for(i=0;i<iomodel->numberofvertices;i++){ 126 126 if(iomodel->my_vertices[i]){ 127 if 127 if(reCast<int,IssmDouble>(vertices_type[i])==NoneApproximationEnum){ 128 128 constraints->AddObject(new SpcStatic(count+1,iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofelements+i+1,1,g*rho_ice*(surface[i]-z[i])/FSreconditioning,StressbalanceAnalysisEnum)); 129 129 count++; … … 135 135 for(i=0;i<iomodel->numberofvertices;i++){ 136 136 if(iomodel->my_vertices[i]){ 137 if 137 if(reCast<int,IssmDouble>(vertices_type[i])==NoneApproximationEnum){ 138 138 constraints->AddObject(new SpcStatic(count+1,iomodel->nodecounter+iomodel->numberofvertices+iomodel->numberofedges+i+1,1,g*rho_ice*(surface[i]-z[i])/FSreconditioning,StressbalanceAnalysisEnum)); 139 139 count++; … … 162 162 if(iomodel->dim==3)iomodel->FetchData(&nodeonFS,NULL,NULL,FlowequationBorderFSEnum); 163 163 if(iomodel->dim==3)iomodel->FetchData(&nodeonbed,NULL,NULL,MeshVertexonbedEnum); 164 if(iomodel->dim==3)iomodel->FetchData(& nodeonicesheet,NULL,NULL,MaskVertexongroundediceEnum);164 if(iomodel->dim==3)iomodel->FetchData(&groundedice_ls,NULL,NULL,MaskGroundediceLevelsetEnum); 165 165 iomodel->FetchData(&vertices_type,NULL,NULL,FlowequationVertexEquationEnum); 166 166 iomodel->FetchData(&surface,NULL,NULL,SurfaceEnum); … … 383 383 384 384 /*Constraint at the bedrock interface (v.n = vz = 0) (Coordinates will be updated according to the bed slope)*/ 385 if (iomodel->dim==3) if( reCast<int,IssmDouble>(nodeonbed[i]) && reCast<int,IssmDouble>(nodeonicesheet[i]) && reCast<int,IssmDouble>(nodeonFS[i])){385 if (iomodel->dim==3) if(nodeonbed[i]>0. && groundedice_ls[i]>=0. && nodeonFS[i]>0.){ 386 386 switch(reCast<int,IssmDouble>(vertices_type[i])){ 387 387 case SSAFSApproximationEnum: … … 411 411 if(iomodel->dim==3)iomodel->DeleteData(nodeonFS,FlowequationBorderFSEnum); 412 412 if(iomodel->dim==3)iomodel->DeleteData(nodeonbed,MeshVertexonbedEnum); 413 if(iomodel->dim==3)iomodel->DeleteData( nodeonicesheet,MaskVertexongroundediceEnum);413 if(iomodel->dim==3)iomodel->DeleteData(groundedice_ls,MaskGroundediceLevelsetEnum); 414 414 iomodel->DeleteData(vertices_type,FlowequationVertexEquationEnum); 415 415 iomodel->DeleteData(surface,SurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateLoadsStressbalance.cpp
r15843 r15986 74 74 if(numrifts){ 75 75 iomodel->FetchData(&riftinfo,&numriftsegments,NULL,RiftsRiftstructEnum); 76 iomodel->FetchData(5,RiftsRiftstructEnum,ThicknessEnum,BedEnum,SurfaceEnum,Mask VertexonfloatingiceEnum);76 iomodel->FetchData(5,RiftsRiftstructEnum,ThicknessEnum,BedEnum,SurfaceEnum,MaskGroundediceLevelsetEnum); 77 77 for(i=0;i<numriftsegments;i++){ 78 78 if(iomodel->my_elements[reCast<int,IssmDouble>(*(riftinfo+RIFTINFOSIZE*i+2))-1]){ … … 81 81 } 82 82 } 83 iomodel->DeleteData(5,RiftsRiftstructEnum,ThicknessEnum,BedEnum,SurfaceEnum,Mask VertexonfloatingiceEnum);83 iomodel->DeleteData(5,RiftsRiftstructEnum,ThicknessEnum,BedEnum,SurfaceEnum,MaskGroundediceLevelsetEnum); 84 84 xDelete<IssmDouble>(riftinfo); 85 85 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateNodesStressbalance.cpp
r15923 r15986 50 50 iomodel->Constant(&finiteelement,FlowequationFeFSEnum); 51 51 } 52 iomodel->FetchData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,53 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);52 iomodel->FetchData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 53 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 54 54 CreateNodes(pnodes,iomodel,StressbalanceAnalysisEnum,finiteelement,approximation); 55 iomodel->DeleteData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,56 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);55 iomodel->DeleteData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 56 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 57 57 } 58 58 else{ … … 64 64 if(!nodes) nodes = new Nodes(); 65 65 66 iomodel->FetchData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,67 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);66 iomodel->FetchData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 67 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 68 68 if(isFS){ 69 69 /*P1+ velocity*/ … … 101 101 } 102 102 } 103 iomodel->DeleteData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,104 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);103 iomodel->DeleteData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 104 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 105 105 106 106 /*Assign output pointer: */ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp
r15952 r15986 94 94 iomodel->FetchDataToInput(elements,FrictionPEnum); 95 95 iomodel->FetchDataToInput(elements,FrictionQEnum); 96 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);97 96 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 98 97 iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/StressbalanceSIA/CreateNodesStressbalanceSIA.cpp
r15771 r15986 25 25 if(!nodes) nodes = new Nodes(); 26 26 27 iomodel->FetchData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,28 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);27 iomodel->FetchData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 28 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 29 29 30 30 for(int i=0;i<iomodel->numberofvertices;i++){ … … 44 44 } 45 45 46 iomodel->DeleteData( 8,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum,47 MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,StressbalanceReferentialEnum);46 iomodel->DeleteData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationBorderSSAEnum,FlowequationBorderFSEnum, 47 FlowequationVertexEquationEnum,StressbalanceReferentialEnum); 48 48 49 49 /*Assign output pointer: */ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/StressbalanceVertical/CreateNodesStressbalanceVertical.cpp
r15771 r15986 14 14 if(iomodel->dim==2) return; 15 15 16 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);16 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 17 17 CreateNodes(pnodes,iomodel,StressbalanceVerticalAnalysisEnum,P1Enum); 18 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);18 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 19 19 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/StressbalanceVertical/UpdateElementsStressbalanceVertical.cpp
r15941 r15986 29 29 iomodel->FetchDataToInput(elements,BedEnum); 30 30 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 31 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);32 31 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 33 32 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateNodesSurfaceSlope.cpp
r15749 r15986 11 11 void CreateNodesSurfaceSlope(Nodes** pnodes, IoModel* iomodel){ 12 12 13 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);13 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 14 14 CreateNodes(pnodes,iomodel,SurfaceSlopeAnalysisEnum,P1Enum); 15 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);15 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 16 16 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateNodesThermal.cpp
r15749 r15986 11 11 void CreateNodesThermal(Nodes** pnodes, IoModel* iomodel){ 12 12 13 iomodel->FetchData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);13 iomodel->FetchData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 14 14 CreateNodes(pnodes,iomodel,ThermalAnalysisEnum,P1Enum); 15 iomodel->DeleteData( 5,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum);15 iomodel->DeleteData(3,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,FlowequationVertexEquationEnum); 16 16 } -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/UpdateElementsThermal.cpp
r15979 r15986 32 32 iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum); 33 33 iomodel->FetchDataToInput(elements,MaskGroundediceLevelsetEnum); 34 iomodel->FetchDataToInput(elements,MaskElementonfloatingiceEnum);35 34 iomodel->FetchDataToInput(elements,MeshElementonbedEnum); 36 35 iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum); -
issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h
r15941 r15986 142 142 InversionVyObsEnum, 143 143 InversionVzObsEnum, 144 MaskElementonfloatingiceEnum,145 MaskElementongroundediceEnum,146 MaskVertexonfloatingiceEnum,147 MaskVertexongroundediceEnum,148 144 MaskIceLevelsetEnum, 149 145 MaterialsBetaEnum, … … 484 480 WatercolumnEnum, 485 481 BasalFrictionEnum, 486 BasalMeltrateEnum,487 482 ViscousHeatingEnum, 488 483 QmuTemperatureEnum, -
issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp
r15941 r15986 150 150 case InversionVyObsEnum : return "InversionVyObs"; 151 151 case InversionVzObsEnum : return "InversionVzObs"; 152 case MaskElementonfloatingiceEnum : return "MaskElementonfloatingice";153 case MaskElementongroundediceEnum : return "MaskElementongroundedice";154 case MaskVertexonfloatingiceEnum : return "MaskVertexonfloatingice";155 case MaskVertexongroundediceEnum : return "MaskVertexongroundedice";156 152 case MaskIceLevelsetEnum : return "MaskIceLevelset"; 157 153 case MaterialsBetaEnum : return "MaterialsBeta"; … … 478 474 case WatercolumnEnum : return "Watercolumn"; 479 475 case BasalFrictionEnum : return "BasalFriction"; 480 case BasalMeltrateEnum : return "BasalMeltrate";481 476 case ViscousHeatingEnum : return "ViscousHeating"; 482 477 case QmuTemperatureEnum : return "QmuTemperature"; -
issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp
r15941 r15986 153 153 else if (strcmp(name,"InversionVyObs")==0) return InversionVyObsEnum; 154 154 else if (strcmp(name,"InversionVzObs")==0) return InversionVzObsEnum; 155 else if (strcmp(name,"MaskElementonfloatingice")==0) return MaskElementonfloatingiceEnum;156 else if (strcmp(name,"MaskElementongroundedice")==0) return MaskElementongroundediceEnum;157 else if (strcmp(name,"MaskVertexonfloatingice")==0) return MaskVertexonfloatingiceEnum;158 else if (strcmp(name,"MaskVertexongroundedice")==0) return MaskVertexongroundediceEnum;159 155 else if (strcmp(name,"MaskIceLevelset")==0) return MaskIceLevelsetEnum; 160 156 else if (strcmp(name,"MaterialsBeta")==0) return MaterialsBetaEnum; … … 260 256 else if (strcmp(name,"GiaMantleViscosity")==0) return GiaMantleViscosityEnum; 261 257 else if (strcmp(name,"GiaLithosphereThickness")==0) return GiaLithosphereThicknessEnum; 258 else if (strcmp(name,"Thickness")==0) return ThicknessEnum; 259 else if (strcmp(name,"TimesteppingStartTime")==0) return TimesteppingStartTimeEnum; 260 else if (strcmp(name,"TimesteppingFinalTime")==0) return TimesteppingFinalTimeEnum; 261 else if (strcmp(name,"TimesteppingCflCoefficient")==0) return TimesteppingCflCoefficientEnum; 262 262 else stage=3; 263 263 } 264 264 if(stage==3){ 265 if (strcmp(name,"Thickness")==0) return ThicknessEnum; 266 else if (strcmp(name,"TimesteppingStartTime")==0) return TimesteppingStartTimeEnum; 267 else if (strcmp(name,"TimesteppingFinalTime")==0) return TimesteppingFinalTimeEnum; 268 else if (strcmp(name,"TimesteppingCflCoefficient")==0) return TimesteppingCflCoefficientEnum; 269 else if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum; 265 if (strcmp(name,"TimesteppingTimeAdapt")==0) return TimesteppingTimeAdaptEnum; 270 266 else if (strcmp(name,"TimesteppingTimeStep")==0) return TimesteppingTimeStepEnum; 271 267 else if (strcmp(name,"TransientIsstressbalance")==0) return TransientIsstressbalanceEnum; … … 383 379 else if (strcmp(name,"Profiler")==0) return ProfilerEnum; 384 380 else if (strcmp(name,"MatrixParam")==0) return MatrixParamEnum; 381 else if (strcmp(name,"NodeSId")==0) return NodeSIdEnum; 382 else if (strcmp(name,"VectorParam")==0) return VectorParamEnum; 383 else if (strcmp(name,"Riftfront")==0) return RiftfrontEnum; 384 else if (strcmp(name,"RiftfrontType")==0) return RiftfrontTypeEnum; 385 385 else stage=4; 386 386 } 387 387 if(stage==4){ 388 if (strcmp(name,"NodeSId")==0) return NodeSIdEnum; 389 else if (strcmp(name,"VectorParam")==0) return VectorParamEnum; 390 else if (strcmp(name,"Riftfront")==0) return RiftfrontEnum; 391 else if (strcmp(name,"RiftfrontType")==0) return RiftfrontTypeEnum; 392 else if (strcmp(name,"Segment")==0) return SegmentEnum; 388 if (strcmp(name,"Segment")==0) return SegmentEnum; 393 389 else if (strcmp(name,"SegmentRiftfront")==0) return SegmentRiftfrontEnum; 394 390 else if (strcmp(name,"SpcDynamic")==0) return SpcDynamicEnum; … … 487 483 else if (strcmp(name,"Watercolumn")==0) return WatercolumnEnum; 488 484 else if (strcmp(name,"BasalFriction")==0) return BasalFrictionEnum; 489 else if (strcmp(name,"BasalMeltrate")==0) return BasalMeltrateEnum;490 485 else if (strcmp(name,"ViscousHeating")==0) return ViscousHeatingEnum; 491 486 else if (strcmp(name,"QmuTemperature")==0) return QmuTemperatureEnum; … … 506 501 else if (strcmp(name,"P1DG")==0) return P1DGEnum; 507 502 else if (strcmp(name,"P1bubble")==0) return P1bubbleEnum; 508 else stage=5; 509 } 510 if(stage==5){ 511 if (strcmp(name,"P1bubblecondensed")==0) return P1bubblecondensedEnum; 503 else if (strcmp(name,"P1bubblecondensed")==0) return P1bubblecondensedEnum; 512 504 else if (strcmp(name,"P2")==0) return P2Enum; 513 505 else if (strcmp(name,"P2xP1")==0) return P2xP1Enum; 514 506 else if (strcmp(name,"P1xP2")==0) return P1xP2Enum; 515 507 else if (strcmp(name,"P1P1")==0) return P1P1Enum; 516 else if (strcmp(name,"P1P1GLS")==0) return P1P1GLSEnum; 508 else stage=5; 509 } 510 if(stage==5){ 511 if (strcmp(name,"P1P1GLS")==0) return P1P1GLSEnum; 517 512 else if (strcmp(name,"MINI")==0) return MINIEnum; 518 513 else if (strcmp(name,"MINIcondensed")==0) return MINIcondensedEnum; -
issm/trunk-jpl/src/m/classes/mask.m
r15951 r15986 52 52 end % }}} 53 53 function marshall(obj,md,fid) % {{{ 54 WriteData(fid,'object',obj,'fieldname','elementonfloatingice','format','BooleanMat','mattype',2);55 WriteData(fid,'object',obj,'fieldname','elementongroundedice','format','BooleanMat','mattype',2);56 WriteData(fid,'object',obj,'fieldname','vertexonfloatingice','format','DoubleMat','mattype',1);57 WriteData(fid,'object',obj,'fieldname','vertexongroundedice','format','DoubleMat','mattype',1);58 54 WriteData(fid,'object',obj,'fieldname','groundedice_levelset','format','DoubleMat','mattype',1); 59 55 WriteData(fid,'object',obj,'fieldname','ice_levelset','format','DoubleMat','mattype',1); -
issm/trunk-jpl/src/m/classes/mask.py
r15961 r15986 54 54 # }}} 55 55 def marshall(self,md,fid): # {{{ 56 WriteData(fid,'object',self,'fieldname','elementonfloatingice','format','BooleanMat','mattype',2)57 WriteData(fid,'object',self,'fieldname','elementongroundedice','format','BooleanMat','mattype',2)58 WriteData(fid,'object',self,'fieldname','vertexonfloatingice','format','DoubleMat','mattype',1)59 WriteData(fid,'object',self,'fieldname','vertexongroundedice','format','DoubleMat','mattype',1)60 56 WriteData(fid,'object',self,'fieldname','groundedice_levelset','format','DoubleMat','mattype',1) 61 57 WriteData(fid,'object',self,'fieldname','ice_levelset','format','DoubleMat','mattype',1) -
issm/trunk-jpl/src/m/enum/EnumDefinitions.py
r15941 r15986 142 142 def InversionVyObsEnum(): return StringToEnum("InversionVyObs")[0] 143 143 def InversionVzObsEnum(): return StringToEnum("InversionVzObs")[0] 144 def MaskElementonfloatingiceEnum(): return StringToEnum("MaskElementonfloatingice")[0]145 def MaskElementongroundediceEnum(): return StringToEnum("MaskElementongroundedice")[0]146 def MaskVertexonfloatingiceEnum(): return StringToEnum("MaskVertexonfloatingice")[0]147 def MaskVertexongroundediceEnum(): return StringToEnum("MaskVertexongroundedice")[0]148 144 def MaskIceLevelsetEnum(): return StringToEnum("MaskIceLevelset")[0] 149 145 def MaterialsBetaEnum(): return StringToEnum("MaterialsBeta")[0] … … 470 466 def WatercolumnEnum(): return StringToEnum("Watercolumn")[0] 471 467 def BasalFrictionEnum(): return StringToEnum("BasalFriction")[0] 472 def BasalMeltrateEnum(): return StringToEnum("BasalMeltrate")[0]473 468 def ViscousHeatingEnum(): return StringToEnum("ViscousHeating")[0] 474 469 def QmuTemperatureEnum(): return StringToEnum("QmuTemperature")[0]
Note:
See TracChangeset
for help on using the changeset viewer.