| 106 | IF ( time4SnapShot ) THEN | IF ( time4SnapShot ) THEN | 
| 107 | C--       switch ON diagnostics of output-stream # n | C--       switch ON diagnostics of output-stream # n | 
| 108 | DO m=1,nActive(n) | DO m=1,nActive(n) | 
| 109 | c             nd = jdiag(m,n) | c             nd = ABS(jdiag(m,n)) | 
| 110 | c             IF (ndiag(nd).NE.0) WRITE(0,*) myIter,nd,ndiag(nd),' ->',0 | c             IF (ndiag(nd).NE.0) WRITE(0,*) myIter,nd,ndiag(nd),' ->',0 | 
| 111 | ip = idiag(m,n) | ip = idiag(m,n) | 
| 112 | IF (ip.GT.0) ndiag(ip,bi,bj) = 0 | IF (ip.GT.0) ndiag(ip,bi,bj) = 0 | 
| 114 | ELSE | ELSE | 
| 115 | C--       switch OFF diagnostics of output-stream # n | C--       switch OFF diagnostics of output-stream # n | 
| 116 | DO m=1,nActive(n) | DO m=1,nActive(n) | 
| 117 | c             nd = jdiag(m,n) | c             nd = ABS(jdiag(m,n)) | 
| 118 | c             IF (ndiag(nd).NE.-1) WRITE(0,*) myIter,nd,ndiag(nd),' ->',-1 | c             IF (ndiag(nd).NE.-1) WRITE(0,*) myIter,nd,ndiag(nd),' ->',-1 | 
| 119 | ip = idiag(m,n) | ip = idiag(m,n) | 
| 120 | IF (ip.GT.0) ndiag(ip,bi,bj) = -1 | IF (ip.GT.0) ndiag(ip,bi,bj) = -1 |