85 |
IF ( fluidIsAir ) THEN |
IF ( fluidIsAir ) THEN |
86 |
tUnit4c= 'K ' |
tUnit4c= 'K ' |
87 |
sUnit4c= 'g/kg' |
sUnit4c= 'g/kg' |
88 |
|
#ifdef ALLOW_FIZHI |
89 |
|
IF(useFIZHI)THEN |
90 |
|
sUnit4c= 'kg/kg' |
91 |
|
ENDIF |
92 |
|
#endif /* ALLOW_FIZHI */ |
93 |
ELSE |
ELSE |
94 |
tUnit4c= 'degC' |
tUnit4c= 'degC' |
95 |
sUnit4c= 'psu ' |
sUnit4c= 'psu ' |
106 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
107 |
|
|
108 |
diagName = 'ETANSQ ' |
diagName = 'ETANSQ ' |
109 |
diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)' |
diagTitle = 'Square of Perturbation of Sfc position (r_unit^2)' |
110 |
diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid ) |
diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid ) |
111 |
diagCode = 'SM M1 ' |
diagCode = 'SM M1 ' |
112 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
126 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
127 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
128 |
|
|
129 |
diagName = 'SST ' |
c diagName = 'SST ' |
130 |
diagTitle = 'Sea Surface Temperature (degC,K)' |
c diagTitle = 'Sea Surface Temperature (degC,K)' |
131 |
diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid ) |
c diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid ) |
132 |
diagCode = 'SM M1 ' |
c diagCode = 'SM M1 ' |
133 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
c CALL DIAGNOSTICS_ADD2LIST( diagNum, |
134 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
c I diagName, diagCode, diagUnits, diagTitle, myThid ) |
135 |
|
|
136 |
diagName = 'SALT ' |
diagName = 'SALT ' |
137 |
diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)' |
diagTitle = 'Salt or Water Vapor mixing ratio' |
138 |
diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
139 |
diagCode = 'SM MR ' |
diagCode = 'SM MR ' |
140 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
141 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
142 |
|
|
143 |
diagName = 'SSS ' |
diagName = 'RELHUM ' |
144 |
diagTitle = 'Sea Surface Salinity (g/kg)' |
diagTitle = 'Relative Humidity ' |
145 |
diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
diagUnits = 'percent ' |
146 |
diagCode = 'SM M1 ' |
diagCode = 'SM MR ' |
147 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
148 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
149 |
|
|
150 |
|
c diagName = 'SSS ' |
151 |
|
c diagTitle = 'Sea Surface Salinity ' |
152 |
|
c diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
153 |
|
c diagCode = 'SM M1 ' |
154 |
|
c CALL DIAGNOSTICS_ADD2LIST( diagNum, |
155 |
|
c I diagName, diagCode, diagUnits, diagTitle, myThid ) |
156 |
|
|
157 |
diagName = 'SALTanom' |
diagName = 'SALTanom' |
158 |
diagTitle = 'Salt anomaly (=SALT-35; g/kg)' |
diagTitle = 'Salt anomaly (=SALT-35; g/kg)' |
159 |
diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) |
175 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
176 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
177 |
|
|
|
diagName = 'UVEL_k2 ' |
|
|
diagTitle = 'Zonal Component of Velocity at level 2 (m/s)' |
|
|
diagUnits = 'm/s ' |
|
|
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'M1 ' |
|
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
|
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
|
|
|
|
|
diagName = 'VVEL_k2 ' |
|
|
diagTitle = 'Meridional Component of Velocity at level 2 (m/s)' |
|
|
diagUnits = 'm/s ' |
|
|
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'M1 ' |
|
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
|
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
|
|
|
|
178 |
diagName = 'WVEL ' |
diagName = 'WVEL ' |
179 |
diagTitle = 'Vertical Component of Velocity (r_units/s)' |
diagTitle = 'Vertical Component of Velocity (r_units/s)' |
180 |
diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid ) |
diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid ) |
295 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
296 |
|
|
297 |
diagName = 'USLTMASS' |
diagName = 'USLTMASS' |
298 |
diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)' |
diagTitle = 'Zonal Mass-Weight Transp of Salt ' |
299 |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) |
300 |
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
301 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
302 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
303 |
|
|
304 |
diagName = 'VSLTMASS' |
diagName = 'VSLTMASS' |
305 |
diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)' |
diagTitle = 'Meridional Mass-Weight Transp of Salt' |
306 |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) |
307 |
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
308 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
309 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
310 |
|
|
311 |
diagName = 'WSLTMASS' |
diagName = 'WSLTMASS' |
312 |
diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)' |
diagTitle = 'Vertical Mass-Weight Transp of Salt ' |
313 |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) |
314 |
diagCode = 'WM LR ' |
diagCode = 'WM LR ' |
315 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
337 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
338 |
|
|
339 |
diagName = 'UVELSLT ' |
diagName = 'UVELSLT ' |
340 |
diagTitle = 'Zonal Transp of Salt (g/kg.m/s)' |
diagTitle = 'Zonal Transp of Salt ' |
341 |
diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) |
diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) |
342 |
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
343 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
344 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
345 |
|
|
346 |
diagName = 'VVELSLT ' |
diagName = 'VVELSLT ' |
347 |
diagTitle = 'Meridional Transp of Salt (g/kg.m/s)' |
diagTitle = 'Meridional Transp of Salt ' |
348 |
diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) |
diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) |
349 |
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
350 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
351 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
352 |
|
|
353 |
diagName = 'WVELSLT ' |
diagName = 'WVELSLT ' |
354 |
diagTitle = 'Vertical Transp of Salt (g/kg.m/s)' |
diagTitle = 'Vertical Transp of Salt ' |
355 |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) |
diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) |
356 |
diagCode = 'WM LR ' |
diagCode = 'WM LR ' |
357 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
358 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
359 |
|
|
360 |
|
diagName = 'UVELPHI ' |
361 |
|
diagTitle = 'Zonal Transp of Geopotential Height (m2/s2.m/s)' |
362 |
|
diagUnits = 'm^3/s^3 ' |
363 |
|
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
364 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
365 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
366 |
|
|
367 |
|
diagName = 'VVELPHI ' |
368 |
|
diagTitle = 'Meridional Transp of Geopotential Height (m2/s2.m/s)' |
369 |
|
diagUnits = 'm^3/s^3 ' |
370 |
|
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
371 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
372 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
373 |
|
|
374 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
375 |
|
|
376 |
diagName = 'RHOAnoma' |
diagName = 'RHOAnoma' |
444 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
445 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
446 |
|
|
447 |
C-- Momentum terms: |
diagName = 'CONVADJ ' |
448 |
diagName = 'VISCA4 ' |
diagTitle = 'Convective Adjustment Index [0-1] ' |
449 |
diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s) ' |
diagUnits = 'fraction ' |
450 |
diagUnits = 'm^4/s ' |
diagCode = 'SM LR ' |
|
diagCode = 'SM MR ' |
|
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
|
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
|
|
|
|
|
diagName = 'VISCAH ' |
|
|
diagTitle = 'Harmonic Viscosity Coefficient in (m2/s) ' |
|
|
diagUnits = 'm^2/s ' |
|
|
diagCode = 'SM MR ' |
|
451 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
452 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
453 |
|
|
501 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
502 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
503 |
|
|
504 |
|
#ifdef SHORTWAVE_HEATING |
505 |
|
diagName = 'SWFLUX ' |
506 |
|
diagTitle = 'net upward SW radiation, >0 increases theta ' |
507 |
|
diagUnits = 'W/m^2 ' |
508 |
|
diagCode = 'SM U1 ' |
509 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
510 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
511 |
|
#endif |
512 |
|
|
513 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
514 |
diagName = 'PRESSURE' |
diagName = 'PRESSURE' |
515 |
diagTitle = 'Atmospheric Pressure (Pa) ' |
diagTitle = 'Atmospheric Pressure (Pa) ' |
518 |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
519 |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
520 |
|
|
521 |
|
diagName = 'RSURF ' |
522 |
|
diagTitle = 'Surface (Pressure, Height) (Pa,m)' |
523 |
|
diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid ) |
524 |
|
diagCode = 'SM M1 ' |
525 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
526 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
527 |
|
|
528 |
|
diagName = 'TOTUTEND' |
529 |
|
diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)' |
530 |
|
diagUnits = 'm/s^2 ' |
531 |
|
WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' |
532 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
533 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
534 |
|
|
535 |
|
diagName = 'TOTVTEND' |
536 |
|
diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)' |
537 |
|
diagUnits = 'm/s^2 ' |
538 |
|
WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' |
539 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
540 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
541 |
|
|
542 |
|
diagName = 'TOTTTEND' |
543 |
|
diagTitle = 'Tendency of Potential Temperature (degC,K)' |
544 |
|
diagUnits = 'deg K/sec' |
545 |
|
diagCode = 'SM MR ' |
546 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
547 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
548 |
|
|
549 |
|
diagName = 'TOTSTEND' |
550 |
|
diagTitle = 'Tendency of Salt or Water Vapor mixing ratio' |
551 |
|
diagUnits = 'kg/kg-sec' |
552 |
|
diagCode = 'SM MR ' |
553 |
|
CALL DIAGNOSTICS_ADD2LIST( diagNum, |
554 |
|
I diagName, diagCode, diagUnits, diagTitle, myThid ) |
555 |
|
|
556 |
|
|
557 |
RETURN |
RETURN |
558 |
END |
END |