model/src/diags_phi_hyd.F

C $Header:  $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: DIAGS_PHI_HYD
C     !INTERFACE:
      SUBROUTINE DIAGS_PHI_HYD( 
     I                       k, bi, bj, iMin,iMax, jMin,jMax,
     I                       phiHydC,
     I                       myTime, myIter, myThid)
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | S/R DIAGS_PHI_HYD                                    
C     | o Diagnose full hydrostatic Potential at cell center ; 
C     |   used for output & with EOS funct. of P
C     *==========================================================*
C     | NOTE: For now, only contains the (total) Potential anomaly
C     |  since phiRef (for Atmos) is not available (not in common)
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SURFACE.h"
#include "DYNVARS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine Arguments ==
C     k, bi,bj      :: level & tile indices
C     iMin,iMax,jMin,jMax :: Loop counters
C     phiHydC    :: hydrostatic potential anomaly at cell center
C                  (atmos: =Geopotential ; ocean-z: =Pressure/rho)
C     myTime :: Current time
C     myIter :: Current iteration number
C     myThid :: Instance number for this call of the routine.
      INTEGER k, bi,bj, iMin,iMax, jMin,jMax
      _RL phiHydC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL myTime
      INTEGER myIter, myThid

#ifdef INCLUDE_PHIHYD_CALCULATION_CODE

C     !LOCAL VARIABLES:
C     == Local variables ==
C     i,j :: Loop counters
      INTEGER i,j
CEOP
        DO j=jMin,jMax
         DO i=iMin,iMax
           totPhiHyd(i,j,k,bi,bj) = phiHydC(i,j)
c    &            - gravity*rC(k)*recip_horiVertRatio
     &            + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
     &            + phi0surf(i,j,bi,bj) 
         ENDDO
        ENDDO

#endif /* INCLUDE_PHIHYD_CALCULATION_CODE */

      RETURN
      END
1	jmc	1.1	C $Header: $
2			C $Name: $
3
4			#include "CPP_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: DIAGS_PHI_HYD
8			C !INTERFACE:
9			SUBROUTINE DIAGS_PHI_HYD(
10			I k, bi, bj, iMin,iMax, jMin,jMax,
11			I phiHydC,
12			I myTime, myIter, myThid)
13			C !DESCRIPTION: \bv
14			C ==========================================================
15			C \| S/R DIAGS_PHI_HYD
16			C \| o Diagnose full hydrostatic Potential at cell center ;
17			C \| used for output & with EOS funct. of P
18			C ==========================================================
19			C \| NOTE: For now, only contains the (total) Potential anomaly
20			C \| since phiRef (for Atmos) is not available (not in common)
21			C ==========================================================
22			C \ev
23
24			C !USES:
25			IMPLICIT NONE
26			C == Global variables ==
27			#include "SIZE.h"
28			#include "EEPARAMS.h"
29			#include "PARAMS.h"
30			#include "GRID.h"
31			#include "SURFACE.h"
32			#include "DYNVARS.h"
33
34			C !INPUT/OUTPUT PARAMETERS:
35			C == Routine Arguments ==
36			C k, bi,bj :: level & tile indices
37			C iMin,iMax,jMin,jMax :: Loop counters
38			C phiHydC :: hydrostatic potential anomaly at cell center
39			C (atmos: =Geopotential ; ocean-z: =Pressure/rho)
40			C myTime :: Current time
41			C myIter :: Current iteration number
42			C myThid :: Instance number for this call of the routine.
43			INTEGER k, bi,bj, iMin,iMax, jMin,jMax
44			_RL phiHydC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45			_RL myTime
46			INTEGER myIter, myThid
47
48			#ifdef INCLUDE_PHIHYD_CALCULATION_CODE
49
50			C !LOCAL VARIABLES:
51			C == Local variables ==
52			C i,j :: Loop counters
53			INTEGER i,j
54			CEOP
55			DO j=jMin,jMax
56			DO i=iMin,iMax
57			totPhiHyd(i,j,k,bi,bj) = phiHydC(i,j)
58			c & - gravityrC(k)recip_horiVertRatio
59			& + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
60			& + phi0surf(i,j,bi,bj)
61			ENDDO
62			ENDDO
63
64			#endif /* INCLUDE_PHIHYD_CALCULATION_CODE */
65
66			RETURN
67			END