1 |
C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_diagnostics_init.F,v 1.7 2010/01/12 21:31:15 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "EXF_OPTIONS.h" |
5 |
|
6 |
CBOP |
7 |
C !ROUTINE: EXF_DIAGNOSTICS_INIT |
8 |
C !INTERFACE: |
9 |
SUBROUTINE EXF_DIAGNOSTICS_INIT( myThid ) |
10 |
|
11 |
C !DESCRIPTION: \bv |
12 |
C *==========================================================* |
13 |
C | SUBROUTINE EXF_DIAGNOSTICS_INIT |
14 |
C | o Routine to initialize list of all available diagnostics |
15 |
C | for EXF package |
16 |
C *==========================================================* |
17 |
C \ev |
18 |
C !USES: |
19 |
IMPLICIT NONE |
20 |
|
21 |
C === Global variables === |
22 |
#include "EEPARAMS.h" |
23 |
#include "EXF_PARAM.h" |
24 |
|
25 |
C !INPUT/OUTPUT PARAMETERS: |
26 |
C === Routine arguments === |
27 |
C myThid :: my Thread Id number |
28 |
INTEGER myThid |
29 |
CEOP |
30 |
|
31 |
#ifdef ALLOW_DIAGNOSTICS |
32 |
C !LOCAL VARIABLES: |
33 |
C === Local variables === |
34 |
C diagNum :: diagnostics number in the (long) list of available diag. |
35 |
C diagName :: local short name (8c) of a diagnostics |
36 |
C diagCode :: local parser field with characteristics of the diagnostics |
37 |
C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT |
38 |
C diagUnits :: local string (16c): physical units of a diagnostic field |
39 |
C diagTitle :: local string (80c): description of field in diagnostic |
40 |
INTEGER diagNum |
41 |
CHARACTER*8 diagName |
42 |
CHARACTER*16 diagCode |
43 |
CHARACTER*16 diagUnits |
44 |
CHARACTER*(80) diagTitle |
45 |
|
46 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
47 |
|
48 |
diagName = 'EXFhs ' |
49 |
diagTitle = 'Sensible heat flux into ocean, >0 increases theta' |
50 |
diagUnits = 'W/m^2 ' |
51 |
diagCode = 'SM U1 ' |
52 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
53 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
54 |
|
55 |
diagName = 'EXFhl ' |
56 |
diagTitle = 'Latent heat flux into ocean, >0 increases theta' |
57 |
diagUnits = 'W/m^2 ' |
58 |
diagCode = 'SM U1 ' |
59 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
60 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
61 |
|
62 |
diagName = 'EXFlwnet' |
63 |
diagTitle = 'Net upward longwave radiation, >0 decreases theta' |
64 |
diagUnits = 'W/m^2 ' |
65 |
diagCode = 'SM U1 ' |
66 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
67 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
68 |
|
69 |
diagName = 'EXFswnet' |
70 |
diagTitle = 'Net upward shortwave radiation, >0 decreases theta' |
71 |
diagUnits = 'W/m^2 ' |
72 |
diagCode = 'SM U1 ' |
73 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
74 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
75 |
|
76 |
diagName = 'EXFlwdn ' |
77 |
diagTitle = 'Downward longwave radiation, >0 increases theta' |
78 |
diagUnits = 'W/m^2 ' |
79 |
diagCode = 'SM U1 ' |
80 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
81 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
82 |
|
83 |
diagName = 'EXFswdn ' |
84 |
diagTitle = 'Downward shortwave radiation, >0 increases theta' |
85 |
diagUnits = 'W/m^2 ' |
86 |
diagCode = 'SM U1 ' |
87 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
88 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
89 |
|
90 |
diagName = 'EXFqnet ' |
91 |
diagTitle = 'Net upward heat flux (turb+rad), >0 decreases theta' |
92 |
diagUnits = 'W/m^2 ' |
93 |
diagCode = 'SM U1 ' |
94 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
95 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
96 |
|
97 |
C-- surface fluxes: |
98 |
diagName = 'EXFtaux ' |
99 |
diagTitle = 'zonal surface wind stress, >0 increases uVel ' |
100 |
diagUnits = 'N/m^2 ' |
101 |
IF ( stressIsOnCgrid ) THEN |
102 |
diagCode = 'UU U1 ' |
103 |
ELSE |
104 |
diagCode = 'UM U1 ' |
105 |
ENDIF |
106 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
107 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
108 |
|
109 |
diagName = 'EXFtauy ' |
110 |
diagTitle = 'meridional surface wind stress, >0 increases vVel' |
111 |
diagUnits = 'N/m^2 ' |
112 |
IF ( stressIsOnCgrid ) THEN |
113 |
diagCode = 'VV U1 ' |
114 |
ELSE |
115 |
diagCode = 'VM U1 ' |
116 |
ENDIF |
117 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
118 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
119 |
|
120 |
diagName = 'EXFuwind' |
121 |
diagTitle = 'zonal 10-m wind speed, >0 increases uVel' |
122 |
diagUnits = 'm/s ' |
123 |
diagCode = 'UM U1 ' |
124 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
125 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
126 |
|
127 |
diagName = 'EXFvwind' |
128 |
diagTitle = 'meridional 10-m wind speed, >0 increases uVel' |
129 |
diagUnits = 'm/s ' |
130 |
diagCode = 'VM U1 ' |
131 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
132 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
133 |
|
134 |
diagName = 'EXFwspee' |
135 |
diagTitle = '10-m wind speed modulus ( >= 0 )' |
136 |
diagUnits = 'm/s ' |
137 |
diagCode = 'SM U1 ' |
138 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
139 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
140 |
|
141 |
diagName = 'EXFatemp' |
142 |
diagTitle = 'surface (2-m) air temperature' |
143 |
diagUnits = 'degK ' |
144 |
diagCode = 'SM U1 ' |
145 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
146 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
147 |
|
148 |
diagName = 'EXFaqh ' |
149 |
diagTitle = 'surface (2-m) specific humidity' |
150 |
diagUnits = 'kg/kg ' |
151 |
diagCode = 'SM U1 ' |
152 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
153 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
154 |
|
155 |
diagName = 'EXFevap ' |
156 |
diagTitle = 'evaporation, > 0 increases salinity' |
157 |
diagUnits = 'm/s ' |
158 |
diagCode = 'SM U1 ' |
159 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
160 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
161 |
|
162 |
diagName = 'EXFpreci' |
163 |
diagTitle = 'precipitation, > 0 decreases salinity' |
164 |
diagUnits = 'm/s ' |
165 |
diagCode = 'SM U1 ' |
166 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
167 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
168 |
|
169 |
diagName = 'EXFsnow ' |
170 |
diagTitle = 'snow precipitation, > 0 decreases salinity' |
171 |
diagUnits = 'm/s ' |
172 |
diagCode = 'SM U1 ' |
173 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
174 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
175 |
|
176 |
diagName = 'EXFempmr' |
177 |
diagTitle = 'net upward freshwater flux, > 0 increases salinity' |
178 |
diagUnits = 'm/s ' |
179 |
diagCode = 'SM U1 ' |
180 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
181 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
182 |
|
183 |
diagName = 'EXFpress' |
184 |
diagTitle = 'atmospheric pressure field' |
185 |
diagUnits = 'N/m^2 ' |
186 |
diagCode = 'SM U1 ' |
187 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
188 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
189 |
|
190 |
diagName = 'EXFroff ' |
191 |
diagTitle = 'river runoff, > 0 decreases salinity' |
192 |
diagUnits = 'm/s ' |
193 |
diagCode = 'SM U1 ' |
194 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
195 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
196 |
|
197 |
diagName = 'EXFroft ' |
198 |
diagTitle = 'river runoff temperature' |
199 |
diagUnits = 'deg C ' |
200 |
diagCode = 'SM U1 ' |
201 |
CALL DIAGNOSTICS_ADDTOLIST( diagNum, |
202 |
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) |
203 |
|
204 |
#endif /* ALLOW_DIAGNOSTICS */ |
205 |
|
206 |
RETURN |
207 |
END |