forums.silverfrost.com

[tara87]

I have a program with error775.please help me to run it.what should I do?

PROGRAM p72
!-------------------------------------------------------------------------
! Program 7.2 Plane or axisymmetric analysis of steady seepage using
! 4-node rectangular quadrilaterals. Mesh numbered
! in x(r)- or y(z)- direction.
!-------------------------------------------------------------------------

IMPLICIT NONE
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::fixed_freedoms,i,iel,k,loaded_nodes,nci,ndim=2,nels,neq,nip=4, &
nod=4,nn,np_types,nxe,nye,determinant
REAL(iwp)::det,one=1.0_iwp,penalty=1.0e20_iwp,zero=0.0_iwp
CHARACTER(LEN=15)::dir,element='quadrilateral',type_2d
!-----------------------dynamic arrays------------------------------------
INTEGER,ALLOCATABLE::etype(,g_num(:,,kdiag(,node(,num(
REAL(iwp),ALLOCATABLE::coord(:,,der(:,,deriv(:,,disps(,fun(, &
gc(,g_coord(:,,jac(:,,kay(:,,kc(:,,kv(,kvh(,loads(, &
points(:,,prop(:,,value(,weights(,x_coords(,y_coords(
!-----------------------input and initialisation--------------------------
OPEN(10,FILE='fe95.dat')
OPEN(11,FILE='fe95.res')
READ(10,*)type_2d,dir,nxe,nye,np_types
CALL mesh_size(element,nod,nels,nn,nxe,nye)
neq=nn
ALLOCATE(points(nip,ndim),g_coord(ndim,nn),coord(nod,ndim), &
jac(ndim,ndim),weights(nip),der(ndim,nod),deriv(ndim,nod), &
kc(nod,nod),num(nod),g_num(nod,nels),kay(ndim,ndim),etype(nels), &
x_coords(nxe+1),y_coords(nye+1),prop(ndim,np_types),gc(ndim),fun(nod), &
kdiag(neq),loads(0:neq),disps(0:neq))
READ(10,*)prop
etype=1
IF(np_types>1)READ(10,*)etype
READ(10,*)x_coords,y_coords
!-----------------------loop the elements to find global arrays sizes-----
kdiag=0
elements_1: DO iel=1,nels
CALL geom_rect(element,iel,x_coords,y_coords,coord,num,dir)
g_num(:,iel)=num
g_coord(:,num)=TRANSPOSE(coord)
CALL fkdiag(kdiag,num)
END DO elements_1
CALL mesh(g_coord,g_num,12)
DO i=2,neq
kdiag(i)=kdiag(i)+kdiag(i-1)
END DO
ALLOCATE(kv(kdiag(neq)),kvh(kdiag(neq)))
WRITE(11,'(2(A,I5))') &
"There are",neq," equations and the skyline storage is",kdiag(neq)
CALL sample(element,points,weights)
!-----------------------global conductivity matrix assembly---------------
kv=zero
gc=one
elements_2: DO iel=1,nels
kay=zero
DO i=1,ndim
kay(i,i)=prop(i,etype(iel))
END DO
num=g_num(:,iel)
coord=TRANSPOSE(g_coord(:,num))
kc=zero
gauss_pts_1: DO i=1,nip
CALL shape_der(der,points,i)
CALL shape_fun(fun,points,i)
jac=MATMUL(der,coord)
det=determinant(jac)
CALL invert(jac)
deriv=MATMUL(jac,der)
IF(type_2d=='axisymmetric')gc=MATMUL(fun,coord)
kc=kc+MATMUL(MATMUL(TRANSPOSE(deriv),kay),deriv)*det*weights(i)*gc(1)
END DO gauss_pts_1
CALL fsparv(kv,kc,num,kdiag)
END DO elements_2
kvh=kv
!-----------------------specify boundary values---------------------------
loads=zero
READ(10,*)loaded_nodes
IF(loaded_nodes/=0)READ(10,*)(k,loads(k),i=1,loaded_nodes)
READ(10,*)fixed_freedoms
IF(fixed_freedoms/=0)THEN
ALLOCATE(node(fixed_freedoms),value(fixed_freedoms))
READ(10,*)(node(i),value(i),i=1,fixed_freedoms)
kv(kdiag(node))=kv(kdiag(node))+penalty
loads(node)=kv(kdiag(node))*value
END IF
!-----------------------equation solution---------------------------------
CALL sparin(kv,kdiag)
CALL spabac(kv,loads,kdiag)
!-----------------------retrieve nodal net flow rates---------------------
CALL linmul_sky(kvh,loads,disps,kdiag)
WRITE(11,'(/A)')" Node Total Head Flow rate"
DO k=1,nn
WRITE(11,'(I5,2E12.4)')k,loads(k),disps(k)
END DO
disps(0)=zero
WRITE(11,'(/A)')" Inflow Outflow"
WRITE(11,'(5X,2E12.4)')

[Robert]

Can you give us a clue: on which line do you get the error?

[jjgermis]

Perhaps you could include the program as code and not as normal test. This makes reading it easier.

[brucebowler]

And since it's not "all there", can you make a trim out the stuff that's irrelevant and post a smaller version that exhibits the problem?

[JohnCampbell]

Are you developing this code, as the error 775 is probably due to using an END rather than an END DO statement.

Your use of ALLOCATE is beyond my experience, although it may be valid.
I would not have used the ALLOCATE statements like:
.... x_coords(nxe+1),y_coords(nye+1) ... or
ALLOCATE(kv(kdiag(neq)),kvh(kdiag(neq)))

and preferred each array in a seperate allocate statement

Your use of arrays appears luxurious by my old practices, but can work for a smallish problem.