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

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

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

Revision 1.2 - (show 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 C $Header: /u/gcmpack/MITgcm/pkg/ggl90/GGL90.h,v 1.1 2004/09/16 11:27:18 mlosch Exp $
2 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 _RL GGL90diffTKEh
61 _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 & GGL90diffTKEh,
69 & 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