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 |