11character(*),
intent(in) :: fext
12integer,
intent(in) :: ik
13real(8),
intent(in) :: evalfv(nstfv,nspnfv)
18fname=
'EVALFV'//trim(fext)
22 call putrd(fname,ik,v1=
vkl(1:3,ik),n1=nstfv,n2=nspnfv,nrv=nstfv*nspnfv, &
28 inquire(iolength=recl)
vkl(1:3,ik),nstfv,nspnfv,evalfv
29 open(200,
file=fname,form=
'UNFORMATTED',access=
'DIRECT',recl=recl)
30 write(200,rec=ik)
vkl(1:3,ik),nstfv,nspnfv,evalfv
real(8), dimension(:,:), allocatable vkl
subroutine putrd(fname, irec, n1, n2, n3, v1, v2, nrv, rva, nzv, zva)