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

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

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


Revision 1.2 - (show annotations) (download)
Wed Dec 15 00:20:37 2004 UTC (19 years, 4 months ago) by jmc
Branch: MAIN
Changes since 1.1: +22 -1 lines
add 3 diagnostics: Total PhiHyd, Bottom Pressure & Sq of Bot.Pressure

1 C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.1 2004/12/13 21:43:54 jmc Exp $
2 C $Name: $
3
4 #include "DIAG_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 0
8 C !ROUTINE: DIAGNOSTICS_MAIN_INIT
9
10 C !INTERFACE:
11 SUBROUTINE DIAGNOSTICS_MAIN_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize available diagnostics list for variables of the main code
15 C (not part of a package): set the following attributes:
16 C name (=cdiag), parsing code (=gdiag), units (=udiag), and title (=tdiag)
17 C Notes: 1) diagnostics defined here do not require any EQUIVALENCE
18 C since they get filled-in with S/R FILL_DIAGNOSTICS
19 C 2) GDIAG is defined as character*16 and can be to character*1
20 C parse(16) with the following codes currently defined:
21
22 C \begin{center}
23 C \begin{tabular}[h]{|c|c|}\hline
24 C \textbf{Positions} & \textbf{Characters}
25 C & \textbf{Meanings} \\\hline
26 C parse(1) & S & scalar \\
27 C & U & U-vector component \\
28 C & V & V-vector component \\
29 C parse(2) & U & C-grid U-Point \\
30 C & V & C-grid V-Point \\
31 C & M & C-grid Mass Point \\
32 C & Z & C-grid Corner Point \\
33 C parse(3) & -- & --Unused-- \\
34 C parse(4) & P & positive definite \\
35 C parse(5 ) & C & with counter array \\
36 C & D & disable an array for output \\
37 C parse(6--8) & '123' & 3-digit mate number \\
38 C parse(9) & U & model-level plus 1/2 \\
39 C & M & model-level middle \\
40 C & L & model-level minus 1/2 \\
41 C parse(10) & 0 & levels = 0 \\
42 C & 1 & levels = 1 \\
43 C & R & levels = Nr \\
44 C & L & levels = MAX(Nr,NrPhys) \\
45 C & M & levels = MAX(Nr,NrPhys) - 1 \\
46 C \end{tabular}
47 C \end{center}
48
49 C !USES:
50 IMPLICIT NONE
51 #include "SIZE.h"
52 #include "EEPARAMS.h"
53 c #include "EESUPPORT.h"
54 #include "PARAMS.h"
55
56 C !INPUT PARAMETERS:
57 INTEGER myThid
58 CEOP
59
60 C !LOCAL VARIABLES:
61 CHARACTER*8 diagName
62 CHARACTER*16 diagCode
63 CHARACTER*16 diagUnits
64 CHARACTER*(80) diagTitle
65
66 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
67 C For each output variable,
68 C specify Name (cdiag, 8c), Descriptions (tdiag, *c), Units (udiag, 16c)
69 C and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
70 C----------------------------------------------------------------------
71
72 diagName = 'PHIHYD '
73 diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
74 diagUnits = 'm2/s2 '
75 diagCode = 'SM MR '
76 CALL DIAGNOSTICS_ADD2LIST(
77 I diagName, diagCode, diagUnits, diagTitle, myThid )
78
79 diagName = 'PHIBOT '
80 diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
81 diagUnits = 'm2/s2 '
82 diagCode = 'SM MR '
83 CALL DIAGNOSTICS_ADD2LIST(
84 I diagName, diagCode, diagUnits, diagTitle, myThid )
85
86 diagName = 'PHIBOTSQ'
87 diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
88 diagUnits = 'm4/s4 '
89 diagCode = 'SM M1 '
90 CALL DIAGNOSTICS_ADD2LIST(
91 I diagName, diagCode, diagUnits, diagTitle, myThid )
92
93 diagName = 'DRHODR '
94 diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit) '
95 diagUnits = 'kg/m4 '
96 diagCode = 'SM LR '
97 CALL DIAGNOSTICS_ADD2LIST(
98 I diagName, diagCode, diagUnits, diagTitle, myThid )
99
100 diagName = 'DETADT2 '
101 diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
102 diagUnits = 'm2 or Pa2 /sec2 '
103 diagCode = 'SM M1 '
104 CALL DIAGNOSTICS_ADD2LIST(
105 I diagName, diagCode, diagUnits, diagTitle, myThid )
106
107 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
108
109 RETURN
110 END

  ViewVC Help
Powered by ViewVC 1.1.22