/[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.2 - (hide annotations) (download)
Mon Jul 2 09:40:58 2012 UTC (13 years ago) by benw
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt63p_20120707
Changes since 1.1: +40 -2 lines 
Allows for MNC output of nitrogen, iron, temperature and light limitation factors
1 benw 1.2 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_mnc_init.F,v 1.1 2011/04/13 18:56:24 jahn Exp $
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     #ifdef DAR_DIAG_ABSORP
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_a',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     enddo
166     #endif
167     #ifdef DAR_DIAG_SCATTER
168     do ii=1,tlam
169     DO jj = 1,MAX_LEN_FNAM
170     name(jj:jj) = ' '
171     ENDDO
172     WRITE(name,'(a,i2.2)') 'dar_bt',ii
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     WRITE(name,'(a,i2.2)') 'dar_bb',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_PART_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_apart',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_btpart',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     WRITE(name,'(a,i2.2)') 'dar_bbpart',ii
196     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
197     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
198     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
199     enddo
200     #endif
201    
202     c ANNA_TAVE
203     #ifdef WAVES_DIAG_PCHL
204     do np=1,npmax
205     ii=iPhy+np-1
206     DO jj = 1,MAX_LEN_FNAM
207     name(jj:jj) = ' '
208     ENDDO
209     ilnb = ILNBLNK(PTRACERS_names(ii))
210     WRITE(name,'(a,a)') 'dar_Pchl_',PTRACERS_names(ii)(1:ilnb)
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     c ANNA end TAVE
217     #ifdef DAR_DIAG_RSTAR
218     do np=1,npmax
219     ii=iPhy+np-1
220     DO jj = 1,MAX_LEN_FNAM
221     name(jj:jj) = ' '
222     ENDDO
223     ilnb = ILNBLNK(PTRACERS_names(ii))
224     WRITE(name,'(a,a)') 'dar_Rstar_',PTRACERS_names(ii)(1:ilnb)
225     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
226     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
227     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
228     enddo
229     do np=1,npmax
230     ii=iPhy+np-1
231     DO jj = 1,MAX_LEN_FNAM
232     name(jj:jj) = ' '
233     ENDDO
234     ilnb = ILNBLNK(PTRACERS_names(ii))
235     WRITE(name,'(a,a)') 'dar_RNstar_',PTRACERS_names(ii)(1:ilnb)
236     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
237     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
238     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
239     enddo
240     #endif
241    
242     #ifdef DAR_DIAG_DIVER
243     CALL MNC_CW_ADD_VNAME(
244     & 'dar_Diver1_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
245     CALL MNC_CW_ADD_VATTR_TEXT(
246     & 'dar_Diver1_ave','units','--', myThid)
247     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver1_ave','description',
248     & '', myThid)
249    
250     CALL MNC_CW_ADD_VNAME(
251     & 'dar_Diver2_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
252     CALL MNC_CW_ADD_VATTR_TEXT(
253     & 'dar_Diver2_ave','units','--', myThid)
254     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver2_ave','description',
255     & '', myThid)
256    
257     CALL MNC_CW_ADD_VNAME(
258     & 'dar_Diver3_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
259     CALL MNC_CW_ADD_VATTR_TEXT(
260     & 'dar_Diver3_ave','units','--', myThid)
261     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver3_ave','description',
262     & '', myThid)
263    
264     CALL MNC_CW_ADD_VNAME(
265     & 'dar_Diver4_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
266     CALL MNC_CW_ADD_VATTR_TEXT(
267     & 'dar_Diver4_ave','units','--', myThid)
268     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver4_ave','description',
269     & '', myThid)
270    
271     #endif
272     #ifdef DAR_DIAG_GROW
273     do np=1,npmax
274     ii=iPhy+np-1
275     DO jj = 1,MAX_LEN_FNAM
276     name(jj:jj) = ' '
277     ENDDO
278     ilnb = ILNBLNK(PTRACERS_names(ii))
279     WRITE(name,'(a,a)') 'dar_Grow_',PTRACERS_names(ii)(1:ilnb)
280     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
281     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
282     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
283     enddo
284     do np=1,npmax
285     ii=iPhy+np-1
286     DO jj = 1,MAX_LEN_FNAM
287     name(jj:jj) = ' '
288     ENDDO
289     ilnb = ILNBLNK(PTRACERS_names(ii))
290     WRITE(name,'(a,a)') 'dar_Growsq_',PTRACERS_names(ii)(1:ilnb)
291     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
292     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
293     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
294     enddo
295     #endif
296    
297     #ifdef ALLOW_DIAZ
298     #ifdef DAR_DIAG_NFIXP
299     do np=1,npmax
300     ii=iPhy+np-1
301     DO jj = 1,MAX_LEN_FNAM
302     name(jj:jj) = ' '
303     ENDDO
304     ilnb = ILNBLNK(PTRACERS_names(ii))
305     WRITE(name,'(a,a)') 'dar_NfixP_',PTRACERS_names(ii)(1:ilnb)
306     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
307     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
308     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
309     enddo
310     #endif
311     #endif
312    
313    
314    
315     c CALL MNC_CW_ADD_VNAME(
316     c & 'dic_pH_ave', 'Cen_xy_Hn__-__t', 3,4, myThid)
317     c CALL MNC_CW_ADD_VATTR_TEXT(
318     c & 'dic_pH_ave','units','--', myThid)
319     c CALL MNC_CW_ADD_VATTR_TEXT('dic_pH_ave','description',
320     c & '', myThid)
321    
322     ENDIF
323    
324     #endif /* ALLOW_MNC */
325     #endif /* ALLOW_DARWIN */
326    
327     RETURN
328     END
  ViewVC Help
Powered by ViewVC 1.1.22