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

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

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

Revision 1.9 - (show 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 C $Header: /u/gcmpack/MITgcm/model/src/integrate_for_w.F,v 1.8 2002/10/07 16:24:45 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: INTEGRATE_FOR_W
8 C !INTERFACE:
9 SUBROUTINE INTEGRATE_FOR_W(
10 I bi,bj,k,uFld,vFld,
11 O wFld,
12 I myThid)
13
14 C !DESCRIPTION: \bv
15 C *==========================================================*
16 C | SUBROUTINE INTEGRATE_FOR_W
17 C | o Integrate for vertical velocity.
18 C *==========================================================*
19 C \ev
20
21 C !USES:
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 "SURFACE.h"
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C == Routine arguments ==
32 C uFld, vFld :: Zonal and meridional flow
33 C wFld :: Vertical flow
34 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 C !LOCAL VARIABLES:
41 C == Local variables ==
42 C uTrans, vTrans :: Temps. for volume transports
43 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 CEOP
47
48 C-- Calculate velocity field "volume transports" through
49 C tracer cell faces.
50 DO j=1,sNy+1
51 DO i=1,sNx+1
52 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 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 DO j=1,sNy
67 DO i=1,sNx
68 wFld(i,j,k,bi,bj) = 0.
69 ENDDO
70 ENDDO
71 ELSEIF (k.eq.Nr) THEN
72 DO j=1,sNy
73 DO i=1,sNx
74 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 DO j=1,sNy
83 DO i=1,sNx
84 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 ENDDO
92 ENDIF
93 #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 ELSE
123 C- o Free Surface case (r-Coordinate):
124 C non zero at surface ; zero under-ground and at r_lower boundary
125 IF (k.eq.Nr) THEN
126 DO j=1,sNy
127 DO i=1,sNx
128 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 ENDDO
135 ELSE
136 DO j=1,sNy
137 DO i=1,sNx
138 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 C- endif - rigid-lid / Free-Surf.
148 ENDIF
149
150 RETURN
151 END
  ViewVC Help
Powered by ViewVC 1.1.22