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

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

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

Revision 1.1 - (hide annotations) (download)
Wed Jan 20 23:33:45 2010 UTC (14 years, 4 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 dimitri 1.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