- 积分
- 0
- 注册时间
- 2010-10-4
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-5-9 10:19:34
|
显示全部楼层
来自 美国
To make the functions of "module" clear, I just created a small example. It may fit someone's interests.
program test
use mydata
IMPLICIT REAL(A-H, O-Z)
call readfile
write (*,*) ndim, npt
do i = 1, npt
write(*,*) ptcoord(1:3,i)
end do
end program test
module mydata
IMPLICIT REAL(A-H, O-Z)
integer :: npt,ndim
real, allocatable :: ptcoord(:,
end module mydata
subroutine readfile
use mydata
IMPLICIT REAL(A-H, O-Z)
open(502, FILE= 'input.dat')
read(502,FMT='(I8)') ndim
read(502,FMT='(I8)') npt
allocate (ptcoord(ndim,npt))
do ipt = 1, npt
read(502,FMT='(3E16.6)') ptcoord(1:ndim,ipt)
end do !i
close(502, STATUS='KEEP')
end subroutine readfile |
|