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

Annotation of /MITgcm/model/src/pressure_for_eos.F

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

Revision 1.1 - (hide annotations) (download)
Tue Feb 18 15:25:09 2003 UTC (21 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint50c_post, checkpoint52d_pre, checkpoint50c_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint54d_post, checkpoint54e_post, checkpoint51l_post, checkpoint48i_post, checkpoint50d_pre, checkpoint52k_post, checkpoint55, checkpoint54, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52f_post, checkpoint50b_pre, checkpoint54f_post, checkpoint51f_post, checkpoint51d_post, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint55c_post, checkpoint51j_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint52b_pre, checkpoint54b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint55g_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52f_pre, checkpoint55d_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, branchpoint-genmake2, checkpoint54a_post, checkpoint51r_post, checkpoint51i_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint55f_post, checkpoint52d_post, checkpoint53g_post, checkpoint50g_post, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint53f_post, checkpoint52j_post, checkpoint50e_post, branch-netcdf, checkpoint52l_post, checkpoint52n_post, checkpoint53b_pre, checkpoint51e_post, checkpoint55a_post, checkpoint49, checkpoint51o_post, checkpoint51f_pre, checkpoint48g_post, checkpoint53b_post, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint50b_post, checkpoint51m_post, checkpoint53d_pre, checkpoint55e_post, checkpoint54c_post, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-genmake2, branch-nonh, tg2-branch, ecco-branch, netcdf-sm0, checkpoint51n_branch 
o compute locally the pressure for use in EOS: UNESCO, JMD95P or MDJWF
o store total Potential in totPhyHyd for diagnostic & EOS funct. of P
o fix restart and overlap Pb when using Z-coord and EOS funct. of P
1 jmc 1.1 C $Header: $
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) = -rhoConst*rC(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
  ViewVC Help
Powered by ViewVC 1.1.22