/[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.1 - (hide annotations) (download)
Wed Apr 13 18:56:24 2011 UTC (14 years, 3 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 jahn 1.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