/[MITgcm]/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_diagnostics_init.F

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


Revision 1.4 - (hide annotations) (download)
Fri May 2 06:10:48 2014 UTC (11 years, 3 months ago) by atn
Branch: MAIN
CVS Tags: HEAD
Changes since 1.3: +37 -22 lines
reduce some global arrays to local

1 atn 1.4 C $Header: /u/gcmpack/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_diagnostics_init.F,v 1.3 2014/04/29 06:49:40 atn Exp $
2 atn 1.1 C $Name: $
3    
4     #include "SALT_PLUME_OPTIONS.h"
5    
6     CBOP 0
7     C !ROUTINE: SALT_PLUME_DIAGNOSTICS_INIT
8     C !INTERFACE:
9     SUBROUTINE SALT_PLUME_DIAGNOSTICS_INIT( myThid )
10    
11     C !DESCRIPTION:
12     C Initialize list of all available diagnostics
13    
14     C !USES:
15     IMPLICIT NONE
16     #include "EEPARAMS.h"
17    
18     C !INPUT/OUTPUT PARAMETERS:
19     C myThid :: my Thread Id number
20     INTEGER myThid
21     CEOP
22    
23     #ifdef ALLOW_SALT_PLUME
24     #ifdef ALLOW_DIAGNOSTICS
25    
26     C !LOCAL VARIABLES:
27     C === Local variables ===
28     C diagNum :: diagnostics number in the (long) list of available diag.
29     C diagName :: local short name (8c) of a diagnostics
30     C diagCode :: local parser field with characteristics of the diagnostics
31     C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
32     C diagUnits :: local string (16c): physical units of a diagnostic field
33     C diagTitle :: local string (80c): description of field in diagnostic
34     INTEGER diagNum
35     CHARACTER*8 diagName
36     CHARACTER*16 diagCode
37     CHARACTER*16 diagUnits
38     CHARACTER*(80) diagTitle
39    
40     diagName = 'PLUMEKB '
41     diagTitle = 'fractional plume: [0-1] (unitless)'
42     diagUnits = ' '
43     diagCode = 'SM MR '
44     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
45     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
46    
47     diagName = 'oceSPtnd'
48     diagTitle = 'salt tendency due to salt plume flux'
49     & //' >0 increases salinity'
50     diagUnits = 'g/m^2/s '
51     diagCode = 'SM MR '
52     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
53     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
54    
55     #ifdef SALT_PLUME_VOLUME
56 atn 1.4 c diagName = 'oceEPtnd'
57     diagName = 'SPtd_TH '
58 atn 1.1 diagTitle = 'Energy tendency associated with salt plume flux'
59     & //' >0 increases temperature'
60     diagUnits = 'W/m^2 '
61     diagCode = 'SM MR '
62     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
63     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
64 atn 1.4 cSPtd_SLT is a duplicate of oceSPtnd
65     diagName = 'SPtd_SLT'
66     diagTitle = 'salt tendency associated with salt plume flux'
67     & //' >0 increases salinity'
68     diagUnits = 'g/m^2/s '
69     diagCode = 'SM MR '
70     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
71     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
72 atn 1.2 c
73 atn 1.4 diagName = 'oceSPvol'
74 atn 1.2 diagTitle = ''
75 atn 1.4 & //'surf brineVol associated w salt plume flux'
76     diagUnits = 'm/s '
77 atn 1.2 diagCode = 'SM U1 '
78     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80 atn 1.4 cbelow is retired
81     c diagName = 'SPkbottm'
82     c diagTitle = ''
83     c & //' range [1:Nr] '
84     c diagUnits = 'unitless '
85     c diagCode = 'SM U1 '
86     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
87     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
88     cbelow is retired
89     c diagName = 'SPbrineS'
90     c diagTitle = 'Salinity of brine, f(SPalpha) or f(SPbrineSconst)'
91     c diagUnits = 'psu '
92     c diagCode = 'SM U1 '
93     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
94     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
95 atn 1.2 c
96     diagName = 'SPVsuf2k'
97     diagTitle = 'salt plume volume flux from surface INTO lev k'
98     & //' >0 increases volume(k)'
99     diagUnits = 'm/s '
100     diagCode = 'SM MR '
101     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
102     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
103 atn 1.4 cbelow is retired
104     c diagName = 'SPVp2k '
105     c diagTitle = 'vol flux from k+1 INTO k due to SPVsuf2k'
106     c & //' >0 increases volume(k)'
107     c diagUnits = 'm/s '
108     c diagCode = 'SM MR '
109     c CALL DIAGNOSTICS_ADDTOLIST( diagNum,
110     c I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
111 atn 1.2 c
112     diagName = 'SPVk2m '
113     diagTitle = 'vol flux LEAVING k to k-1 due to SPVsuf2k&SPVp2k'
114     & //' <0 decreases volume(k)'
115     diagUnits = 'm/s '
116     diagCode = 'SM MR '
117     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
118     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
119 atn 1.1 #endif /* SALT_PLUME_VOLUME */
120    
121     diagName = 'oceSPflx'
122     diagTitle = 'net surface Salt flux rejected into the ocean during'
123     & //' freezing, (+=down),'
124     c & //' >0 increases salinity'
125     diagUnits = 'g/m^2/s '
126     diagCode = 'SM U1 '
127     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
128     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
129    
130     diagName = 'oceSPDep'
131     diagTitle = 'Salt plume depth based on density criterion (>0)'
132     diagUnits = 'm '
133     diagCode = 'SM U1 '
134     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
135     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
136    
137     #endif /* ALLOW_DIAGNOSTICS */
138     #endif /* ALLOW_SALT_PLUME */
139    
140     RETURN
141     END

  ViewVC Help
Powered by ViewVC 1.1.22