/[MITgcm]/MITgcm/model/inc/CG3D.h
ViewVC logotype

Contents of /MITgcm/model/inc/CG3D.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.9 - (show annotations) (download)
Tue Dec 20 20:08:34 2005 UTC (18 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58b_post, checkpoint62c, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint62a, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint60, checkpoint61, checkpoint62, checkpoint58w_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint58r_post, checkpoint58n_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint58l_post, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint58j_post, checkpoint61q, checkpoint61e, checkpoint58i_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.8: +9 -11 lines
File MIME type: text/plain 
store solver main diagonal element in common block.
1 C $Header: /u/gcmpack/MITgcm/model/inc/CG3D.h,v 1.8 2005/11/08 06:03:21 jmc Exp $
2 C $Name: $
3
4 #ifdef ALLOW_NONHYDROSTATIC
5 CBOP
6 C !ROUTINE: CG3D.h
7 C !INTERFACE:
8 C include "CG3D.h"
9 C !DESCRIPTION: \bv
10 C *==========================================================*
11 C | CG3D.h
12 C | o Three-dimensional conjugate gradient solver header.
13 C *==========================================================*
14 C | The common blocks set up here are used in the elliptic
15 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
17 C | solver set the appropriate array below. To read the
18 C | solution read from the appropriate array below.
19 C *==========================================================*
20 C \ev
21 CEOP
22
23 C-- COMMON /CG3D_R/ DEL**2 Laplacian operators
24 C aW3d :: East-west operator.
25 C aS3d :: North-south operator.
26 C aV3d :: Vertical operator.
27 C aC3d :: 3D operator main diagonal term.
28 C zMC, zML, zMU :: preconditioner 3D solver
29 C cg3dNorm - A matrix normalisation factor.
30 COMMON /CG3D_R/
31 & aW3d, aS3d, aV3d,
32 & aC3d,
33 & zMC, zML, zMU,
34 & cg3dNorm
35 _RS aW3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
36 _RS aS3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
37 _RS aV3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
38 _RS aC3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
39 _RS zMC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40 _RS zML (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41 _RS zMU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
42 _RL cg3dNorm
43
44 C-- COMMON /CG3D_WK_R/ Work array common block
45 C cg3d_q - Intermediate matrix-vector product term
46 C cg3d_r - "
47 C cg3d_s - "
48 COMMON /CG3D_WK_R/
49 & cg3d_q, cg3d_r, cg3d_s
50 _RL cg3d_q(1-1 :sNx+1 ,1-1 :sNy+1 ,Nr,nSx,nSy)
51 _RL cg3d_r(1-1 :sNx+1 ,1-1 :sNy+1 ,Nr,nSx,nSy)
52 _RL cg3d_s(1-1 :sNx+1 ,1-1 :sNy+1 ,Nr,nSx,nSy)
53
54 #endif /* ALLOW_NONHYDROSTATIC */
  ViewVC Help
Powered by ViewVC 1.1.22