/[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.2 - (hide annotations) (download)
Mon Sep 27 08:02:05 2004 UTC (19 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint57m_post, checkpoint57s_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint57g_pre, checkpoint56c_post, checkpoint57y_pre, checkpoint57f_pre, checkpoint57a_post, checkpoint55g_post, checkpoint55f_post, checkpoint57r_post, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint57q_post, checkpoint57z_post, checkpoint57c_post, checkpoint55e_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.1: +3 -1 lines
File MIME type: text/plain 
  - add horizontal diffusion of TKE, requires exchanges in
    do_fields_blocking_exchanges, horizontal diffusivity is zero by
    default. In OPA there is no horizontal diffusion of TKE but the
    mixing coefficients are computed from a horizontal average of TKE of
    6 points or so. I think that diffusion has a little more physical
    justification.
  - clean up ggl90_calc in the hope of reducing memory usage (this hope
    was in vain)
  - mask tke-variable in ggl90_init
1 mlosch 1.2 C $Header: /u/gcmpack/MITgcm/pkg/ggl90/GGL90.h,v 1.1 2004/09/16 11:27:18 mlosch 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     C GGL90TKEmin - minimum kinetic energy, leads to minimum mixing if TKE=0.
31     C GGL90TKEbottom - bottom bounardy condition for kinetic energy
32     C (default=GGL90TKEmin)
33     C GGL90TKEFile - File with initial field of TKE
34     C
35     C GGL90dumpFreq - analogue of dumpFreq (= default)
36     C GGL90taveFreq - analogue of taveFreq (= default)
37     C GGL90mixingMaps - output to standard out (default = .FALSE.)
38     C GGL90writeState - output to files (default = .FALSE.)
39     C
40     C Time varying parameters computed by subroutine ggl90_calc
41     C GGL90TKE - prognostic variable stepped forward in time ((m/s)^2)
42     C GGL90viscAr - Vertical eddy viscosity coefficient (m^2/s)
43     C GGL90diffKr - Vertical diffusion coefficient for heat,
44     C salt and tracers (m^2/s)
45     C
46     C-----------------------------------------------------------------------
47     C \ev
48     CEOP
49     _RL SQRTTWO
50     PARAMETER ( SQRTTWO = 1.41421356237310D0 )
51     _RL GGL90eps
52     PARAMETER ( GGL90eps = 2.23D-16 )
53    
54     CHARACTER*(MAX_LEN_FNAM) GGL90TKEFile
55     COMMON /GGL90_PARMS_C/ GGL90TKEFile
56    
57    
58     _RL GGL90ck, GGL90ceps
59     _RL GGL90alpha, GGL90m2
60 mlosch 1.2 _RL GGL90diffTKEh
61 mlosch 1.1 _RL GGL90mixingLengthMin
62     _RL GGL90TKEmin, GGL90TKEbottom
63     _RL GGL90viscMax, GGL90diffMax
64     _RL GGL90dumpFreq, GGL90taveFreq
65     COMMON /GGL90_PARMS_R/
66     & GGL90ck, GGL90ceps,
67     & GGL90alpha, GGL90m2,
68 mlosch 1.2 & GGL90diffTKEh,
69 mlosch 1.1 & GGL90mixingLengthMin,
70     & GGL90TKEmin, GGL90TKEbottom,
71     & GGL90viscMax, GGL90diffMax,
72     & GGL90dumpFreq, GGL90taveFreq
73    
74     _RL GGL90TKE (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
75     _RL GGL90viscAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
76     _RL GGL90diffKr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
77     _RL recip_hFacI(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
78     COMMON /GGL90_FIELDS/ GGL90TKE,
79     & GGL90viscAr, GGL90diffKr,
80     & recip_hFacI
81    
82     LOGICAL GGL90isOn, GGL90mixingMaps, GGL90writeState
83     COMMON /GGL90_PARMS_L/
84     & GGL90isOn, GGL90mixingMaps, GGL90writeState
85    
86     #endif /* ALLOW_GGL90 */
  ViewVC Help
Powered by ViewVC 1.1.22