v4_3deg/code/exf_diagnostics_init.F

C $Header: /u/gcmpack/MITgcm_contrib/ecco_darwin/v4_3deg/code/exf_diagnostics_init.F,v 1.1 2020/01/23 01:27:38 dimitri Exp $
C $Name:  $

#include "EXF_OPTIONS.h"

CBOP
C     !ROUTINE: EXF_DIAGNOSTICS_INIT
C     !INTERFACE:
      SUBROUTINE EXF_DIAGNOSTICS_INIT( myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE EXF_DIAGNOSTICS_INIT
C     | o Routine to initialize list of all available diagnostics
C     |   for EXF package
C     *==========================================================*
C     \ev
C     !USES:
      IMPLICIT NONE

C     === Global variables ===
#include "EEPARAMS.h"
#include "EXF_PARAM.h"

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid ::  my Thread Id number
      INTEGER myThid
CEOP

#ifdef ALLOW_DIAGNOSTICS
C     !LOCAL VARIABLES:
C     === Local variables ===
C     diagNum   :: diagnostics number in the (long) list of available diag.
C     diagName  :: local short name (8c) of a diagnostics
C     diagCode  :: local parser field with characteristics of the diagnostics
C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
C     diagUnits :: local string (16c): physical units of a diagnostic field
C     diagTitle :: local string (80c): description of field in diagnostic
      INTEGER       diagNum
      CHARACTER*8   diagName
      CHARACTER*16  diagCode
      CHARACTER*16  diagUnits
      CHARACTER*(80) diagTitle

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      diagName  = 'EXFhs   '
      diagTitle = 'Sensible heat flux into ocean, >0 increases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFhl   '
      diagTitle = 'Latent heat flux into ocean, >0 increases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFlwnet'
      diagTitle = 'Net upward longwave radiation, >0 decreases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFswnet'
      diagTitle = 'Net upward shortwave radiation, >0 decreases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFlwdn '
      diagTitle = 'Downward longwave radiation, >0 increases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFswdn '
      diagTitle = 'Downward shortwave radiation, >0 increases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFqnet '
      diagTitle = 'Net upward heat flux (turb+rad), >0 decreases theta'
      diagUnits = 'W/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

C--   surface fluxes:
      diagName  = 'EXFtaux '
      diagTitle = 'zonal surface wind stress, >0 increases uVel    '
      diagUnits = 'N/m^2           '
      IF ( stressIsOnCgrid ) THEN
       diagCode = 'UU      U1      '
      ELSE
       diagCode = 'UM      U1      '
      ENDIF
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFtauy '
      diagTitle = 'meridional surface wind stress, >0 increases vVel'
      diagUnits = 'N/m^2           '
      IF ( stressIsOnCgrid ) THEN
       diagCode = 'VV      U1      '
      ELSE
       diagCode = 'VM      U1      '
      ENDIF
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFuwind'
      diagTitle = 'zonal 10-m wind speed, >0 eastward'
      diagUnits = 'm/s             '
      diagCode  = 'UM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFvwind'
      diagTitle = 'meridional 10-m wind speed, >0 northward'
      diagUnits = 'm/s             '
      diagCode  = 'VM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFwspee'
      diagTitle = '10-m wind speed modulus ( >= 0 )'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFatemp'
      diagTitle = 'surface (2-m) air temperature'
      diagUnits = 'degK            '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFaqh  '
      diagTitle = 'surface (2-m) specific humidity'
      diagUnits = 'kg/kg           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFevap '
      diagTitle = 'evaporation, > 0 increases salinity'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFpreci'
      diagTitle = 'precipitation, > 0 decreases salinity'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFsnow '
      diagTitle = 'snow precipitation, > 0 decreases salinity'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFempmr'
      diagTitle = 'net upward freshwater flux, > 0 increases salinity'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFpress'
      diagTitle = 'atmospheric pressure field'
      diagUnits = 'N/m^2           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFroff '
      diagTitle = 'river runoff, > 0 decreases salinity'
      diagUnits = 'm/s             '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFroft '
      diagTitle = 'river runoff temperature'
      diagUnits = 'deg C           '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFsalfx'
      diagTitle = 'upward salt flux, > 0 decreases salinity'
      diagUnits = 'psu.kg/m^2/s    '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'EXFapco2 '
      diagTitle = 'Atmospheric partial pressure of CO2'
      diagUnits = '10^-6 ppm       '
      diagCode  = 'SM      U1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )

#endif /* ALLOW_DIAGNOSTICS */

      RETURN
      END
1	dimitri	1.2	C $Header: /u/gcmpack/MITgcm_contrib/ecco_darwin/v4_3deg/code/exf_diagnostics_init.F,v 1.1 2020/01/23 01:27:38 dimitri Exp $
2	dimitri	1.1	C $Name: $
3
4			#include "EXF_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: EXF_DIAGNOSTICS_INIT
8			C !INTERFACE:
9			SUBROUTINE EXF_DIAGNOSTICS_INIT( myThid )
10
11			C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| SUBROUTINE EXF_DIAGNOSTICS_INIT
14			C \| o Routine to initialize list of all available diagnostics
15			C \| for EXF package
16			C ==========================================================
17			C \ev
18			C !USES:
19			IMPLICIT NONE
20
21			C === Global variables ===
22			#include "EEPARAMS.h"
23			#include "EXF_PARAM.h"
24
25			C !INPUT/OUTPUT PARAMETERS:
26			C === Routine arguments ===
27			C myThid :: my Thread Id number
28			INTEGER myThid
29			CEOP
30
31			#ifdef ALLOW_DIAGNOSTICS
32			C !LOCAL VARIABLES:
33			C === Local variables ===
34			C diagNum :: diagnostics number in the (long) list of available diag.
35			C diagName :: local short name (8c) of a diagnostics
36			C diagCode :: local parser field with characteristics of the diagnostics
37			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
39			C diagTitle :: local string (80c): description of field in diagnostic
40			INTEGER diagNum
41			CHARACTER*8 diagName
42			CHARACTER*16 diagCode
43			CHARACTER*16 diagUnits
44			CHARACTER*(80) diagTitle
45
46			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
47
48			diagName = 'EXFhs '
49			diagTitle = 'Sensible heat flux into ocean, >0 increases theta'
50			diagUnits = 'W/m^2 '
51			diagCode = 'SM U1 '
52			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
53			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
54
55			diagName = 'EXFhl '
56			diagTitle = 'Latent heat flux into ocean, >0 increases theta'
57			diagUnits = 'W/m^2 '
58			diagCode = 'SM U1 '
59			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
60			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
61
62			diagName = 'EXFlwnet'
63			diagTitle = 'Net upward longwave radiation, >0 decreases theta'
64			diagUnits = 'W/m^2 '
65			diagCode = 'SM U1 '
66			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
67			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
68
69			diagName = 'EXFswnet'
70			diagTitle = 'Net upward shortwave radiation, >0 decreases theta'
71			diagUnits = 'W/m^2 '
72			diagCode = 'SM U1 '
73			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
74			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
75
76			diagName = 'EXFlwdn '
77			diagTitle = 'Downward longwave radiation, >0 increases theta'
78			diagUnits = 'W/m^2 '
79			diagCode = 'SM U1 '
80			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
81			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
82
83			diagName = 'EXFswdn '
84			diagTitle = 'Downward shortwave radiation, >0 increases theta'
85			diagUnits = 'W/m^2 '
86			diagCode = 'SM U1 '
87			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
88			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
89
90			diagName = 'EXFqnet '
91			diagTitle = 'Net upward heat flux (turb+rad), >0 decreases theta'
92			diagUnits = 'W/m^2 '
93			diagCode = 'SM U1 '
94			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
95			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
96
97			C-- surface fluxes:
98			diagName = 'EXFtaux '
99			diagTitle = 'zonal surface wind stress, >0 increases uVel '
100			diagUnits = 'N/m^2 '
101			IF ( stressIsOnCgrid ) THEN
102			diagCode = 'UU U1 '
103			ELSE
104			diagCode = 'UM U1 '
105			ENDIF
106			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
107			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
108
109			diagName = 'EXFtauy '
110			diagTitle = 'meridional surface wind stress, >0 increases vVel'
111			diagUnits = 'N/m^2 '
112			IF ( stressIsOnCgrid ) THEN
113			diagCode = 'VV U1 '
114			ELSE
115			diagCode = 'VM U1 '
116			ENDIF
117			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
118			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
119
120			diagName = 'EXFuwind'
121			diagTitle = 'zonal 10-m wind speed, >0 eastward'
122			diagUnits = 'm/s '
123			diagCode = 'UM U1 '
124			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
125			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
126
127			diagName = 'EXFvwind'
128			diagTitle = 'meridional 10-m wind speed, >0 northward'
129			diagUnits = 'm/s '
130			diagCode = 'VM U1 '
131			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
132			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
133
134			diagName = 'EXFwspee'
135			diagTitle = '10-m wind speed modulus ( >= 0 )'
136			diagUnits = 'm/s '
137			diagCode = 'SM U1 '
138			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
139			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
140
141			diagName = 'EXFatemp'
142			diagTitle = 'surface (2-m) air temperature'
143			diagUnits = 'degK '
144			diagCode = 'SM U1 '
145			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
146			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
147
148			diagName = 'EXFaqh '
149			diagTitle = 'surface (2-m) specific humidity'
150			diagUnits = 'kg/kg '
151			diagCode = 'SM U1 '
152			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
153			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
154
155			diagName = 'EXFevap '
156			diagTitle = 'evaporation, > 0 increases salinity'
157			diagUnits = 'm/s '
158			diagCode = 'SM U1 '
159			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
160			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
161
162			diagName = 'EXFpreci'
163			diagTitle = 'precipitation, > 0 decreases salinity'
164			diagUnits = 'm/s '
165			diagCode = 'SM U1 '
166			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
167			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
168
169			diagName = 'EXFsnow '
170			diagTitle = 'snow precipitation, > 0 decreases salinity'
171			diagUnits = 'm/s '
172			diagCode = 'SM U1 '
173			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
174			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
175
176			diagName = 'EXFempmr'
177			diagTitle = 'net upward freshwater flux, > 0 increases salinity'
178			diagUnits = 'm/s '
179			diagCode = 'SM U1 '
180			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
181			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
182
183			diagName = 'EXFpress'
184			diagTitle = 'atmospheric pressure field'
185			diagUnits = 'N/m^2 '
186			diagCode = 'SM U1 '
187			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
188			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
189
190			diagName = 'EXFroff '
191			diagTitle = 'river runoff, > 0 decreases salinity'
192			diagUnits = 'm/s '
193			diagCode = 'SM U1 '
194			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
195			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
196
197			diagName = 'EXFroft '
198			diagTitle = 'river runoff temperature'
199			diagUnits = 'deg C '
200			diagCode = 'SM U1 '
201			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
202			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
203
204			diagName = 'EXFsalfx'
205			diagTitle = 'upward salt flux, > 0 decreases salinity'
206			diagUnits = 'psu.kg/m^2/s '
207			diagCode = 'SM U1 '
208			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
209			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
210
211			diagName = 'EXFapco2 '
212			diagTitle = 'Atmospheric partial pressure of CO2'
213			diagUnits = '10^-6 ppm '
214			diagCode = 'SM U1 '
215			CALL DIAGNOSTICS_ADDTOLIST( diagNum,
216			I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
217
218			#endif /* ALLOW_DIAGNOSTICS */
219
220			RETURN
221			END