/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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

revision 1.35 by jmc, Sat Nov 15 19:52:32 2008 UTC revision 1.43 by jmc, Wed Nov 30 20:58:35 2011 UTC
# Line 31  C         parse(2)  &  U  &  C-grid U-Po Line 31  C         parse(2)  &  U  &  C-grid U-Po
31  C                   &  V  &  C-grid V-Point  \\  C                   &  V  &  C-grid V-Point  \\
32  C                   &  M  &  C-grid Mass Point  \\  C                   &  M  &  C-grid Mass Point  \\
33  C                   &  Z  &  C-grid Corner Point  \\  C                   &  Z  &  C-grid Corner Point  \\
34  C         parse(3)  &  --  &  --Unused--  \\  C         parse(3)  &     &  Used for Level Integrated output: cumulate levels \\
35    C                   &  r  &  same but cumulate product by model level thickness \\
36    C                   &  R  &  same but cumulate product by hFac & level thickness \\
37  C         parse(4)  &  P  &  positive definite  \\  C         parse(4)  &  P  &  positive definite  \\
38  C         parse(5 ) &  C  &  with counter array  \\  C         parse(5 ) &  C  &  with counter array  \\
39    C                   &  P  &  post-processed (not filled up) from other diags  \\
40  C                   &  D  &  disable an array for output  \\  C                   &  D  &  disable an array for output  \\
41  C         parse(6--8) & '123'  &  3-digit mate number \\  C         parse(6--8) & '123'  &  retired, formerly: 3-digit mate number \\
42  C         parse(9)  &  U  &  model-level plus 1/2  \\  C         parse(9)  &  U  &  model-level plus 1/2  \\
43  C                   &  M  &  model-level middle  \\  C                   &  M  &  model-level middle  \\
44  C                   &  L  &  model-level minus 1/2  \\  C                   &  L  &  model-level minus 1/2  \\
# Line 46  C                   &  L  &  levels = MA Line 49  C                   &  L  &  levels = MA
49  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
50  C                   &  G  &  levels = Ground_level Number \\  C                   &  G  &  levels = Ground_level Number \\
51  C                   &  I  &  levels = sea-Ice_level Number \\  C                   &  I  &  levels = sea-Ice_level Number \\
52    C                   &  X  &  free levels option (need to be set explicitly) \\
53  C       \end{tabular}  C       \end{tabular}
54  C     \end{center}  C     \end{center}
55    
# Line 152  c    &diagTitle = 'Sea Surface Elevation Line 156  c    &diagTitle = 'Sea Surface Elevation
156        diagName  = 'THETA   '        diagName  = 'THETA   '
157        diagTitle = 'Potential Temperature'        diagTitle = 'Potential Temperature'
158        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
159        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
160        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
162    
# Line 166  c    I          diagName, diagCode, diag Line 170  c    I          diagName, diagCode, diag
170        diagName  = 'SALT    '        diagName  = 'SALT    '
171        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
172        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
173        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
174        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
175       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
176    
177        diagName  = 'RELHUM  '        diagName  = 'RELHUM  '
178        diagTitle = 'Relative Humidity'        diagTitle = 'Relative Humidity'
179        diagUnits = 'percent         '        diagUnits = 'percent         '
180        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
181        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
182       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
183    
# Line 188  c    I          diagName, diagCode, diag Line 192  c    I          diagName, diagCode, diag
192        diagName  = 'SALTanom'        diagName  = 'SALTanom'
193        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
194        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
195        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
196        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
197       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
198        ENDIF        ENDIF
# Line 196  c    I          diagName, diagCode, diag Line 200  c    I          diagName, diagCode, diag
200        diagName  = 'UVEL    '        diagName  = 'UVEL    '
201        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
202        diagUnits = 'm/s             '        diagUnits = 'm/s             '
203        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
204        diagMate  = diagNum + 2        diagMate  = diagNum + 2
205        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
206       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 204  c    I          diagName, diagCode, diag Line 208  c    I          diagName, diagCode, diag
208        diagName  = 'VVEL    '        diagName  = 'VVEL    '
209        diagTitle = 'Meridional Component of Velocity (m/s)'        diagTitle = 'Meridional Component of Velocity (m/s)'
210        diagUnits = 'm/s             '        diagUnits = 'm/s             '
211        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
212        diagMate  = diagNum        diagMate  = diagNum
213        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
214       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 219  c    I          diagName, diagCode, diag Line 223  c    I          diagName, diagCode, diag
223        diagName  = 'THETASQ '        diagName  = 'THETASQ '
224        diagTitle = 'Square of Potential Temperature'        diagTitle = 'Square of Potential Temperature'
225        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
226        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
227        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
229    
230        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
231        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
232        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
233        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
234        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
235       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
236    
# Line 234  c    I          diagName, diagCode, diag Line 238  c    I          diagName, diagCode, diag
238        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
239        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
240        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
241        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
242        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
243       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
244        ENDIF        ENDIF
# Line 242  c    I          diagName, diagCode, diag Line 246  c    I          diagName, diagCode, diag
246        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
247        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
248        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
249        diagCode  = 'UU P    MR      '        diagCode  = 'UURP    MR      '
250        diagMate  = diagNum + 2        diagMate  = diagNum + 2
251        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
252       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 250  c    I          diagName, diagCode, diag Line 254  c    I          diagName, diagCode, diag
254        diagName  = 'VVELSQ  '        diagName  = 'VVELSQ  '
255        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
256        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
257        diagCode  = 'VV P    MR      '        diagCode  = 'VVRP    MR      '
258        diagMate  = diagNum        diagMate  = diagNum
259        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
260       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 265  c    I          diagName, diagCode, diag Line 269  c    I          diagName, diagCode, diag
269        diagName  = 'UE_VEL_C'        diagName  = 'UE_VEL_C'
270        diagTitle = 'Eastward Velocity (m/s) (cell center)'        diagTitle = 'Eastward Velocity (m/s) (cell center)'
271        diagUnits = 'm/s             '        diagUnits = 'm/s             '
272        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
273        diagMate  = diagNum + 2        diagMate  = diagNum + 2
274        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
275       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 273  c    I          diagName, diagCode, diag Line 277  c    I          diagName, diagCode, diag
277        diagName  = 'VN_VEL_C'        diagName  = 'VN_VEL_C'
278        diagTitle = 'Northward Velocity (m/s) (cell center)'        diagTitle = 'Northward Velocity (m/s) (cell center)'
279        diagUnits = 'm/s             '        diagUnits = 'm/s             '
280        diagCode  = 'VM      MR      '        diagCode  = 'VMR     MR      '
281        diagMate  = diagNum        diagMate  = diagNum
282        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
283       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 281  c    I          diagName, diagCode, diag Line 285  c    I          diagName, diagCode, diag
285        diagName  = 'UV_VEL_C'        diagName  = 'UV_VEL_C'
286        diagTitle ='Product of horizontal Comp of velocity (cell center)'        diagTitle ='Product of horizontal Comp of velocity (cell center)'
287        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
288        diagCode  = 'UM      MR      '        diagCode  = 'UMR     MR      '
289        diagMate  = diagNum + 1        diagMate  = diagNum + 1
290        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
291       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 289  c    I          diagName, diagCode, diag Line 293  c    I          diagName, diagCode, diag
293        diagName  = 'UV_VEL_Z'        diagName  = 'UV_VEL_Z'
294        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'        diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
295        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
296        diagCode  = 'UZ      MR      '        diagCode  = 'UZR     MR      '
297        diagMate  = diagNum + 1        diagMate  = diagNum + 1
298        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
299       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 311  c    I          diagName, diagCode, diag Line 315  c    I          diagName, diagCode, diag
315        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
316        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
317        diagUnits = 'm/s             '        diagUnits = 'm/s             '
318        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
319        diagMate  = diagNum + 2        diagMate  = diagNum + 2
320        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
321       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 319  c    I          diagName, diagCode, diag Line 323  c    I          diagName, diagCode, diag
323        diagName  = 'VVELMASS'        diagName  = 'VVELMASS'
324        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
325        diagUnits = 'm/s             '        diagUnits = 'm/s             '
326        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
327        diagMate  = diagNum        diagMate  = diagNum
328        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
329       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 331  c    I          diagName, diagCode, diag Line 335  c    I          diagName, diagCode, diag
335        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
336       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
337    
338          diagName  = 'PhiVEL  '
339          diagTitle = 'Horizontal Velocity Potential (m^2/s)'
340          diagUnits = 'm^2/s           '
341          diagCode  = 'SMR P   MR      '
342    C-    use 'UVELMASS' as mate.
343          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
344         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
345    
346          diagName  = 'PsiVEL  '
347          diagTitle = 'Horizontal Velocity Stream-Function'
348          diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
349          diagCode  = 'SZ  P   MR      '
350    C-    use 'PhiVEL' as mate.
351          diagMate  = diagNum
352          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
353         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
354    
355        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
356        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
357        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
358        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
359        diagMate  = diagNum + 2        diagMate  = diagNum + 2
360        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
361       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 342  c    I          diagName, diagCode, diag Line 363  c    I          diagName, diagCode, diag
363        diagName  = 'VTHMASS '        diagName  = 'VTHMASS '
364        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'        diagTitle = 'Meridional Mass-Weight Transp of Pot Temp'
365        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
366        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
367        diagMate  = diagNum        diagMate  = diagNum
368        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
369       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 358  c    I          diagName, diagCode, diag Line 379  c    I          diagName, diagCode, diag
379        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
380       I                           //sTitle, myThid )       I                           //sTitle, myThid )
381        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
382        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
383        diagMate  = diagNum + 2        diagMate  = diagNum + 2
384        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
385       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 367  c    I          diagName, diagCode, diag Line 388  c    I          diagName, diagCode, diag
388        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
389       I                           //sTitle, myThid )       I                           //sTitle, myThid )
390        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
391        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
392        diagMate  = diagNum        diagMate  = diagNum
393        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
394       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 383  c    I          diagName, diagCode, diag Line 404  c    I          diagName, diagCode, diag
404        diagName  = 'UVELTH  '        diagName  = 'UVELTH  '
405        diagTitle = 'Zonal Transport of Pot Temp'        diagTitle = 'Zonal Transport of Pot Temp'
406        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
407        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
408        diagMate  = diagNum + 2        diagMate  = diagNum + 2
409        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
410       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 391  c    I          diagName, diagCode, diag Line 412  c    I          diagName, diagCode, diag
412        diagName  = 'VVELTH  '        diagName  = 'VVELTH  '
413        diagTitle = 'Meridional Transport of Pot Temp'        diagTitle = 'Meridional Transport of Pot Temp'
414        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
415        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
416        diagMate  = diagNum        diagMate  = diagNum
417        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
418       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 407  c    I          diagName, diagCode, diag Line 428  c    I          diagName, diagCode, diag
428        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
429       I                          //sTitle, myThid )       I                          //sTitle, myThid )
430        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
431        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
432        diagMate  = diagNum + 2        diagMate  = diagNum + 2
433        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
434       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 416  c    I          diagName, diagCode, diag Line 437  c    I          diagName, diagCode, diag
437        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
438       I                          //sTitle, myThid )       I                          //sTitle, myThid )
439        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
440        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
441        diagMate  = diagNum        diagMate  = diagNum
442        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
443       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 433  c    I          diagName, diagCode, diag Line 454  c    I          diagName, diagCode, diag
454        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
455       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
456        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
457        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
458        diagMate  = diagNum + 2        diagMate  = diagNum + 2
459        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
460       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 442  c    I          diagName, diagCode, diag Line 463  c    I          diagName, diagCode, diag
463        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
464       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
465        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
466        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
467        diagMate  = diagNum        diagMate  = diagNum
468        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
469       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 452  C---+----1----+----2----+----3----+----4 Line 473  C---+----1----+----2----+----3----+----4
473        diagName  = 'RHOAnoma'        diagName  = 'RHOAnoma'
474        diagTitle = 'Density Anomaly (=Rho-rhoConst)'        diagTitle = 'Density Anomaly (=Rho-rhoConst)'
475        diagUnits = 'kg/m^3          '        diagUnits = 'kg/m^3          '
476        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
477        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
478       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
479    
480        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
481        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
482        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6          '
483        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
484        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
485       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
486    
487        diagName  = 'URHOMASS'        diagName  = 'URHOMASS'
488        diagTitle = 'Zonal Transport of Density'        diagTitle = 'Zonal Transport of Density'
489        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
490        diagCode  = 'UU      MR      '        diagCode  = 'UUr     MR      '
491        diagMate  = diagNum + 2        diagMate  = diagNum + 2
492        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
493       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 474  C---+----1----+----2----+----3----+----4 Line 495  C---+----1----+----2----+----3----+----4
495        diagName  = 'VRHOMASS'        diagName  = 'VRHOMASS'
496        diagTitle = 'Meridional Transport of Density'        diagTitle = 'Meridional Transport of Density'
497        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
498        diagCode  = 'VV      MR      '        diagCode  = 'VVr     MR      '
499        diagMate  = diagNum        diagMate  = diagNum
500        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
501       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
502    
503        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
504        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Density'
505          diagUnits = 'kg/m^2/s        '
506          diagCode  = 'WM      LR      '
507          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
508         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
509    
510          diagName  = 'WdRHO_P '
511          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
512          diagUnits = 'kg/m^2/s        '
513          diagCode  = 'WM      LR      '
514          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
515         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
516    
517          diagName  = 'WdRHOdP '
518          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
519        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
520        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
521        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 490  C---+----1----+----2----+----3----+----4 Line 525  C---+----1----+----2----+----3----+----4
525        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '        diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
526       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
527        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
528        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
529        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
530       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
531    
# Line 498  C---+----1----+----2----+----3----+----4 Line 533  C---+----1----+----2----+----3----+----4
533        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '        diagTitle = DIAGS_MK_TITLE( 'Square of Hyd. '
534       I                           //pTitle//' Anomaly', myThid )       I                           //pTitle//' Anomaly', myThid )
535        diagUnits = 'm^4/s^4         '        diagUnits = 'm^4/s^4         '
536        diagCode  = 'SM P    MR      '        diagCode  = 'SMRP    MR      '
537        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
538       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
539    
# Line 520  c     diagTitle = 'Square of ocean botto Line 555  c     diagTitle = 'Square of ocean botto
555        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
556       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
557    
558    #ifdef NONLIN_FRSURF
559          diagName  = 'PHIHYDcR'
560          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
561         I                       //pTitle//' Anomaly @ const r', myThid )
562          diagUnits = 'm^2/s^2         '
563          diagCode  = 'SMR     MR      '
564          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
565         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
566    #endif
567    
568  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
569        diagName  = 'PHI_NH  '        diagName  = 'PHI_NH  '
570        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
571        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
572        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
573        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
574       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
575  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
# Line 547  c     diagTitle = 'Square of ocean botto Line 592  c     diagTitle = 'Square of ocean botto
592        diagName  = 'CONVADJ '        diagName  = 'CONVADJ '
593        diagTitle = 'Convective Adjustment Index [0-1] '        diagTitle = 'Convective Adjustment Index [0-1] '
594        diagUnits = 'fraction        '        diagUnits = 'fraction        '
595        diagCode  = 'SM      LR      '        diagCode  = 'SMR     LR      '
596        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
597       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
598    
# Line 688  c     diagTitle = 'Free-Surface r-Positi Line 733  c     diagTitle = 'Free-Surface r-Positi
733        diagName  = 'TOTUTEND'        diagName  = 'TOTUTEND'
734        diagTitle = 'Tendency of Zonal Component of Velocity'        diagTitle = 'Tendency of Zonal Component of Velocity'
735        diagUnits = 'm/s/day         '        diagUnits = 'm/s/day         '
736        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
737        diagMate  = diagNum + 2        diagMate  = diagNum + 2
738        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
739       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 696  c     diagTitle = 'Free-Surface r-Positi Line 741  c     diagTitle = 'Free-Surface r-Positi
741        diagName  = 'TOTVTEND'        diagName  = 'TOTVTEND'
742        diagTitle = 'Tendency of Meridional Component of Velocity'        diagTitle = 'Tendency of Meridional Component of Velocity'
743        diagUnits = 'm/s/day         '        diagUnits = 'm/s/day         '
744        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
745        diagMate  = diagNum        diagMate  = diagNum
746        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
747       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 704  c     diagTitle = 'Free-Surface r-Positi Line 749  c     diagTitle = 'Free-Surface r-Positi
749        diagName  = 'TOTTTEND'        diagName  = 'TOTTTEND'
750        diagTitle = 'Tendency of Potential Temperature'        diagTitle = 'Tendency of Potential Temperature'
751        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/day', myThid )
752        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
753        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
754       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
755    
756        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
757        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
758        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )
759        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
760          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
761         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
762    
763          diagName  = 'AB_gU   '
764          diagTitle = 'U momentum tendency from Adams-Bashforth'
765          diagUnits = 'm/s^2           '
766          diagCode  = 'UUR     MR      '
767          diagMate  = diagNum + 2
768          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
769         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
770    
771          diagName  = 'AB_gV   '
772          diagTitle = 'V momentum tendency from Adams-Bashforth'
773          diagUnits = 'm/s^2           '
774          diagCode  = 'VVR     MR      '
775          diagMate  = diagNum
776          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
777         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
778    
779    #ifdef ALLOW_NONHYDROSTATIC
780          diagName  = 'AB_gW   '
781          diagTitle = 'W momentum tendency from Adams-Bashforth'
782          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
783          diagCode  = 'WM      LR      '
784          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
785         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
786    #endif /* ALLOW_NONHYDROSTATIC */
787    
788    #ifdef ALLOW_GENERIC_ADVDIFF
789          diagName  = 'AB_gT   '
790          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
791          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
792          diagCode  = 'SMR     MR      '
793          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
794         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
795    
796          diagName  = 'AB_gS   '
797          diagTitle = DIAGS_MK_TITLE(
798         &            sTitle//'tendency from Adams-Bashforth', myThid )
799          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
800          diagCode  = 'SMR     MR      '
801        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
802       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
803    
804          diagName  = 'gTinAB  '
805          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
806          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
807          diagCode  = 'SMR     MR      '
808          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
809         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
810    
811          diagName  = 'gSinAB  '
812          diagTitle = DIAGS_MK_TITLE(
813         &            sTitle//'tendency going in Adams-Bashforth', myThid )
814          diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )
815          diagCode  = 'SMR     MR      '
816          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
817         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
818    #endif /* ALLOW_GENERIC_ADVDIFF */
819    
820        RETURN        RETURN
821        END        END
Legend:
Removed from v.1.35  
changed lines
  Added in v.1.43
  ViewVC Help
Powered by ViewVC 1.1.22