64 |
_RL smallNo |
_RL smallNo |
65 |
c _RL Rjjm,Rjjp |
c _RL Rjjm,Rjjp |
66 |
|
|
67 |
IF (inAdMode) THEN |
c IF (inAdMode) THEN |
68 |
|
c smallNo = 1.0D-20 |
69 |
|
c ELSE |
70 |
smallNo = 1.0D-20 |
smallNo = 1.0D-20 |
71 |
ELSE |
c ENDIF |
|
smallNo = 1.0D-20 |
|
|
ENDIF |
|
72 |
#endif |
#endif |
73 |
|
|
74 |
DO j=1-Oly,sNy+Oly |
DO j=1-OLy,sNy+OLy |
75 |
uT(1-Olx,j)=0. |
uT(1-OLx,j)=0. |
76 |
uT(2-Olx,j)=0. |
uT(2-OLx,j)=0. |
77 |
uT(sNx+Olx,j)=0. |
uT(sNx+OLx,j)=0. |
78 |
ENDDO |
ENDDO |
79 |
DO j=1-Oly,sNy+Oly |
DO j=1-OLy,sNy+OLy |
80 |
DO i=1-Olx+2,sNx+Olx-1 |
DO i=1-OLx+2,sNx+OLx-1 |
81 |
Rjp=(tracer(i+1,j)-tracer( i ,j))*maskLocW(i+1,j) |
Rjp=(tracer(i+1,j)-tracer( i ,j))*maskLocW(i+1,j) |
82 |
Rj =(tracer( i ,j)-tracer(i-1,j))*maskLocW( i ,j) |
Rj =(tracer( i ,j)-tracer(i-1,j))*maskLocW( i ,j) |
83 |
Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskLocW(i-1,j) |
Rjm=(tracer(i-1,j)-tracer(i-2,j))*maskLocW(i-1,j) |