/[MITgcm]/MITgcm/pkg/ggl90/GGL90.h
ViewVC logotype

Annotation of /MITgcm/pkg/ggl90/GGL90.h

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

Revision 1.5 - (hide annotations) (download)
Fri Aug 6 18:37:05 2010 UTC (13 years, 10 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint63p, checkpoint63q, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint63, checkpoint62k, checkpoint62j, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x
Changes since 1.4: +7 -7 lines
File MIME type: text/plain 
Minor changes in pkg/ggl90:
 o GGL90diffKrS was removed --> always use GGL90diffKr
 o GGL90viscAr was removed --> replaced with GGL90viscArU, GGL90viscArV
 o hack of mxlMaxFlag=2 --> ensure mixing between first and second level (commented out for now)
 o change in max/min operations to ensure that smoothing is ok
 o smoothing of GGL90viscAr was moved to ggl90_calc.F (as done for GGL90diffKr)
 o always use diffKrNrT as background profile (i.e. never use diffKr field)
1 gforget 1.5 C $Header: /u/gcmpack/MITgcm/pkg/ggl90/GGL90.h,v 1.4 2009/01/30 02:23:56 dfer Exp $
2 mlosch 1.1 C $Name: $
3    
4     #ifdef ALLOW_GGL90
5    
6     CBOP
7     C !ROUTINE: GGL90.h
8    
9     C !DESCRIPTION: \bv
10     C /==========================================================\
11     C | GGL90.h |
12     C | o Basic header for Garpar et al. (1990) |
13     C | vertical mixing parameterization. Implementation |
14     C | follows OPA implementation of Blanke+Delecuse (1993) |
15     C | Contains all GGL90 field declarations. |
16     C \==========================================================/
17    
18     C-----------------------------------------------------------------------
19     C
20     C Parameters that can be set in data.ggl90
21     C unless otherwise noted, equation numbers refer to Gaspar et al.
22     C (1990), JGR, 95(C9) pp. 16,179ff.
23     C
24     C GGL90ck - constant in viscosity coefficient (eq.10)
25     C GGL90ceps - dissipation constant according to Kolmogorov (1942)
26     C GGL90m2 - constant relating wind stress to vertical stress of TKE
27     C in K [d(TKE)/dz](0)=m2*ustar^3
28     C GGL90alpha - constant relating viscosity to GGL90 diffusivity
29     C (default=1 in Gaspar etal.)
30 mlosch 1.3 C GGL90TKEsurfmin - minimum of surface kinetic energy boundary condition
31     C (default=GGL90TKEmin)
32 mlosch 1.1 C GGL90TKEmin - minimum kinetic energy, leads to minimum mixing if TKE=0.
33     C GGL90TKEbottom - bottom bounardy condition for kinetic energy
34     C (default=GGL90TKEmin)
35     C GGL90TKEFile - File with initial field of TKE
36     C
37     C GGL90dumpFreq - analogue of dumpFreq (= default)
38     C GGL90taveFreq - analogue of taveFreq (= default)
39     C GGL90mixingMaps - output to standard out (default = .FALSE.)
40     C GGL90writeState - output to files (default = .FALSE.)
41     C
42     C Time varying parameters computed by subroutine ggl90_calc
43     C GGL90TKE - prognostic variable stepped forward in time ((m/s)^2)
44     C GGL90viscAr - Vertical eddy viscosity coefficient (m^2/s)
45     C GGL90diffKr - Vertical diffusion coefficient for heat,
46     C salt and tracers (m^2/s)
47     C
48     C-----------------------------------------------------------------------
49     C \ev
50     CEOP
51     _RL SQRTTWO
52     PARAMETER ( SQRTTWO = 1.41421356237310D0 )
53     _RL GGL90eps
54     PARAMETER ( GGL90eps = 2.23D-16 )
55    
56     CHARACTER*(MAX_LEN_FNAM) GGL90TKEFile
57     COMMON /GGL90_PARMS_C/ GGL90TKEFile
58    
59    
60     _RL GGL90ck, GGL90ceps
61     _RL GGL90alpha, GGL90m2
62 mlosch 1.2 _RL GGL90diffTKEh
63 mlosch 1.1 _RL GGL90mixingLengthMin
64 mlosch 1.3 _RL GGL90TKEmin, GGL90TKEsurfMin, GGL90TKEbottom
65 mlosch 1.1 _RL GGL90viscMax, GGL90diffMax
66     _RL GGL90dumpFreq, GGL90taveFreq
67 dfer 1.4 INTEGER mxlMaxFlag
68 mlosch 1.1 COMMON /GGL90_PARMS_R/
69     & GGL90ck, GGL90ceps,
70     & GGL90alpha, GGL90m2,
71 mlosch 1.2 & GGL90diffTKEh,
72 mlosch 1.1 & GGL90mixingLengthMin,
73 mlosch 1.3 & GGL90TKEmin, GGL90TKEsurfMin, GGL90TKEbottom,
74 mlosch 1.1 & GGL90viscMax, GGL90diffMax,
75 dfer 1.4 & GGL90dumpFreq, GGL90taveFreq, mxlMaxFlag
76 mlosch 1.1
77 gforget 1.5 _RL GGL90TKE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
78     _RL GGL90viscArU(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
79     _RL GGL90viscArV(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
80     _RL GGL90diffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
81 dfer 1.4 c _RL recip_hFacI(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
82 mlosch 1.1 COMMON /GGL90_FIELDS/ GGL90TKE,
83 gforget 1.5 & GGL90viscArU, GGL90viscArV, GGL90diffKr
84 dfer 1.4 c & ,recip_hFacI
85 mlosch 1.1
86     LOGICAL GGL90isOn, GGL90mixingMaps, GGL90writeState
87     COMMON /GGL90_PARMS_L/
88     & GGL90isOn, GGL90mixingMaps, GGL90writeState
89    
90 dfer 1.4 #ifdef ALLOW_GGL90_SMOOTH
91 gforget 1.5 COMMON /GGL90_CORNER/ mskCor
92 dfer 1.4 _RL mskCor(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93     #endif
94    
95 mlosch 1.1 #endif /* ALLOW_GGL90 */
  ViewVC Help
Powered by ViewVC 1.1.22