4 |
#ifdef ALLOW_NONHYDROSTATIC |
#ifdef ALLOW_NONHYDROSTATIC |
5 |
CBOP |
CBOP |
6 |
C !ROUTINE: CG3D.h |
C !ROUTINE: CG3D.h |
7 |
C !INTERFACE: |
C !INTERFACE: |
8 |
C include "CG3D.h" |
C include "CG3D.h" |
9 |
C !DESCRIPTION: \bv |
C !DESCRIPTION: \bv |
10 |
C *==========================================================* |
C *==========================================================* |
11 |
C | CG3D.h |
C | CG3D.h |
12 |
C | o Three-dimensional conjugate gradient solver header. |
C | o Three-dimensional conjugate gradient solver header. |
13 |
C *==========================================================* |
C *==========================================================* |
14 |
C | The common blocks set up here are used in the elliptic |
C | The common blocks set up here are used in the elliptic |
15 |
C | equation inversion. They are also used as the interface |
C | equation inversion. They are also used as the interface |
16 |
C | to the rest of the model. To set the source term for the |
C | to the rest of the model. To set the source term for the |
17 |
C | solver set the appropriate array below. To read the |
C | solver set the appropriate array below. To read the |
18 |
C | solution read from the appropriate array below. |
C | solution read from the appropriate array below. |
19 |
C *==========================================================* |
C *==========================================================* |
20 |
C \ev |
C \ev |
21 |
CEOP |
CEOP |
43 |
|
|
44 |
C-- COMMON /CG3D_WK_R/ Work array common block |
C-- COMMON /CG3D_WK_R/ Work array common block |
45 |
C cg3d_q - Intermediate matrix-vector product term |
C cg3d_q - Intermediate matrix-vector product term |
46 |
C cg3d_r - " |
C cg3d_r - idem |
47 |
C cg3d_s - " |
C cg3d_s - idem |
48 |
COMMON /CG3D_WK_R/ |
COMMON /CG3D_WK_R/ |
49 |
& cg3d_q, cg3d_r, cg3d_s |
& cg3d_q, cg3d_r, cg3d_s |
50 |
_RL cg3d_q(1-1 :sNx+1 ,1-1 :sNy+1 ,Nr,nSx,nSy) |
_RL cg3d_q(1-1 :sNx+1 ,1-1 :sNy+1 ,Nr,nSx,nSy) |