/[MITgcm]/MITgcm/pkg/gmredi/gmredi_calc_uvflow.F
ViewVC logotype

Contents of /MITgcm/pkg/gmredi/gmredi_calc_uvflow.F

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

Revision 1.3 - (show annotations) (download)
Mon Mar 29 03:33:51 2004 UTC (20 years, 2 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint54d_post, checkpoint54e_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint53, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint54f_post, checkpoint55a_post, checkpoint55i_post, checkpoint57l_post, checkpoint57h_post, checkpoint57t_post, checkpoint55c_post, checkpoint57v_post, checkpoint57f_post, checkpoint53d_post, checkpoint57a_post, checkpoint57h_pre, checkpoint57x_post, checkpoint54b_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint57c_post, checkpoint55d_post, checkpoint58e_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint55j_post, checkpoint54a_post, checkpoint55h_post, checkpoint57e_post, checkpoint55b_post, checkpoint53a_post, checkpoint55f_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint56a_post, checkpoint53f_post, checkpoint53b_pre, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint58g_post, checkpoint52n_post, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint57a_pre, checkpoint57o_post, checkpoint57k_post, checkpoint53b_post, checkpoint57w_post, checkpoint58i_post, checkpoint58c_post, checkpoint53d_pre, checkpoint55e_post, checkpoint54c_post
Changes since 1.2: +4 -8 lines 
 o new "poster children" for the API reference:
   - generic_advdiff
   - mnc
1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_calc_uvflow.F,v 1.2 2002/11/14 22:43:49 heimbach Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: GMREDI_CALC_UVFLOW
8 C !INTERFACE:
9 SUBROUTINE GMREDI_CALC_UVFLOW(
10 U uTrans, vTrans,
11 I bi, bj, k, myThid)
12 C !DESCRIPTION:
13 C Add GM-bolus flow to Eulerian horizontal transport.
14
15 C !USES:
16 IMPLICIT NONE
17
18 C == GLobal variables ==
19 #include "SIZE.h"
20 #include "EEPARAMS.h"
21 #include "PARAMS.h"
22 #include "GRID.h"
23 #include "GMREDI.h"
24
25 C !INPUT/OUTPUT PARAMETERS:
26 C == Routine arguments ==
27 C bi,bj :: tile indices
28 C k :: level index
29 C myThid :: thread number
30 C uTrans :: zonal volume transport (updated)
31 C vTrans :: meridional volume transport (updated)
32 INTEGER bi, bj, k, myThid
33 _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
34 _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
35
36 #ifdef ALLOW_GMREDI
37 #ifdef GM_BOLUS_ADVEC
38
39 C !LOCAL VARIABLES:
40 C == Local variables ==
41 C i, j :: loop indices
42 INTEGER i, j
43 INTEGER kp1
44 _RL maskp1
45 CEOP
46
47 IF (GM_AdvForm .AND. .NOT.GM_AdvSeparate) THEN
48
49 kp1 = MIN(k+1,Nr)
50 maskp1 = 1.
51 IF (k.GE.Nr) maskp1 = 0.
52 DO j=1-OLy,sNy+OLy
53 DO i=1-OLx,sNx+OLx
54 uTrans(i,j) = uTrans(i,j)
55 & + dyG(i,j,bi,bj)*( GM_PsiX(i,j,kp1,bi,bj)*maskp1
56 & -GM_PsiX(i,j,k,bi,bj) )
57 & *maskW(i,j,k,bi,bj)
58 vTrans(i,j) = vTrans(i,j)
59 & + dxG(i,j,bi,bj)*( GM_PsiY(i,j,kp1,bi,bj)*maskp1
60 & -GM_PsiY(i,j,k,bi,bj) )
61 & *maskS(i,j,k,bi,bj)
62 ENDDO
63 ENDDO
64
65 ENDIF
66 #endif /* GM_BOLUS_ADVEC */
67 #endif /* ALLOW_GMREDI */
68
69 RETURN
70 END
  ViewVC Help
Powered by ViewVC 1.1.22