66 |
REAL rhsMax |
REAL rhsMax |
67 |
REAL rhsNorm |
REAL rhsNorm |
68 |
|
|
69 |
|
CcnhDebugStarts |
70 |
|
CHARACTER*(MAX_LEN_FNAM) suff |
71 |
|
CcnhDebugEnds |
72 |
|
|
73 |
|
|
74 |
C-- Initialise inverter |
C-- Initialise inverter |
75 |
etaNBuf(1,myThid) = 0. D0 |
etaNBuf(1,myThid) = 0. D0 |
76 |
errBuf(1,myThid) = 0. D0 |
errBuf(1,myThid) = 0. D0 |
80 |
CcnhDebugStarts |
CcnhDebugStarts |
81 |
C _EXCH_XY_R8( cg2d_b, myThid ) |
C _EXCH_XY_R8( cg2d_b, myThid ) |
82 |
C CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.0 CG2D_B' , 1, myThid ) |
C CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.0 CG2D_B' , 1, myThid ) |
83 |
|
C suff = 'unnormalised' |
84 |
|
C CALL WRITE_FLD_XY_RL ( 'cg2d_b.',suff, cg2d_b, 1, myThid) |
85 |
CcnhDebugEnds |
CcnhDebugEnds |
86 |
|
|
87 |
C-- Normalise RHS |
C-- Normalise RHS |
117 |
_EXCH_XY_R8( cg2d_x, myThid ) |
_EXCH_XY_R8( cg2d_x, myThid ) |
118 |
CcnhDebugStarts |
CcnhDebugStarts |
119 |
C CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.1 CG2D_B' , 1, myThid ) |
C CALL PLOT_FIELD_XYRL( cg2d_b, 'CG2D.1 CG2D_B' , 1, myThid ) |
120 |
|
C suff = 'normalised' |
121 |
|
C CALL WRITE_FLD_XY_RL ( 'cg2d_b.',suff, cg2d_b, 1, myThid) |
122 |
CcnhDebugEnds |
CcnhDebugEnds |
123 |
|
|
124 |
C-- Initial residual calculation |
C-- Initial residual calculation |