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

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

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

Revision 1.1 - (show 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 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