69 |
factor = FLOAT(ndiag(ip,bi,bj)) |
factor = FLOAT(ndiag(ip,bi,bj)) |
70 |
IF (ndiag(ip,bi,bj).NE.0) factor = 1. _d 0 / factor |
IF (ndiag(ip,bi,bj).NE.0) factor = 1. _d 0 / factor |
71 |
|
|
72 |
|
#ifdef ALLOW_FIZHI |
73 |
DO j = 1,sNy+1 |
DO j = 1,sNy+1 |
74 |
DO i = 1,sNx+1 |
DO i = 1,sNx+1 |
75 |
IF ( qdiag(i,j,ipnt,bi,bj) .LE. undef ) THEN |
IF ( qdiag(i,j,ipnt,bi,bj) .LE. undef ) THEN |
79 |
ENDIF |
ENDIF |
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
|
#else /* ALLOW_FIZHI */ |
83 |
|
DO j = 1,sNy+1 |
84 |
|
DO i = 1,sNx+1 |
85 |
|
qtmp(i,j) = qdiag(i,j,ipnt,bi,bj)*factor |
86 |
|
ENDDO |
87 |
|
ENDDO |
88 |
|
#endif /* ALLOW_FIZHI */ |
89 |
|
|
90 |
ELSE |
ELSE |
91 |
C- With counter diagnostics => average = Sum / counter: |
C- With counter diagnostics => average = Sum / counter: |