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

Annotation 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 - (hide annotations) (download)
Wed Jan 6 00:36:02 2016 UTC (8 years, 4 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 jmc 1.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