Tek-Tips is the largest IT community on the Internet today!

Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!

  • Congratulations cowski on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

nblock and coordNp in vdisp

Status
Not open for further replies.

liq001

Aerospace
Joined
Feb 25, 2022
Messages
56
Location
US
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
 
Anyone has experience of using Vdisp? Need help. Thanks
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor

Back
Top