186 |
ENDDO |
ENDDO |
187 |
|
|
188 |
C Z is special since it varies |
C Z is special since it varies |
189 |
WRITE(dn(1),'(a,i6.6)') 'Zd', nlevels(listId) |
WRITE(dn(1),'(a,i6.6)') 'Zd', kdiag(ndId) |
190 |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
191 |
& .AND. (gdiag(ndId)(9:9) .EQ. 'M') ) THEN |
& .AND. (gdiag(ndId)(9:9) .EQ. 'M') ) THEN |
192 |
WRITE(dn(1),'(a,i6.6)') 'Zmd', nlevels(listId) |
WRITE(dn(1),'(a,i6.6)') 'Zmd', kdiag(ndId) |
193 |
ENDIF |
ENDIF |
194 |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
195 |
& .AND. (gdiag(ndId)(9:9) .EQ. 'L') ) THEN |
& .AND. (gdiag(ndId)(9:9) .EQ. 'L') ) THEN |
196 |
WRITE(dn(1),'(a,i6.6)') 'Zld', nlevels(listId) |
WRITE(dn(1),'(a,i6.6)') 'Zld', kdiag(ndId) |
197 |
ENDIF |
ENDIF |
198 |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
IF ( (gdiag(ndId)(10:10) .EQ. 'R') |
199 |
& .AND. (gdiag(ndId)(9:9) .EQ. 'U') ) THEN |
& .AND. (gdiag(ndId)(9:9) .EQ. 'U') ) THEN |
200 |
WRITE(dn(1),'(a,i6.6)') 'Zud', nlevels(listId) |
WRITE(dn(1),'(a,i6.6)') 'Zud', kdiag(ndId) |
201 |
ENDIF |
ENDIF |
202 |
dim(1) = Nr+Nrphys+1 |
dim(1) = Nr+Nrphys+1 |
203 |
ib(1) = 1 |
ib(1) = 1 |