/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/monod/monod_diagnostics_init.F

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

Revision 1.2 - (hide annotations) (download)
Wed Nov 9 23:34:17 2011 UTC (13 years, 8 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt64f_20130405, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt64a_20121116, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt64e_20130305, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63s_20120908, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt63r_20120817, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt64c_20130120, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63p_20120707, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt63q_20120731, ctrb_darwin2_ckpt64b_20121224, ctrb_darwin2_ckpt64d_20130219, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt64_20121012, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt65h_20141217, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.1: +18 -2 lines 
add Shannon and Simpson diversity indices
1 jahn 1.2 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/monod/monod_diagnostics_init.F,v 1.1 2011/04/13 18:56:25 jahn Exp $
2     C $Name: $
3 jahn 1.1
4     #include "CPP_OPTIONS.h"
5     #include "PTRACERS_OPTIONS.h"
6     #include "DARWIN_OPTIONS.h"
7    
8     #ifdef ALLOW_PTRACERS
9     #ifdef ALLOW_MONOD
10    
11     CBOP
12     C !ROUTINE: MONOD_DIAGNOSTICS_INIT
13     C !INTERFACE:
14     SUBROUTINE MONOD_DIAGNOSTICS_INIT( myThid )
15     C !DESCRIPTION:
16     C Routine to initialize Darwin diagnostics
17    
18     C !USES:
19     IMPLICIT NONE
20     C === Global variables ===
21     #include "SIZE.h"
22     #include "PTRACERS_SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "GCHEM.h"
26     #include "MONOD_SIZE.h"
27    
28     C !INPUT/OUTPUT PARAMETERS:
29     C === Routine arguments ===
30     C myThid - Number of this instance of MONOD_DIAGNOSTICS_INIT
31     INTEGER myThid
32     CEOP
33    
34     #ifdef ALLOW_DIAGNOSTICS
35     C !LOCAL VARIABLES:
36     C === Local variables ===
37     C msgBuf - Informational/error meesage buffer
38     c CHARACTER*(MAX_LEN_MBUF) msgBuf
39    
40     INTEGER diagNum
41     CHARACTER*8 diagName
42     CHARACTER*16 diagCode
43     CHARACTER*16 diagUnits
44     CHARACTER*(80) diagTitle
45    
46     INTEGER np
47    
48     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
49    
50     C-- Add diagnostics to the (long) list
51    
52     C-- Add diagnostic PAR, PP and Nfix
53     diagName = 'PAR '
54     WRITE(diagTitle, '(A)') 'PAR'
55     diagUnits = 'uEin/m^2/s '
56     diagCode = 'SM P MR '
57     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
58     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
59    
60     diagName = 'PP '
61     WRITE(diagTitle, '(A)') 'Primary Production'
62     diagUnits = 'mmol P/m^3/s '
63     diagCode = 'SM P MR '
64     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
65     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
66    
67     diagName = 'Nfix '
68     WRITE(diagTitle, '(A)') 'N fixation'
69     diagUnits = 'mmol N/m^3/s '
70     diagCode = 'SM P MR '
71     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
72     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
73    
74     diagName = 'Denit '
75     WRITE(diagTitle, '(A)') 'Denitrification'
76     diagUnits = 'mmol N/m^3/s '
77     diagCode = 'SM P MR '
78     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
79     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
80    
81     c ANNA_TAVE
82     #ifdef WAVES_DIAG_PCHL
83     DO np=1,MIN(99,npmax)
84     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
85     C-- Add diagnostics of Pchl
86     c ANNA_Q format to allow more Ptracers I2.3 ??
87     WRITE(diagName,'(A,I2.2,A)') 'Pchl',np,' '
88     WRITE(diagTitle,'(A,I2.2)') 'Chl of Phy', np
89     diagUnits = 'mg chl/m^3 '
90     diagCode = 'SM P MR '
91     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
92     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
93    
94     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
95     ENDDO /* np */
96     #endif /* WAVES_DIAG_PCHL */
97     c ANNA end TAVE
98    
99     #ifdef DAR_DIAG_RSTAR
100     DO np=1,MIN(99,npmax)
101     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
102     C-- Add diagnostics of Rstar
103     WRITE(diagName,'(A,I2.2,A)') 'Rstar',np,' '
104     WRITE(diagTitle,'(A,I2.2)') 'Rstar^PO4 of Phy', np
105     diagUnits = 'mmol P/m^3 '
106     diagCode = 'SM P MR '
107     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
108     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
109    
110     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
111     ENDDO /* np */
112     #endif /* DAR_DIAG_RSTAR */
113    
114     #ifdef DAR_DIAG_DIVER
115     diagName = 'Diver1 '
116     WRITE(diagTitle, '(A)')
117     & '# Phy species over absolute threshold'
118     diagUnits = '1 '
119     diagCode = 'SM P MR '
120     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
121     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
122    
123     diagName = 'Diver2 '
124     WRITE(diagTitle, '(A)')
125     & '# Phy species over threshold for proportion of biomass'
126     diagUnits = '1 '
127     diagCode = 'SM P MR '
128     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
129     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
130    
131     diagName = 'Diver3 '
132     WRITE(diagTitle, '(A)')
133     & '# Phy species constituting bulk of biomass'
134     diagUnits = '1 '
135     diagCode = 'SM P MR '
136     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
137     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
138    
139     diagName = 'Diver4 '
140     WRITE(diagTitle, '(A)')
141     & '# Phy species over threshold for ratio to maximum species'
142     diagUnits = '1 '
143     diagCode = 'SM P MR '
144     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
145     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
146 jahn 1.2
147     diagName = 'Shannon '
148     WRITE(diagTitle, '(A)')
149     & '# Shannon index'
150     diagUnits = '1 '
151     diagCode = 'SM P MR '
152     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
153     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
154    
155     diagName = 'Simpson '
156     WRITE(diagTitle, '(A)')
157     & '# Simpson diversity index 1/sum p*p'
158     diagUnits = '1 '
159     diagCode = 'SM P MR '
160     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
162 jahn 1.1 #endif /* DAR_DIAG_DIVER */
163    
164     #ifdef ALLOW_DIAZ
165     #ifdef DAR_DIAG_NFIXP
166     DO np=1,MIN(99,npmax)
167     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|C-- Add diagnostics of NFix
168     WRITE(diagName,'(A,I2.2,A)') 'NfixP',np,' '
169     WRITE(diagTitle,'(A,I2.2)') 'Nfix of Phy', np
170     diagUnits = 'mmol N/m^3 '
171     diagCode = 'SM P MR '
172     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
173     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
174    
175     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| ENDDO /* np */
176     ENDDO /* np */
177     #endif /* DAR_DIAG_NFIXP */
178     #endif
179    
180     #ifdef ALLOW_PAR_DAY
181     diagName = 'PARday '
182     WRITE(diagTitle, '(A)') 'daily-averaged PAR'
183     diagUnits = 'uEin/m^2/s '
184     diagCode = 'SM P MR '
185     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
186     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
187     #endif /* ALLOW_PAR_DAY */
188    
189     #ifdef DAR_DIAG_CHL
190     diagName = 'ChlGeide'
191     WRITE(diagTitle, '(A)') 'Chlorophyll a concentration a la Geider'
192     diagUnits = 'mg Chla/m^3 '
193     diagCode = 'SM P MR '
194     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
195     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
196     diagName = 'Chl2CGei'
197     WRITE(diagTitle, '(A)')'Chlorophyll a to carbon ratio a la Geider'
198     diagUnits = 'mg Chla/mmol C '
199     diagCode = 'SM P MR '
200     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
201     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
202    
203     diagName = 'ChlDoney'
204     WRITE(diagTitle, '(A)') 'Chlorophyll a concentration a la Doney'
205     diagUnits = 'mg Chla/m^3 '
206     diagCode = 'SM P MR '
207     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
208     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
209    
210     diagName = 'Chl2CDon'
211     WRITE(diagTitle, '(A)') 'Chlorophyll a to carbon ratio a la Doney'
212     diagUnits = 'mg Chla/mmol C '
213     diagCode = 'SM P MR '
214     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
215     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
216    
217     diagName = 'ChlCloer'
218     WRITE(diagTitle, '(A)') 'Chlorophyll a concentration a la Cloern'
219     diagUnits = 'mg Chla/m^3 '
220     diagCode = 'SM P MR '
221     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
222     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
223     diagName = 'Chl2CClo'
224     WRITE(diagTitle, '(A)')'Chlorophyll a to carbon ratio a la Cloern'
225     diagUnits = 'mg Chla/mmol C '
226     diagCode = 'SM P MR '
227     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
228     I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
229     #endif /* DAR_DIAG_CHL */
230    
231     #endif /* ALLOW_DIAGNOSTICS */
232    
233     RETURN
234     END
235     C============================================================================
236     #endif
237     #endif
238    
  ViewVC Help
Powered by ViewVC 1.1.22