/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_u_metric_nh.F
ViewVC logotype

Annotation of /MITgcm/pkg/mom_fluxform/mom_u_metric_nh.F

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

Revision 1.4 - (hide annotations) (download)
Tue Nov 5 19:20:42 2002 UTC (21 years, 6 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint47j_post, checkpoint48d_pre, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint48i_post, checkpoint47f_post, checkpoint48d_post, checkpoint47c_post, checkpoint50e_post, checkpoint50c_post, checkpoint47d_post, checkpoint47a_post, checkpoint48a_post, checkpoint46n_post, checkpoint51f_pre, checkpoint48e_post, checkpoint48h_post, checkpoint50c_pre, branchpoint-genmake2, checkpoint50d_pre, checkpoint47i_post, checkpoint47h_post, checkpoint48c_post, checkpoint50b_pre, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint48g_post, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint46m_post, checkpoint51b_pre, checkpoint47g_post, checkpoint51h_pre, checkpoint50g_post, checkpoint51g_post, checkpoint51f_post, checkpoint48b_post, checkpoint50b_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47d_pre, checkpoint51d_post, checkpoint48c_pre, checkpoint47, checkpoint50h_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre
Branch point for: branch-genmake2, branch-exfmods-curt
Changes since 1.3: +2 -2 lines 
Fixed sign of non-hydrostatic metric terms of form u*w/r and v*w/r.
 - now depends on sign of gravitySign which in P coordinates should be
   -1 and otherwise +1.
1 adcroft 1.4 C $Header: /u/gcmpack/MITgcm/pkg/mom_fluxform/mom_u_metric_nh.F,v 1.3 2001/09/26 19:05:21 adcroft Exp $
2 adcroft 1.3 C $Name: $
3 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6 adcroft 1.3 CBOP
7     C !ROUTINE: MOM_U_METRIC_NH
8    
9     C !INTERFACE: ==========================================================
10 adcroft 1.2 SUBROUTINE MOM_U_METRIC_NH(
11     I bi,bj,k,
12     I uFld, wFld,
13     O uMetricTerms,
14     I myThid)
15 adcroft 1.3
16     C !DESCRIPTION:
17     C Calculates the zonal metric term due to non-hydrostaticity on the sphere:
18     C \begin{equation}
19     C \frac{u}{a} \overline{w}^{ik}
20     C \end{equation}
21    
22     C !USES: ===============================================================
23 adcroft 1.2 IMPLICIT NONE
24     #include "SIZE.h"
25     #include "EEPARAMS.h"
26     #include "PARAMS.h"
27     #include "GRID.h"
28    
29 adcroft 1.3 C !INPUT PARAMETERS: ===================================================
30     C bi,bj :: tile indices
31     C k :: vertical level
32     C uFld :: zonal flow
33     C wFld :: vertical flow
34     C myThid :: thread number
35 adcroft 1.2 INTEGER bi,bj,k
36     _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37     _RL wFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
38 adcroft 1.3 INTEGER myThid
39    
40     C !OUTPUT PARAMETERS: ==================================================
41     C uMetricTerms :: metric term
42 adcroft 1.2 _RL uMetricTerms(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43    
44 adcroft 1.3 C !LOCAL VARIABLES: ====================================================
45     C i,j :: loop indices
46     INTEGER i,j,Kp1
47 adcroft 1.2 _RL wVelBottomOverride
48 adcroft 1.3 CEOP
49 adcroft 1.2
50     Kp1=min(k+1,Nr)
51     wVelBottomOverride=1.
52     IF (k.EQ.Nr) wVelBottomOverride=0.
53    
54     DO j=1-Olx,sNy+Oly
55     DO i=1-Olx+1,sNx+Olx
56     uMetricTerms(i,j) =
57     & -1.* uFld(i,j)*recip_RSphere
58     & *0.25*(wVelBottomOverride*
59     & (wFld(i-1,j,kp1,bi,bj)+wFld(i,j,kp1,bi,bj))
60     & +wFld(i-1,j, k ,bi,bj)+wFld(i,j, k ,bi,bj)
61 adcroft 1.4 & )*rkFac*recip_horiVertRatio*gravitySign
62 adcroft 1.2 ENDDO
63     ENDDO
64    
65     RETURN
66     END
  ViewVC Help
Powered by ViewVC 1.1.22