83 |
ENDDO |
ENDDO |
84 |
ENDDO |
ENDDO |
85 |
|
|
|
C- jmc: Patrick, will this work for TAF ? |
|
86 |
IF (K.EQ.1) THEN |
IF (K.EQ.1) THEN |
87 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
88 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
93 |
DO j=1-Oly,sNy+Oly |
DO j=1-Oly,sNy+Oly |
94 |
DO i=1-Olx,sNx+Olx |
DO i=1-Olx,sNx+Olx |
95 |
sigmaR(i,j,k)= maskC(i,j,k,bi,bj) |
sigmaR(i,j,k)= maskC(i,j,k,bi,bj) |
96 |
& *recip_drC(k)*rkFac |
& *recip_drC(k)*rkSign |
97 |
& *(sigKm1(i,j)-sigKp1(i,j)) |
& *(sigKp1(i,j)-sigKm1(i,j)) |
98 |
ENDDO |
ENDDO |
99 |
ENDDO |
ENDDO |
100 |
ENDIF |
ENDIF |
|
C- jmc: leave the old code commented (in case ...) |
|
|
c DO j=1-Oly,sNy+Oly |
|
|
c DO i=1-Olx,sNx+Olx |
|
|
c IF (K.NE.1 .AND. hfacC(i,j,k,bi,bj).NE.0.) THEN |
|
|
c sigmaR(i,j,k)=recip_drC(k)*rkFac*(sigKm1(i,j)-sigKp1(i,j)) |
|
|
c ELSE |
|
|
c sigmaR(i,j,k)=0. |
|
|
c ENDIF |
|
|
c ENDDO |
|
|
c ENDDO |
|
101 |
|
|
102 |
RETURN |
RETURN |
103 |
END |
END |