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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
SUBROUTINE mxpeak(Sxx,Tpeak,Domfqt,Ntfreq,Speak,Domfqs,Nsfreq,
& Maxsxx,Nform,Spclab)
IMPLICIT NONE
c-----------------------------------------------------------------------
INCLUDE 'notset.prm'
c-----------------------------------------------------------------------
CHARACTER labvec*(2),domfrq*(2),Spclab*(*)
DOUBLE PRECISION Sxx,Maxsxx
INTEGER i,i2,Domfqt,Domfqs,Nform,Pkidx,frq1,Tpeak,Ntfreq,Speak,
& Nsfreq
DIMENSION Sxx(*),Tpeak(*),Speak(*),labvec(11)
c-----------------------------------------------------------------------
LOGICAL dpeq
EXTERNAL dpeq
c-----------------------------------------------------------------------
DATA labvec/'t1','t2','t3','t4','t5','s1','s2','s3','s4','s5',
& 's6'/
c-----------------------------------------------------------------------
domfrq = 'no'
i2 = 0
c-----------------------------------------------------------------------
IF(Domfqt.eq.NOTSET.and.Domfqs.eq.NOTSET)THEN
WRITE(Nform,1010)Spclab,domfrq
RETURN
ELSE IF(Domfqt.eq.NOTSET)THEN
frq1=Domfqs
i2=5
ELSE IF(Domfqs.eq.NOTSET)THEN
frq1=Domfqt
ELSE IF(Sxx(Domfqt).gt.Sxx(Domfqs))THEN
frq1=Domfqt
ELSE
frq1=Domfqs
i2=5
END IF
c-----------------------------------------------------------------------
IF(dpeq(Maxsxx,Sxx(frq1)))THEN
IF(i2.eq.0)THEN
i=1
DO WHILE (i.le.Ntfreq)
IF(Tpeak(i).eq.frq1)THEN
i2=i
i=Ntfreq
END IF
i=i+1
END DO
ELSE
i=1
DO WHILE (i.le.Nsfreq)
IF(Speak(i).eq.frq1)THEN
i2=i+i2
i=Nsfreq
END IF
i=i+1
END DO
END IF
domfrq=labvec(i2)
END IF
c-----------------------------------------------------------------------
WRITE(Nform,1010)Spclab,domfrq
c-----------------------------------------------------------------------
1010 FORMAT(a,'.dom: ',a)
RETURN
END
|
