model/src/solve_for_pressure.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/solve_for_pressure.F,v 1.4 1998/06/08 21:43:02 cnh Exp $

#include "CPP_OPTIONS.h"

CStartOfInterface
      SUBROUTINE SOLVE_FOR_PRESSURE( myThid )
C     /==========================================================\
C     | SUBROUTINE SOLVE_FOR_PRESSURE                            |
C     | o Controls inversion of two and/or three-dimensional     |
C     |   elliptic problems for the pressure field.              |
C     \==========================================================/

C     == Global variables
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "CG2D.h"

C     == Routine arguments ==
C     myThid - Number of this instance of SOLVE_FOR_PRESSURE
      INTEGER myThid
CEndOfInterface

C     Local variables
      INTEGER i,j,bi,bj

#ifdef ALLOW_CD
C--   Save previous solution.
      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          cg2d_xNM1(i,j,bi,bj) = cg2d_x(i,j,bi,bj)
         ENDDO
        ENDDO
       ENDDO
      ENDDO
#endif

C--   Find the surface pressure using a two-dimensional conjugate
C--   gradient solver.
C     see CG2D.h for the interface to this routine.
      CALL CG2D(
     I           myThid )

C--   Solve for a three-dimensional pressure term (NH or IGW or both ).
C     see CG3D.h for the interface to this routine.
C*P*  CALL CG3D(
C*P* I           myThid )

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/solve_for_pressure.F,v 1.4 1998/06/08 21:43:02 cnh Exp $
2
3	#include "CPP_OPTIONS.h"
4
5	CStartOfInterface
6	SUBROUTINE SOLVE_FOR_PRESSURE( myThid )
7	C /==========================================================\
8	C \| SUBROUTINE SOLVE_FOR_PRESSURE \|
9	C \| o Controls inversion of two and/or three-dimensional \|
10	C \| elliptic problems for the pressure field. \|
11	C \==========================================================/
12
13	C == Global variables
14	#include "SIZE.h"
15	#include "EEPARAMS.h"
16	#include "PARAMS.h"
17	#include "DYNVARS.h"
18	#include "CG2D.h"
19
20	C == Routine arguments ==
21	C myThid - Number of this instance of SOLVE_FOR_PRESSURE
22	INTEGER myThid
23	CEndOfInterface
24
25	C Local variables
26	INTEGER i,j,bi,bj
27
28	#ifdef ALLOW_CD
29	C-- Save previous solution.
30	DO bj=myByLo(myThid),myByHi(myThid)
31	DO bi=myBxLo(myThid),myBxHi(myThid)
32	DO j=1-OLy,sNy+OLy
33	DO i=1-OLx,sNx+OLx
34	cg2d_xNM1(i,j,bi,bj) = cg2d_x(i,j,bi,bj)
35	ENDDO
36	ENDDO
37	ENDDO
38	ENDDO
39	#endif
40
41	C-- Find the surface pressure using a two-dimensional conjugate
42	C-- gradient solver.
43	C see CG2D.h for the interface to this routine.
44	CALL CG2D(
45	I myThid )
46
47	C-- Solve for a three-dimensional pressure term (NH or IGW or both ).
48	C see CG3D.h for the interface to this routine.
49	CP CALL CG3D(
50	CP I myThid )
51
52	RETURN
53	END