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

Annotation 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.3 - (hide annotations) (download)
Mon Jul 30 15:21:51 2012 UTC (12 years, 11 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt63r_20120817, ctrb_darwin2_ckpt63q_20120731
Changes since 1.2: +25 -1 lines 
add timeave diagnostics rmud and c1,c2 for irradiance amplitudes
1 jahn 1.3 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_mnc_init.F,v 1.2 2012/07/02 09:40:58 benw Exp $
2 benw 1.2 C $Name: $
3 jahn 1.1
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 benw 1.2 cbenw limitation diagnostics
64     #ifdef QUOTA_DIAG_LIMIT
65     do np=1,npmax
66     c nitrogen
67     DO jj = 1,MAX_LEN_FNAM
68     name(jj:jj) = ' '
69     ENDDO
70     WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
71     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
72     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
73     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
74     c iron
75     DO jj = 1,MAX_LEN_FNAM
76     name(jj:jj) = ' '
77     ENDDO
78     WRITE(name,'(a,i3.3)') 'dar_Flim_',np
79     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
80     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
81     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
82     c light
83     DO jj = 1,MAX_LEN_FNAM
84     name(jj:jj) = ' '
85     ENDDO
86     WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
87     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
88     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
89     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
90     enddo
91     c temperature
92     DO jj = 1,MAX_LEN_FNAM
93     name(jj:jj) = ' '
94     ENDDO
95     WRITE(name,'(a)') 'dar_Tlim_ave'
96     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
97     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
98     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
99     #endif
100     cbenw limitation diagnostics
101 jahn 1.1 #ifndef ALLOW_QUOTA
102     CALL MNC_CW_ADD_VNAME(
103     & 'dar_Chl_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
104     CALL MNC_CW_ADD_VATTR_TEXT(
105     & 'dar_Chl_ave','units','--', myThid)
106     CALL MNC_CW_ADD_VATTR_TEXT('dar_Chl_ave','description',
107     & '', myThid)
108     #endif
109     CALL MNC_CW_ADD_VNAME(
110     & 'dar_Nfix_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
111     CALL MNC_CW_ADD_VATTR_TEXT(
112     & 'dar_Nfix_ave','units','--', myThid)
113     CALL MNC_CW_ADD_VATTR_TEXT('dar_Nfix_ave','description',
114     & '', myThid)
115    
116     CALL MNC_CW_ADD_VNAME(
117     & 'dar_Denit_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
118     CALL MNC_CW_ADD_VATTR_TEXT(
119     & 'dar_Denit_ave','units','--', myThid)
120     CALL MNC_CW_ADD_VATTR_TEXT('dar_Denit_ave','description',
121     & '', myThid)
122    
123     #ifdef DAR_DIAG_ACDOM
124     CALL MNC_CW_ADD_VNAME(
125     & 'dar_aCDOM_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
126     CALL MNC_CW_ADD_VATTR_TEXT(
127     & 'dar_aCDOM_ave','units','--', myThid)
128     CALL MNC_CW_ADD_VATTR_TEXT('dar_aCDOM_ave','description',
129     & '', myThid)
130     #endif
131    
132     #ifdef DAR_DIAG_IRR
133     do ii=1,tlam
134     DO jj = 1,MAX_LEN_FNAM
135     name(jj:jj) = ' '
136     ENDDO
137     WRITE(name,'(a,i2.2)') 'dar_Ed',ii
138     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
139     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
140     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
141     WRITE(name,'(a,i2.2)') 'dar_Es',ii
142     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
143     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
144     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
145     WRITE(name,'(a,i2.2)') 'dar_Eu',ii
146     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
147     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
148     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
149     WRITE(name,'(a,i2.2)') 'dar_Eutop',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     enddo
154     #endif
155    
156 jahn 1.3 #ifdef DAR_DIAG_IRR_AMPS
157     do ii=1,tlam
158     DO jj = 1,MAX_LEN_FNAM
159     name(jj:jj) = ' '
160     ENDDO
161     WRITE(name,'(a,i2.2)') 'dar_c1_',ii
162     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
163     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
164     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
165     WRITE(name,'(a,i2.2)') 'dar_c2_',ii
166     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
167     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
168     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
169     enddo
170     #endif
171 jahn 1.1 #ifdef DAR_DIAG_ABSORP
172     do ii=1,tlam
173     DO jj = 1,MAX_LEN_FNAM
174     name(jj:jj) = ' '
175     ENDDO
176     WRITE(name,'(a,i2.2)') 'dar_a',ii
177     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
178     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
179     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
180     enddo
181     #endif
182     #ifdef DAR_DIAG_SCATTER
183     do ii=1,tlam
184     DO jj = 1,MAX_LEN_FNAM
185     name(jj:jj) = ' '
186     ENDDO
187     WRITE(name,'(a,i2.2)') 'dar_bt',ii
188     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
189     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
190     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
191     WRITE(name,'(a,i2.2)') 'dar_bb',ii
192     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
193     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
194     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
195     enddo
196     #endif
197     #ifdef DAR_DIAG_PART_SCATTER
198     do ii=1,tlam
199     DO jj = 1,MAX_LEN_FNAM
200     name(jj:jj) = ' '
201     ENDDO
202     WRITE(name,'(a,i2.2)') 'dar_apart',ii
203     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
204     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
205     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
206     WRITE(name,'(a,i2.2)') 'dar_btpart',ii
207     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
208     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
209     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
210     WRITE(name,'(a,i2.2)') 'dar_bbpart',ii
211     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
212     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
213     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
214     enddo
215     #endif
216 jahn 1.3 #ifdef DAR_RADTRANS
217     DO jj = 1,MAX_LEN_FNAM
218     name(jj:jj) = ' '
219     ENDDO
220     WRITE(name,'(a)') 'dar_rmud'
221     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)
222     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
223     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
224     #endif
225 jahn 1.1
226     c ANNA_TAVE
227     #ifdef WAVES_DIAG_PCHL
228     do np=1,npmax
229     ii=iPhy+np-1
230     DO jj = 1,MAX_LEN_FNAM
231     name(jj:jj) = ' '
232     ENDDO
233     ilnb = ILNBLNK(PTRACERS_names(ii))
234     WRITE(name,'(a,a)') 'dar_Pchl_',PTRACERS_names(ii)(1:ilnb)
235     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
236     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
237     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
238     enddo
239     #endif
240     c ANNA end TAVE
241     #ifdef DAR_DIAG_RSTAR
242     do np=1,npmax
243     ii=iPhy+np-1
244     DO jj = 1,MAX_LEN_FNAM
245     name(jj:jj) = ' '
246     ENDDO
247     ilnb = ILNBLNK(PTRACERS_names(ii))
248     WRITE(name,'(a,a)') 'dar_Rstar_',PTRACERS_names(ii)(1:ilnb)
249     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
250     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
251     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
252     enddo
253     do np=1,npmax
254     ii=iPhy+np-1
255     DO jj = 1,MAX_LEN_FNAM
256     name(jj:jj) = ' '
257     ENDDO
258     ilnb = ILNBLNK(PTRACERS_names(ii))
259     WRITE(name,'(a,a)') 'dar_RNstar_',PTRACERS_names(ii)(1:ilnb)
260     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
261     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
262     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
263     enddo
264     #endif
265    
266     #ifdef DAR_DIAG_DIVER
267     CALL MNC_CW_ADD_VNAME(
268     & 'dar_Diver1_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
269     CALL MNC_CW_ADD_VATTR_TEXT(
270     & 'dar_Diver1_ave','units','--', myThid)
271     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver1_ave','description',
272     & '', myThid)
273    
274     CALL MNC_CW_ADD_VNAME(
275     & 'dar_Diver2_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
276     CALL MNC_CW_ADD_VATTR_TEXT(
277     & 'dar_Diver2_ave','units','--', myThid)
278     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver2_ave','description',
279     & '', myThid)
280    
281     CALL MNC_CW_ADD_VNAME(
282     & 'dar_Diver3_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
283     CALL MNC_CW_ADD_VATTR_TEXT(
284     & 'dar_Diver3_ave','units','--', myThid)
285     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver3_ave','description',
286     & '', myThid)
287    
288     CALL MNC_CW_ADD_VNAME(
289     & 'dar_Diver4_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
290     CALL MNC_CW_ADD_VATTR_TEXT(
291     & 'dar_Diver4_ave','units','--', myThid)
292     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver4_ave','description',
293     & '', myThid)
294    
295     #endif
296     #ifdef DAR_DIAG_GROW
297     do np=1,npmax
298     ii=iPhy+np-1
299     DO jj = 1,MAX_LEN_FNAM
300     name(jj:jj) = ' '
301     ENDDO
302     ilnb = ILNBLNK(PTRACERS_names(ii))
303     WRITE(name,'(a,a)') 'dar_Grow_',PTRACERS_names(ii)(1:ilnb)
304     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
305     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
306     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
307     enddo
308     do np=1,npmax
309     ii=iPhy+np-1
310     DO jj = 1,MAX_LEN_FNAM
311     name(jj:jj) = ' '
312     ENDDO
313     ilnb = ILNBLNK(PTRACERS_names(ii))
314     WRITE(name,'(a,a)') 'dar_Growsq_',PTRACERS_names(ii)(1:ilnb)
315     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
316     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
317     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
318     enddo
319     #endif
320    
321     #ifdef ALLOW_DIAZ
322     #ifdef DAR_DIAG_NFIXP
323     do np=1,npmax
324     ii=iPhy+np-1
325     DO jj = 1,MAX_LEN_FNAM
326     name(jj:jj) = ' '
327     ENDDO
328     ilnb = ILNBLNK(PTRACERS_names(ii))
329     WRITE(name,'(a,a)') 'dar_NfixP_',PTRACERS_names(ii)(1:ilnb)
330     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
331     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
332     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
333     enddo
334     #endif
335     #endif
336    
337    
338    
339     c CALL MNC_CW_ADD_VNAME(
340     c & 'dic_pH_ave', 'Cen_xy_Hn__-__t', 3,4, myThid)
341     c CALL MNC_CW_ADD_VATTR_TEXT(
342     c & 'dic_pH_ave','units','--', myThid)
343     c CALL MNC_CW_ADD_VATTR_TEXT('dic_pH_ave','description',
344     c & '', myThid)
345    
346     ENDIF
347    
348     #endif /* ALLOW_MNC */
349     #endif /* ALLOW_DARWIN */
350    
351     RETURN
352     END
  ViewVC Help
Powered by ViewVC 1.1.22