I cant compile below subroutine ,
subroutine vdload (
C Read only (unmodifiable)variables -
1 nblock, ndim, stepTime, totalTime,
2 amplitude, curCoords, velocity, dirCos, jltyp, sname,
C Write only (modifiable) variable -
1 value )
C
include 'vaba_param.inc'
C
dimension curCoords(nblock,ndim), velocity(nblock,ndim),
1 dirCos(nblock,ndim,ndim), value(nblock)
character*80 sname
Real R,r,i,j,P0
t=steptime
P0=2.3000000
R=0.1
x=(2*R**2)*t
C
do 100 km = 1, nblock
r=sqrt((curCoords(km,1))**2+(curCoords(km,2))**2)
P=P0*t*exp(-(r/x))
value=P
100 continue
return
end
several times i debog it and retry to run, I encountered this error
anyone can help me solve this problem ?
hi Reza Amini
I had the same problem
the solution is in the Fortran fixed format, for me the problem was in the number of spaces befor the starting word of each line.
try this code
SUBROUTINE vdload (nBlock, ndim, stepTime, totalTime,
1 amplitude, curCoords, velocity, dirCos, jltyp, sname,
1 value )
INCLUDE 'VABA_PARAM.INC'
DIMENSION curCoords(nBlock,ndim), velocity(nBlock,ndim),
1 dirCos(nBlock,ndim,ndim),
1 value(nBlock)
CHARACTER*80 sname
parameter(Pressure = 5000.d0, Radi = 4.d0)
parameter(iX =1, iY =2, iZ =3)
xini= 10.d0
yini= 10.d0
xvelo=200.d0
yvelo=200.d0
xpos= xini + (xvelo * stepTime)
ypos= yini + (yvelo * stepTime)
DO i = 1, nBlock
CurPos = (curCoords(i,iX)-xpos)**2+(curCoords(i,iY)-ypos)**2
if (CurPos < Radi**2) then
value(i) = Pressure
else
value(i) = 0
endif
END DO
RETURN
END
- Badges
- Science topic
- Similar topics
- Biological Science
- Kinesiology
- Running