/[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.7 - (hide annotations) (download)
Tue May 19 14:30:15 2015 UTC (10 years, 2 months ago) by benw
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.6: +46 -12 lines
Ben Ward - some superficial structural changes allowing runs with no pfts
         - more significant structural and parameter changes to follow later

1 benw 1.7 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_mnc_init.F,v 1.6 2012/10/23 17:55:30 stephd 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 benw 1.7 c carbon
67 benw 1.2 DO jj = 1,MAX_LEN_FNAM
68 benw 1.7 name(jj:jj) = ' '
69 benw 1.2 ENDDO
70 benw 1.7 WRITE(name,'(a,i3.3)') 'dar_AP_C_',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     WRITE(name,'(a,i3.3)') 'dar_HP_C_',np
75     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
76     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
77     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
78     c nitrogen
79     WRITE(name,'(a,i3.3)') 'dar_AP_N_',np
80     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
81     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
82     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
83     WRITE(name,'(a,i3.3)') 'dar_HP_N_',np
84     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
85     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
86     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
87 benw 1.2 WRITE(name,'(a,i3.3)') 'dar_Nlim_',np
88     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
89     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
90     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
91 benw 1.7 c phosphorus
92     #ifdef PQUOTA
93     WRITE(name,'(a,i3.3)') 'dar_AP_P_',np
94     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
95     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
96     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
97     WRITE(name,'(a,i3.3)') 'dar_HP_P_',np
98     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
99     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
100     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
101     WRITE(name,'(a,i3.3)') 'dar_Plim_',np
102     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
103     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
104     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
105     #endif
106 benw 1.2 c iron
107 benw 1.7 #ifdef FQUOTA
108     WRITE(name,'(a,i3.3)') 'dar_AP_F_',np
109     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
110     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
111     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
112     WRITE(name,'(a,i3.3)') 'dar_HP_F_',np
113     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
114     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
115     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
116 benw 1.2 WRITE(name,'(a,i3.3)') 'dar_Flim_',np
117     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
118     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
119     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
120 benw 1.7 #endif
121 benw 1.2 c light
122     WRITE(name,'(a,i3.3)') 'dar_Ilim_',np
123     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
124     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
125     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
126 benw 1.7 !
127 benw 1.2 enddo
128     c temperature
129     WRITE(name,'(a)') 'dar_Tlim_ave'
130     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
131     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
132     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
133     #endif
134     cbenw limitation diagnostics
135 jahn 1.1 #ifndef ALLOW_QUOTA
136     CALL MNC_CW_ADD_VNAME(
137     & 'dar_Chl_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
138     CALL MNC_CW_ADD_VATTR_TEXT(
139     & 'dar_Chl_ave','units','--', myThid)
140     CALL MNC_CW_ADD_VATTR_TEXT('dar_Chl_ave','description',
141     & '', myThid)
142     #endif
143     CALL MNC_CW_ADD_VNAME(
144     & 'dar_Nfix_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
145     CALL MNC_CW_ADD_VATTR_TEXT(
146     & 'dar_Nfix_ave','units','--', myThid)
147     CALL MNC_CW_ADD_VATTR_TEXT('dar_Nfix_ave','description',
148     & '', myThid)
149    
150     CALL MNC_CW_ADD_VNAME(
151     & 'dar_Denit_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
152     CALL MNC_CW_ADD_VATTR_TEXT(
153     & 'dar_Denit_ave','units','--', myThid)
154     CALL MNC_CW_ADD_VATTR_TEXT('dar_Denit_ave','description',
155     & '', myThid)
156    
157     #ifdef DAR_DIAG_ACDOM
158     CALL MNC_CW_ADD_VNAME(
159     & 'dar_aCDOM_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
160     CALL MNC_CW_ADD_VATTR_TEXT(
161     & 'dar_aCDOM_ave','units','--', myThid)
162     CALL MNC_CW_ADD_VATTR_TEXT('dar_aCDOM_ave','description',
163     & '', myThid)
164     #endif
165    
166 stephd 1.6 #ifdef DAR_DIAG_PARW
167     do ii=1,tlam
168     DO jj = 1,MAX_LEN_FNAM
169     name(jj:jj) = ' '
170     ENDDO
171     WRITE(name,'(a,i2.2)') 'dar_PARw',ii
172     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
173     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','uEin/m2/s', myThid)
174     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
175     & 'par at centre grid cell', myThid)
176     enddo
177     do np=1,npmax
178     DO jj = 1,MAX_LEN_FNAM
179     name(jj:jj) = ' '
180     ENDDO
181     WRITE(name,'(a,i2.2)') 'dar_chl2c',np
182     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
183     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','mg/mmol', myThid)
184     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
185     & 'local chl to carbon ratio', myThid)
186     enddo
187     #endif
188    
189 jahn 1.1 #ifdef DAR_DIAG_IRR
190     do ii=1,tlam
191     DO jj = 1,MAX_LEN_FNAM
192     name(jj:jj) = ' '
193     ENDDO
194     WRITE(name,'(a,i2.2)') 'dar_Ed',ii
195     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
196 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
197     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
198     & 'direct downwelling irradiance at top of layer', myThid)
199 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_Es',ii
200     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
201 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
202     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
203     & 'diffuse downwelling irradiance at top of layer', myThid)
204 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_Eu',ii
205     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
206 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
207     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
208     & 'upwelling irradiance just above top of layer', myThid)
209 jahn 1.5 WRITE(name,'(a,i2.2)') 'dar_Estop',ii
210     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
211     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
212     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
213     & 'diffuse downwelling irradiance just below top of layer',
214     & myThid)
215 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_Eutop',ii
216     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
217 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
218     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
219     & 'upwelling irradiance just below top of layer', myThid)
220 jahn 1.1 enddo
221     #endif
222    
223 jahn 1.3 #ifdef DAR_DIAG_IRR_AMPS
224     do ii=1,tlam
225     DO jj = 1,MAX_LEN_FNAM
226     name(jj:jj) = ' '
227     ENDDO
228 jahn 1.5 WRITE(name,'(a,i2.2)') 'dar_amp1_',ii
229 jahn 1.3 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
230 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
231     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
232     & 'amplitude of downward increasing mode', myThid)
233 jahn 1.5 WRITE(name,'(a,i2.2)') 'dar_amp2_',ii
234 jahn 1.3 CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
235 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','W/m^2', myThid)
236     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
237     & 'amplitude of downward decreasing mode', myThid)
238 jahn 1.3 enddo
239     #endif
240 jahn 1.1 #ifdef DAR_DIAG_ABSORP
241     do ii=1,tlam
242     DO jj = 1,MAX_LEN_FNAM
243     name(jj:jj) = ' '
244     ENDDO
245     WRITE(name,'(a,i2.2)') 'dar_a',ii
246     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
247 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
248     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
249     & 'total absorption coefficient', myThid)
250 jahn 1.1 enddo
251     #endif
252     #ifdef DAR_DIAG_SCATTER
253     do ii=1,tlam
254     DO jj = 1,MAX_LEN_FNAM
255     name(jj:jj) = ' '
256     ENDDO
257     WRITE(name,'(a,i2.2)') 'dar_bt',ii
258     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
259 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
260     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
261     & 'total scattering coefficient', myThid)
262 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_bb',ii
263     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
264 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
265     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
266     & 'total backscattering coefficient', myThid)
267 jahn 1.1 enddo
268     #endif
269     #ifdef DAR_DIAG_PART_SCATTER
270     do ii=1,tlam
271     DO jj = 1,MAX_LEN_FNAM
272     name(jj:jj) = ' '
273     ENDDO
274     WRITE(name,'(a,i2.2)') 'dar_apart',ii
275     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
276 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
277     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
278     & 'absorption coefficient due to particles', myThid)
279 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_btpart',ii
280     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
281 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
282     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
283     & 'scattering coefficient due to particles', myThid)
284 jahn 1.1 WRITE(name,'(a,i2.2)') 'dar_bbpart',ii
285     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
286 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'units','1/m', myThid)
287     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
288     & 'backscattering coefficient due to particles', myThid)
289 jahn 1.1 enddo
290     #endif
291 jahn 1.3 #ifdef DAR_RADTRANS
292     DO jj = 1,MAX_LEN_FNAM
293     name(jj:jj) = ' '
294     ENDDO
295     WRITE(name,'(a)') 'dar_rmud'
296     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__-__t',3,4,myThid)
297     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
298 jahn 1.4 CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
299     & 'inverse cosine of under-water solar zenith angle', myThid)
300 jahn 1.3 #endif
301 stephd 1.6 #ifdef DAR_DIAG_EK
302     do np=1,npmax
303     DO jj = 1,MAX_LEN_FNAM
304     name(jj:jj) = ' '
305     ENDDO
306     WRITE(name,'(a,i2.2)') 'dar_Ek',np
307     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
308     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','uEin/m2/s', myThid)
309     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
310     & 'Ek', myThid)
311     WRITE(name,'(a,i2.2)') 'dar_EkoverE',np
312     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
313     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
314     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
315     & 'Ek over E', myThid)
316     WRITE(name,'(a,i2.2)') 'dar_acclim',np
317     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
318     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','mg/mmol', myThid)
319     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
320     & 'acclimated Chl to C ratio', myThid)
321     do ii=1,tlam
322     DO jj = 1,MAX_LEN_FNAM
323     name(jj:jj) = ' '
324     ENDDO
325     WRITE(name,'(a,i2.2,a,i2.2)') 'dar_Ek_nl',np,'_',ii
326     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
327     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','uEin/m2/s', myThid)
328     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
329     & 'spectral Ek', myThid)
330     WRITE(name,'(a,i2.2,a,i2.2)') 'dar_EkoverE_nl',np,'_',ii
331     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
332     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
333     CALL MNC_CW_ADD_VATTR_TEXT(name,'description',
334     & 'spectral Ek over E', myThid)
335     enddo
336     enddo
337     #endif
338 jahn 1.1 c ANNA_TAVE
339     #ifdef WAVES_DIAG_PCHL
340     do np=1,npmax
341     ii=iPhy+np-1
342     DO jj = 1,MAX_LEN_FNAM
343     name(jj:jj) = ' '
344     ENDDO
345     ilnb = ILNBLNK(PTRACERS_names(ii))
346     WRITE(name,'(a,a)') 'dar_Pchl_',PTRACERS_names(ii)(1:ilnb)
347     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
348     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
349     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
350     enddo
351     #endif
352     c ANNA end TAVE
353     #ifdef DAR_DIAG_RSTAR
354     do np=1,npmax
355     ii=iPhy+np-1
356     DO jj = 1,MAX_LEN_FNAM
357     name(jj:jj) = ' '
358     ENDDO
359     ilnb = ILNBLNK(PTRACERS_names(ii))
360     WRITE(name,'(a,a)') 'dar_Rstar_',PTRACERS_names(ii)(1:ilnb)
361     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
362     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
363     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
364     enddo
365     do np=1,npmax
366     ii=iPhy+np-1
367     DO jj = 1,MAX_LEN_FNAM
368     name(jj:jj) = ' '
369     ENDDO
370     ilnb = ILNBLNK(PTRACERS_names(ii))
371     WRITE(name,'(a,a)') 'dar_RNstar_',PTRACERS_names(ii)(1:ilnb)
372     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
373     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
374     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
375     enddo
376     #endif
377    
378     #ifdef DAR_DIAG_DIVER
379     CALL MNC_CW_ADD_VNAME(
380     & 'dar_Diver1_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
381     CALL MNC_CW_ADD_VATTR_TEXT(
382     & 'dar_Diver1_ave','units','--', myThid)
383     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver1_ave','description',
384     & '', myThid)
385    
386     CALL MNC_CW_ADD_VNAME(
387     & 'dar_Diver2_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
388     CALL MNC_CW_ADD_VATTR_TEXT(
389     & 'dar_Diver2_ave','units','--', myThid)
390     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver2_ave','description',
391     & '', myThid)
392    
393     CALL MNC_CW_ADD_VNAME(
394     & 'dar_Diver3_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
395     CALL MNC_CW_ADD_VATTR_TEXT(
396     & 'dar_Diver3_ave','units','--', myThid)
397     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver3_ave','description',
398     & '', myThid)
399    
400     CALL MNC_CW_ADD_VNAME(
401     & 'dar_Diver4_ave', 'Cen_xy_Hn__C__t', 4,5, myThid)
402     CALL MNC_CW_ADD_VATTR_TEXT(
403     & 'dar_Diver4_ave','units','--', myThid)
404     CALL MNC_CW_ADD_VATTR_TEXT('dar_Diver4_ave','description',
405     & '', myThid)
406    
407     #endif
408     #ifdef DAR_DIAG_GROW
409     do np=1,npmax
410     ii=iPhy+np-1
411     DO jj = 1,MAX_LEN_FNAM
412     name(jj:jj) = ' '
413     ENDDO
414     ilnb = ILNBLNK(PTRACERS_names(ii))
415     WRITE(name,'(a,a)') 'dar_Grow_',PTRACERS_names(ii)(1:ilnb)
416     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
417     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
418     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
419     enddo
420     do np=1,npmax
421     ii=iPhy+np-1
422     DO jj = 1,MAX_LEN_FNAM
423     name(jj:jj) = ' '
424     ENDDO
425     ilnb = ILNBLNK(PTRACERS_names(ii))
426     WRITE(name,'(a,a)') 'dar_Growsq_',PTRACERS_names(ii)(1:ilnb)
427     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
428     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
429     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
430     enddo
431     #endif
432    
433     #ifdef ALLOW_DIAZ
434     #ifdef DAR_DIAG_NFIXP
435     do np=1,npmax
436     ii=iPhy+np-1
437     DO jj = 1,MAX_LEN_FNAM
438     name(jj:jj) = ' '
439     ENDDO
440     ilnb = ILNBLNK(PTRACERS_names(ii))
441     WRITE(name,'(a,a)') 'dar_NfixP_',PTRACERS_names(ii)(1:ilnb)
442     CALL MNC_CW_ADD_VNAME(name,'Cen_xy_Hn__C__t',4,5,myThid)
443     CALL MNC_CW_ADD_VATTR_TEXT(name,'units','--', myThid)
444     CALL MNC_CW_ADD_VATTR_TEXT(name,'description','', myThid)
445     enddo
446     #endif
447     #endif
448    
449    
450    
451     c CALL MNC_CW_ADD_VNAME(
452     c & 'dic_pH_ave', 'Cen_xy_Hn__-__t', 3,4, myThid)
453     c CALL MNC_CW_ADD_VATTR_TEXT(
454     c & 'dic_pH_ave','units','--', myThid)
455     c CALL MNC_CW_ADD_VATTR_TEXT('dic_pH_ave','description',
456     c & '', myThid)
457    
458     ENDIF
459    
460     #endif /* ALLOW_MNC */
461     #endif /* ALLOW_DARWIN */
462    
463     RETURN
464     END

  ViewVC Help
Powered by ViewVC 1.1.22