37 |
INTEGER bi, bj |
INTEGER bi, bj |
38 |
INTEGER I, J, K |
INTEGER I, J, K |
39 |
_RL faceArea |
_RL faceArea |
40 |
_RL myNorm |
_RS myNorm |
41 |
_RL aC, aCw, aCs |
_RL aC, aCw, aCs |
42 |
|
|
43 |
C-- Initialise laplace operator |
C-- Initialise laplace operator |
101 |
|
|
102 |
C-- Update overlap regions |
C-- Update overlap regions |
103 |
CcnhDebugStarts |
CcnhDebugStarts |
104 |
C CALL PLOT_FIELD_XYR8( aW2d, 'AW2D INI_CG2D.1' , 1, myThid ) |
C CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.1' , 1, myThid ) |
105 |
C CALL PLOT_FIELD_XYR8( aS2d, 'AS2D INI_CG2D.1' , 1, myThid ) |
C CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.1' , 1, myThid ) |
106 |
CcnhDebugEnds |
CcnhDebugEnds |
107 |
_EXCH_XY_R4(aW2d, myThid) |
_EXCH_XY_R4(aW2d, myThid) |
108 |
_EXCH_XY_R4(aS2d, myThid) |
_EXCH_XY_R4(aS2d, myThid) |
109 |
CcnhDebugStarts |
CcnhDebugStarts |
110 |
C CALL PLOT_FIELD_XYR8( aW2d, 'AW2D INI_CG2D.2' , 1, myThid ) |
C CALL PLOT_FIELD_XYRS( aW2d, 'AW2D INI_CG2D.2' , 1, myThid ) |
111 |
C CALL PLOT_FIELD_XYR8( aS2d, 'AS2D INI_CG2D.2' , 1, myThid ) |
C CALL PLOT_FIELD_XYRS( aS2d, 'AS2D INI_CG2D.2' , 1, myThid ) |
112 |
CcnhDebugEnds |
CcnhDebugEnds |
113 |
|
|
114 |
C-- Initialise preconditioner |
C-- Initialise preconditioner |