1 |
cnh |
1.1 |
C $Id$ |
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 |
31 |
|
|
C $Id: $ |