23 |
|
|
24 |
C == GLobal variables == |
C == GLobal variables == |
25 |
#include "SIZE.h" |
#include "SIZE.h" |
26 |
#include "GRID.h" |
#ifdef OLD_DST3_FORMULATION |
27 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
28 |
#include "PARAMS.h" |
#include "PARAMS.h" |
29 |
|
#endif |
30 |
|
#include "GRID.h" |
31 |
#include "GAD.h" |
#include "GAD.h" |
32 |
|
|
33 |
C == Routine arguments == |
C == Routine arguments == |
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 |
79 |
|
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) |