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

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

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

revision 1.16 by jmc, Wed Sep 19 13:58:08 2001 UTC revision 1.17 by cnh, Wed Sep 26 18:09:13 2001 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    CBOP
6  C     /==========================================================\  C     !ROUTINE: CALC_DIV_GHAT
7  C     | S/R CALC_DIV_GHAT                                        |  C     !INTERFACE:
8  C     | o Form the right hand-side of the surface pressure eqn.  |        SUBROUTINE CALC_DIV_GHAT(
 C     |==========================================================|  
 C     \==========================================================/  
       SUBROUTINE CALC_DIV_GHAT(  
9       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
10       I        xA,yA,       I        xA,yA,
11       U        cg2d_b,       U        cg2d_b,
12       I        myThid)       I        myThid)
13    C     !DESCRIPTION: \bv
14    C     *==========================================================*
15    C     | S/R CALC_DIV_GHAT                                        
16    C     | o Form the right hand-side of the surface pressure eqn.  
17    C     *==========================================================*
18    C     | Right hand side of pressure equation is divergence
19    C     | of veclocity tendency (GHAT) term along with a relaxation
20    C     | term equal to the barotropic flow field divergence.
21    C     *==========================================================*
22    C     \ev
23    
24    C     !USES:
25        IMPLICIT NONE        IMPLICIT NONE
   
26  C     == Global variables ==  C     == Global variables ==
27  #include "SIZE.h"  #include "SIZE.h"
28  #include "DYNVARS.h"  #include "DYNVARS.h"
# Line 25  C     == Global variables == Line 32  C     == Global variables ==
32  #include "GRID.h"  #include "GRID.h"
33  #include "SOLVE_FOR_PRESSURE3D.h"  #include "SOLVE_FOR_PRESSURE3D.h"
34    
35    C     !INPUT/OUTPUT PARAMETERS:
36  C     == Routine arguments ==  C     == Routine arguments ==
37  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
38  C                                      results will be set.  C                                      results will be set.
# Line 39  C     myThid - Instance number for this Line 47  C     myThid - Instance number for this
47        _RL cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL cg2d_b(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48        INTEGER myThid        INTEGER myThid
49    
50    C     !LOCAL VARIABLES:
51  C     == Local variables ==  C     == Local variables ==
52    C     i,j :: Loop counters
53    C     pf  :: Intermediate array for building RHS source term.
54        INTEGER i,j        INTEGER i,j
55        _RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL pf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
56    CEOP
57    
58  C--   Pressure equation source term  C--   Pressure equation source term
59  C     Term is the vertical integral of the divergence of the  C     Term is the vertical integral of the divergence of the
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22