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

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

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

Revision 1.9 - (hide annotations) (download)
Tue Dec 20 20:08:34 2005 UTC (18 years, 5 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 jmc 1.9 C $Header: /u/gcmpack/MITgcm/model/inc/CG3D.h,v 1.8 2005/11/08 06:03:21 jmc Exp $
2 adcroft 1.4 C $Name: $
3 adcroft 1.1
4     #ifdef ALLOW_NONHYDROSTATIC
5 cnh 1.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 adcroft 1.1
23     C-- COMMON /CG3D_R/ DEL**2 Laplacian operators
24 jmc 1.9 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 adcroft 1.1 C cg3dNorm - A matrix normalisation factor.
30     COMMON /CG3D_R/
31 jmc 1.9 & aW3d, aS3d, aV3d,
32     & aC3d,
33 adcroft 1.1 & zMC, zML, zMU,
34 adcroft 1.2 & cg3dNorm
35 adcroft 1.1 _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 jmc 1.9 _RS aC3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
39 adcroft 1.1 _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 adcroft 1.4 & cg3d_q, cg3d_r, cg3d_s
50 jmc 1.8 _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 adcroft 1.1
54     #endif /* ALLOW_NONHYDROSTATIC */
  ViewVC Help
Powered by ViewVC 1.1.22