/[MITgcm]/MITgcm/pkg/aim_v23/aim_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/aim_diagnostics_init.F

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

revision 1.7 by dfer, Fri Aug 10 03:09:07 2007 UTC revision 1.8 by jmc, Tue Feb 5 15:37:58 2008 UTC
# Line 31  CEOP Line 31  CEOP
31  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
32  C     === Local variables ===  C     === Local variables ===
33  C     diagNum   :: diagnostics number in the (long) list of available diag.  C     diagNum   :: diagnostics number in the (long) list of available diag.
34    C     diagMate  :: diag. mate number in the (long) list of available diag.
35  C     diagName  :: local short name (8c) of a diagnostics  C     diagName  :: local short name (8c) of a diagnostics
36  C     diagCode  :: local parser field with characteristics of the diagnostics  C     diagCode  :: local parser field with characteristics of the diagnostics
37  C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT  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  C     diagUnits :: local string (16c): physical units of a diagnostic field
39  C     diagTitle :: local string (80c): description of field in diagnostic  C     diagTitle :: local string (80c): description of field in diagnostic
40        INTEGER       diagNum        INTEGER       diagNum
41          INTEGER       diagMate
42        CHARACTER*8   diagName        CHARACTER*8   diagName
43        CHARACTER*16  diagCode        CHARACTER*16  diagCode
44        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
# Line 49  C---+----1----+----2----+----3----+----4 Line 51  C---+----1----+----2----+----3----+----4
51       &  'Pot. Temp.  Tendency (Mass-Weighted) from Diabatic Processes'       &  'Pot. Temp.  Tendency (Mass-Weighted) from Diabatic Processes'
52        diagUnits = 'K/s             '        diagUnits = 'K/s             '
53        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
54        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
55       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
56    
57        diagName  = 'DIABQ   '        diagName  = 'DIABQ   '
58        diagTitle =        diagTitle =
59       &  'Spec.Humid. Tendency (Mass-Weighted) from Diabatic Processes'       &  'Spec.Humid. Tendency (Mass-Weighted) from Diabatic Processes'
60        diagUnits = 'g/kg/s          '        diagUnits = 'g/kg/s          '
61        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
62        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
63       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
64    
65        diagName  = 'RADSW   '        diagName  = 'RADSW   '
66        diagTitle =        diagTitle =
67       &  'Temperature Tendency due to Shortwave Radiation (TT_RSW)'       &  'Temperature Tendency due to Shortwave Radiation (TT_RSW)'
68        diagUnits = 'K/s             '        diagUnits = 'K/s             '
69        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
70        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
71       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
72    
73        diagName  = 'RADLW   '        diagName  = 'RADLW   '
74        diagTitle =        diagTitle =
75       &  'Temperature Tendency due to Longwave  Radiation (TT_RLW)'       &  'Temperature Tendency due to Longwave  Radiation (TT_RLW)'
76        diagUnits = 'K/s             '        diagUnits = 'K/s             '
77        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
78        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80    
81        diagName  = 'DTCONV  '        diagName  = 'DTCONV  '
82        diagTitle = 'Temperature Tendency due to Convection (TT_CNV)'        diagTitle = 'Temperature Tendency due to Convection (TT_CNV)'
83        diagUnits = 'K/s             '        diagUnits = 'K/s             '
84        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
85        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
86       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
87    
88        diagName  = 'TURBT   '        diagName  = 'TURBT   '
89        diagTitle =        diagTitle =
90       &  'Temperature Tendency due to Turbulence in PBL (TT_PBL)'       &  'Temperature Tendency due to Turbulence in PBL (TT_PBL)'
91        diagUnits = 'K/s             '        diagUnits = 'K/s             '
92        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
93        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
94       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
95    
96        diagName  = 'DTLS    '        diagName  = 'DTLS    '
97        diagTitle =        diagTitle =
98       &  'Temperature Tendency due to Large-scale condens. (TT_LSC)'       &  'Temperature Tendency due to Large-scale condens. (TT_LSC)'
99        diagUnits = 'K/s             '        diagUnits = 'K/s             '
100        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
101        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
102       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
103    
104        diagName  = 'DQCONV  '        diagName  = 'DQCONV  '
105        diagTitle = 'Spec. Humidity Tendency due to Convection (QT_CNV)'        diagTitle = 'Spec. Humidity Tendency due to Convection (QT_CNV)'
106        diagUnits = 'g/kg/s          '        diagUnits = 'g/kg/s          '
107        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
108        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
109       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
110    
111        diagName  = 'TURBQ   '        diagName  = 'TURBQ   '
112        diagTitle =        diagTitle =
113       &  'Spec. Humidity Tendency due to Turbulence in PBL (QT_PBL)'       &  'Spec. Humidity Tendency due to Turbulence in PBL (QT_PBL)'
114        diagUnits = 'g/kg/s          '        diagUnits = 'g/kg/s          '
115        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
116        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
117       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
118    
119        diagName  = 'DQLS    '        diagName  = 'DQLS    '
120        diagTitle =        diagTitle =
121       &  'Spec. Humidity Tendency due to Large-Scale Condens. (QT_LSC)'       &  'Spec. Humidity Tendency due to Large-Scale Condens. (QT_LSC)'
122        diagUnits = 'g/kg/s          '        diagUnits = 'g/kg/s          '
123        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
124        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
125       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
126    
127        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
128        diagTitle = 'Relative Humidity (fraction)'        diagTitle = 'Relative Humidity (fraction)'
129        diagUnits = '0-1             '        diagUnits = '0-1             '
130        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
131        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
132       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
133    
134        diagName  = 'TSR     '        diagName  = 'TSR     '
135        diagTitle = 'Top-of-atm. net Shortwave Radiation (+=dw)'        diagTitle = 'Top-of-atm. net Shortwave Radiation (+=dw)'
136        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
137        diagCode  = 'SM P    U1      '        diagCode  = 'SM P    U1      '
138        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
139       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
140    
141        diagName  = 'OLR     '        diagName  = 'OLR     '
142        diagTitle = 'Outgoing Longwave  Radiation (+=up)'        diagTitle = 'Outgoing Longwave  Radiation (+=up)'
143        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
144        diagCode  = 'SM P    U1      '        diagCode  = 'SM P    U1      '
145        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
146       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
147    
148        diagName  = 'RADSWG  '        diagName  = 'RADSWG  '
149        diagTitle = 'Net Shortwave Radiation at the Ground (+=dw)'        diagTitle = 'Net Shortwave Radiation at the Ground (+=dw)'
150        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
151        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
152        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
153       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
154    
155        diagName  = 'RADLWG  '        diagName  = 'RADLWG  '
156        diagTitle = 'Net Longwave  Radiation at the Ground (+=up)'        diagTitle = 'Net Longwave  Radiation at the Ground (+=up)'
157        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
158        diagCode  = 'SM      L1      '        diagCode  = 'SM      L1      '
159        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
160       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
161    
162        diagName  = 'HFLUX   '        diagName  = 'HFLUX   '
163        diagTitle = 'Sensible Heat Flux (+=up)'        diagTitle = 'Sensible Heat Flux (+=up)'
164        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
165        diagCode  = 'SM      L1      '        diagCode  = 'SM      L1      '
166        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
167       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
168    
169        diagName  = 'EVAP    '        diagName  = 'EVAP    '
170        diagTitle = 'Surface Evaporation (g/m2/s)'        diagTitle = 'Surface Evaporation (g/m2/s)'
171        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
172        diagCode  = 'SM      L1      '        diagCode  = 'SM      L1      '
173        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
174       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
175    
176        diagName  = 'PRECON  '        diagName  = 'PRECON  '
177        diagTitle = 'Convective  Precipitation (g/m2/s)'        diagTitle = 'Convective  Precipitation (g/m2/s)'
178        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
179        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
180        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
181       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
182    
183        diagName  = 'PRECLS  '        diagName  = 'PRECLS  '
184        diagTitle = 'Large Scale Precipitation (g/m2/s)'        diagTitle = 'Large Scale Precipitation (g/m2/s)'
185        diagUnits = 'g/m^2/s         '        diagUnits = 'g/m^2/s         '
186        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
187        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
188       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
   
       diagName  = 'CLDFRC  '  
       diagTitle = 'Total Cloud Fraction (0-1)'  
       diagUnits = '0-1             '  
       diagCode  = 'SM P    M1      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
