1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
!
! CalculiX - A 3-dimensional finite element program
! Copyright (C) 1998-2015 Guido Dhondt
!
! This program is free software; you can redistribute it and/or
! modify it under the terms of the GNU General Public License as
! published by the Free Software Foundation(version 2);
!
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! You should have received a copy of the GNU General Public License
! along with this program; if not, write to the Free Software
! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
!
subroutine writepf(d,bjr,bji,freq,nev,mode,nherm)
!
! writes the participation factors to unit 5
!
implicit none
!
integer j,nev,mode,nherm
real*8 d(*),bjr(*),bji(*),freq,pi
!
pi=4.d0*datan(1.d0)
!
write(5,*)
if(mode.eq.0) then
write(5,100) freq
else
write(5,101) mode
endif
!
100 format('P A R T I C I P A T I O N F A C T O R S F O R',
&' F R E Q U E N C Y ',e20.13,' (CYCLES/TIME)')
101 format('P A R T I C I P A T I O N F A C T O R S F O R',
&' M O D E ',i5)
!
if(nherm.eq.1) then
write(5,*)
write(5,*) 'MODE NO FREQUENCY FACTOR'
write(5,*) ' (CYCLES/TIME) REAL IMAGINARY'
write(5,*)
do j=1,nev
write(5,'(i7,3(2x,e14.7))') j,d(j)/(2.d0*pi),bjr(j),bji(j)
enddo
else
write(5,*)
write(5,*)
&'MODE NO FREQ. (REAL) FREQ. (IMAG) FACTOR'
write(5,*)
&' (CYCLES/TIME) (RAD/TIME) REAL IMAGINARY'
write(5,*)
do j=1,nev
write(5,'(i7,4(2x,e14.7))') j,d(2*j-1)/(2.d0*pi),d(2*j),
& bjr(j),bji(j)
enddo
endif
!
return
end
|
