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

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

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

Revision 1.4 - (show annotations) (download)
Sun Dec 16 18:54:49 2001 UTC (22 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint44e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, checkpoint46l_pre, chkpt44d_post, checkpoint44e_pre, checkpoint46d_pre, release1-branch_tutorials, checkpoint45d_post, checkpoint46j_pre, chkpt44a_post, checkpoint44h_pre, checkpoint46a_post, checkpoint46j_post, checkpoint46k_post, chkpt44c_pre, checkpoint45a_post, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, chkpt44a_pre, checkpoint46i_post, checkpoint46c_post, checkpoint46e_post, checkpoint44b_pre, checkpoint44, checkpoint45, checkpoint46h_post, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post, release1-branch_branchpoint
Branch point for: release1_final, release1-branch
Changes since 1.3: +5 -15 lines 
Modification to the GMREDI package :
 change units of tensor-K arrays, scale now like diffusivity
 initialise all common block arrays in S/R gmredi_init
 add option to use different isopycnal(Redi) & GM diffusivity
 add option to use the advective GM form or the skew-flux form (=default)
 bug in non_unity_diagonal part fixed.
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/gmredi/gmredi_calc_diff.F,v 1.3 2001/02/04 14:38:49 cnh Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 SUBROUTINE GMREDI_CALC_DIFF(
7 I bi,bj,iMin,iMax,jMin,jMax,K,
8 I maskUp,
9 U KappaRT,KappaRS,
10 I myThid)
11
12 C /==========================================================\
13 C | SUBROUTINE GMREDI_CALC_DIFF |
14 C | o Add contrubution to net diffusivity from GM/Redi |
15 C \==========================================================/
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
24 C == Routine arguments ==
25 C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
26 C maskUp - land/water mask for Wvel points (above tracer level)
27 C myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
28 C
29 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
30 _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31 _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
32 _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
33 INTEGER myThid
34
35 #ifdef ALLOW_GMREDI
36 #include "GMREDI.h"
37
38 C == Local variables ==
39 C I, J - Loop counters
40 INTEGER i,j
41
42 DO j=jMin,jMax
43 DO i=iMin,iMax
44 KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj)
45 ENDDO
46 ENDDO
47
48 DO j=jMin,jMax
49 DO i=iMin,iMax
50 KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj)
51 ENDDO
52 ENDDO
53
54 #endif /* ALLOW_GMREDI */
55
56 RETURN
57 END
  ViewVC Help
Powered by ViewVC 1.1.22