File indexing completed on 2025-12-02 20:42:55 UTC
view on githubraw file Latest commit a092808e on 2025-12-02 20:09:45 UTC
8fbfd1f382 Oliv*0001 #include "DARWIN_OPTIONS.h"
0002
0003
0004
0005
0006
0007 SUBROUTINE DARWIN_DIAGNOSTICS_INIT( myThid )
0008
0009
0010
0011
0012
0013 IMPLICIT NONE
0014 #include "SIZE.h"
0015 #include "EEPARAMS.h"
0016 #include "PARAMS.h"
0017 #include "PTRACERS_SIZE.h"
0018 #include "PTRACERS_PARAMS.h"
0019 #ifdef ALLOW_RADTRANS
0020 #include "RADTRANS_SIZE.h"
0021 #endif
0022 #include "DARWIN_SIZE.h"
0023 #include "DARWIN_INDICES.h"
0024
0025
0026
0027 INTEGER myThid
0028
0029
0030 #ifdef ALLOW_DARWIN
0031 #ifdef ALLOW_DIAGNOSTICS
0032
0033
0034 INTEGER ILNBLNK
0035 EXTERNAL ILNBLNK
0036 CHARACTER*16 DIAGS_MK_UNITS
0037 EXTERNAL DIAGS_MK_UNITS
0038
0039
0040
0041
0042 INTEGER j
0043 INTEGER diagNum
0044 CHARACTER*8 diagName
0045 CHARACTER*16 diagCode
0046 CHARACTER*16 diagUnits
0047 CHARACTER*(80) diagTitle
0048 INTEGER ilnb
0049 CHARACTER*7 trcUnits
0050 CHARACTER*30 locName
0051 CHARACTER*4 diagSufx
0052
0053
0054
2834b1a102 Oliv*0055 WRITE(diagName,'(A)')'plankC'
865e47606a Oliv*0056 WRITE(diagTitle,'(A)')'Total plankton carbon biomass'
0057 diagUnits = 'mmol C/m^3 '
0058 diagCode = 'SMR MR '
0059 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0060 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0061
2834b1a102 Oliv*0062 WRITE(diagName,'(A)')'Chl'
865e47606a Oliv*0063 WRITE(diagTitle,'(A)')'Total Chlorophyll a'
0064 diagUnits = 'mg Chl a/m^3 '
0065 diagCode = 'SMR MR '
0066 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0067 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0068
8fbfd1f382 Oliv*0069 WRITE(diagName,'(A)')'PP'
0070 WRITE(diagTitle,'(A)')'Primary Production'
0071 diagUnits = 'mmol C/m^3/s '
0072 diagCode = 'SMRP MR '
0073 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0074 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0075
0076 WRITE(diagName,'(A)')'Nfix'
0077 WRITE(diagTitle,'(A)')'N fixation'
0078 diagUnits = 'mmol N/m^3/s '
0079 diagCode = 'SMRP MR '
0080 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0081 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0082
0083 WRITE(diagName,'(A)')'Denit'
0084 WRITE(diagTitle,'(A)')'Denitrification'
0085 diagUnits = 'mmol N/m^3/s '
0086 diagCode = 'SMRP MR '
0087 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0088 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
3179881a2b Oliv*0089
0090 WRITE(diagName,'(A)')'DenitN'
0091 WRITE(diagTitle,'(A)')'Nitrogen loss due to denitrification'
0092 diagUnits = 'mmol N/m^3/s '
0093 diagCode = 'SMRP MR '
0094 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0095 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
8fbfd1f382 Oliv*0096
c7b6c66d45 Oliv*0097 #ifdef DARWIN_ALLOW_CSTORE
0098 WRITE(diagName,'(A)')'EXU'
0099 WRITE(diagTitle,'(A)')'Carbon exudation'
0100 diagUnits = 'mmol C/m^3/s '
0101 diagCode = 'SM P MR '
0102 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0103 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0104
0105 WRITE(diagName,'(A)')'BioSyn'
0106 WRITE(diagTitle,'(A)')'Biosynthesis rate'
0107 diagUnits = 'mmol C/m^3/s '
0108 diagCode = 'SM P MR '
0109 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0110 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0111
0112 WRITE(diagName,'(A)')'DmdN'
0113 WRITE(diagTitle,'(A)')'Carbon demand from N limit'
0114 diagUnits = 'mmol C/m^3/s '
0115 diagCode = 'SM P MR '
0116 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0117 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0118
0119 WRITE(diagName,'(A)')'DmdP'
0120 WRITE(diagTitle,'(A)')'Carbon demand from P limit'
0121 diagUnits = 'mmol C/m^3/s '
0122 diagCode = 'SM P MR '
0123 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0124 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0125
0126 WRITE(diagName,'(A)')'DmdFe'
0127 WRITE(diagTitle,'(A)')'Carbon demand from Fe limit'
0128 diagUnits = 'mmol C/m^3/s '
0129 diagCode = 'SM P MR '
0130 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0131 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0132
0133 WRITE(diagName,'(A)')'DmdSi'
0134 WRITE(diagTitle,'(A)')'Carbon demand from Si limit'
0135 diagUnits = 'mmol C/m^3/s '
0136 diagCode = 'SM P MR '
0137 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0138 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0139
0140 WRITE(diagName,'(A)')'Dmdmin'
0141 WRITE(diagTitle,'(A)')'Minimum carbon demand'
0142 diagUnits = 'mmol C/m^3/s '
0143 diagCode = 'SM P MR '
0144 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0145 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0146 #endif
0147
8fbfd1f382 Oliv*0148 #ifdef ALLOW_RADTRANS
0149 DO j = 1, nlam
0150 WRITE(diagName,'(A,I3.3)')'PAR',j
0151 WRITE(diagTitle,'(A,I3)')'PAR waveband ',j
0152 diagUnits = 'uEin/m^2/s '
0153 diagCode = 'SMRP MR '
0154 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0155 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0156
0157 WRITE(diagName,'(A,I3.3)')'PARF',j
0158 WRITE(diagTitle,'(A,I3)')'PAR at W point, waveband ',j
0159 diagUnits = 'uEin/m^2/s '
55aff0bb70 Oliv*0160 diagCode = 'SM P LR '
8fbfd1f382 Oliv*0161 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0162 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0163
0164 #ifdef DARWIN_DIAG_IOP
0165 WRITE(diagName,'(A,I3.3)')'a',j
0166 WRITE(diagTitle,'(A,I3)')
0167 & 'total absorption for waveband ',j
0168 diagUnits = '1/m '
0169 diagCode = 'SMRP MR '
0170 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0171 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0172
0173 WRITE(diagName,'(A,I3.3)')'bt',j
0174 WRITE(diagTitle,'(A,I3)')
0175 & 'total scattering for waveband ',j
0176 diagUnits = '1/m '
0177 diagCode = 'SMRP MR '
0178 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0179 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0180
0181 WRITE(diagName,'(A,I3.3)')'bb',j
0182 WRITE(diagTitle,'(A,I3)')
0183 & 'total backscattering for waveband ',j
0184 diagUnits = '1/m '
0185 diagCode = 'SMRP MR '
0186 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0187 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0188
0189 WRITE(diagName,'(A,I3.3)')'aplk',j
0190 WRITE(diagTitle,'(A,I3)')
0191 & 'absorption by plankton for waveband ',j
0192 diagUnits = '1/m '
0193 diagCode = 'SMRP MR '
0194 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0195 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0196
0197 WRITE(diagName,'(A,I3.3)')'btplk',j
0198 WRITE(diagTitle,'(A,I3)')
0199 & 'scattering by plankton for waveband ',j
0200 diagUnits = '1/m '
0201 diagCode = 'SMRP MR '
0202 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0203 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0204
0205 WRITE(diagName,'(A,I3.3)')'bbplk',j
0206 WRITE(diagTitle,'(A,I3)')
0207 & 'backscattering by plankton for waveband ',j
0208 diagUnits = '1/m '
0209 diagCode = 'SMRP MR '
0210 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0211 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0212
0213 WRITE(diagName,'(A,I3.3)')'aprt',j
0214 WRITE(diagTitle,'(A,I3)')
0215 & 'absorption by particles for waveband ',j
0216 diagUnits = '1/m '
0217 diagCode = 'SMRP MR '
0218 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0219 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0220
0221 WRITE(diagName,'(A,I3.3)')'btprt',j
0222 WRITE(diagTitle,'(A,I3)')
0223 & 'scattering by particles for waveband ',j
0224 diagUnits = '1/m '
0225 diagCode = 'SMRP MR '
0226 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0227 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0228
0229 WRITE(diagName,'(A,I3.3)')'bbprt',j
0230 WRITE(diagTitle,'(A,I3)')
0231 & 'backscattering by particles for waveband ',j
0232 diagUnits = '1/m '
0233 diagCode = 'SMRP MR '
0234 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0235 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0236
0237 WRITE(diagName,'(A,I3.3)')'aCDOM',j
0238 WRITE(diagTitle,'(A,I3)')
0239 & 'absorption by CDOM for waveband ',j
0240 diagUnits = '1/m '
0241 diagCode = 'SMRP MR '
0242 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0243 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0244 #endif
0245
0246 ENDDO
0247 #endif /* ALLOW_RADTRANS */
0248
0249 WRITE(diagName,'(A)')'atten'
0250 WRITE(diagTitle,'(A)')'attenuation in layer'
0251 diagUnits = '1 '
0252 diagCode = 'SMRP MR '
0253 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0254 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0255
0256 WRITE(diagName,'(A)')'PARF'
0257 WRITE(diagTitle,'(A)')'PAR at top of layer'
0258 diagUnits = 'uEin/m^2/s '
55aff0bb70 Oliv*0259 diagCode = 'SM P LR '
8fbfd1f382 Oliv*0260 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0261 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0262
0263 WRITE(diagName,'(A)')'PAR'
0264 WRITE(diagTitle,'(A,I3)')'total PAR at layer center'
0265 diagUnits = 'uEin/m^2/s '
0266 diagCode = 'SMRP MR '
0267 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0268 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0269
0270
0271 WRITE(diagName,'(A)')'C_DIN'
0272 WRITE(diagTitle,'(A)')'consumption of DIN'
4891738544 Oliv*0273 diagUnits = 'mmol N/m3/s '
2834b1a102 Oliv*0274 diagCode = 'SMRP MR '
8fbfd1f382 Oliv*0275 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0276 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0277
b61dca872c Oliv*0278 WRITE(diagName,'(A)')'C_NO3'
0279 WRITE(diagTitle,'(A)')'consumption of NO3'
4891738544 Oliv*0280 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0281 diagCode = 'SMRP MR '
0282 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0283 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0284
0285 WRITE(diagName,'(A)')'C_NO2'
0286 WRITE(diagTitle,'(A)')'consumption of NO2'
4891738544 Oliv*0287 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0288 diagCode = 'SMRP MR '
0289 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0290 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0291
0292 WRITE(diagName,'(A)')'C_NH4'
0293 WRITE(diagTitle,'(A)')'consumption of NH4'
4891738544 Oliv*0294 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0295 diagCode = 'SMRP MR '
0296 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0297 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0298
8fbfd1f382 Oliv*0299 WRITE(diagName,'(A)')'C_PO4'
0300 WRITE(diagTitle,'(A)')'consumption of PO4'
4891738544 Oliv*0301 diagUnits = 'mmol P/m3/s '
2834b1a102 Oliv*0302 diagCode = 'SMRP MR '
8fbfd1f382 Oliv*0303 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0304 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0305
0306 WRITE(diagName,'(A)')'C_Si'
0307 WRITE(diagTitle,'(A)')'consumption of Si'
4891738544 Oliv*0308 diagUnits = 'mmol Si/m3/s '
2834b1a102 Oliv*0309 diagCode = 'SMRP MR '
8fbfd1f382 Oliv*0310 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0311 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0312
0313 WRITE(diagName,'(A)')'C_Fe'
0314 WRITE(diagTitle,'(A)')'consumption of Fe'
4891738544 Oliv*0315 diagUnits = 'mmol Fe/m3/s '
2834b1a102 Oliv*0316 diagCode = 'SMRP MR '
8fbfd1f382 Oliv*0317 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0318 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0319
0320
0321 WRITE(diagName,'(A)')'S_DIN'
0322 WRITE(diagTitle,'(A)')'non-transport source of DIN'
4891738544 Oliv*0323 diagUnits = 'mmol N/m3/s '
2834b1a102 Oliv*0324 diagCode = 'SMR MR '
8fbfd1f382 Oliv*0325 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0326 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0327
b61dca872c Oliv*0328 WRITE(diagName,'(A)')'S_NO3'
0329 WRITE(diagTitle,'(A)')'non-transport source of NO3'
4891738544 Oliv*0330 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0331 diagCode = 'SMR MR '
0332 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0333 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0334
0335 WRITE(diagName,'(A)')'S_NO2'
0336 WRITE(diagTitle,'(A)')'non-transport source of NO2'
4891738544 Oliv*0337 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0338 diagCode = 'SMR MR '
0339 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0340 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0341
0342 WRITE(diagName,'(A)')'S_NH4'
0343 WRITE(diagTitle,'(A)')'non-transport source of NH4'
4891738544 Oliv*0344 diagUnits = 'mmol N/m3/s '
b61dca872c Oliv*0345 diagCode = 'SMR MR '
0346 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0347 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0348
8fbfd1f382 Oliv*0349 WRITE(diagName,'(A)')'S_PO4'
0350 WRITE(diagTitle,'(A)')'non-transport source of PO4'
4891738544 Oliv*0351 diagUnits = 'mmol P/m3/s '
2834b1a102 Oliv*0352 diagCode = 'SMR MR '
8fbfd1f382 Oliv*0353 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0354 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0355
0356 WRITE(diagName,'(A)')'S_Si'
0357 WRITE(diagTitle,'(A)')'non-transport source of Si'
4891738544 Oliv*0358 diagUnits = 'mmol Si/m3/s '
2834b1a102 Oliv*0359 diagCode = 'SMR MR '
8fbfd1f382 Oliv*0360 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0361 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0362
0363 WRITE(diagName,'(A)')'S_Fe'
0364 WRITE(diagTitle,'(A)')'non-transport source of Fe'
4891738544 Oliv*0365 diagUnits = 'mmol Fe/m3/s '
2834b1a102 Oliv*0366 diagCode = 'SMR MR '
8fbfd1f382 Oliv*0367 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0368 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0369
0370 DO j = 1, nDarwin
0371 WRITE(diagName, '(2A)') 'gDAR', PTRACERS_ioLabel(j)
0372
0373 WRITE(locName,'(A,I4.4,A)') 'Tracer ', j
0374 trcUnits = 'kg/kg '
0375
0376 ilnb = ILNBLNK(PTRACERS_units(j))
0377 IF ( ilnb.GE.1 ) THEN
0378 ilnb = LEN(trcUnits)
0379 trcUnits = PTRACERS_units(j)(1:ilnb)
0380 ENDIF
0381 ilnb = ILNBLNK(PTRACERS_names(j))
0382 IF ( ilnb.GE.1 ) THEN
0383 ilnb = MIN(LEN(locName),ilnb)
0384 WRITE(locName,'(A)') PTRACERS_names(j)(1:ilnb)
0385 ENDIF
0386 ilnb = MAX(ILNBLNK(locName),1)
0387 WRITE(diagTitle,'(2A)') locName(1:ilnb), ' tendency from DARWIN'
0388 diagUnits = DIAGS_MK_UNITS( trcUnits//'/s', myThid )
0389 diagCode = 'SMR MR '
0390 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0391 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
5a97f77379 Oliv*0392
0393 WRITE(diagName, '(2A)') 'gECO', PTRACERS_ioLabel(j)
0394 ilnb = MAX(ILNBLNK(locName),1)
0395 WRITE(diagTitle,'(2A)') locName(1:ilnb),
0396 & ' tendency from DARWIN w/o sink/swim'
faa67d1773 Oliv*0397 diagUnits = DIAGS_MK_UNITS( trcUnits//'/s', myThid )
0398 diagCode = 'SMR MR '
5a97f77379 Oliv*0399 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0400 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
8fbfd1f382 Oliv*0401 ENDDO
0402
ee4178dfbe Oliv*0403 #ifdef DARWIN_DIAG_PERTYPE
0404 DO j = 1, nplank
8fbfd1f382 Oliv*0405 WRITE(diagName,'(A,I4.4)')'PP',j
0406 WRITE(diagTitle,'(A,I4)')'Primary Production plankton ',j
0407 diagUnits = 'mmol C/m^3/s '
0408 diagCode = 'SMRP MR '
0409 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0410 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0411
d5f8f1c735 Oliv*0412 WRITE(diagName,'(A,I4.4)')'PC',j
a93256cee4 Oliv*0413 WRITE(diagTitle,'(A,I4)')
0414 & 'Carbon-specific phototrophic growth rate plankton ',j
0415 diagUnits = '1/s '
0416 diagCode = 'SMRP MR '
0417 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0418 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0419
0420 WRITE(diagName,'(A,I4.4)')'HP',j
0421 WRITE(diagTitle,'(A,I4)')'Heterotrophic production plankton ',j
0422 diagUnits = 'mmol C/m^3/s '
0423 diagCode = 'SMRP MR '
0424 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0425 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0426
0427 WRITE(diagName,'(A,I4.4)')'HC',j
0428 WRITE(diagTitle,'(A,I4)')
0429 & 'Carbon-specific heterotrophic growth rate plankton ' ,j
d5f8f1c735 Oliv*0430 diagUnits = '1/s '
0431 diagCode = 'SMRP MR '
0432 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0433 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0434
8fbfd1f382 Oliv*0435 WRITE(diagName,'(A,I4.4)')'GR',j
0436 WRITE(diagTitle,'(A,I4)')'Grazing loss of plankton ',j
0437 diagUnits = 'mmol C/m^3/s '
0438 diagCode = 'SMRP MR '
0439 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0440 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0441
73aabcca42 Oliv*0442 WRITE(diagName,'(A,I4.4)')'GrGn',j
0443 WRITE(diagTitle,'(A,I4)')'Grazing gain of plankton ',j
55aff0bb70 Oliv*0444 diagUnits = 'mmol C/m^3/s '
a93256cee4 Oliv*0445 diagCode = 'SMRP MR '
0446 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0447 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0448
0449 WRITE(diagName,'(A,I4.4)')'GrGC',j
0450 WRITE(diagTitle,'(A,I4)')
0451 & 'Carbon-specific grazing rate plankton ',j
0452 diagUnits = '1/s '
73aabcca42 Oliv*0453 diagCode = 'SMRP MR '
0454 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0455 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
fd20fa810c Oliv*0456
0457 WRITE(diagName,'(A,I4.4)')'Mort',j
0458 WRITE(diagTitle,'(A,I4)')'Mortality loss plankton ',j
0459 diagUnits = 'mmol C/m^3/s '
0460 diagCode = 'SMRP MR '
0461 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0462 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0463
0464 WRITE(diagName,'(A,I4.4)')'Resp',j
0465 WRITE(diagTitle,'(A,I4)')'Respiration loss plankton ',j
0466 diagUnits = 'mmol C/m^3/s '
0467 diagCode = 'SMRP MR '
0468 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0469 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0470
0471 WRITE(diagName,'(A,I4.4)')'limN',j
0472 WRITE(diagTitle,'(A,I4)')'N growth limitation plankton ',j
0473 diagUnits = '[0/1] '
0474 diagCode = 'SMRP MR '
0475 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0476 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0477
0478 WRITE(diagName,'(A,I4.4)')'limP',j
0479 WRITE(diagTitle,'(A,I4)')'P growth limitation plankton ',j
0480 diagUnits = '[0/1] '
0481 diagCode = 'SMRP MR '
0482 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0483 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0484
0485 WRITE(diagName,'(A,I4.4)')'limF',j
0486 WRITE(diagTitle,'(A,I4)')'Fe growth limitation plankton ',j
0487 diagUnits = '[0/1] '
0488 diagCode = 'SMRP MR '
0489 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0490 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0491
0492 WRITE(diagName,'(A,I4.4)')'limS',j
0493 WRITE(diagTitle,'(A,I4)')'Si growth limitation plankton ',j
0494 diagUnits = '[0/1] '
0495 diagCode = 'SMRP MR '
0496 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0497 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0498
0499 WRITE(diagName,'(A,I4.4)')'fnut',j
0500 WRITE(diagTitle,'(A,I4)')'Nutrient growth factor plankton ',j
0501 diagUnits = '[0/1] '
0502 diagCode = 'SMRP MR '
0503 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0504 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0505
0506 WRITE(diagName,'(A,I4.4)')'fIph',j
0507 WRITE(diagTitle,'(A,I4)')'Light growth factor plankton ',j
0508 diagUnits = '[0/1] '
0509 diagCode = 'SMRP MR '
0510 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0511 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0512
0513 WRITE(diagName,'(A,I4.4)')'fTph',j
0514 WRITE(diagTitle,'(A,I4)')'Temperature growth factor plankton ',j
0515 diagUnits = '[0/1] '
0516 diagCode = 'SMRP MR '
0517 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0518 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
73aabcca42 Oliv*0519 ENDDO
a092808e6b shlo*0520
0521 # ifdef DARWIN_MACROMOLECULAR_GROWTH
0522 DO j = 1,nPhoto
0523 WRITE(diagName,'(A,I4.4)')'PChl', j
0524 WRITE(diagTitle,'(A,I4)')'Chl-spec carbohydrat fix rate plank ',j
0525 diagUnits = 'mol C/mol Chl/s '
0526 diagCode = 'SM P MR '
0527 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0528 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0529
0530 WRITE(diagName,'(A,I4.4)')'VN', j
0531 WRITE(diagTitle,'(A,I4)')
0532 & 'Carbon-spec nitrogen uptake rate plank ',j
0533 diagUnits = 'mmol N/mmol C/s '
0534 diagCode = 'SM P MR '
0535 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0536 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0537
0538 WRITE(diagName,'(A,I4.4)')'VP', j
0539 WRITE(diagTitle,'(A,I4)')
0540 & 'Carbon-spec phosphorus uptake rate plank ',j
0541 diagUnits = 'mmol N/mmol C/s '
0542 diagCode = 'SM P MR '
0543 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0544 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0545
0546 WRITE(diagName,'(A,I4.4)')'MODE', j
0547 WRITE(diagTitle,'(A,I4)')'Growth limitation enum plankton ',j
0548 diagUnits = 'dimensionless '
0549 diagCode = 'SM P MR '
0550 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0551 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0552
0553 WRITE(diagName,'(A,I4.4)')'Fe_C', j
0554 WRITE(diagTitle,'(A,I4)')'Fe/C plankton ',j
0555 diagUnits = 'mol Fe/mol C '
0556 diagCode = 'SM P MR '
0557 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0558 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0559
0560 WRITE(diagName,'(A,I4.4)')'exQc', j
0561 WRITE(diagTitle,'(A,I4)')'exQc plankton ',j
0562 diagUnits = 's-1 '
0563 diagCode = 'SM P MR '
0564 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0565 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0566
0567 WRITE(diagName,'(A,I4.4)')'CChl', j
0568 WRITE(diagTitle,'(A,I4)')'C quota in Chlorophyll plankton ',j
0569 diagUnits = 'molC molC-1 '
0570 diagCode = 'SM P MR '
0571 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0572 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0573
0574 WRITE(diagName,'(A,I4.4)')'NChl', j
0575 WRITE(diagTitle,'(A,I4)')'N quota in Chlorophyll plankton ',j
0576 diagUnits = 'molN molC-1 '
0577 diagCode = 'SM P MR '
0578 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0579 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0580
0581 WRITE(diagName,'(A,I4.4)')'NPho', j
0582 WRITE(diagTitle,'(A,I4)')'N quota related to photosynth plank ',j
0583 diagUnits = 'molN molC-1 '
0584 diagCode = 'SM P MR '
0585 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0586 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0587
0588 WRITE(diagName,'(A,I4.4)')'NSyn', j
0589 WRITE(diagTitle,'(A,I4)')'N quota related to biosynth plank ',j
0590 diagUnits = 'molN molC-1 '
0591 diagCode = 'SM P MR '
0592 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0593 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0594
0595 WRITE(diagName,'(A,I4.4)')'NPrn', j
0596 WRITE(diagTitle,'(A,I4)')'N quota in proteins plankton ',j
0597 diagUnits = 'molN molC-1 '
0598 diagCode = 'SM P MR '
0599 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0600 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0601
0602 WRITE(diagName,'(A,I4.4)')'NRNA', j
0603 WRITE(diagTitle,'(A,I4)')'N quota in RNA plankton ',j
0604 diagUnits = 'molN molC-1 '
0605 diagCode = 'SM P MR '
0606 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0607 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0608
0609 WRITE(diagName,'(A,I4.4)')'NDNA', j
0610 WRITE(diagTitle,'(A,I4)')'N quota in DNA plankton ',j
0611 diagUnits = 'molN molC-1 '
0612 diagCode = 'SM P MR '
0613 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0614 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0615
0616 WRITE(diagName,'(A,I4.4)')'NSTO', j
0617 WRITE(diagTitle,'(A,I4)')'N quota in QN_store plankton ',j
0618 diagUnits = 'molN molC-1 '
0619 diagCode = 'SM P MR '
0620 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0621 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0622
0623 WRITE(diagName,'(A,I4.4)')'NEXC', j
0624 WRITE(diagTitle,'(A,I4)')'N quota in QN_excess plankton ',j
0625 diagUnits = 'molN molC-1 '
0626 diagCode = 'SM P MR '
0627 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0628 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0629
0630 WRITE(diagName,'(A,I4.4)')'PRNA', j
0631 WRITE(diagTitle,'(A,I4)')'P quota in RNA plankton ',j
0632 diagUnits = 'molP molC-1 '
0633 diagCode = 'SM P MR '
0634 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0635 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0636
0637 WRITE(diagName,'(A,I4.4)')'PDNA', j
0638 WRITE(diagTitle,'(A,I4)')'P quota in DNA plankton ',j
0639 diagUnits = 'molP molC-1 '
0640 diagCode = 'SM P MR '
0641 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0642 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0643
0644 WRITE(diagName,'(A,I4.4)')'PTHY', j
0645 WRITE(diagTitle,'(A,I4)')'P quota in Thylakoid membrane plank ',j
0646 diagUnits = 'molP molC-1 '
0647 diagCode = 'SM P MR '
0648 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0649 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0650
0651 WRITE(diagName,'(A,I4.4)')'PCON', j
0652 WRITE(diagTitle,'(A,I4)')'P quota a constant part plankton ',j
0653 diagUnits = 'molP molC-1 '
0654 diagCode = 'SM P MR '
0655 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0656 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0657
0658 WRITE(diagName,'(A,I4.4)')'PSTO', j
0659 WRITE(diagTitle,'(A,I4)')'P quota in P storage plankton ',j
0660 diagUnits = 'molP molC-1 '
0661 diagCode = 'SM P MR '
0662 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0663 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0664
0665 WRITE(diagName,'(A,I4.4)')'PEXC', j
0666 WRITE(diagTitle,'(A,I4)')'P quota in QP_excess plankton ',j
0667 diagUnits = 'mol P molC-1 '
0668 diagCode = 'SM P MR '
0669 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0670 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0671
0672 WRITE(diagName,'(A,I4.4)')'FPHO', j
0673 WRITE(diagTitle,'(A,I4)')'Fe quota in Photosystem plankton ',j
0674 diagUnits = 'mol Fe molC-1 '
0675 diagCode = 'SM P MR '
0676 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0677 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0678
0679 WRITE(diagName,'(A,I4.4)')'FSTO', j
0680 WRITE(diagTitle,'(A,I4)')'Fe quota in QFe_store plankton ',j
0681 diagUnits = 'mol Fe molC-1 '
0682 diagCode = 'SM P MR '
0683 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0684 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0685
0686 WRITE(diagName,'(A,I4.4)')'FEXC', j
0687 WRITE(diagTitle,'(A,I4)')'Fe quota in QFe_excess plankton ',j
0688 diagUnits = 'mol Fe molC-1 '
0689 diagCode = 'SM P MR '
0690 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0691 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0692
0693 WRITE(diagName,'(A,I4.4)')'Y_RQ', j
0694 WRITE(diagTitle,'(A,I4)')'O2:C ratio for respiration plankton ',j
0695 diagUnits = 'mol O2 molC-1 '
0696 diagCode = 'SM P MR '
0697 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0698 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0699
0700 WRITE(diagName,'(A,I4.4)')'limC', j
0701 WRITE(diagTitle,'(A,I4)')'C growth limitation plankton ',j
0702 diagUnits = '[0/1] '
0703 diagCode = 'SM P MR '
0704 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0705 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0706
0707 WRITE(diagName,'(A,I4.4)')'limL', j
0708 WRITE(diagTitle,'(A,I4)')'Light growth limitation plankton ',j
0709 diagUnits = '[0/1] '
0710 diagCode = 'SM P MR '
0711 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0712 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0713 ENDDO
0714 # endif /* DARWIN_MACROMOLECULAR_GROWTH */
0715 #endif /* DARWIN_DIAG_PERTYPE */
73aabcca42 Oliv*0716
c7b6c66d45 Oliv*0717 #ifdef DARWIN_ALLOW_CSTORE
0718 #ifdef DARWIN_ALLOW_CSTORE_DIAGS
0719 DO j = 1, nPhoto
0720 WRITE(diagName,'(A,I4.4)')'EXU',j
0721 WRITE(diagTitle,'(A,I4)')'Exudation plankton ',j
0722 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0723 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0724 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0725 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0726 ENDDO
0727
0728 DO j = 1, nPhoto
0729 WRITE(diagName,'(A,I4.4)')'BS',j
0730 WRITE(diagTitle,'(A,I4)')'Biosynthesis rate plankton ',j
0731 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0732 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0733 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0734 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0735 ENDDO
0736
0737 DO j = 1, nPhoto
0738 WRITE(diagName,'(A,I4.4)')'DN',j
0739 WRITE(diagTitle,'(A,I4)')
0740 & 'Carbon demand from N limit plankton ',j
0741 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0742 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0743 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0744 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0745 ENDDO
0746
0747 DO j = 1, nPhoto
0748 WRITE(diagName,'(A,I4.4)')'DP',j
0749 WRITE(diagTitle,'(A,I4)')
0750 & 'Carbon demand from P limit plankton ',j
0751 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0752 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0753 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0754 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0755 ENDDO
0756
0757 DO j = 1, nPhoto
0758 WRITE(diagName,'(A,I4.4)')'DFe',j
0759 WRITE(diagTitle,'(A,I4)')
0760 & 'Carbon demand from Fe limit plankton ',j
0761 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0762 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0763 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0764 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0765 ENDDO
0766
0767 DO j = 1, nPhoto
0768 WRITE(diagName,'(A,I4.4)')'DSi',j
0769 WRITE(diagTitle,'(A,I4)')
0770 & 'Carbon demand from Si limit plankton ',j
0771 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0772 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0773 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0774 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0775 ENDDO
0776
0777 DO j = 1, nPhoto
0778 WRITE(diagName,'(A,I4.4)')'Dmin',j
0779 WRITE(diagTitle,'(A,I4)')'Minimum carbon demand plankton ',j
0780 diagUnits = 'mmol C/m^3/s '
55aff0bb70 Oliv*0781 diagCode = 'SMRP MR '
c7b6c66d45 Oliv*0782 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0783 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0784 ENDDO
0785 #endif
0786 #endif
0787
3997029bbb Oliv*0788 WRITE(diagName,'(A)')'sfcSolFe'
0789 WRITE(diagTitle,'(A)')'Soluble iron input at sea surface'
0790 diagUnits = 'mmol Fe/m2/s '
0791 diagCode = 'SM P L1 '
0792 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0793 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0794
fd32b5ab01 Oliv*0795 WRITE(diagName,'(A)')'scvLosFe'
0796 WRITE(diagTitle,'(A)')'Iron loss from scavenging'
0797 diagUnits = 'mmol Fe/m3/s '
0798 diagCode = 'SMRP MR '
0799 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0800 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0801
82fd1e0266 Oliv*0802 WRITE(diagName,'(A)')'scavRate'
0803 WRITE(diagTitle,'(A)')'Iron scavenging rate'
0804 diagUnits = '1/s '
0805 diagCode = 'SMRP MR '
0806 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0807 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0808
b61dca872c Oliv*0809 WRITE(diagName,'(A)')'sedFe'
0810 WRITE(diagTitle,'(A)')'Iron input from sediment'
a8f4b41a08 Oliv*0811 diagUnits = 'mmol Fe/m3/s '
0812 diagCode = 'SMRP MR '
63ffa46034 Oliv*0813 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0814 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0815
0816 WRITE(diagName,'(A)')'freeFeLs'
0817 WRITE(diagTitle,'(A)')'Iron loss due to free iron limit'
0818 diagUnits = 'mmol Fe/m3/s '
0819 diagCode = 'SMRP MR '
b61dca872c Oliv*0820 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0821 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0822
c8be5f4655 Oliv*0823 WRITE(diagName,'(A)')'sedFlxFe'
0824 WRITE(diagTitle,'(A)')'Iron sediment flux'
0825 diagUnits = 'mmol Fe/m2/s '
0826 diagCode = 'SM P U1 '
0827 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0828 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0829
8fbfd1f382 Oliv*0830 #ifdef DARWIN_ALLOW_CARBON
90ddb10039 Oliv*0831 WRITE(diagName,'(A)')'apCO2'
0832 WRITE(diagTitle,'(A)')'atmospheric pCO2'
0833 diagUnits = 'atm '
0834 diagCode = 'SM P L1 '
0835 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0836 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0837
0838 WRITE(diagName,'(A)')'apCO2sat'
0839 WRITE(diagTitle,'(A)')'atmospheric pCO2 sat'
0840 diagUnits = 'atm '
0841 diagCode = 'SM P L1 '
0842 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0843 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0844
8fbfd1f382 Oliv*0845 WRITE(diagName,'(A)')'pH'
0846 WRITE(diagTitle,'(A)')'pH'
0847 diagUnits = '1 '
0848 diagCode = 'SMR MR '
0849 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0850 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0851
0852 WRITE(diagName,'(A)')'pCO2'
90ddb10039 Oliv*0853 WRITE(diagTitle,'(A)')'Partial Pressure of CO2 in the ocean'
8fbfd1f382 Oliv*0854 diagUnits = 'atm '
0855 diagCode = 'SMRP MR '
0856 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0857 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0858
2e7c9eb2a6 Oliv*0859 WRITE(diagName,'(A)')'CO3'
0860 WRITE(diagTitle,'(A)')'Concentration of carbonate'
0861 diagUnits = 'mol/kg '
0862 diagCode = 'SMR MR '
0863 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0864 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0865
0866 WRITE(diagName,'(A)')'KspTPClc'
0867 WRITE(diagTitle,'(A)')'Solubility product for calcite'
0868 diagUnits = 'mol^2/kg^2 '
0869 diagCode = 'SMR MR '
0870 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0871 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0872
0873 WRITE(diagName,'(A)')'OmegaC'
0874 WRITE(diagTitle,'(A)')'Calcite saturation ratio'
0875 diagUnits = '1 '
0876 diagCode = 'SMR MR '
0877 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0878 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0879
875e37ed0a Oliv*0880 WRITE(diagName,'(A)')'fugfCO2'
cdf461a04e Oliv*0881 WRITE(diagTitle,'(A)')'Fugacity factor of CO2 at surface'
0882 diagUnits = '1 '
875e37ed0a Oliv*0883 diagCode = 'SM P L1 '
0884 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0885 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0886
a97de47f46 Oliv*0887 WRITE(diagName,'(A)')'fCO2'
0888 WRITE(diagTitle,'(A)')'Fugacity of CO2 (atm)'
0889 diagUnits = 'atm '
0890 diagCode = 'SM P L1 '
0891 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0892 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0893
8fbfd1f382 Oliv*0894 WRITE(diagName,'(A)')'fluxCO2'
0895 WRITE(diagTitle,'(A)')'Flux of CO2 - air-sea exch'
0896 diagUnits = 'mmol C/m2/s '
0897 diagCode = 'SM L1 '
0898 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0899 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0900
1b73f5d93e Oliv*0901 WRITE(diagName,'(A)')'fluxO2'
0902 WRITE(diagTitle,'(A)')'Flux of O2 - air-sea exch'
0903 diagUnits = 'mmol O2/m2/s '
0904 diagCode = 'SM L1 '
0905 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0906 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0907
0908 WRITE(diagName,'(A)')'VfluxCO2'
0909 WRITE(diagTitle,'(A)')'Virtual flux of CO2'
0910 diagUnits = 'mmol C/m2/s '
0911 diagCode = 'SM L1 '
0912 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0913 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0914
0915 WRITE(diagName,'(A)')'VfluxAlk'
0916 WRITE(diagTitle,'(A)')'Virtual flux of alkalinity'
e8cfd117fd Oliv*0917 diagUnits = 'meq/m2/s '
1b73f5d93e Oliv*0918 diagCode = 'SM L1 '
0919 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0920 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0921
8fbfd1f382 Oliv*0922 WRITE(diagName,'(A)')'gDICsurf'
0923 WRITE(diagTitle,'(A)')
0924 & 'Tendency of DIC due to air-sea exch + oldvirt.flux'
0925 diagUnits = 'mmol C/m3/s '
0926 diagCode = 'SM L1 '
0927 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0928 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0929
0930 WRITE(diagName,'(A)')'gALKsurf'
0931 WRITE(diagTitle,'(A)')'Tendency of ALK due to oldvirt.flux'
e8cfd117fd Oliv*0932 diagUnits = 'meq/m3/s '
8fbfd1f382 Oliv*0933 diagCode = 'SM L1 '
0934 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0935 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0936
0937 WRITE(diagName,'(A)')'gO2surf'
0938 WRITE(diagTitle,'(A)')'Tendency of O2 due to air-sea exch'
623d207dd8 Oliv*0939 diagUnits = 'mmol O2/m3/s '
8fbfd1f382 Oliv*0940 diagCode = 'SM L1 '
0941 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0942 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
b61dca872c Oliv*0943
0944 WRITE(diagName,'(A)')'C_DIC'
0945 WRITE(diagTitle,'(A)')'Consumption of DIC'
4891738544 Oliv*0946 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0947 diagCode = 'SMR MR '
0948 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0949 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0950
0951 WRITE(diagName,'(A)')'C_DICPIC'
0952 WRITE(diagTitle,'(A)')'Consumption of DIC due to PIC'
4891738544 Oliv*0953 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0954 diagCode = 'SMR MR '
0955 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0956 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0957
0958 WRITE(diagName,'(A)')'respDIC'
0959 WRITE(diagTitle,'(A)')'DIC due to respiration'
4891738544 Oliv*0960 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0961 diagCode = 'SMR MR '
0962 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0963 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0964
0965 WRITE(diagName,'(A)')'rDIC_DOC'
0966 WRITE(diagTitle,'(A)')'DIC due to remineralization of DOC'
4891738544 Oliv*0967 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0968 diagCode = 'SMR MR '
0969 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0970 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0971
0972 WRITE(diagName,'(A)')'rDIC_POC'
0973 WRITE(diagTitle,'(A)')'DIC due to remineralization of POC'
4891738544 Oliv*0974 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0975 diagCode = 'SMR MR '
0976 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0977 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0978
0979 WRITE(diagName,'(A)')'dDIC_PIC'
0980 WRITE(diagTitle,'(A)')'DIC due to dissolution of PIC'
4891738544 Oliv*0981 diagUnits = 'mmol C/m3/s '
b61dca872c Oliv*0982 diagCode = 'SMR MR '
0983 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0984 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0985
0986 WRITE(diagName,'(A)')'C_ALK'
0987 WRITE(diagTitle,'(A)')'Consumption of alkalinity'
4891738544 Oliv*0988 diagUnits = 'meq/m3/s '
b61dca872c Oliv*0989 diagCode = 'SMR MR '
0990 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0991 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0992
0993 WRITE(diagName,'(A)')'S_ALK'
0994 WRITE(diagTitle,'(A)')'Non-transport source of alkalinity'
4891738544 Oliv*0995 diagUnits = 'meq/m3/s '
b61dca872c Oliv*0996 diagCode = 'SMR MR '
0997 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0998 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0999
1000 WRITE(diagName,'(A)')'C_O2'
1001 WRITE(diagTitle,'(A)')'Consumption of oxygen'
4891738544 Oliv*1002 diagUnits = 'mmol O2/m3/s '
b61dca872c Oliv*1003 diagCode = 'SMR MR '
1004 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1005 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1006
1007 WRITE(diagName,'(A)')'S_O2'
1008 WRITE(diagTitle,'(A)')'Non-transport source of oxygen'
4891738544 Oliv*1009 diagUnits = 'mmol O2/m3/s '
b61dca872c Oliv*1010 diagCode = 'SMR MR '
1011 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1012 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1013
8a10163480 Oliv*1014 WRITE(diagName,'(A)')'OmegCbot'
1015 WRITE(diagTitle,'(A)')'Calcite saturation ratio above the seabed'
1016 diagUnits = '1 '
1017 diagCode = 'SM P M1 '
1018 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1019 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1020
b6fafcd98a Oliv*1021 WRITE(diagName,'(A)')'fIph_avg'
1022 WRITE(diagTitle,'(A)')'Community-averaged light growth factor'
1023 diagUnits = '1 '
1024 diagCode = 'SMR MR '
1025 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1026 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1027
1028 WRITE(diagName,'(A)')'fTph_avg'
1029 WRITE(diagTitle,'(A)')'Community-averaged temp growth factor'
1030 diagUnits = '1 '
1031 diagCode = 'SMR MR '
1032 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1033 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1034
1035 WRITE(diagName,'(A)')'fnut_avg'
1036 WRITE(diagTitle,'(A)')'Community-averaged nutrients growth factor'
1037 diagUnits = '1 '
1038 diagCode = 'SMR MR '
1039 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1040 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1041
8a10163480 Oliv*1042 # ifdef DARWIN_ALLOW_RADIv1
6aa5674e10 Oliv*1043 WRITE(diagName,'(A)')'DICFsed '
1044 WRITE(diagTitle,'(A)')'DIC sediment flux'
1045 diagUnits = 'mmol C/m2/s '
1046 diagCode = 'SM U1 '
1047 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1048 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1049
1050 WRITE(diagName,'(A)')'ALKFsed '
1051 WRITE(diagTitle,'(A)')'ALK sediment flux'
1052 diagUnits = 'meq/m2/s '
1053 diagCode = 'SM U1 '
1054 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1055 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1056
1057 WRITE(diagName,'(A)')'O2Fsed '
1058 WRITE(diagTitle,'(A)')'O2 sediment flux'
1059 diagUnits = 'mmol O2/m2/s '
1060 diagCode = 'SM U1 '
1061 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1062 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1063
1064 WRITE(diagName,'(A)')'POCFbur '
1065 WRITE(diagTitle,'(A)')'POC buried flux'
1066 diagUnits = 'mmol C/m2/s '
1067 diagCode = 'SM U1 '
1068 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1069 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1070
1071 WRITE(diagName,'(A)')'CALFbur '
1072 WRITE(diagTitle,'(A)')'Calcite buried flux'
4891738544 Oliv*1073 diagUnits = 'mmol C/m2/s '
6aa5674e10 Oliv*1074 diagCode = 'SM U1 '
1075 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1076 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1077
1078 WRITE(diagName,'(A)')'sedFlPOC'
1079 WRITE(diagTitle,'(A)')'POC flux from water column to sediment'
1080 diagUnits = 'mmol C/m2/s '
1081 diagCode = 'SM P U1 '
1082 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1083 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1084
1085 WRITE(diagName,'(A)')'sedFlPIC'
1086 WRITE(diagTitle,'(A)')'PIC flux from water column to sediment'
1087 diagUnits = 'mmol C/m2/s '
1088 diagCode = 'SM P U1 '
1089 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1090 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1091
1092 WRITE(diagName,'(A)')'TcorrO2 '
1093 WRITE(diagTitle,'(A)')
1094 & 'Temperature correction factor for diffusive O2 flux'
1095 diagUnits = '1 '
1096 diagCode = 'SM P M1 '
1097 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1098 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1099
1100 WRITE(diagName,'(A)')'TcorrALK'
1101 WRITE(diagTitle,'(A)')
1102 & 'Temperature correction factor for diffusive ALK flux'
1103 diagUnits = '1 '
1104 diagCode = 'SM P M1 '
1105 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1106 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1107
1108 WRITE(diagName,'(A)')'TcorrDIC'
1109 WRITE(diagTitle,'(A)')
1110 & 'Temperature correction factor for diffusive DIC flux'
1111 diagUnits = '1 '
1112 diagCode = 'SM P M1 '
1113 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1114 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1115
1116 WRITE(diagName,'(A)')'OmegCbot'
1117 WRITE(diagTitle,'(A)')'Calcite saturation ratio above the seabed'
1118 diagUnits = '1 '
1119 diagCode = 'SM P M1 '
1120 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1121 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
8a10163480 Oliv*1122 # endif /* DARWIN_ALLOW_RADIv1 */
1123
1124 # ifdef DARWIN_ALLOW_RADIv2
1125 WRITE(diagName,'(A)')'DICFsed '
1126 WRITE(diagTitle,'(A)')'DIC sediment flux'
1127 diagUnits = 'mmol C/m2/s '
1128 diagCode = 'SM U1 '
1129 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1130 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1131
1132 WRITE(diagName,'(A)')'ALKFsed '
1133 WRITE(diagTitle,'(A)')'ALK sediment flux'
1134 diagUnits = 'meq/m2/s '
1135 diagCode = 'SM U1 '
1136 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1137 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1138
1139 WRITE(diagName,'(A)')'O2Fsed '
1140 WRITE(diagTitle,'(A)')'O2 sediment flux'
1141 diagUnits = 'mmol O2/m2/s '
1142 diagCode = 'SM U1 '
1143 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1144 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1145
1146 WRITE(diagName,'(A)')'NO3Fsed '
1147 WRITE(diagTitle,'(A)')'NO3 sediment flux'
1148 diagUnits = 'mmol NO3/m2/s '
1149 diagCode = 'SM U1 '
1150 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1151 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1152
1153 WRITE(diagName,'(A)')'PO4Fsed '
1154 WRITE(diagTitle,'(A)')'PO4 sediment flux'
1155 diagUnits = 'mmol PO4/m2/s '
1156 diagCode = 'SM U1 '
1157 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1158 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1159
1160 WRITE(diagName,'(A)')'NH4Fsed '
1161 WRITE(diagTitle,'(A)')'NH4 sediment flux'
1162 diagUnits = 'mmol NH4/m2/s '
1163 diagCode = 'SM U1 '
1164 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1165 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1166
1167 WRITE(diagName,'(A)')'sedFlPOC'
1168 WRITE(diagTitle,'(A)')'POC flux from water column to sediment'
1169 diagUnits = 'mmol C/m2/s '
1170 diagCode = 'SM P U1 '
1171 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1172 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1173
1174 WRITE(diagName,'(A)')'sedFlPIC'
1175 WRITE(diagTitle,'(A)')'PIC flux from water column to sediment'
1176 diagUnits = 'mmol C/m2/s '
1177 diagCode = 'SM P U1 '
1178 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1179 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1180
1181 WRITE(diagName,'(A)')'sedFlPOP'
1182 WRITE(diagTitle,'(A)')'POP flux from water column to sediment'
1183 diagUnits = 'mmol P/m2/s '
1184 diagCode = 'SM P U1 '
1185 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1186 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1187
1188 WRITE(diagName,'(A)')'sedFlPON'
1189 WRITE(diagTitle,'(A)')'PON flux from water column to sediment'
1190 diagUnits = 'mmol N/m2/s '
1191 diagCode = 'SM P U1 '
1192 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1193 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1194
1195 WRITE(diagName,'(A)')'sedFplkC'
1196 WRITE(diagTitle,'(A)')'C flux from bottom sink plankton to sed'
1197 diagUnits = 'mmol C/m2/s '
1198 diagCode = 'SM P U1 '
1199 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1200 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1201
1202 WRITE(diagName,'(A)')'sedFplkN'
1203 WRITE(diagTitle,'(A)')'N flux from bottom sink plankton to sed'
1204 diagUnits = 'mmol N/m2/s '
1205 diagCode = 'SM P U1 '
1206 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1207 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1208
1209 WRITE(diagName,'(A)')'sedFplkP'
1210 WRITE(diagTitle,'(A)')'P flux from bottom sink plankton to sed'
1211 diagUnits = 'mmol P/m2/s '
1212 diagCode = 'SM P U1 '
1213 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1214 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1215
1216 # endif /* DARWIN_ALLOW_RADIv2 */
8fbfd1f382 Oliv*1217 #endif /* DARWIN_ALLOW_CARBON */
1218
90ddb10039 Oliv*1219 #ifdef DARWIN_DIAG_TENDENCIES
1220 WRITE(diagName,'(A)')'gDICEpr'
1221 WRITE(diagTitle,'(A)')
1222 & 'Tendency of DIC due to E/P/runoff'
1223 diagUnits = 'mmol C/m3/s '
1224 diagCode = 'SM L1 '
1225 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1226 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1227
1228 WRITE(diagName,'(A)')'gNO3Epr'
1229 WRITE(diagTitle,'(A)')
1230 & 'Tendency of DIC due to E/P/runoff'
1231 diagUnits = 'mmol N/m3/s '
1232 diagCode = 'SM L1 '
1233 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1234 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1235
1236 WRITE(diagName,'(A)')'gNO2Epr'
1237 WRITE(diagTitle,'(A)')
1238 & 'Tendency of DIC due to E/P/runoff'
1239 diagUnits = 'mmol N/m3/s '
1240 diagCode = 'SM L1 '
1241 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1242 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1243
1244 WRITE(diagName,'(A)')'gNH4Epr'
1245 WRITE(diagTitle,'(A)')
1246 & 'Tendency of DIC due to E/P/runoff'
1247 diagUnits = 'mmol N/m3/s '
1248 diagCode = 'SM L1 '
1249 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1250 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1251
1252 WRITE(diagName,'(A)')'gPO4Epr'
1253 WRITE(diagTitle,'(A)')
1254 & 'Tendency of PO4 due to E/P/runoff'
1255 diagUnits = 'mmol P/m3/s '
1256 diagCode = 'SM L1 '
1257 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1258 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1259
1260 WRITE(diagName,'(A)')'gFeTEpr'
1261 WRITE(diagTitle,'(A)')
1262 & 'Tendency of FeT due to E/P/runoff'
1263 diagUnits = 'mmol Fe/m3/s '
1264 diagCode = 'SM L1 '
1265 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1266 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1267
1268 WRITE(diagName,'(A)')'gSiO2Epr'
1269 WRITE(diagTitle,'(A)')
1270 & 'Tendency of SiO2 due to E/P/runoff'
1271 diagUnits = 'mmol Si/m3/s '
1272 diagCode = 'SM L1 '
1273 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1274 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1275
1276 # ifdef DARWIN_ALLOW_CARBON
1277 WRITE(diagName,'(A)')'gALKEpr'
1278 WRITE(diagTitle,'(A)')
1279 & 'Tendency of ALK due to E/P/runoff'
1280 diagUnits = 'meq/m3/s '
1281 diagCode = 'SM L1 '
1282 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1283 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1284
1285 WRITE(diagName,'(A)')'gO2Epr'
1286 WRITE(diagTitle,'(A)')
1287 & 'Tendency of O2 due to E/P/runoff'
1288 diagUnits = 'mmol O2/m3/s '
1289 diagCode = 'SM L1 '
1290 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1291 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1292 # endif
1293 #endif /* DARWIN_DIAG_TENDENCIES */
1294
8fbfd1f382 Oliv*1295 #ifdef ALLOW_EXF
1296 WRITE(diagName,'(A)')'surfPAR'
1297 WRITE(diagTitle,'(A)')'PAR forcing at surface'
1298 diagUnits = 'uEin/m^2/s '
1299 diagCode = 'SM P L1 '
1300 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1301 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1302
1303 WRITE(diagName,'(A)')'surfiron'
1304 WRITE(diagTitle,'(A)')'iron forcing at surface'
1305 diagUnits = 'mmol Fe/m^2/s '
1306 diagCode = 'SM P L1 '
1307 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1308 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1309
1310 WRITE(diagName,'(A)')'DARice'
1311 WRITE(diagTitle,'(A)')'ice area fraction'
1312 diagUnits = 'm^2/m^2 '
1313 diagCode = 'SM P L1 '
1314 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1315 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1316
1317 WRITE(diagName,'(A)')'DARwind'
1318 WRITE(diagTitle,'(A)')'wind speed used for carbon exchange'
1319 diagUnits = 'm/s '
1320 diagCode = 'SM P L1 '
1321 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1322 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1323
1324 WRITE(diagName,'(A)')'surfpCO2'
2834b1a102 Oliv*1325 WRITE(diagTitle,'(A)')'atmospheric surface pCO2'
8fbfd1f382 Oliv*1326 diagUnits = 'mol/mol '
1327 diagCode = 'SM P L1 '
1328 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
1329 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
1330 #endif
1331
1332 #endif /* ALLOW_DIAGNOSTICS */
1333 #endif /* ALLOW_DARWIN */
1334
1335 RETURN
1336 END
1337
