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

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

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

Revision 1.9 - (hide annotations) (download)
Sun Jan 26 21:06:11 2003 UTC (21 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint57m_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint57g_pre, checkpoint50c_post, checkpoint57s_post, checkpoint57b_post, checkpoint52d_pre, checkpoint57g_post, checkpoint48e_post, checkpoint56b_post, checkpoint50c_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint54d_post, checkpoint54e_post, checkpoint51l_post, checkpoint48i_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint52l_post, checkpoint52k_post, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52f_post, checkpoint57n_post, checkpoint50b_pre, checkpoint54f_post, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint51t_post, checkpoint51n_post, checkpoint55i_post, checkpoint57l_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint48d_pre, checkpoint51j_post, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint57a_post, checkpoint48d_post, checkpoint57h_pre, checkpoint48f_post, checkpoint52b_pre, checkpoint54b_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52f_pre, checkpoint55d_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint55j_post, branchpoint-genmake2, checkpoint54a_post, checkpoint55h_post, checkpoint51r_post, checkpoint51i_post, checkpoint57e_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint55f_post, checkpoint52d_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint50g_post, checkpoint57q_post, eckpoint57e_pre, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint56a_post, checkpoint53f_post, checkpoint57h_done, checkpoint52j_post, checkpoint50e_post, checkpoint57j_post, checkpoint57f_pre, branch-netcdf, checkpoint50d_pre, checkpoint52n_post, checkpoint53b_pre, checkpoint56c_post, checkpoint51e_post, checkpoint57a_pre, checkpoint55a_post, checkpoint49, checkpoint57o_post, checkpoint51o_post, checkpoint57k_post, checkpoint51f_pre, checkpoint48g_post, checkpoint53b_post, checkpoint52a_post, checkpoint57w_post, checkpoint51g_post, ecco_c52_e35, checkpoint57x_post, 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, netcdf-sm0, checkpoint51n_branch
Changes since 1.8: +32 -2 lines 
r* coordinate added in #ifdef NONLIN_FRSURF block.
 (modification to pressure gradient not yet implemented)
1 jmc 1.9 C $Header: /u/gcmpack/MITgcm/model/src/integrate_for_w.F,v 1.8 2002/10/07 16:24:45 jmc Exp $
2 cnh 1.4 C $Name: $
3 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6 cnh 1.7 CBOP
7     C !ROUTINE: INTEGRATE_FOR_W
8     C !INTERFACE:
9 adcroft 1.2 SUBROUTINE INTEGRATE_FOR_W(
10     I bi,bj,k,uFld,vFld,
11     O wFld,
12     I myThid)
13    
14 cnh 1.7 C !DESCRIPTION: \bv
15     C *==========================================================*
16 jmc 1.8 C | SUBROUTINE INTEGRATE_FOR_W
17 cnh 1.7 C | o Integrate for vertical velocity.
18     C *==========================================================*
19     C \ev
20    
21     C !USES:
22 adcroft 1.2 IMPLICIT NONE
23     C == GLobal variables ==
24     #include "SIZE.h"
25     #include "EEPARAMS.h"
26     #include "PARAMS.h"
27     #include "GRID.h"
28 jmc 1.9 #include "SURFACE.h"
29 adcroft 1.2
30 cnh 1.7 C !INPUT/OUTPUT PARAMETERS:
31 adcroft 1.2 C == Routine arguments ==
32 cnh 1.7 C uFld, vFld :: Zonal and meridional flow
33     C wFld :: Vertical flow
34 adcroft 1.2 INTEGER bi,bj,k
35     _RL uFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
36     _RL vFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
37     _RL wFld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
38     INTEGER myThid
39    
40 cnh 1.7 C !LOCAL VARIABLES:
41 adcroft 1.2 C == Local variables ==
42 cnh 1.7 C uTrans, vTrans :: Temps. for volume transports
43 adcroft 1.2 INTEGER i,j
44     _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45     _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46 cnh 1.7 CEOP
47 adcroft 1.2
48     C-- Calculate velocity field "volume transports" through
49     C tracer cell faces.
50 jmc 1.8 DO j=1,sNy+1
51     DO i=1,sNx+1
52 adcroft 1.2 uTrans(i,j) = uFld(i,j,k,bi,bj)*
53     & _dyG(i,j,bi,bj)
54     & *drF(k)*_hFacW(i,j,k,bi,bj)
55     vTrans(i,j) = vFld(i,j,k,bi,bj)*
56     & _dxG(i,j,bi,bj)
57     & *drF(k)*_hFacS(i,j,k,bi,bj)
58     ENDDO
59     ENDDO
60    
61 adcroft 1.6 C-- Calculate vertical "volume transport" through face k
62     C between tracer cell k-1 & k
63     IF (rigidLid) THEN
64     C- o Rigid-Lid case: zero at lower and upper boundaries
65     IF (k.eq.1) THEN
66 jmc 1.8 DO j=1,sNy
67     DO i=1,sNx
68 adcroft 1.6 wFld(i,j,k,bi,bj) = 0.
69     ENDDO
70 adcroft 1.2 ENDDO
71 adcroft 1.6 ELSEIF (k.eq.Nr) THEN
72 jmc 1.8 DO j=1,sNy
73     DO i=1,sNx
74 adcroft 1.6 wFld(i,j,k,bi,bj) =
75     & -( uTrans(i+1,j)-uTrans(i,j)
76     & +vTrans(i,j+1)-vTrans(i,j)
77     & )*recip_rA(i,j,bi,bj)
78     & *maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
79     ENDDO
80     ENDDO
81     ELSE
82 jmc 1.8 DO j=1,sNy
83     DO i=1,sNx
84 adcroft 1.6 wFld(i,j,k,bi,bj) =
85     & ( wFld(i,j,k+1,bi,bj)
86     & -( uTrans(i+1,j)-uTrans(i,j)
87     & +vTrans(i,j+1)-vTrans(i,j)
88     & )*recip_rA(i,j,bi,bj)
89     & )*maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
90     ENDDO
91 adcroft 1.2 ENDDO
92 adcroft 1.6 ENDIF
93 jmc 1.9 #ifdef NONLIN_FRSURF
94     ELSEIF (select_rStar .NE. 0) THEN
95     C- o rStar case: zero under-ground and at r_lower boundary
96     C can be non-zero at surface (useRealFreshWaterFlux)
97     IF (k.eq.Nr) THEN
98     DO j=1,sNy
99     DO i=1,sNx
100     wFld(i,j,k,bi,bj) =
101     & ( -( uTrans(i+1,j)-uTrans(i,j)
102     & +vTrans(i,j+1)-vTrans(i,j)
103     & )*recip_rA(i,j,bi,bj)
104     & - rStarDhCDt(i,j,bi,bj)*drF(k)*h0FacC(i,j,k,bi,bj)
105     & )*maskC(i,j,k,bi,bj)
106     ENDDO
107     ENDDO
108     ELSE
109     DO j=1,sNy
110     DO i=1,sNx
111     wFld(i,j,k,bi,bj) =
112     & ( wFld(i,j,k+1,bi,bj)
113     & -( uTrans(i+1,j)-uTrans(i,j)
114     & +vTrans(i,j+1)-vTrans(i,j)
115     & )*recip_rA(i,j,bi,bj)
116     & - rStarDhCDt(i,j,bi,bj)*drF(k)*h0FacC(i,j,k,bi,bj)
117     & )*maskC(i,j,k,bi,bj)
118     ENDDO
119     ENDDO
120     ENDIF
121     #endif /* NONLIN_FRSURF */
122 adcroft 1.2 ELSE
123 jmc 1.9 C- o Free Surface case (r-Coordinate):
124 adcroft 1.6 C non zero at surface ; zero under-ground and at r_lower boundary
125     IF (k.eq.Nr) THEN
126 jmc 1.8 DO j=1,sNy
127     DO i=1,sNx
128 adcroft 1.6 wFld(i,j,k,bi,bj) =
129     & -( uTrans(i+1,j)-uTrans(i,j)
130     & +vTrans(i,j+1)-vTrans(i,j)
131     & )*recip_rA(i,j,bi,bj)
132     & *maskC(i,j,k,bi,bj)
133     ENDDO
134 adcroft 1.2 ENDDO
135 adcroft 1.6 ELSE
136 jmc 1.8 DO j=1,sNy
137     DO i=1,sNx
138 adcroft 1.6 wFld(i,j,k,bi,bj) =
139     & ( wFld(i,j,k+1,bi,bj)
140     & -( uTrans(i+1,j)-uTrans(i,j)
141     & +vTrans(i,j+1)-vTrans(i,j)
142     & )*recip_rA(i,j,bi,bj)
143     & )*maskC(i,j,k,bi,bj)
144     ENDDO
145     ENDDO
146     ENDIF
147 jmc 1.8 C- endif - rigid-lid / Free-Surf.
148 adcroft 1.2 ENDIF
149    
150     RETURN
151     END
  ViewVC Help
Powered by ViewVC 1.1.22