/[MITgcm]/MITgcm/pkg/gmredi/GMREDI.h
ViewVC logotype

Contents of /MITgcm/pkg/gmredi/GMREDI.h

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

Revision 1.6 - (show annotations) (download)
Sun Dec 16 18:54:49 2001 UTC (22 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint44f_post, checkpoint43a-release1mods, chkpt44d_post, checkpoint44e_pre, release1-branch_tutorials, chkpt44a_post, chkpt44c_pre, checkpoint44g_post, release1-branch-end, release1_final_v1, checkpoint44b_post, chkpt44a_pre, checkpoint44b_pre, checkpoint44, chkpt44c_post, checkpoint44f_pre, release1-branch_branchpoint
Branch point for: release1_final, release1-branch
Changes since 1.5: +35 -8 lines
File MIME type: text/plain 
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.h,v 1.5 2001/02/04 14:38:49 cnh Exp $
2 C $Name: $
3
4 #ifdef ALLOW_GMREDI
5
6 C Package flag
7 logical GMREDIisON
8 logical GM_AdvForm
9 logical GM_ExtraDiag
10 COMMON /GM_PACKAGE/ GMREDIisON, GM_AdvForm, GM_ExtraDiag
11
12 C GM/Redi parameters
13 _RL GM_isopycK
14 _RL GM_background_K
15 _RL GM_maxSlope
16 _RL GM_Kmin_horiz
17 _RL GM_Visbeck_alpha
18 _RL GM_Visbeck_length
19 _RL GM_Visbeck_depth
20 _RL GM_Visbeck_maxval_K
21 CHARACTER*(MAX_LEN_FNAM) GM_taper_scheme
22 _RL GM_Scrit
23 _RL GM_Sd
24 COMMON /GM_PARAMS/ GM_isopycK, GM_background_K,
25 & GM_maxSlope,
26 & GM_Kmin_horiz,
27 & GM_Visbeck_alpha,
28 & GM_Visbeck_length,
29 & GM_Visbeck_depth,
30 & GM_Visbeck_maxval_K,
31 & GM_taper_scheme,
32 & GM_Scrit,
33 & GM_Sd
34
35
36 C More GM/Redi parameters diagnosed from previous block
37 C (not directly user configured)
38 _RL GM_rMaxSlope
39 _RL GM_skewflx
40 _RL GM_advect
41 COMMON /GM_PARAMS2/ GM_rMaxSlope,
42 & GM_skewflx, GM_advect
43
44 C GM/Redi tensor elements
45
46 C Bottom row of tensor corresponds to W points
47 C Kwx is K_31 element, X direction at W point
48 C Kwy is K_32 element, Y direction at W point
49 C Kwz is K_33 element, Z direction at W point
50 _RL Kwx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
51 _RL Kwy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
52 _RL Kwz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
53 COMMON /GM_Wtensor/ Kwx,Kwy,Kwz
54
55 #ifdef GM_NON_UNITY_DIAGONAL
56 C Horizontal part of the tensor
57 C Kux is K_11 element, X direction at U point
58 C Kvy is K_22 element, Y direction at V point
59 _RL Kux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
60 _RL Kvy(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
61 COMMON /GM_HorTensor/ Kux,Kvy
62 #else
63 _RL Kux,Kvy
64 PARAMETER(Kux=1.,Kvy=1.)
65 #endif
66
67 #ifdef GM_EXTRA_DIAGONAL
68 C First/second rows of tensor corresponds to U/V points
69 C Kuz is K_13 element, Z direction at U point
70 C Kvz is K_23 element, Z direction at V point
71 _RL Kuz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
72 _RL Kvz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
73 COMMON /GM_UVtensor/ Kuz,Kvz
74 #else
75 _RL Kuz,Kvz
76 PARAMETER(Kuz=1.,Kvz=1.)
77 #endif
78
79 #ifdef GM_BOLUS_ADVEC
80 C GM advection formulation: bolus velocities are derived from 2
81 C streamfunctions PsiX and PsiY :
82 _RL GM_PsiX(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
83 _RL GM_PsiY(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
84 COMMON /GM_BOLUS/ GM_PsiX,GM_PsiY
85 #endif
86
87 #ifdef GM_VISBECK_VARIABLE_K
88 C GM mixing/stirring coefficient (spatially variable in horizontal
89 C for Visbeck et al. parameterization)
90 _RL VisbeckK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91 COMMON /GM_Visbeck/ VisbeckK
92 #endif
93
94 #endif /* ALLOW_GMREDI */
  ViewVC Help
Powered by ViewVC 1.1.22