1 |
C $Header: /u/gcmpack/models/MITgcmUV/model/src/solve_for_pressure.F,v 1.2 1998/04/24 02:05:42 cnh Exp $ |
2 |
|
3 |
#include "CPP_EEOPTIONS.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 == Routine arguments == |
14 |
C myThid - Number of this instance of SOLVE_FOR_PRESSURE |
15 |
INTEGER myThid |
16 |
CEndOfInterface |
17 |
|
18 |
C-- Find the surface pressure using a two-dimensional conjugate |
19 |
C-- gradient solver. |
20 |
C see CG2D.h for the interface to this routine. |
21 |
CALL CG2D( |
22 |
I myThid ) |
23 |
|
24 |
C-- Solve for a three-dimensional pressure term (NH or IGW or both ). |
25 |
C see CG3D.h for the interface to this routine. |
26 |
C*P* CALL CG3D( |
27 |
C*P* I myThid ) |
28 |
|
29 |
RETURN |
30 |
END |