/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ResultOutputSolve/ |
H A D | DXOutputSolver.F90 | 503 NodeIndex(1,:3) = Elem % NodeIndexes(1:3) 506 NodeIndex(1,:3) = Elem % NodeIndexes(1:3) 507 NodeIndex(2,:3) = Elem % NodeIndexes((/ 1,3,4 /)) 510 NodeIndex(1,:3) = Elem % NodeIndexes((/ 1,5,9 /)) 511 NodeIndex(2,:3) = Elem % NodeIndexes((/ 5,2,9 /)) 512 NodeIndex(3,:3) = Elem % NodeIndexes((/ 2,6,9 /)) 513 NodeIndex(4,:3) = Elem % NodeIndexes((/ 6,3,9 /)) 514 NodeIndex(5,:3) = Elem % NodeIndexes((/ 3,7,9 /)) 515 NodeIndex(6,:3) = Elem % NodeIndexes((/ 7,4,9 /)) 516 NodeIndex(7,:3) = Elem % NodeIndexes((/ 4,8,9 /)) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/ |
H A D | Differentials.F90 | 70 INTEGER, POINTER :: NodeIndexes(:) local 92 NodeIndexes => Element % NodeIndexes 113 n,NodeIndexes, Gotit ) 116 n,NodeIndexes, Gotit ) 293 INTEGER, POINTER :: NodeIndexes(:) local 347 NodeIndexes => Element % NodeIndexes 399 n, NodeIndexes ) 411 n,NodeIndexes, Gotit ) 589 NodeIndexes => Element % NodeIndexes 601 q = Reorder(NodeIndexes(p)) [all …]
|
H A D | Adaptive.F90 | 1110 IF ( Edge % NodeIndexes(1) == RefMesh % Elements(i) % NodeIndexes(1) .AND. & 1111 Edge % NodeIndexes(2) == RefMesh % Elements(i) % NodeIndexes(2) .OR. & 1112 Edge % NodeIndexes(2) == RefMesh % Elements(i) % NodeIndexes(1) .AND. & 1113 Edge % NodeIndexes(1) == RefMesh % Elements(i) % NodeIndexes(2) ) EXIT 1439 IF ( Edge % NodeIndexes(1) == RefMesh % Elements(j) % NodeIndexes(1) .AND. & 1440 Edge % NodeIndexes(2) == RefMesh % Elements(j) % NodeIndexes(2) .OR. & 1441 Edge % NodeIndexes(2) == RefMesh % Elements(j) % NodeIndexes(1) .AND. & 1442 Edge % NodeIndexes(1) == RefMesh % Elements(j) % NodeIndexes(2) ) EXIT 1533 RefElement % NodeIndexes 1724 IF ( Child % NodeIndexes(k0) == Element % NodeIndexes(j0) ) THEN [all …]
|
H A D | MeshUtils.F90 | 3620 NodeIndexes => Element % NodeIndexes 4424 NodeIndexes => Element % NodeIndexes 7007 NodeIndexes => Element % NodeIndexes 11113 NodeIndexes => Element % NodeIndexes 18693 NodeIndexes => Element % NodeIndexes 18925 NodeIndexes => Element % NodeIndexes 19262 NodeIndexes => Element % NodeIndexes 20295 NodeIndexes => Element % NodeIndexes 20448 NodeIndexes => Element % NodeIndexes 20489 NodeIndexes => Element % NodeIndexes [all …]
|
H A D | SaveUtils.F90 | 131 INTEGER :: NodeIndexes(:) local 163 IF(Element % NodeIndexes(j) == Parent % NodeIndexes(k)) THEN 185 UseIndexes => Element % NodeIndexes 212 NodeIndexes(1:n) = UseIndexes( Order(1:n) ) 214 NodeIndexes(1:n) = UseIndexes(1:n) 481 Indexes => Element % NodeIndexes 526 NodePerm( Element % NodeIndexes ) = 1 557 INTEGER, POINTER :: NodeIndexes(:) local 593 NodeIndexes => Element % NodeIndexes 606 InvNodePerm(l) = NodeIndexes(k) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/AdjointSSA/ |
H A D | AdjointSSA_GradientSolver.F90 | 141 NodeIndexes(:) local 356 NodeIndexes => Element % NodeIndexes 445 … DJDBeta(DJDBetaPerm(NodeIndexes(1:n)))=DJDBeta(DJDBetaPerm(NodeIndexes(1:n)))+Nodalbetab(1:n) 447 DJDZs(DJDZsPerm(NodeIndexes(1:n)))=DJDZs(DJDZsPerm(NodeIndexes(1:n)))+Nodalzsb(1:n) 449 DJDZb(DJDZbPerm(NodeIndexes(1:n)))=DJDZb(DJDZbPerm(NodeIndexes(1:n)))+Nodalzbb(1:n) 451 DJDRho(DJDRhoPerm(NodeIndexes(1:n)))=DJDRho(DJDRhoPerm(NodeIndexes(1:n)))+Nodalrhob(1:n) 453 DJDEta(DJDEtaPerm(NodeIndexes(1:n)))=DJDEta(DJDEtaPerm(NodeIndexes(1:n)))+Nodaletab(1:n) 467 NodeIndexes => BoundaryElement % NodeIndexes 536 DJDZs(DJDZsPerm(NodeIndexes(1:n)))=DJDZs(DJDZsPerm(NodeIndexes(1:n)))+Nodalzsb(1:n) 538 DJDZb(DJDZbPerm(NodeIndexes(1:n)))=DJDZb(DJDZbPerm(NodeIndexes(1:n)))+Nodalzbb(1:n) [all …]
|
H A D | AdjointSSA_AdjointSolver.F90 | 88 INTEGER, POINTER :: NodeIndexes(:),Perm(:) local 272 NodeIndexes => Element % NodeIndexes 281 IF ( ALL(BoundaryReorder(NodeIndexes(1:n))<1) ) CheckNT = .FALSE. 292 Sol%Values(NSDOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp 307 k = Perm(NodeIndexes(j)) 316 CALL RotateNTSystem( RotVec, NodeIndexes(j) ) 326 Sol%Values(NSDOFs*(Perm(NodeIndexes(j))-1) + kmax)=0._dp 329 Sol%Values(NSDOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp 348 NodeIndexes => Element % NodeIndexes 358 Sol%Values(NSDOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp [all …]
|
H A D | AdjointSSA_CostFluxDivSolver.F90 | 126 INTEGER, POINTER :: NodeIndexes(:) local 270 NodeIndexes => Element % NodeIndexes 274 ElementNodes % x(1:n) = Solver % Mesh % Nodes % x(NodeIndexes) 279 ElementNodes % y(1:n) = Solver % Mesh % Nodes % y(NodeIndexes) 318 Velo(i,1:n)=Velocity(DIM*(VeloPerm(NodeIndexes(1:n))-1)+i) 321 NodeH(1:n)=Zs(ZsPerm(NodeIndexes(1:n)))-Zb(ZbPerm(NodeIndexes(1:n))) 378 DJDZb(DJDZbPerm(NodeIndexes(l)))=DJDZb(DJDZbPerm(NodeIndexes(l))) - & ! - car h=Zs-Zb 381 DJDZs(DJDZsPerm(NodeIndexes(l)))=DJDZs(DJDZsPerm(NodeIndexes(l))) + & ! + car h=Zs-Zb 384 k=(VbPerm(NodeIndexes(l))-1)*c+j 389 … DJDZb(DJDZbPerm(NodeIndexes(l)))=DJDZb(DJDZbPerm(NodeIndexes(l))) - & ! - car h=Zs-Zb [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/src/modules/ |
H A D | FDiffusion.F90 | 76 INTEGER, POINTER :: NodeIndexes(:) local 240 NodeIndexes => CurrentElement % NodeIndexes 258 'Frequency', n, NodeIndexes, GotIt ) 273 'Electric Conductivity',n,NodeIndexes) 277 'Magnetic Permeability',n,NodeIndexes) 306 SomeQuantityPerm(NodeIndexes) ) 329 NodeIndexes => CurrentElement % NodeIndexes 655 INTEGER, POINTER :: NodeIndexes(:),Visited(:) local 682 NodeIndexes => Element % NodeIndexes 689 'Magnetic Permeability',n,NodeIndexes) [all …]
|
H A D | TransportEquation.F90 | 257 NodeIndexes => CurrentElement % NodeIndexes 270 NodeIndexes) 275 NodeIndexes) 277 NodeIndexes) 281 NodeIndexes) 283 NodeIndexes) 285 NodeIndexes) 334 NodeIndexes => CurrentElement % NodeIndexes 369 NodeIndexes) 371 NodeIndexes) [all …]
|
H A D | DCRComplexSolve.F90 | 64 INTEGER, POINTER :: NodeIndexes(:) local 238 NodeIndexes => CurrentElement % NodeIndexes 326 NodeIndexes => CurrentElement % NodeIndexes 429 INTEGER :: n, NodeIndexes(:) local 481 INTEGER :: n, NodeIndexes(:) local 866 IF ( Edge % NodeIndexes(l) == Element % NodeIndexes(k) ) THEN 1070 INTEGER :: n, NodeIndexes(:) local 1253 IF ( Edge % NodeIndexes(j) == Element % NodeIndexes(k) ) THEN 1370 INTEGER :: n, NodeIndexes(:) local 1422 INTEGER :: n, NodeIndexes(:) local [all …]
|
H A D | MovingElstatSolver.F90 | 209 NodeIndexes => CurrentElement % NodeIndexes 212 ElementNodes % x(1:n) = xorig(NodeIndexes) 213 ElementNodes % y(1:n) = yorig(NodeIndexes) 214 ElementNodes % z(1:n) = zorig(NodeIndexes) 296 NodeIndexes => CurrentElement % NodeIndexes 494 NodeIndexes => CurrentElement % NodeIndexes 621 NodeIndexes => CurrentElement % NodeIndexes 646 NodeIndexes => CurrentElement % NodeIndexes 704 NodeIndexes => CurrentElement % NodeIndexes 747 NodeIndexes => CurrentElement % NodeIndexes [all …]
|
H A D | FDiffusion3D.F90 | 76 INTEGER, POINTER :: NodeIndexes(:) local 243 NodeIndexes => CurrentElement % NodeIndexes 261 'Frequency', n, NodeIndexes, GotIt ) 277 'Electric Conductivity',n,NodeIndexes) 280 'Magnetic Permeability',n,NodeIndexes) 322 SomeQuantityPerm(NodeIndexes) ) 656 INTEGER, POINTER :: NodeIndexes(:),Visited(:) local 686 NodeIndexes => Element % NodeIndexes 693 'Magnetic Permeability',n,NodeIndexes) 728 IF ( MINVAL(Reorder(NodeIndexes)) > 0 ) THEN [all …]
|
H A D | PhaseChangeSolve.F90 | 330 NodeIndexes => CurrentElement % NodeIndexes 381 NodeIndexes => CurrentElement % NodeIndexes 456 NodeIndexes => CurrentElement % NodeIndexes 587 NodeIndexes => CurrentElement % NodeIndexes 683 NodeIndexes => CurrentElement % NodeIndexes 701 k = SurfPerm(NodeIndexes(nn)) 1045 IF( NodeIndexes(j) == Parent % NodeIndexes(k) ) THEN 1079 IF(NodeIndexes(p) /= Trip_node) THEN 1218 NodeIndexes => CurrentElement % NodeIndexes 1289 NodeIndexes => CurrentElement % NodeIndexes [all …]
|
H A D | PoissonBoltzmannSolve.F90 | 89 INTEGER, POINTER :: NodeIndexes(:) local 324 NodeIndexes => CurrentElement % NodeIndexes 339 'Charge Density', n, NodeIndexes, GotIt ) 384 NodeIndexes => CurrentElement % NodeIndexes 391 n,NodeIndexes,gotIt ) 493 'Relative Permittivity',n, NodeIndexes ) 742 INTEGER, POINTER :: NodeIndexes(:) local 775 NodeIndexes => Element % NodeIndexes 780 IF ( ANY(Reorder(NodeIndexes) == 0) ) CYCLE 789 'Charge Density', n, NodeIndexes, GotIt ) [all …]
|
H A D | StatElecBoundary.F90 | 58 INTEGER, POINTER :: NodeIndexes(:) local 68 SAVE :: AllocationsDone, n, NodeIndexes, List, CurrentElement, & 86 NodeIndexes => CurrentElement % NodeIndexes 95 'Potential Difference',n,NodeIndexes,gotIt ) 101 'Permittivity',n,NodeIndexes,GotIt) 103 'Relative Permittivity',n,NodeIndexes) 108 'Layer Thickness',n,NodeIndexes,LayerExists) 112 'Layer Permittivity',n,NodeIndexes) 119 HoleSize(1:n) = ListGetReal(List,'Hole Size', n, NodeIndexes) 124 ElemAperture(1:n) = ListGetReal(List,'Gap Height',n,NodeIndexes) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/ |
H A D | DJDBeta_Robin.F90 | 65 INTEGER, POINTER :: Permutation(:), VeloNPerm(:),VeloDPerm(:),BetaPerm(:),NodeIndexes(:) local 180 NodeIndexes => Element % NodeIndexes 184 VisitedNode(NodeIndexes(i))=VisitedNode(NodeIndexes(i))+1.0_dp 188 vn(1) = VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(i))-1)+1) 189 vn(2) = VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(i))-1)+2) 190 vd(1) = VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(i))-1)+1) 191 vd(2) = VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(i))-1)+2) 193 vn(3)=VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(i))-1)+3) 194 vd(3)=VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(i))-1)+3) 202 NodeDJ(i)=NodeDJ(i)*2.0_dp*BetaValues(BetaPerm(NodeIndexes(i))) [all …]
|
H A D | AIFlowSolve_nlD2.F90 | 342 NodeIndexes => CurrentElement % NodeIndexes 402 k = TempPerm(NodeIndexes(i)) 450 NodeIndexes => CurrentElement % NodeIndexes 576 NodeIndexes => CurrentElement % NodeIndexes 611 k = TempPerm(NodeIndexes(i)) 666 RefD(2*dim*(SRPerm(NodeIndexes(i))-1)+1 : & 668 =RefD(2*dim*(SRPerm(NodeIndexes(i))-1)+1 : & 672 RefS(2*dim*(DSPerm(NodeIndexes(i))-1)+1 : & 674 =RefS(2*dim*(DSPerm(NodeIndexes(i))-1)+1 : & 1107 INTEGER, POINTER :: NodeIndexes(:) local [all …]
|
H A D | AIFlowSolve_nlS2.F90 | 342 NodeIndexes => CurrentElement % NodeIndexes 403 k = TempPerm(NodeIndexes(i)) 451 NodeIndexes => CurrentElement % NodeIndexes 577 NodeIndexes => CurrentElement % NodeIndexes 612 k = TempPerm(NodeIndexes(i)) 667 RefD(2*dim*(SRPerm(NodeIndexes(i))-1)+1 : & 669 =RefD(2*dim*(SRPerm(NodeIndexes(i))-1)+1 : & 673 RefS(2*dim*(DSPerm(NodeIndexes(i))-1)+1 : & 675 =RefS(2*dim*(DSPerm(NodeIndexes(i))-1)+1 : & 1142 INTEGER, POINTER :: NodeIndexes(:) local [all …]
|
H A D | Flotation.F90 | 124 INTEGER, POINTER,SAVE :: NodeIndexes(:) local 216 NodeIndexes => Element % NodeIndexes 220 Density(1:n) = ListGetReal( Material, 'SSA Mean Density',n, NodeIndexes,UnFoundFatal=.TRUE.) 226 H=HVar%Values(HVar%Perm(NodeIndexes(i))) 231 bedrock=BedVar%Values(BedVar%Perm(NodeIndexes(i))) 240 ZbVar%Values(ZbVar%Perm(NodeIndexes(i)))=zb 241 ZsVar%Values(ZsVar%Perm(NodeIndexes(i)))=zs 245 DZbDt%Values(DZbDt%Perm(NodeIndexes(i)))=(zb-ZbPrev(ZbVar%Perm(NodeIndexes(i))))/dt 247 DZsDt%Values(DZsDt%Perm(NodeIndexes(i)))=(zs-ZsPrev(ZsVar%Perm(NodeIndexes(i))))/dt 251 topnode=TopPointer(NodeIndexes(i)) [all …]
|
H A D | DJDmu_Adjoint.F90 | 84 INTEGER, POINTER :: NodeIndexes(:) local 204 NodeIndexes => Element % NodeIndexes 206 VisitedNode(NodeIndexes(1:n))=VisitedNode(NodeIndexes(1:n))+1.0_dp 211 Ux(1:n)=VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(1:n))-1)+1) 212 Uy(1:n)=VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(1:n))-1)+2) 237 (- VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(q))-1)+i) * & 238 VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(p))-1)+i)) 241 (- VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(q))-1)+j) * & 242 & VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(p))-1)+i)) 264 c2n = ListGetReal( Material, 'Viscosity Exponent', n, NodeIndexes ) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/fem/tests/HelmholtzFace/ |
H A D | Helmholtz.f90 | 158 IF ( Edge % NodeIndexes(l) == Element % NodeIndexes(k) ) n = n + 1 193 IF( Edge % NodeIndexes(1) > Edge % NodeIndexes(2)) Normal2 = -Normal2 200 Edgeh = (Model % Mesh % Nodes % x(Edge % NodeIndexes(2)) & 307 k = Element % NodeIndexes(j) 445 IF ( Edge % NodeIndexes(l) == Element % NodeIndexes(k) ) n = n + 1 466 IF( Edge % NodeIndexes(2) > Edge % NodeIndexes(1) ) THEN 474 IF( Edge % NodeIndexes(2) > Edge % NodeIndexes(1) ) THEN 494 IF ( Face % NodeIndexes(l) == Element % NodeIndexes(k) ) n = n + 1 557 IF( Element % NodeIndexes(2) < Element % NodeIndexes(1) ) THEN 562 IF( Element % NodeIndexes(3) < Element % NodeIndexes(2) ) THEN [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/Adjoint/ |
H A D | Adjoint_LinearSolver.F90 | 83 INTEGER, POINTER :: NodeIndexes(:) local 267 NodeIndexes => Element % NodeIndexes 276 IF ( ALL(BoundaryReorder(NodeIndexes(1:n))<1) ) CheckNT = .FALSE. 287 Sol%Values(DOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp 302 k = Perm(NodeIndexes(j)) 311 CALL RotateNTSystem( RotVec, NodeIndexes(j) ) 321 Sol%Values(DOFs*(Perm(NodeIndexes(j))-1) + kmax)=0._dp 324 Sol%Values(DOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp 343 NodeIndexes => Element % NodeIndexes 353 Sol%Values(DOFs*(Perm(NodeIndexes(j))-1)+i)=0._dp [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/Solvers/AdjointStokes/ |
H A D | AdjointStokes_GradientMu.F90 | 84 INTEGER, POINTER :: NodeIndexes(:) local 204 NodeIndexes => Element % NodeIndexes 206 VisitedNode(NodeIndexes(1:n))=VisitedNode(NodeIndexes(1:n))+1.0_dp 211 Ux(1:n)=VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(1:n))-1)+1) 212 Uy(1:n)=VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(1:n))-1)+2) 237 (- VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(q))-1)+i) * & 238 VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(p))-1)+i)) 241 (- VelocityN((DIM+1)*(VeloNPerm(NodeIndexes(q))-1)+j) * & 242 & VelocityD((DIM+1)*(VeloDPerm(NodeIndexes(p))-1)+i)) 264 c2n = ListGetReal( Material, 'Viscosity Exponent', n, NodeIndexes ) [all …]
|
/dports/science/elmerfem/elmerfem-release-9.0/elmerice/IceSheet/Greenland/SSA/ |
H A D | Scalar_OUTPUT.F90 | 289 Indexes => Element % NodeIndexes 331 INTEGER,POINTER :: NodeIndexes(:) local 358 NodeIndexes => Element % NodeIndexes 368 NodalH(1:n) = HVar%Values(HVar%Perm(NodeIndexes(1:n))) 383 NodalGM(1:n) = GMVar%Values(GMVar%Perm(NodeIndexes(1:n))) 440 … HRVar%Values(HRVar%Perm(NodeIndexes(i)))*LocalArea(i)/NodeArea(Permutation(NodeIndexes(i))) 469 INTEGER,POINTER :: NodeIndexes(:) local 488 NodeIndexes => Element % NodeIndexes 491 NodalH(1:n) = HVar%Values(HVar%Perm(NodeIndexes(1:n))) 524 INTEGER, POINTER :: NodeIndexes(:) local [all …]
|