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

Annotation of /MITgcm/pkg/gmredi/gmredi_calc_wflow.F

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

Revision 1.2 - (hide annotations) (download)
Thu Nov 14 22:43:49 2002 UTC (21 years, 7 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint47e_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint47c_post, checkpoint50c_post, checkpoint52d_pre, checkpoint48e_post, checkpoint50g_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint48c_post, checkpoint51l_post, checkpoint48i_post, checkpoint50d_pre, checkpoint52k_post, checkpoint51, checkpoint50, checkpoint52, checkpoint50d_post, checkpoint52f_post, checkpoint50b_pre, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint52b_pre, checkpoint51l_pre, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint52h_pre, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, hrcube_1, checkpoint52f_pre, checkpoint47j_post, checkpoint47a_post, branch-exfmods-tag, branchpoint-genmake2, checkpoint51r_post, checkpoint51i_post, checkpoint51b_post, checkpoint51c_post, checkpoint47b_post, checkpoint52d_post, checkpoint50c_pre, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint48a_post, checkpoint52j_post, checkpoint47f_post, checkpoint50e_post, branch-netcdf, checkpoint52l_post, checkpoint51e_post, checkpoint47d_pre, checkpoint47, checkpoint48, checkpoint49, checkpoint51o_post, checkpoint51f_pre, checkpoint47h_post, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint50b_post, checkpoint51m_post, checkpoint51a_post, checkpoint51p_post, checkpoint48g_post, checkpoint51u_post
Branch point for: branch-exfmods-curt, branch-genmake2, branch-nonh, tg2-branch, ecco-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.1: +2 -2 lines 
o * "clean" adjoint code (in terms of extensive recomputations)
    can now be obtained for all GMREDI options (i.e. for
    - GM_VISBECK_VARIABLE_K
    - GM_NON_UNITY_DIAGONAL
    - GM_EXTRA_DIAGONAL
    - GM_BOLUS_ADVEC )
  * However, wrong gradient check problem remains unsolved.
  * New CPP options have been introduced for different
    tapering schemes
1 heimbach 1.2 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_calc_wflow.F,v 1.1 2002/03/06 01:42:37 jmc Exp $
2     C $Name: checkpoint46 $
3 jmc 1.1
4     #include "GMREDI_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: GMREDI_CALC_WFLOW
8     C !INTERFACE:
9     SUBROUTINE GMREDI_CALC_WFLOW(
10     U rTrans,
11     I bi, bj, k, myThid)
12     C !DESCRIPTION: \bv
13     C *==========================================================*
14     C | o SUBROUTINE GMREDI_CALC_WFLOW
15     C | Add GM-bolus flow to Eulerian vertical transport
16     C *==========================================================*
17     C \ev
18    
19     C !USES:
20     IMPLICIT NONE
21    
22     C == GLobal variables ==
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26     #include "GRID.h"
27     #include "GMREDI.h"
28    
29     C !INPUT/OUTPUT PARAMETERS:
30     C == Routine arguments ==
31     C bi,bj :: tile indices
32     C k :: level index
33     C myThid :: thread number
34     C rTrans :: vertical volume transport (updated)
35     INTEGER bi, bj, k, myThid
36     _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37    
38     #ifdef ALLOW_GMREDI
39     #ifdef GM_BOLUS_ADVEC
40    
41     C !LOCAL VARIABLES:
42     C == Local variables ==
43     C i, j :: loop indices
44     INTEGER i, j
45     CEOP
46    
47     IF (GM_AdvForm .AND. .NOT.GM_AdvSeparate) THEN
48    
49     DO j=1-Oly,sNy+Oly-1
50     DO i=1-Olx,sNx+Olx-1
51     rTrans(i,j) = rTrans(i,j)
52     & +( dyG(i+1,j,bi,bj)*GM_PsiX(i+1,j,k,bi,bj)
53     & -dyG( i ,j,bi,bj)*GM_PsiX( i ,j,k,bi,bj)
54     & +dxG(i,j+1,bi,bj)*GM_PsiY(i,j+1,k,bi,bj)
55     & -dxG(i, j ,bi,bj)*GM_PsiY(i, j ,k,bi,bj)
56     & )*maskC(i,j,k,bi,bj)
57     ENDDO
58     ENDDO
59    
60     ENDIF
61     #endif /* GM_BOLUS_ADVEC */
62     #endif /* ALLOW_GMREDI */
63    
64     RETURN
65     END
  ViewVC Help
Powered by ViewVC 1.1.22