189    
190  #ifdef ALLOW_CLOUD_3D  #ifdef ALLOW_CLOUD_3D
191        diagName  = 'CLDCLW  '        diagName  = 'CLDCLW  '
192        diagTitle = 'Local Cloud Fraction (0-1)'        diagTitle = 'Local Cloud Fraction (0-1)'
193        diagUnits = '0-1             '        diagUnits = '0-1             '
194        diagCode  = 'SM      ML      '        diagCode  = 'SM P    ML      '
195        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
196       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
197  #endif  #endif
198    
199          diagName  = 'CLDFRC  '
200          diagTitle = 'Total Cloud Fraction (0-1)'
201          diagUnits = '0-1             '
202          diagCode  = 'SM P    M1      '
203          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
204         I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
205    
206        diagName  = 'CLDPRS  '        diagName  = 'CLDPRS  '
207        diagTitle = 'Cloud Top Pressure (normalized)'        diagTitle = 'Cloud Top Pressure (normalized)'
208        diagUnits = '0-1             '        diagUnits = '0-1             '
209        WRITE(diagCode,'(A,I3.3,A)') 'SM PC',diagNum,'M1      '        diagCode  = 'SM PC   M1      '
210        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
211       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
212         I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
213    
214        diagName  = 'CLDMAS  '        diagName  = 'CLDMAS  '
215        diagTitle = 'Cloud-base Mass Flux  (kg/m^2/s)'        diagTitle = 'Cloud-base Mass Flux  (kg/m^2/s)'
216        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
217        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
218        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
219       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
220    
221        diagName  = 'DRAG    '        diagName  = 'DRAG    '
222        diagTitle = 'Surface Drag Coefficient (kg/m^2/s)'        diagTitle = 'Surface Drag Coefficient (kg/m^2/s)'
223        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
224        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
225        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
226       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
227    
228        diagName  = 'WINDS   '        diagName  = 'WINDS   '
229        diagTitle = 'Surface Wind Speed  (m/s)'        diagTitle = 'Surface Wind Speed  (m/s)'
230        diagUnits = 'm/s             '        diagUnits = 'm/s             '
231        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
232        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
233       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
234    
235        diagName  = 'TS      '        diagName  = 'TS      '
236        diagTitle = 'near Surface Air Temperature  (K)'        diagTitle = 'near Surface Air Temperature  (K)'
237        diagUnits = 'K               '        diagUnits = 'K               '
238        diagCode  = 'SM      L1      '        diagCode  = 'SM      L1      '
239        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
240       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
241    
242        diagName  = 'QS      '        diagName  = 'QS      '
243        diagTitle = 'near Surface Specific Humidity  (g/kg)'        diagTitle = 'near Surface Specific Humidity  (g/kg)'
244        diagUnits = 'g/kg            '        diagUnits = 'g/kg            '
245        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
246        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
247       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
248    
249        diagName  = 'ENPREC  '        diagName  = 'ENPREC  '
250        diagTitle ='Energy flux associated with precip. (snow, rain Temp)'        diagTitle ='Energy flux associated with precip. (snow, rain Temp)'
251        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
252        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
253        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
254       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
255    
256        diagName  = 'ALBVISDF'        diagName  = 'ALBVISDF'
257        diagTitle = 'Surface Albedo (Visible band) (0-1)'        diagTitle = 'Surface Albedo (Visible band) (0-1)'
258        diagUnits = '0-1             '        diagUnits = '0-1             '
259        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
260        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
261       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
262    
263        diagName  = 'DWNLWG  '        diagName  = 'DWNLWG  '
264        diagTitle =        diagTitle =
265       &   'Downward Component of Longwave Flux at the Ground (+=dw)'       &   'Downward Component of Longwave Flux at the Ground (+=dw)'
266        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
267        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
268        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
269       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
270    
271        diagName  = 'SWCLR   '        diagName  = 'SWCLR   '
272        diagTitle = 'Clear Sky Temp. Tendency due to Shortwave Radiation'        diagTitle = 'Clear Sky Temp. Tendency due to Shortwave Radiation'
273        diagUnits = 'K/s             '        diagUnits = 'K/s             '
274        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
275        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
276       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
277    
278        diagName  = 'LWCLR   '        diagName  = 'LWCLR   '
279        diagTitle = 'Clear Sky Temp. Tendency due to Longwave  Radiation'        diagTitle = 'Clear Sky Temp. Tendency due to Longwave  Radiation'
280        diagUnits = 'K/s             '        diagUnits = 'K/s             '
281        diagCode  = 'SM      ML      '        diagCode  = 'SM      ML      '
282        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
283       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
284    
285        diagName  = 'TSRCLR  '        diagName  = 'TSRCLR  '
286        diagTitle = 'Clear Sky Top-of-atm. net Shortwave Radiation (+=dw)'        diagTitle = 'Clear Sky Top-of-atm. net Shortwave Radiation (+=dw)'
287        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
288        diagCode  = 'SM P    U1      '        diagCode  = 'SM P    U1      '
289        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
290       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
291    
292        diagName  = 'OLRCLR  '        diagName  = 'OLRCLR  '
293        diagTitle = 'Clear Sky Outgoing Longwave  Radiation  (+=up)'        diagTitle = 'Clear Sky Outgoing Longwave  Radiation  (+=up)'
294        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
295        diagCode  = 'SM P    U1      '        diagCode  = 'SM P    U1      '
296        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
297       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
298    
299        diagName  = 'SWGCLR  '        diagName  = 'SWGCLR  '
300        diagTitle =        diagTitle =
301       &    'Clear Sky Net Shortwave Radiation at the Ground (+=dw)'       &    'Clear Sky Net Shortwave Radiation at the Ground (+=dw)'
302        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
303        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
304        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
305       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
306    
307        diagName  = 'LWGCLR  '        diagName  = 'LWGCLR  '
308        diagTitle =        diagTitle =
309       &    'Clear Sky Net Longwave  Radiation at the Ground (+=up)'       &    'Clear Sky Net Longwave  Radiation at the Ground (+=up)'
310        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
311        diagCode  = 'SM      L1      '        diagCode  = 'SM      L1      '
312        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
313       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
314    
315        diagName  = 'UFLUX   '        diagName  = 'UFLUX   '
316        diagTitle = 'Zonal Wind Surface Stress  (N/m^2)'        diagTitle = 'Zonal Wind Surface Stress  (N/m^2)'
317        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
318        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'L1      '        diagCode  = 'UU      L1      '
319        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
320       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
321         I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
322    
323        diagName  = 'VFLUX   '        diagName  = 'VFLUX   '
324        diagTitle = 'Meridional Wind Surface Stress  (N/m^2)'        diagTitle = 'Meridional Wind Surface Stress  (N/m^2)'
325        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
326        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'L1      '        diagCode  = 'VV      L1      '
327        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
328       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
329         I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
330    
331        diagName  = 'dKE_Ustr'        diagName  = 'dKE_Ustr'
332        diagTitle =        diagTitle =
333       &    'KE sink due to Zonal Surf. Stress and Strat. drag (<0)'       &    'KE sink due to Zonal Surf. Stress and Strat. drag (<0)'
334        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
335        WRITE(diagCode,'(A,I3.3,A)') 'SU   ',diagNum+2,'L1      '        diagCode  = 'SU      L1      '
336        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum + 2
337       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
338         I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
339    
340        diagName  = 'dKE_Vstr'        diagName  = 'dKE_Vstr'
341        diagTitle =        diagTitle =
342       &    'KE sink due to Meridional Surf. Stress and Strat. drag (<0)'       &    'KE sink due to Merid. Surf. Stress and Strat. drag (<0)'
343        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
344        WRITE(diagCode,'(A,I3.3,A)') 'SV   ', diagNum ,'L1      '        diagCode  = 'SV      L1      '
345        CALL DIAGNOSTICS_ADD2LIST( diagNum,        diagMate  = diagNum
346       I          diagName, diagCode, diagUnits, diagTitle, myThid )        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
347         I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
348    
349        diagName  = 'DTSIMPL '        diagName  = 'DTSIMPL '
350        diagTitle = 'Surf. Temp Change after 1 implicit time step'        diagTitle = 'Surf. Temp Change after 1 implicit time step'
351        diagUnits = 'K               '        diagUnits = 'K               '
352        diagCode  = 'SM P    L1      '        diagCode  = 'SM P    L1      '
353        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
354       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I            diagName, diagCode, diagUnits, diagTitle, 0, myThid )
355    
356  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
357    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22