model/src/calc_grad_phi_surf.F

C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_grad_eta_surf.F,v 1.6 2001/03/06 16:51:02 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

C     /==========================================================\
C     | S/R CALC_GRAD_PHI_SURF                                   |
C     | o Calculate the gradient of the surface Potential anomaly|
C     \==========================================================/
      SUBROUTINE CALC_GRAD_PHI_SURF( bi, bj, iMin, iMax, jMin, jMax,
     I                       etaFld,
     O                       phiSurfX, phiSurfY,
     I                       myThid )
      IMPLICIT NONE
C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SURFACE.h"

C     == Routine Arguments ==
C     etaFld             - free-surface r-anomaly (r unit).
C     phiSurfX, phiSurfY - Gradient in the X and Y directions of surface
C       Potentiel anomaly (atmos: =Geopotential ; ocean: =Pressure/rho) 
C     bi,bj,iMin,iMax,jMin,jMax - Loop counters
C     myThid - Instance number for this call of the routine.
      INTEGER bi,bj,iMin,iMax,jMin,jMax
      _RL etaFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL phiSurfX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      _RL phiSurfY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      INTEGER myThid

C     == Local variables ==
      INTEGER i,j

C     Zonal term
      DO j=jMin,jMax
       DO i=iMin,iMax
        phiSurfX(i,j)=_recip_dxC(i,j,bi,bj)*
     &   ( Bo_surf(i,j,bi,bj)*etaFld(i,j,bi,bj)
     &   - Bo_surf(i-1,j,bi,bj)*etaFld(i-1,j,bi,bj) )      
       ENDDO
      ENDDO

C     Meridional term
      DO j=jMin,jMax
       DO i=iMin,iMax
        phiSurfY(i,j)=_recip_dyC(i,j,bi,bj)*
     &   ( Bo_surf(i,j,bi,bj)*etaFld(i,j,bi,bj)
     &   - Bo_surf(i,j-1,bi,bj)*etaFld(i,j-1,bi,bj) ) 
       ENDDO
      ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_grad_eta_surf.F,v 1.6 2001/03/06 16:51:02 jmc Exp $
2	C $Name: $
3
4	#include "CPP_OPTIONS.h"
5
6	C /==========================================================\
7	C \| S/R CALC_GRAD_PHI_SURF \|
8	C \| o Calculate the gradient of the surface Potential anomaly\|
9	C \==========================================================/
10	SUBROUTINE CALC_GRAD_PHI_SURF( bi, bj, iMin, iMax, jMin, jMax,
11	I etaFld,
12	O phiSurfX, phiSurfY,
13	I myThid )
14	IMPLICIT NONE
15	C == Global variables ==
16	#include "SIZE.h"
17	#include "EEPARAMS.h"
18	#include "PARAMS.h"
19	#include "GRID.h"
20	#include "SURFACE.h"
21
22	C == Routine Arguments ==
23	C etaFld - free-surface r-anomaly (r unit).
24	C phiSurfX, phiSurfY - Gradient in the X and Y directions of surface
25	C Potentiel anomaly (atmos: =Geopotential ; ocean: =Pressure/rho)
26	C bi,bj,iMin,iMax,jMin,jMax - Loop counters
27	C myThid - Instance number for this call of the routine.
28	INTEGER bi,bj,iMin,iMax,jMin,jMax
29	_RL etaFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
30	_RL phiSurfX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31	_RL phiSurfY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
32	INTEGER myThid
33
34	C == Local variables ==
35	INTEGER i,j
36
37	C Zonal term
38	DO j=jMin,jMax
39	DO i=iMin,iMax
40	phiSurfX(i,j)=_recip_dxC(i,j,bi,bj)*
41	& ( Bo_surf(i,j,bi,bj)*etaFld(i,j,bi,bj)
42	& - Bo_surf(i-1,j,bi,bj)*etaFld(i-1,j,bi,bj) )
43	ENDDO
44	ENDDO
45
46	C Meridional term
47	DO j=jMin,jMax
48	DO i=iMin,iMax
49	phiSurfY(i,j)=_recip_dyC(i,j,bi,bj)*
50	& ( Bo_surf(i,j,bi,bj)*etaFld(i,j,bi,bj)
51	& - Bo_surf(i,j-1,bi,bj)*etaFld(i,j-1,bi,bj) )
52	ENDDO
53	ENDDO
54
55	RETURN
56	END