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	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