/[MITgcm]/MITgcm/pkg/atm_compon_interf/cpl_diagnostics_init.F
ViewVC logotype

Contents of /MITgcm/pkg/atm_compon_interf/cpl_diagnostics_init.F

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

Revision 1.1 - (show annotations) (download)
Wed Jan 6 00:36:02 2016 UTC (8 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65s, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, HEAD 
- add diagnostics for all imported fields from coupler (in both OCN & ATM);
  this is entended to replace the OCN cplFldtave.* output files (not updated
  to contain new additional exchanged fields).
1 C $Header: /u/gcmpack/MITgcm/pkg/mypackage/mypackage_diagnostics_init.F,v 1.5 2010/01/15 00:28:14 jmc Exp $
2 C $Name: $
3
4 #include "ATM_CPL_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 0
8 C !ROUTINE: CPL_DIAGNOSTICS_INIT
9
10 C !INTERFACE:
11 SUBROUTINE CPL_DIAGNOSTICS_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize list of available diagnostics
15 C for Coupled Atmospheric Component
16
17 C !USES:
18 IMPLICIT NONE
19 #include "EEPARAMS.h"
20 #include "SIZE.h"
21 #include "CPL_PARAMS.h"
22
23 C !INPUT/OUTPUT PARAMETERS:
24 C myThid :: my Thread Id number
25 INTEGER myThid
26 CEOP
27
28 #ifdef ALLOW_DIAGNOSTICS
29 C !LOCAL VARIABLES:
30 C === Local variables ===
31 C diagNum :: diagnostics number in the (long) list of available diag.
32 C diagMate :: diag. mate number in the (long) list of available diag.
33 C diagName :: local short name (8c) of a diagnostics
34 C diagCode :: local parser field with characteristics of the diagnostics
35 C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
36 C diagUnits :: local string (16c): physical units of a diagnostic field
37 C diagTitle :: local string (80c): description of field in diagnostic
38 INTEGER diagNum
39 c INTEGER diagMate
40 CHARACTER*8 diagName
41 CHARACTER*16 diagCode
42 CHARACTER*16 diagUnits
43 CHARACTER*(80) diagTitle
44 #ifdef ALLOW_THSICE
45 INTEGER numFract, numThick
46 #endif /* ALLOW_THSICE */
47
48 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
49
50 c IF ( useDiagnotics ) THEN
51
52 diagName = 'CPL_MxlD'
53 diagTitle = 'Ocean mixed-layer depth [m] (recv from Coupler)'
54 diagUnits = 'm '
55 diagCode = 'SM M1 '
56 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
57 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
58
59 diagName = 'CPL_SST '
60 diagTitle = 'Ocean surface temperature (recv from Coupler)'
61 diagUnits = 'degC '
62 diagCode = 'SM M1 '
63 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
64 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
65
66 diagName = 'CPL_SSS '
67 diagTitle = 'Ocean surface salinity (recv from Coupler)'
68 diagUnits = 'g/kg '
69 diagCode = 'SM M1 '
70 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
71 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
72
73 diagName = 'CPL_vSq '
74 diagName = 'CPL_Vel2'
75 diagTitle = 'Ocean surface velocity square (recv from Coupler)'
76 diagUnits = 'm^2/s^2 '
77 diagCode = 'SM M1 '
78 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80
81 #ifdef ALLOW_AIM
82 diagName = 'CPL_fCO2'
83 diagTitle ='Ocean air-sea flux of CO2 (+=down)(recv fr Coupler)'
84 diagUnits = 'mol/m^2/s '
85 diagCode = 'SM M1 '
86 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
87 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
88 #endif /* ALLOW_AIM */
89
90 #ifdef ALLOW_THSICE
91 diagName = 'CPL_icFr'
92 diagTitle = 'Seaice fraction (Ocean, recv from Coupler)'
93 diagUnits = '0-1 '
94 diagCode = 'SM P M1 '
95 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
96 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
97 numFract = diagNum
98
99 diagName = 'CPL_iceH'
100 diagTitle = 'Seaice thickness (Ocean, recv from Coupler)'
101 diagUnits = 'm '
102 diagCode = 'SM PC M1 '
103 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
104 I diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
105 numThick = diagNum
106
107 diagName = 'CPL_snwH'
108 diagTitle = 'Snow thickness over seaice (OCN, recv fr Coupler)'
109 diagUnits = 'm '
110 diagCode = 'SM PC M1 '
111 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
112 I diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
113
114 diagName = 'CPL_Qic1'
115 diagTitle ='Enthalpy of seaice layer 1 (Ocean, recv fr Coupler)'
116 diagUnits = 'J/kg '
117 diagCode = 'SM C M1 '
118 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
119 I diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
120
121 diagName = 'CPL_Qic2'
122 diagTitle ='Enthalpy of seaice layer 2 (Ocean, recv fr Coupler)'
123 diagUnits = 'J/kg '
124 diagCode = 'SM C M1 '
125 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
126 I diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
127 #endif /* ALLOW_THSICE */
128
129 c ENDIF
130
131 #endif /* ALLOW_DIAGNOSTICS */
132
133 RETURN
134 END
  ViewVC Help
Powered by ViewVC 1.1.22