/[MITgcm]/MITgcm/pkg/longstep/longstep_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/longstep/longstep_diagnostics_init.F

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

Revision 1.3 - (hide annotations) (download)
Tue Jan 12 23:55:48 2010 UTC (14 years, 5 months ago) by jahn
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, HEAD
Changes since 1.2: +8 -1 lines 
longstep: fix tracer conservation with rstar/nonlinear free surface
1 jahn 1.3 C $Header: /u/gcmpack/MITgcm/pkg/longstep/longstep_diagnostics_init.F,v 1.2 2010/01/12 21:34:09 jmc Exp $
2 jmc 1.2 C $Name: $
3 jahn 1.1
4     #include "LONGSTEP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: LONGSTEP_DIAGNOSTICS_INIT
8     C !INTERFACE:
9     SUBROUTINE LONGSTEP_DIAGNOSTICS_INIT( myThid )
10     C !DESCRIPTION:
11     C Routine to initialize longstep diagnostics
12    
13     C !USES:
14     IMPLICIT NONE
15     C === Global variables ===
16     #include "SIZE.h"
17     #include "EEPARAMS.h"
18     #include "PARAMS.h"
19     #include "LONGSTEP_PARAMS.h"
20    
21     C !INPUT/OUTPUT PARAMETERS:
22     C === Routine arguments ===
23     C myThid - Number of this instance of LONGSTEP_DIAGNOSTICS_INIT
24     INTEGER myThid
25     CEOP
26    
27     #ifdef ALLOW_DIAGNOSTICS
28     C !LOCAL VARIABLES:
29     C === Local variables ===
30 jmc 1.2 C msgBuf - Informational/error message buffer
31 jahn 1.1 c CHARACTER*(MAX_LEN_MBUF) msgBuf
32    
33     INTEGER diagNum
34     INTEGER diagMate
35     CHARACTER*8 diagName
36     CHARACTER*16 diagCode
37     CHARACTER*16 diagUnits
38     CHARACTER*(80) diagTitle
39    
40     C-- Add diagnostics to the (long) list
41    
42     diagName = 'LSwVel '
43     diagTitle = 'Vertical Component of Velocity (m/s)'
44     diagUnits = 'm/s '
45     diagCode = 'WM LR '
46     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
47     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
48    
49     diagName = 'LSuVel '
50     diagTitle = 'Zonal Component of Velocity (m/s)'
51     diagUnits = 'm/s '
52 jmc 1.2 diagCode = 'UUR MR '
53 jahn 1.1 diagMate = diagNum + 2
54     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
55     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
56    
57     diagName = 'LSvVel '
58     diagTitle = 'Meridional Component of Velocity (m/s)'
59     diagUnits = 'm/s '
60 jmc 1.2 diagCode = 'VVR MR '
61 jahn 1.1 diagMate = diagNum
62     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
63     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
64    
65     diagName = 'LStheta '
66     diagTitle = 'Potential Temperature'
67     diagUnits = 'degC '
68 jmc 1.2 diagCode = 'SMR MR '
69 jahn 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
70     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
71    
72     diagName = 'LSsalt '
73     diagTitle = 'Salinity'
74     diagUnits = 'psu '
75 jmc 1.2 diagCode = 'SMR MR '
76 jahn 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
77     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
78    
79     IF ( ivdc_kappa .NE. 0. _d 0 ) THEN
80     diagName = 'LScnvAdj'
81     diagTitle = 'Convective Adjustment Index [0-1] '
82     diagUnits = 'fraction '
83     diagCode = 'SM LR '
84     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
85     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
86     ENDIF
87    
88     #ifdef SHORTWAVE_HEATING
89     diagName = 'LSqsw '
90     diagTitle = 'net Short-Wave radiation (+=up)'
91     diagUnits = 'W/m^2 '
92     diagCode = 'SM U1 '
93     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
94     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
95     #endif
96    
97 jahn 1.3 diagName = 'LSfwFlux'
98     diagTitle = 'net surface Fresh-Water flux into the ocean'
99     diagUnits = 'kg/m^2/s '
100     diagCode = 'SM U1 '
101     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
102     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
103    
104 jahn 1.1 #ifdef ALLOW_GMREDI
105     IF ( useGMRedi ) THEN
106     diagName = 'LSkwx '
107     diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'
108     diagUnits = 'm^2/s '
109 jmc 1.2 diagCode = 'UM LR '
110     diagMate = diagNum + 2
111     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
112     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
113 jahn 1.1
114     diagName = 'LSkwy '
115     diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'
116     diagUnits = 'm^2/s '
117 jmc 1.2 diagCode = 'VM LR '
118     diagMate = diagNum
119     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
120     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
121 jahn 1.1
122     diagName = 'LSkwz '
123     diagTitle = 'K_33 element (W.point, Z.dir) of GM-Redi tensor'
124     diagUnits = 'm^2/s '
125     diagCode = 'WM P LR '
126 jmc 1.2 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
127     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
128 jahn 1.1 ENDIF
129     #endif
130    
131     #ifdef ALLOW_KPP
132     IF ( useKPP ) THEN
133     diagName = 'LSKPPdfS'
134     diagTitle = 'Vertical diffusion coefficient for salt & tracers'
135     diagUnits = 'm^2/s '
136 jmc 1.2 diagCode = 'SM P LR '
137     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
138     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
139 jahn 1.1
140     diagName = 'LSKPPght'
141     diagTitle = 'Nonlocal transport coefficient'
142     diagUnits = 's/m^2 '
143     diagCode = 'SM P LR '
144 jmc 1.2 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
145     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
146 jahn 1.1 ENDIF
147     #endif /* ALLOW_KPP */
148    
149     #endif /* ALLOW_DIAGNOSTICS */
150    
151     RETURN
152     END
153    
  ViewVC Help
Powered by ViewVC 1.1.22