/[MITgcm]/MITgcm/model/src/solve_for_pressure.F
ViewVC logotype

Diff of /MITgcm/model/src/solve_for_pressure.F

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

revision 1.3 by cnh, Fri Apr 24 02:10:21 1998 UTC revision 1.4 by cnh, Mon Jun 8 21:43:02 1998 UTC
# Line 10  C     | o Controls inversion of two and/ Line 10  C     | o Controls inversion of two and/
10  C     |   elliptic problems for the pressure field.              |  C     |   elliptic problems for the pressure field.              |
11  C     \==========================================================/  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 ==  C     == Routine arguments ==
21  C     myThid - Number of this instance of SOLVE_FOR_PRESSURE  C     myThid - Number of this instance of SOLVE_FOR_PRESSURE
22        INTEGER myThid        INTEGER myThid
23  CEndOfInterface  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  C--   Find the surface pressure using a two-dimensional conjugate
42  C--   gradient solver.  C--   gradient solver.
43  C     see CG2D.h for the interface to this routine.  C     see CG2D.h for the interface to this routine.
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22