!ES = ES_H(:, N:N+1)
!ES = reshape([ES_H(:, N), ES_H(:, N+1)], [2*N, 2])
!call ES_TRIPOD(theta, Phi0)
! the output file
open(1, file = 'eb.dat' , status = 'replace')
!do i = 1, 2*N
write(1,*)EV_H(1)
! write(*,*)i
!enddo
close(1)
9.86 KB, 下载次数: 0
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |