/[MITgcm]/MITgcm/pkg/icefront/icefront_diagnostics_init.F
ViewVC logotype

Contents of /MITgcm/pkg/icefront/icefront_diagnostics_init.F

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

Revision 1.1 - (show annotations) (download)
Wed Jan 20 23:33:45 2010 UTC (14 years, 3 months ago) by dimitri
Branch: MAIN 
Adding the shell of, and the hooks to, a new package that will be
used to model melting and freezing of vertical glacier ice fronts:
 Modified Files:
 	doc/tag-index model/inc/PARAMS.h model/src/do_oceanic_phys.F
 	model/src/external_forcing.F model/src/packages_boot.F
 	model/src/packages_check.F model/src/packages_init_fixed.F
 	model/src/packages_init_variables.F
 	model/src/packages_readparms.F
 Added Files:
 	pkg/icefront/ICEFRONT.h pkg/icefront/ICEFRONT_OPTIONS.h
 	pkg/icefront/icefront_check.F
 	pkg/icefront/icefront_description.tex
 	pkg/icefront/icefront_diagnostics_init.F
 	pkg/icefront/icefront_init_fixed.F
 	pkg/icefront/icefront_init_varia.F
 	pkg/icefront/icefront_readparms.F
 	pkg/icefront/icefront_tendency_apply.F
 	pkg/icefront/icefront_thermodynamics.F
1 C $Header: $
2 C $Name: $
3
4 #include "ICEFRONT_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP 0
8 C !ROUTINE: ICEFRONT_DIAGNOSTICS_INIT
9
10 C !INTERFACE:
11 SUBROUTINE ICEFRONT_DIAGNOSTICS_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize list of all available diagnostics
15
16 C !USES:
17 IMPLICIT NONE
18 #include "EEPARAMS.h"
19 #include "SIZE.h"
20 #include "ICEFRONT.h"
21
22 C !INPUT/OUTPUT PARAMETERS:
23 C myThid :: my Thread Id number
24 INTEGER myThid
25 CEOP
26
27 #ifdef ALLOW_DIAGNOSTICS
28 C !LOCAL VARIABLES:
29 C === Local variables ===
30 C diagNum :: diagnostics number in the (long) list of available diag.
31 C diagMate :: diag. mate number in the (long) list of available diag.
32 C diagName :: local short name (8c) of a diagnostics
33 C diagCode :: local parser field with characteristics of the diagnostics
34 C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
35 C diagUnits :: local string (16c): physical units of a diagnostic field
36 C diagTitle :: local string (80c): description of field in diagnostic
37 INTEGER diagNum
38 INTEGER diagMate
39 CHARACTER*8 diagName
40 CHARACTER*16 diagCode
41 CHARACTER*16 diagUnits
42 CHARACTER*(80) diagTitle
43
44 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
45
46 c IF ( useDiagnotics ) THEN
47
48 diagName = 'MYPaSta1'
49 diagTitle = 'Scalar state variable 1 for package ICEFRONT'
50 diagUnits = 'unit '
51 diagCode = 'SM MR '
52 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
53 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
54
55 diagName = 'MYPaSta2'
56 diagTitle = 'Scalar state variable 2 for package ICEFRONT'
57 diagUnits = 'unit '
58 diagCode = 'SM MR '
59 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
60 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
61
62 diagName = 'MYPaStaU'
63 diagTitle = 'uVel-type state variable for package ICEFRONT'
64 diagUnits = 'unit '
65 diagCode = 'UM MR '
66 diagMate = diagNum + 2
67 IF ( myPa_StaV_Cgrid ) diagCode(2:2) = 'U'
68 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
69 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
70
71 diagName = 'MYPaStaV'
72 diagTitle = 'vVel-type state variable for package ICEFRONT'
73 diagUnits = 'unit '
74 diagCode = 'VM MR '
75 diagMate = diagNum
76 IF ( myPa_StaV_Cgrid ) diagCode(2:2) = 'V'
77 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
78 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
79
80 diagName = 'MYPaSur1'
81 diagTitle = 'Surface field 1 for package ICEFRONT'
82 diagUnits = 'unit '
83 diagCode = 'SM L1 '
84 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
85 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
86
87 diagName = 'MYPaSur2'
88 diagTitle = 'Surface field 2 for package ICEFRONT'
89 diagUnits = 'unit '
90 diagCode = 'SM L1 '
91 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
92 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
93
94 diagName = 'MYPadTdt'
95 diagTitle = ' Pot.Temp. tendency from package ICEFRONT'
96 diagUnits = 'K/s '
97 diagCode = 'SM MR '
98 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
99 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
100
101 diagName = 'MYPadSdt'
102 diagTitle = ' Salinity tendency from package ICEFRONT'
103 diagUnits = 'psu/s '
104 diagCode = 'SM MR '
105 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
106 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
107
108 diagName = 'MYPadUdt'
109 diagTitle = 'u-Velocity tendency from package ICEFRONT'
110 diagUnits = 'm/s^2 '
111 diagCode = 'UM MR '
112 diagMate = diagNum + 2
113 IF ( myPa_Tend_Cgrid ) diagCode(2:2) = 'U'
114 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
115 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
116
117 diagName = 'MYPadVdt'
118 diagTitle = 'v-Velocity tendency from package ICEFRONT'
119 diagUnits = 'm/s^2 '
120 diagCode = 'VM MR '
121 diagMate = diagNum
122 IF ( myPa_Tend_Cgrid ) diagCode(2:2) = 'V'
123 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
124 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
125
126 C-- an example of how to set "free levels option" diagnostics
127 diagName = 'MYPaLev3'
128 diagTitle = '3-levels field for package ICEFRONT'
129 diagUnits = 'unit '
130 diagCode = 'SM MX '
131 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
132 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
133 CALL DIAGNOSTICS_SETKLEV( diagName, 3, myThid )
134
135 c ENDIF
136
137 #endif /* ALLOW_DIAGNOSTICS */
138
139 RETURN
140 END
  ViewVC Help
Powered by ViewVC 1.1.22