nblock and coordNp in vdisp

[liq001]

Hello everyone,

I am writing user subroutine using vdisp. My I ask why the return value of nblock is always 1 and the coordNP is always 0? the coordNp should return current nodal coordinates at 1,2 & 3 direction right? Please help if My understanding is not correct. Part of my code is below used to print the coordNp

User subroutine VDISP to prescribe boundary conditions
c
subroutine vdisp(
c Read only -
* nblock, nDof, nCoord, kstep, kinc,
* stepTime, totalTime, dtNext, dt,
* cbname, jBCType, jDof, jNodeUid, amp,
* coordNp, u, v, a, rf, rmass, rotaryI,
c Write only -
* rval )
c
include 'vaba_param.inc'
parameter( zero = 0.d0, half = 0.5d0 )
dimension jDof(nDof), jNodeUid(nblock),
* amp(nblock), coordNp(nCoord, nblock),
* u(nDof,nblock), v(nDof,nblock), a(nDof,nblock),
* rf(nDof,nblock), rmass(nblock), rotaryI(3,3,nblock),
* rval(nDof,nblock)
c
character*80 cbname
print*, 'nblock', nblock

do k = 1, nblock
do j = 1, nCoord
currCoords = coordNp(j,k)
end do
end do

print*, 'currCoords', currCoords


Thanks
lee

 

[liq001]

Anyone has experience of using Vdisp? Need help. Thanks