50 |
LOGICAL flag |
LOGICAL flag |
51 |
|
|
52 |
CHARACTER*8 parms1 |
CHARACTER*8 parms1 |
|
CHARACTER*3 mate_index |
|
53 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
54 |
|
|
55 |
|
|
65 |
flag = .TRUE. |
flag = .TRUE. |
66 |
DO l=1,listId |
DO l=1,listId |
67 |
IF (flag .AND. freq(l) .EQ. freq(listId) |
IF (flag .AND. freq(l) .EQ. freq(listId) |
68 |
& .AND. phase(l).EQ.phase(listId) ) THEN |
& .AND. phase(l).EQ.phase(listId) |
69 |
|
& .AND. averageFreq(l) .EQ.averageFreq(listId) |
70 |
|
& .AND. averagePhase(l).EQ.averagePhase(listId) |
71 |
|
& .AND. averageCycle(l).EQ.averageCycle(listId) ) THEN |
72 |
DO k=1,MIN(nActive(l),numperlist) |
DO k=1,MIN(nActive(l),numperlist) |
73 |
IF (flag .AND. jdiag(k,l).GT.0) THEN |
IF (flag .AND. jdiag(k,l).GT.0) THEN |
74 |
IF ( cdiag(ndId).EQ.cdiag(jdiag(k,l)) ) THEN |
IF ( cdiag(ndId).EQ.cdiag(jdiag(k,l)) ) THEN |
85 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
86 |
|
|
87 |
IF ( flag ) THEN |
IF ( flag ) THEN |
88 |
IF (ndiagmx+kdiag(ndId).GT.numdiags) THEN |
idiag(mId,listId) = ndiagmx + 1 |
89 |
|
ndiagmx = ndiagmx + kdiag(ndId)*averageCycle(listId) |
90 |
|
IF ( ndiagmx.GT.numdiags ) THEN |
91 |
WRITE(msgBuf,'(A,I4,1X,A)') |
WRITE(msgBuf,'(A,I4,1X,A)') |
92 |
& 'SETDIAG: Not enough space for Diagnostic #',ndId,cdiag(ndId) |
& 'SETDIAG: Not enough space for Diagnostic #',ndId,cdiag(ndId) |
93 |
CALL PRINT_MESSAGE(msgBuf, errUnit, SQUEEZE_RIGHT, myThid) |
CALL PRINT_MESSAGE(msgBuf, errUnit, SQUEEZE_RIGHT, myThid) |
94 |
ELSE |
ELSE |
95 |
WRITE(msgBuf,'(A,I3,A,I4,1X,A)') 'SETDIAG: Allocate', |
WRITE(msgBuf,'(A,2(I3,A),I4,1X,A)') 'SETDIAG: Allocate', |
96 |
& kdiag(ndId), ' Levels for Diagnostic #', ndId, cdiag(ndId) |
& kdiag(ndId), ' x', averageCycle(listId), |
97 |
|
& ' Levels for Diagnostic #', ndId, cdiag(ndId) |
98 |
CALL PRINT_MESSAGE(msgBuf, stdUnit, SQUEEZE_RIGHT, myThid) |
CALL PRINT_MESSAGE(msgBuf, stdUnit, SQUEEZE_RIGHT, myThid) |
99 |
ENDIF |
ENDIF |
|
idiag(mId,listId) = ndiagmx + 1 |
|
|
ndiagmx = ndiagmx + kdiag(ndId) |
|
100 |
ELSE |
ELSE |
101 |
WRITE(msgBuf,'(A,I4,1X,2A)') |
WRITE(msgBuf,'(A,I4,1X,2A)') |
102 |
& '- WARNING - SETDIAG: Diagnostic #', ndId, cdiag(ndId), |
& '- WARNING - SETDIAG: Diagnostic #', ndId, cdiag(ndId), |
110 |
c ---------------------------- |
c ---------------------------- |
111 |
mate = 0 |
mate = 0 |
112 |
IF ( parms1(5:5).EQ.'C') THEN |
IF ( parms1(5:5).EQ.'C') THEN |
113 |
mate_index = parms1(6:8) |
READ(parms1,'(5X,I3)') mate |
|
READ (mate_index,'(I3)') mate |
|
114 |
|
|
115 |
C-- Seach for the same diag (with same freq) to see if already set |
C-- Seach for the same diag (with same freq) to see if already set |
116 |
flag = .TRUE. |
flag = .TRUE. |
117 |
DO l=1,listId |
DO l=1,listId |
118 |
IF (flag .AND. freq(l).EQ.freq(listId) |
IF (flag .AND. freq(l) .EQ.freq(listId) |
119 |
& .AND. phase(l).EQ.phase(listId) ) THEN |
& .AND. phase(l).EQ.phase(listId) |
120 |
|
& .AND. averageFreq(l) .EQ.averageFreq(listId) |
121 |
|
& .AND. averagePhase(l).EQ.averagePhase(listId) |
122 |
|
& .AND. averageCycle(l).EQ.averageCycle(listId) ) THEN |
123 |
DO k=1,MIN(nActive(l),numperlist) |
DO k=1,MIN(nActive(l),numperlist) |
124 |
IF (flag .AND. jdiag(k,l).GT.0) THEN |
IF (flag .AND. jdiag(k,l).GT.0) THEN |
125 |
IF (cdiag(mate).EQ.cdiag(jdiag(k,l)) ) THEN |
IF (cdiag(mate).EQ.cdiag(jdiag(k,l)) ) THEN |
133 |
ENDDO |
ENDDO |
134 |
|
|
135 |
IF ( flag ) THEN |
IF ( flag ) THEN |
136 |
IF (ndiagmx+kdiag(mate).GT.numdiags) THEN |
mdiag(mId,listId) = ndiagmx + 1 |
137 |
|
ndiagmx = ndiagmx + kdiag(mate)*averageCycle(listId) |
138 |
|
IF ( ndiagmx.GT.numdiags ) THEN |
139 |
WRITE(msgBuf,'(A,I4,1X,A)') |
WRITE(msgBuf,'(A,I4,1X,A)') |
140 |
& 'SETDIAG: Not enough space for Counter Diagnostic #', |
& 'SETDIAG: Not enough space for Counter Diagnostic #', |
141 |
& mate, cdiag(mate) |
& mate, cdiag(mate) |
142 |
CALL PRINT_MESSAGE(msgBuf, errUnit, SQUEEZE_RIGHT, myThid) |
CALL PRINT_MESSAGE(msgBuf, errUnit, SQUEEZE_RIGHT, myThid) |
143 |
ELSE |
ELSE |
144 |
WRITE(msgBuf,'(A,I3,A,I4,1X,A)') 'SETDIAG: Allocate', |
WRITE(msgBuf,'(A,2(I3,A),I4,1X,A)') 'SETDIAG: Allocate', |
145 |
& kdiag(mate), ' Levels for Count.Diag #', mate, cdiag(mate) |
& kdiag(mate), ' x', averageCycle(listId), |
146 |
|
& ' Levels for Count.Diag #', mate, cdiag(mate) |
147 |
CALL PRINT_MESSAGE(msgBuf, stdUnit, SQUEEZE_RIGHT, myThid) |
CALL PRINT_MESSAGE(msgBuf, stdUnit, SQUEEZE_RIGHT, myThid) |
148 |
ENDIF |
ENDIF |
|
mdiag(mId,listId) = ndiagmx + 1 |
|
|
ndiagmx = ndiagmx + kdiag(mate) |
|
149 |
ELSE |
ELSE |
150 |
WRITE(msgBuf,'(A,I4,1X,2A)') |
WRITE(msgBuf,'(A,I4,1X,2A)') |
151 |
& '- NOTE - SETDIAG: Counter Diagnostic #', mate, cdiag(mate), |
& '- NOTE - SETDIAG: Counter Diagnostic #', mate, cdiag(mate), |