-
Notifications
You must be signed in to change notification settings - Fork 0
/
prepro.f
75 lines (67 loc) · 4.46 KB
/
prepro.f
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
66
67
68
69
70
71
72
73
74
75
c *******************************************************
c **************** pre-process prepro ******************
c
subroutine prepro
c
include 'common.inc'
c
c ********* if double precision ******************
c implicit real*8(d)
c
c ---------- set harmonic matrix dn( , , ) ----------
do j=1,nmd
dn(j,1,1)=0.0
dn(j,1,2)=-2.0d0*float(j)+1.0
dn(j,2,1)= 2.0d0*float(j)-1.0
dn(j,2,2)=0.0
enddo
c ---------- clear matrix ----------
do j=1, nmd
do k=1,nom
dkk(j,k)=0.0
enddo
enddo
c
ddh=dpi/float(ndi)
do n=0,2*(2*nmd-1)
do nt=0,ndi
dsine(n,nt)=sin(float(n*nt)*ddh)
dcosi(n,nt)=cos(float(n*nt)*ddh)
enddo
enddo
c
return
end
c ******************************************
c ********** initialize matrix init ************
c
subroutine init
c
include 'common.inc'
c
c ********* if double precision ******************
c implicit real*8(d)
c
c nom1 nom1=mmu*npo1
c nom2 nom2=mmu*(ncn+nbr)
C nbb nbb =4*nb-1
do i=1, nom1
do j=1, nbb
dh(i,j) =0.0
enddo
do j=1, nom2
dhc(i,j)=0.0
dhl(j,i)=0.0
enddo
dk(i)=0.0
enddo
c....
do i=1,nom2
do j=1,nom2
dhcl(i,j)=0.0
enddo
dk(nom1+i)=0.0
enddo
c
return
end