elements and nodes do not need to start at 1 necessarily

This commit is contained in:
f.basile 2021-02-02 17:12:03 +01:00 committed by Martin Diehl
parent ee40884916
commit 47f52b9095
2 changed files with 3 additions and 3 deletions

@ -1 +1 @@
Subproject commit f7297868d932edbadef170c584ac5349f2c0a6d5
Subproject commit f3cd8f4c4b2d154e179d9bb5f396a77535f88556

View File

@ -370,8 +370,8 @@ subroutine uedinc(inc,incsub)
if (inc > inc_written) then
allocate(d_n(3,size(mesh_FEM2DAMASK_node)))
do n = 1, size(d_n,2)
allocate(d_n(3,count(mesh_FEM2DAMASK_node /= -1)))
do n = lbound(mesh_FEM2DAMASK_node,1), ubound(mesh_FEM2DAMASK_node,1)
if (mesh_FEM2DAMASK_node(n) /= -1) then
call nodvar(1,n,d_n(1:3,mesh_FEM2DAMASK_node(n)),nqncomp,nqdatatype)
if(nqncomp == 2) d_n(3,mesh_FEM2DAMASK_node(n)) = 0.0_pReal