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

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

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

Revision 1.4 - (hide 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 jmc 1.4 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 adcroft 1.1
4     #include "GMREDI_OPTIONS.h"
5    
6     SUBROUTINE GMREDI_CALC_DIFF(
7     I bi,bj,iMin,iMax,jMin,jMax,K,
8     I maskUp,
9 heimbach 1.2 U KappaRT,KappaRS,
10 adcroft 1.1 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 jmc 1.4 C I, J - Loop counters
40 adcroft 1.1 INTEGER i,j
41    
42     DO j=jMin,jMax
43     DO i=iMin,iMax
44 jmc 1.4 KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj)
45 adcroft 1.1 ENDDO
46     ENDDO
47    
48     DO j=jMin,jMax
49     DO i=iMin,iMax
50 jmc 1.4 KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*Kwz(i,j,k,bi,bj)
51 adcroft 1.1 ENDDO
52     ENDDO
53    
54     #endif /* ALLOW_GMREDI */
55    
56     RETURN
57     END
  ViewVC Help
Powered by ViewVC 1.1.22