39 |
INTEGER i,j |
INTEGER i,j |
40 |
_RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
41 |
|
|
|
CcnhDEbugStarts |
|
|
INTEGER i1, j1 |
|
|
CcnhDebugEnds |
|
|
|
|
42 |
C-- Pressure equation source term |
C-- Pressure equation source term |
43 |
C Term is the vertical integral of the divergence of the |
C Term is the vertical integral of the divergence of the |
44 |
C time tendency terms along with a relaxation term that |
C time tendency terms along with a relaxation term that |
78 |
|
|
79 |
DO j=1,sNy |
DO j=1,sNy |
80 |
DO i=1,sNx |
DO i=1,sNx |
81 |
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) + |
cg2d_b(i,j,bi,bj) = cg2d_b(i,j,bi,bj) + |
82 |
& pf(i,j+1)-pf(i,j) |
& pf(i,j+1)-pf(i,j) |
83 |
ENDDO |
ENDDO |
84 |
ENDDO |
ENDDO |
85 |
|
|
|
CcnhDebugSTarts |
|
|
C IF ( K .EQ. 10 .OR. K .EQ. 1 ) THEN |
|
|
C i1 = 54 |
|
|
C j1 = 26 |
|
|
C WRITE(0,*) ' @ I = ', i1, ' , J = ', j1, ' K = ', K |
|
|
C WRITE(0,*) ' cg2d_b = ', cg2d_b(i1,j1,bi,bj) |
|
|
C ENDIF |
|
|
CcnhDebugEnds |
|
|
|
|
86 |
|
|
87 |
RETURN |
RETURN |
88 |
END |
END |