/[MITgcm]/MITgcm_contrib/submesoscale/code/gmredi_mnc_init.F
ViewVC logotype

Contents of /MITgcm_contrib/submesoscale/code/gmredi_mnc_init.F

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


Revision 1.1 - (show annotations) (download)
Fri May 30 21:51:24 2008 UTC (17 years, 2 months ago) by dimitri
Branch: MAIN
This is the code from pkg/gmredi (checkpoint59r) on which Baylor based
his submesoscale parameterization.

1 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_mnc_init.F,v 1.2 2007/10/25 12:08:00 mlosch Exp $
2 C $Name: $
3
4 #include "GMREDI_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 0
8 C !ROUTINE: GMREDI_MNC_INIT
9
10 C !INTERFACE:
11 SUBROUTINE GMREDI_MNC_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize MNC lookup tables for GMRedi.
15
16 C !USES:
17 IMPLICIT NONE
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GMREDI.h"
22
23 C !INPUT PARAMETERS:
24 INTEGER myThid
25 CEOP
26
27 #ifdef ALLOW_GMREDI
28 #ifdef ALLOW_MNC
29
30 IF (useMNC) THEN
31
32 CALL MNC_CW_ADD_VNAME('Kux','U_xy_Hn__C__t', 4,5,myThid)
33 CALL MNC_CW_ADD_VATTR_TEXT('Kux','units','m^2/s', myThid)
34 CALL MNC_CW_ADD_VATTR_TEXT('Kux','description',
35 & 'GM Tensor Kux is K_11 element, X direction at U point',
36 & myThid)
37 C CALL MNC_CW_ADD_VATTR_TEXT('Kux',
38 C & 'coordinates','XC YC RC iter', myThid)
39
40 CALL MNC_CW_ADD_VNAME('Kvy','V_xy_Hn__C__t', 4,5,myThid)
41 CALL MNC_CW_ADD_VATTR_TEXT('Kvy','units','m^2/s', myThid)
42 CALL MNC_CW_ADD_VATTR_TEXT('Kvy','description',
43 & 'GM Tensor Kvy is K_22 element, Y direction at V point',
44 & myThid)
45
46 CALL MNC_CW_ADD_VNAME('Kuz','U_xy_Hn__C__t', 4,5,myThid)
47 CALL MNC_CW_ADD_VATTR_TEXT('Kuz','units','m^2/s', myThid)
48 CALL MNC_CW_ADD_VATTR_TEXT('Kuz','description',
49 & 'GM Tensor Kuz is K_13 element, Z direction at U point',
50 & myThid)
51
52 CALL MNC_CW_ADD_VNAME('Kvz','V_xy_Hn__C__t', 4,5,myThid)
53 CALL MNC_CW_ADD_VATTR_TEXT('Kvz','units','m^2/s', myThid)
54 CALL MNC_CW_ADD_VATTR_TEXT('Kvz','description',
55 & 'GM Tensor Kvz is K_23 element, Z direction at V point',
56 & myThid)
57
58 CALL MNC_CW_ADD_VNAME('Kwx','Cen_xy_Hn__L__t', 4,5,myThid)
59 CALL MNC_CW_ADD_VATTR_TEXT('Kwx','units','m^2/s', myThid)
60 CALL MNC_CW_ADD_VATTR_TEXT('Kwx','description',
61 & 'GM Tensor Kwx is K_31 element, X direction at W point',
62 & myThid)
63
64 CALL MNC_CW_ADD_VNAME('Kwy','Cen_xy_Hn__L__t', 4,5,myThid)
65 CALL MNC_CW_ADD_VATTR_TEXT('Kwy','units','m^2/s', myThid)
66 CALL MNC_CW_ADD_VATTR_TEXT('Kwy','description',
67 & 'GM Tensor Kwy is K_32 element, Y direction at W point',
68 & myThid)
69
70 CALL MNC_CW_ADD_VNAME('Kwz','Cen_xy_Hn__L__t', 4,5,myThid)
71 CALL MNC_CW_ADD_VATTR_TEXT('Kwz','units','m^2/s', myThid)
72 CALL MNC_CW_ADD_VATTR_TEXT('Kwz','description',
73 & 'GM Tensor Kwz is K_33 element, Z direction at W point',
74 & myThid)
75
76 CALL MNC_CW_ADD_VNAME('VisbK','Cen_xy_Hn__-__t',3,4,myThid)
77 CALL MNC_CW_ADD_VATTR_TEXT('VisbK','units','m^2/s', myThid)
78 CALL MNC_CW_ADD_VATTR_TEXT('VisbK','description',
79 & 'GM mixing/stirring coef after Visbeck et al.',
80 & myThid)
81
82 CALL MNC_CW_ADD_VNAME('PsiX','Cen_xy_Hn__C__t',4,5,myThid)
83 CALL MNC_CW_ADD_VATTR_TEXT('PsiX','units','m^2/s', myThid)
84 CALL MNC_CW_ADD_VATTR_TEXT('PsiX','description',
85 & 'GM bolus velocity streamfunction in X',
86 & myThid)
87
88 CALL MNC_CW_ADD_VNAME('PsiY','Cen_xy_Hn__C__t',4,5,myThid)
89 CALL MNC_CW_ADD_VATTR_TEXT('PsiY','units','m^2/s', myThid)
90 CALL MNC_CW_ADD_VATTR_TEXT('PsiY','description',
91 & 'GM bolus velocity streamfunction in Y',
92 & myThid)
93
94 ENDIF
95
96 #endif /* ALLOW_MNC */
97 #endif /* ALLOW_GMREDI */
98
99 RETURN
100 END
101
102 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

  ViewVC Help
Powered by ViewVC 1.1.22