/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/darwin_mnc_init.F
ViewVC logotype

Contents of /MITgcm_contrib/darwin2/pkg/darwin/darwin_mnc_init.F

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

Revision 1.1 - (show annotations) (download)
Wed Apr 13 18:56:24 2011 UTC (14 years, 4 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt62v_20110413, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt62y_20110526, ctrb_darwin2_ckpt62x_20110513, ctrb_darwin2_ckpt62w_20110426, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt63c_20111011, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63e_20111107, ctrb_darwin2_ckpt63b_20110830, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt63a_20110804, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63d_20111107, ctrb_darwin2_ckpt63_20110728, ctrb_darwin2_baseline, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317, ctrb_darwin2_ckpt62z_20110622 
darwin2 initial checkin
1 C $Header$
2 C $Name$
3
4 #include "DARWIN_OPTIONS.h"
5
6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7 CBOP
8 C !ROUTINE: DARWIN_MNC_INIT
9
10 C !INTERFACE:
11 SUBROUTINE DARWIN_MNC_INIT( myThid )
12
13 C !DESCRIPTION:
14 C Initialize MNC for DARWIN
15
16 C !USES:
17 IMPLICIT NONE
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "PTRACERS_SIZE.h"
22 #include "PTRACERS_PARAMS.h"
23 #include "DARWIN_SIZE.h"
24 #include "SPECTRAL_SIZE.h"
25
26 C !INPUT PARAMETERS:
27 C myThid :: thread number
28 INTEGER myThid
29 CEOP
30 INTEGER np
31 INTEGER ii,jj,ilnb
32 CHARACTER*(MAX_LEN_MBUF) suff
33 CHARACTER*(1) pf
34 CHARACTER*(MAX_LEN_FNAM) name
35 #ifdef ALLOW_QUOTA
36 INTEGER iPhy
37 #endif
38 C Functions
39 integer ILNBLNK
40
41 #ifdef ALLOW_DARWIN
42 #ifdef ALLOW_MNC
43
44 IF ( useMNC ) THEN
45
46 #ifdef ALLOW_QUOTA
47 iPhy=ibiomass
48 #endif
49
50 CALL MNC_CW_ADD_VNAME(
51 & 'dar_PAR_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
52 CALL MNC_CW_ADD_VATTR_TEXT(
53 & 'dar_PAR_ave','units','--', myThid)
54 CALL MNC_CW_ADD_VATTR_TEXT('dar_PAR_ave','description',
55 & '', myThid)
56
57 CALL MNC_CW_ADD_VNAME(
58 & 'dar_PP_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
59 CALL MNC_CW_ADD_VATTR_TEXT(
60 & 'dar_PP_ave','units','--', myThid)
61 CALL MNC_CW_ADD_VATTR_TEXT('dar_PP_ave','description',
62 & '', myThid)
63 #ifndef ALLOW_QUOTA
64 CALL MNC_CW_ADD_VNAME(
65 & 'dar_Chl_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
66 CALL MNC_CW_ADD_VATTR_TEXT(
67 & 'dar_Chl_ave','units','--', myThid)
68 CALL MNC_CW_ADD_VATTR_TEXT('dar_Chl_ave','description',
69 & '', myThid)
70 #endif
71 CALL MNC_CW_ADD_VNAME(
72 & 'dar_Nfix_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
73 CALL MNC_CW_ADD_VATTR_TEXT(
74 & 'dar_Nfix_ave','units','--', myThid)
75 CALL MNC_CW_ADD_VATTR_TEXT('dar_Nfix_ave','description',
76 & '', myThid)
77
78 CALL MNC_CW_ADD_VNAME(
79 & 'dar_Denit_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
80 CALL MNC_CW_ADD_VATTR_TEXT(
81 & 'dar_Denit_ave','units','--', myThid)
82 CALL MNC_CW_ADD_VATTR_TEXT('dar_Denit_ave','description',
83 & '', myThid)
84
85 #ifdef DAR_DIAG_ACDOM
86 CALL MNC_CW_ADD_VNAME(
87 & 'dar_aCDOM_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
88 CALL MNC_CW_ADD_VATTR_TEXT(
89 & 'dar_aCDOM_ave','units','--', myThid)
90 CALL MNC_CW_ADD_VATTR_TEXT('dar_aCDOM_ave','description',
91 & '', myThid)
92 #endif
93
94 #ifdef DAR_DIAG_IRR
95 do ii=1,tlam
96 DO jj = 1,MAX_LEN_FNAM
97 name(jj:jj) = ' '
98 ENDDO
99 WRITE(name,'(a,i2.2)') 'dar_Ed',ii
100 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
101 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
102 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
103 WRITE(name,'(a,i2.2)') 'dar_Es',ii
104 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
105 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
106 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
107 WRITE(name,'(a,i2.2)') 'dar_Eu',ii
108 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
109 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
110 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
111 WRITE(name,'(a,i2.2)') 'dar_Eutop',ii
112 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
113 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
114 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
115 enddo
116 #endif
117
118 #ifdef DAR_DIAG_ABSORP
119 do ii=1,tlam
120 DO jj = 1,MAX_LEN_FNAM
121 name(jj:jj) = ' '
122 ENDDO
123 WRITE(name,'(a,i2.2)') 'dar_a',ii
124 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
125 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
126 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
127 enddo
128 #endif
129 #ifdef DAR_DIAG_SCATTER
130 do ii=1,tlam
131 DO jj = 1,MAX_LEN_FNAM
132 name(jj:jj) = ' '
133 ENDDO
134 WRITE(name,'(a,i2.2)') 'dar_bt',ii
135 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
136 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
137 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
138 WRITE(name,'(a,i2.2)') 'dar_bb',ii
139 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
140 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
141 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
142 enddo
143 #endif
144 #ifdef DAR_DIAG_PART_SCATTER
145 do ii=1,tlam
146 DO jj = 1,MAX_LEN_FNAM
147 name(jj:jj) = ' '
148 ENDDO
149 WRITE(name,'(a,i2.2)') 'dar_apart',ii
150 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
151 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
152 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
153 WRITE(name,'(a,i2.2)') 'dar_btpart',ii
154 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
155 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
156 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
157 WRITE(name,'(a,i2.2)') 'dar_bbpart',ii
158 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
159 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
160 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
161 enddo
162 #endif
163
164 c ANNA_TAVE
165 #ifdef WAVES_DIAG_PCHL
166 do np=1,npmax
167 ii=iPhy+np-1
168 DO jj = 1,MAX_LEN_FNAM
169 name(jj:jj) = ' '
170 ENDDO
171 ilnb = ILNBLNK(PTRACERS_names(ii))
172 WRITE(name,'(a,a)') 'dar_Pchl_',PTRACERS_names(ii)(1:ilnb)
173 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
174 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
175 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
176 enddo
177 #endif
178 c ANNA end TAVE
179 #ifdef DAR_DIAG_RSTAR
180 do np=1,npmax
181 ii=iPhy+np-1
182 DO jj = 1,MAX_LEN_FNAM
183 name(jj:jj) = ' '
184 ENDDO
185 ilnb = ILNBLNK(PTRACERS_names(ii))
186 WRITE(name,'(a,a)') 'dar_Rstar_',PTRACERS_names(ii)(1:ilnb)
187 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
188 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
189 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
190 enddo
191 do np=1,npmax
192 ii=iPhy+np-1
193 DO jj = 1,MAX_LEN_FNAM
194 name(jj:jj) = ' '
195 ENDDO
196 ilnb = ILNBLNK(PTRACERS_names(ii))
197 WRITE(name,'(a,a)') 'dar_RNstar_',PTRACERS_names(ii)(1:ilnb)
198 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
199 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
200 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
201 enddo
202 #endif
203
204 #ifdef DAR_DIAG_DIVER
205 CALL MNC_CW_ADD_VNAME(
206 & 'dar_Diver1_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
207 CALL MNC_CW_ADD_VATTR_TEXT(
208 & 'dar_Diver1_ave','units','--', myThid)
209 CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver1_ave','description',
210 & '', myThid)
211
212 CALL MNC_CW_ADD_VNAME(
213 & 'dar_Diver2_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
214 CALL MNC_CW_ADD_VATTR_TEXT(
215 & 'dar_Diver2_ave','units','--', myThid)
216 CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver2_ave','description',
217 & '', myThid)
218
219 CALL MNC_CW_ADD_VNAME(
220 & 'dar_Diver3_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
221 CALL MNC_CW_ADD_VATTR_TEXT(
222 & 'dar_Diver3_ave','units','--', myThid)
223 CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver3_ave','description',
224 & '', myThid)
225
226 CALL MNC_CW_ADD_VNAME(
227 & 'dar_Diver4_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
228 CALL MNC_CW_ADD_VATTR_TEXT(
229 & 'dar_Diver4_ave','units','--', myThid)
230 CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver4_ave','description',
231 & '', myThid)
232
233 #endif
234 #ifdef DAR_DIAG_GROW
235 do np=1,npmax
236 ii=iPhy+np-1
237 DO jj = 1,MAX_LEN_FNAM
238 name(jj:jj) = ' '
239 ENDDO
240 ilnb = ILNBLNK(PTRACERS_names(ii))
241 WRITE(name,'(a,a)') 'dar_Grow_',PTRACERS_names(ii)(1:ilnb)
242 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
243 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
244 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
245 enddo
246 do np=1,npmax
247 ii=iPhy+np-1
248 DO jj = 1,MAX_LEN_FNAM
249 name(jj:jj) = ' '
250 ENDDO
251 ilnb = ILNBLNK(PTRACERS_names(ii))
252 WRITE(name,'(a,a)') 'dar_Growsq_',PTRACERS_names(ii)(1:ilnb)
253 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
254 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
255 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
256 enddo
257 #endif
258
259 #ifdef ALLOW_DIAZ
260 #ifdef DAR_DIAG_NFIXP
261 do np=1,npmax
262 ii=iPhy+np-1
263 DO jj = 1,MAX_LEN_FNAM
264 name(jj:jj) = ' '
265 ENDDO
266 ilnb = ILNBLNK(PTRACERS_names(ii))
267 WRITE(name,'(a,a)') 'dar_NfixP_',PTRACERS_names(ii)(1:ilnb)
268 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
269 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
270 CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
271 enddo
272 #endif
273 #endif
274
275
276
277 c CALL MNC_CW_ADD_VNAME(
278 c & 'dic_pH_ave', 'Cen_xy_Hn__-__t', 3,4, myThid)
279 c CALL MNC_CW_ADD_VATTR_TEXT(
280 c & 'dic_pH_ave','units','--', myThid)
281 c CALL MNC_CW_ADD_VATTR_TEXT('dic_pH_ave','description',
282 c & '', myThid)
283
284 ENDIF
285
286 #endif /* ALLOW_MNC */
287 #endif /* ALLOW_DARWIN */
288
289 RETURN
290 END
  ViewVC Help
Powered by ViewVC 1.1.22