model/src/pressure_for_eos.F

C $Header: /u/gcmpack/MITgcm/model/src/pressure_for_eos.F,v 1.1 2003/02/18 15:25:09 jmc Exp $
C $Name:  $
#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: PRESSURE_FOR_EOS
C     !INTERFACE:
      SUBROUTINE PRESSURE_FOR_EOS( 
     I        bi, bj, iMin, iMax, jMin, jMax,  k,
     O        locPres,
     I        myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE PRESSURE_FOR_EOS                                      
C     | o Provide a local copy of the total pressure 
C     |   at cell center (level k) for use in EOS funct. of P
C     *==========================================================*
C     \ev

C     !USES:

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

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     bi,bj, k  :: tile and level indices
C     iMin,iMax,jMin,jMax :: computational domain
C     myThid - Thread number for this instance of the routine.
      INTEGER bi, bj, k
      INTEGER iMin,iMax,jMin,jMax
      _RL locPres(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER myThid

C     !LOCAL VARIABLES:
C     == Local variables ==
C     i,j :: loop index
      INTEGER  i,j
CEOP

C
C     provide the pressure for use in the equation of state
C
      IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN
C     in Z coordinates the pressure is rho0 * (hydrostatic) Potential
        IF ( useDynP_inEos_Zc ) THEN  
C----------
C     NOTE: For now, totPhiHyd only contains the Potential anomaly
C           since PhiRef is not available for Atmos and has not (yet) 
C           been added in S/R DIAGS_PHI_HYD 
C----------
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
            locPres(i,j) = rhoConst*(
     &                   totPhiHyd(i,j,k,bi,bj)
     &                  -rC(k)*gravity 
     &                            )*maskC(i,j,k,bi,bj)
          ENDDO
         ENDDO
       ELSE
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
            locPres(i,j) = -rhoConst*rC(k)*gravity 
     &                              *maskC(i,j,k,bi,bj)
          ENDDO
         ENDDO
       ENDIF
      ELSEIF ( buoyancyRelation .EQ. 'ATMOSPHERIC'
     &    .OR. buoyancyRelation .EQ. 'OCEANICP' ) THEN
C     in P coordinates the pressure is just the coordinate of
C     the tracer point
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx
            locPres(i,j) = rC(k)
     &                   * maskC(i,j,k,bi,bj)
          ENDDO
         ENDDO
      ENDIF

      RETURN 
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/pressure_for_eos.F,v 1.1 2003/02/18 15:25:09 jmc Exp $
2	C $Name: $
3	#include "CPP_OPTIONS.h"
4
5	CBOP
6	C !ROUTINE: PRESSURE_FOR_EOS
7	C !INTERFACE:
8	SUBROUTINE PRESSURE_FOR_EOS(
9	I bi, bj, iMin, iMax, jMin, jMax, k,
10	O locPres,
11	I myThid )
12	C !DESCRIPTION: \bv
13	C ==========================================================
14	C \| SUBROUTINE PRESSURE_FOR_EOS
15	C \| o Provide a local copy of the total pressure
16	C \| at cell center (level k) for use in EOS funct. of P
17	C ==========================================================
18	C \ev
19
20	C !USES:
21
22	IMPLICIT NONE
23	C == Global variables ==
24	#include "SIZE.h"
25	#include "EEPARAMS.h"
26	#include "PARAMS.h"
27	#include "GRID.h"
28	#include "DYNVARS.h"
29	#include "SURFACE.h"
30
31	C !INPUT/OUTPUT PARAMETERS:
32	C == Routine arguments ==
33	C bi,bj, k :: tile and level indices
34	C iMin,iMax,jMin,jMax :: computational domain
35	C myThid - Thread number for this instance of the routine.
36	INTEGER bi, bj, k
37	INTEGER iMin,iMax,jMin,jMax
38	_RL locPres(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39	INTEGER myThid
40
41	C !LOCAL VARIABLES:
42	C == Local variables ==
43	C i,j :: loop index
44	INTEGER i,j
45	CEOP
46
47	C
48	C provide the pressure for use in the equation of state
49	C
50	IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN
51	C in Z coordinates the pressure is rho0 * (hydrostatic) Potential
52	IF ( useDynP_inEos_Zc ) THEN
53	C----------
54	C NOTE: For now, totPhiHyd only contains the Potential anomaly
55	C since PhiRef is not available for Atmos and has not (yet)
56	C been added in S/R DIAGS_PHI_HYD
57	C----------
58	DO j=1-Oly,sNy+Oly
59	DO i=1-Olx,sNx+Olx
60	locPres(i,j) = rhoConst*(
61	& totPhiHyd(i,j,k,bi,bj)
62	& -rC(k)*gravity
63	& )*maskC(i,j,k,bi,bj)
64	ENDDO
65	ENDDO
66	ELSE
67	DO j=1-Oly,sNy+Oly
68	DO i=1-Olx,sNx+Olx
69	locPres(i,j) = -rhoConstrC(k)gravity
70	& *maskC(i,j,k,bi,bj)
71	ENDDO
72	ENDDO
73	ENDIF
74	ELSEIF ( buoyancyRelation .EQ. 'ATMOSPHERIC'
75	& .OR. buoyancyRelation .EQ. 'OCEANICP' ) THEN
76	C in P coordinates the pressure is just the coordinate of
77	C the tracer point
78	DO j=1-Oly,sNy+Oly
79	DO i=1-Olx,sNx+Olx
80	locPres(i,j) = rC(k)
81	& * maskC(i,j,k,bi,bj)
82	ENDDO
83	ENDDO
84	ENDIF
85
86	RETURN
87	END