/[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.7 - (show annotations) (download)
Fri Oct 31 20:35:32 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint57g_pre, checkpoint57s_post, checkpoint57b_post, checkpoint52d_pre, checkpoint57g_post, checkpoint56b_post, checkpoint52j_pre, checkpoint54d_post, checkpoint54e_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint52l_post, checkpoint52k_post, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52f_post, checkpoint57n_post, checkpoint54f_post, checkpoint51t_post, checkpoint55i_post, checkpoint57l_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint53d_post, checkpoint57a_post, checkpoint57h_pre, checkpoint52b_pre, checkpoint54b_post, checkpoint57h_post, checkpoint52m_post, checkpoint55g_post, checkpoint52b_post, checkpoint52c_post, checkpoint57c_post, checkpoint52f_pre, checkpoint55d_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint55j_post, checkpoint54a_post, checkpoint55h_post, checkpoint51r_post, checkpoint57e_post, checkpoint55b_post, checkpoint53a_post, checkpoint55f_post, checkpoint52d_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint52a_pre, checkpoint52i_post, checkpoint52h_pre, checkpoint56a_post, checkpoint53f_post, checkpoint57h_done, checkpoint52j_post, checkpoint57j_post, checkpoint57f_pre, branch-netcdf, checkpoint52n_post, checkpoint53b_pre, checkpoint56c_post, checkpoint57a_pre, checkpoint55a_post, checkpoint57o_post, checkpoint57k_post, checkpoint53b_post, checkpoint52a_post, checkpoint57w_post, ecco_c52_e35, checkpoint53d_pre, checkpoint55e_post, checkpoint54c_post, checkpoint51u_post
Branch point for: branch-nonh, netcdf-sm0
Changes since 1.6: +1 -3 lines
File MIME type: text/plain 
 o remove all '#include "PACACKAGES_CONFIG.h"' from model/inc/* and cleanup
   the verification tests that this breaks
 o this was confirmed to work for the basic tests ("testreport -ieee") on
   shelley
1 C $Header: /u/u3/gcmpack/MITgcm/model/inc/CG3D.h,v 1.6 2003/10/09 04:19:18 edhill 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 etaNBuf - Shared buffers for accumulating con. grad vector product.
28 C alphaBuf
29 C errBuf
30 C nrmBuf
31 C cg3dNorm - A matrix normalisation factor.
32 COMMON /CG3D_R/
33 & aW3d,
34 & aS3d,
35 & aV3d,
36 & zMC, zML, zMU,
37 & cg3dNorm
38 _RS aW3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
39 _RS aS3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40 _RS aV3d (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41 _RS zMC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
42 _RS zML (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43 _RS zMU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
44 _RL cg3dNorm
45
46 C-- COMMON /CG3D_WK_R/ Work array common block
47 C cg3d_q - Intermediate matrix-vector product term
48 C cg3d_r - "
49 C cg3d_s - "
50 C cg3d_x Solution vector
51 C cg3d_b Right-hand side vector
52 COMMON /CG3D_WK_R/
53 & cg3d_q, cg3d_r, cg3d_s
54 _RL cg3d_q(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
55 _RL cg3d_r(1-1 :sNx+1 ,1-1 :sNy+1 ,nR,nSx,nSy)
56 _RL cg3d_s(1-1 :sNx+1 ,1-1 :sNy+1 ,nR,nSx,nSy)
57
58 #endif /* ALLOW_NONHYDROSTATIC */
  ViewVC Help
Powered by ViewVC 1.1.22