pkg/bling/bling_diagnostics_init.F

e0f9a7ba0b Matt*0001 #include "BLING_OPTIONS.h"
c0d1c06c15 Matt*0002
                0003 CBOP
                0004       subroutine BLING_DIAGNOSTICS_INIT(myThid )
                0005
                0006 C     =================================================================
                0007 C     | subroutine bling_diagnostics_init
4ac06494d5 Matt*0008 C     | o Define diagnostics for bling pkg
c0d1c06c15 Matt*0009 C     =================================================================
                0010
e0f9a7ba0b Matt*0011       IMPLICIT NONE
                0012
c0d1c06c15 Matt*0013 C     === Global variables ===
                0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016
                0017 C  myThid               :: thread number
                0018       INTEGER myThid
                0019 CEOP
                0020
                0021 #ifdef ALLOW_DIAGNOSTICS
                0022
                0023 C     !LOCAL VARIABLES:
                0024       INTEGER       diagNum
                0025       CHARACTER*8   diagName
                0026       CHARACTER*16  diagCode
                0027       CHARACTER*16  diagUnits
                0028       CHARACTER*(80) diagTitle
                0029
                0030 C     Define diagnostics Names :
                0031
cab777d667 Matt*0032 c biology
c0d1c06c15 Matt*0033
                0034       diagName  = 'BLGPSM  '
                0035       diagTitle = 'Small phytoplankton biomass (mol C/m3/s)'
                0036       diagUnits = 'mol C/m3        '
                0037       diagCode  = 'SM P    MR      '
                0038       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0039      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0040
                0041       diagName  = 'BLGPLG  '
                0042       diagTitle = 'Large phytoplankton biomass (mol C/m3/s)'
                0043       diagUnits = 'mol C/m3        '
                0044       diagCode  = 'SM P    MR      '
                0045       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0046      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0047
e0f9a7ba0b Matt*0048 #ifndef USE_BLING_V1
c0d1c06c15 Matt*0049       diagName  = 'BLGPDIA '
                0050       diagTitle = 'Diazotroph phytoplankton biomass (mol C/m3/s)'
                0051       diagUnits = 'mol C/m3        '
                0052       diagCode  = 'SM P    MR      '
                0053       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0054      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e0f9a7ba0b Matt*0055 #endif
c0d1c06c15 Matt*0056
                0057       diagName  = 'BLGCHL  '
                0058       diagTitle = 'Chlorophyll concentration'
fe1e55ba7d Matt*0059       diagUnits = 'mg/m3           '
                0060       diagCode  = 'SM P    MR      '
                0061       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0062      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0063
82e538d851 aver*0064       diagName  = 'BLGCHLSA'
                0065       diagTitle = 'Chlorophyll concentration as seen by satellite'
                0066       diagUnits = 'mg/m3           '
                0067       diagCode  = 'SM P    M1      '
                0068       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0069      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0070
fe1e55ba7d Matt*0071       diagName  = 'BLGCHL2C'
                0072       diagTitle = 'Chlorophyll to Carbon ratio'
                0073       diagUnits = 'mol Chl/mol C   '
                0074       diagCode  = 'SM P    MR      '
                0075       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0076      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0077
                0078       diagName  = 'BLGC2CHL'
                0079       diagTitle = 'Carbon to Chlorophyll ratio'
                0080       diagUnits = 'g C/g Chl       '
c0d1c06c15 Matt*0081       diagCode  = 'SM P    MR      '
                0082       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0083      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0084
3b75ee2a42 Matt*0085       diagName  = 'BLGPOC  '
e0f9a7ba0b Matt*0086       diagTitle = 'Particulate Organic Carbon concentration'
fe1e55ba7d Matt*0087       diagUnits = 'mol C/m3        '
3b75ee2a42 Matt*0088       diagCode  = 'SM P    MR      '
                0089       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0090      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0091
c0d1c06c15 Matt*0092       diagName  = 'BLGPOCF '
                0093       diagTitle = 'Particulate Organic Carbon Flux'
                0094       diagUnits = 'mol C/m3/s      '
                0095       diagCode  = 'SM P    MR      '
                0096       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0097      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0098
                0099       diagName  = 'BLGNPP  '
                0100       diagTitle = 'Net Primary Production'
                0101       diagUnits = 'mol C/m3/s      '
                0102       diagCode  = 'SM P    MR      '
                0103       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0104      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0105
                0106       diagName  = 'BLGNCP  '
                0107       diagTitle = 'Net Community Production'
                0108       diagUnits = 'mol C/m3/s      '
                0109       diagCode  = 'SM P    MR      '
                0110       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0111      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0112
cab777d667 Matt*0113       diagName  = 'BLGPCM  '
                0114       diagTitle = 'Maximum photosynthesis rate'
e0f9a7ba0b Matt*0115       diagUnits = '1/s             '
cab777d667 Matt*0116       diagCode  = 'SM P    MR      '
                0117       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0118      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0119
                0120       diagName  = 'BLGMU   '
                0121       diagTitle = 'Phytoplankton growth rate'
e0f9a7ba0b Matt*0122       diagUnits = '1/s             '
cab777d667 Matt*0123       diagCode  = 'SM P    MR      '
                0124       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0125      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0126
e0f9a7ba0b Matt*0127 #ifndef USE_BLING_V1
cab777d667 Matt*0128       diagName  = 'BLGMUDIA'
                0129       diagTitle = 'Diazotroph growth rate'
e0f9a7ba0b Matt*0130       diagUnits = '1/s             '
cab777d667 Matt*0131       diagCode  = 'SM P    MR      '
                0132       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0133      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e0f9a7ba0b Matt*0134 #endif
cab777d667 Matt*0135
                0136       diagName  = 'BLGIRRK '
                0137       diagTitle = 'Nutrient-limited photosynthesis efficiency'
e0f9a7ba0b Matt*0138       diagUnits = 'W/m2            '
cab777d667 Matt*0139       diagCode  = 'SM P    MR      '
                0140       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0141      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0142
                0143       diagName  = 'BLGLLIM '
                0144       diagTitle = 'Light limitation'
e0f9a7ba0b Matt*0145       diagUnits = '-               '
cab777d667 Matt*0146       diagCode  = 'SM P    MR      '
                0147       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0148      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0149
e0f9a7ba0b Matt*0150 #ifndef USE_BLING_V1
cab777d667 Matt*0151       diagName  = 'BLGNLIM '
                0152       diagTitle = 'N limitation'
e0f9a7ba0b Matt*0153       diagUnits = '-               '
cab777d667 Matt*0154       diagCode  = 'SM P    MR      '
                0155       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0156      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e0f9a7ba0b Matt*0157 #endif
cab777d667 Matt*0158
                0159       diagName  = 'BLGPLIM '
                0160       diagTitle = 'P limitation'
e0f9a7ba0b Matt*0161       diagUnits = '-               '
cab777d667 Matt*0162       diagCode  = 'SM P    MR      '
                0163       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0164      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0165
                0166       diagName  = 'BLGFELIM'
                0167       diagTitle = 'Fe limitation'
e0f9a7ba0b Matt*0168       diagUnits = '-               '
                0169       diagCode  = 'SM P    MR      '
                0170       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0171      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0172
                0173 #ifdef USE_SIBLING
                0174       diagName  = 'BLGSILIM'
                0175       diagTitle = 'Si limitation'
                0176       diagUnits = '-               '
                0177       diagCode  = 'SM P    MR      '
                0178       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0179      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
00fa2d4ddd mmaz*0180
                0181       diagName  = 'BLGSIUP '
                0182       diagTitle = 'Silica uptake by bio'
a7e4405ea4 Jean*0183       diagUnits = 'mol Si/m3/s     '
00fa2d4ddd mmaz*0184       diagCode  = 'SM P    MR      '
                0185       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0186      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0187
                0188       diagName  = 'BLGSIREM'
                0189       diagTitle = 'Silica remin of sinking org matter'
a7e4405ea4 Jean*0190       diagUnits = 'mol Si/m3/s     '
00fa2d4ddd mmaz*0191       diagCode  = 'SM P    MR      '
                0192       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0193      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0194
                0195       diagName  = 'BLGSIREC'
                0196       diagTitle = 'Silica recycling'
a7e4405ea4 Jean*0197       diagUnits = 'mol Si/m3/s     '
00fa2d4ddd mmaz*0198       diagCode  = 'SM P    MR      '
                0199       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0200      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0201
                0202       diagName  = 'BLGSIDIS'
                0203       diagTitle = 'Silica dissolution fraction'
a7e4405ea4 Jean*0204       diagUnits = 'mol Si/m3/s     '
00fa2d4ddd mmaz*0205       diagCode  = 'SM P    MR      '
                0206       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0207      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0208
                0209       diagName  = 'BLGSI2N  '
                0210       diagTitle = 'Si to N stoichiometric ratio'
a7e4405ea4 Jean*0211       diagUnits = '-               '
00fa2d4ddd mmaz*0212       diagCode  = 'SM P    MR      '
                0213       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0214      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0215
e0f9a7ba0b Matt*0216 #endif
                0217
                0218 #ifndef USE_BLING_V1
                0219       diagName  = 'BLGP2N  '
                0220       diagTitle = 'P to N stoichiometric ratio'
                0221       diagUnits = '-               '
cab777d667 Matt*0222       diagCode  = 'SM P    MR      '
                0223       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0224      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0225
e0f9a7ba0b Matt*0226       diagName  = 'BLGFE2N '
                0227       diagTitle = 'Fe to N stoichiometric ratio'
                0228       diagUnits = '-               '
                0229       diagCode  = 'SM P    MR      '
                0230       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0231      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0232 #endif
                0233
cab777d667 Matt*0234 c carbon system
                0235
                0236       diagName  = 'BLGPCO2 '
                0237       diagTitle = 'Surface ocean partial pressure of CO2'
                0238       diagUnits = 'atm             '
                0239       diagCode  = 'SM P    M1      '
                0240       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0241      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0242
                0243       diagName  = 'BLGPH3D '
                0244       diagTitle = 'pH'
                0245       diagUnits = 'dimensionless   '
                0246       diagCode  = 'SM P    MR      '
                0247       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0248      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0249
                0250       diagName  = 'BLGOMAR '
                0251       diagTitle = 'Aragonite saturation'
                0252       diagUnits = 'dimensionless   '
                0253       diagCode  = 'SM P    MR      '
                0254       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0255      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0256
                0257       diagName  = 'BLGCFLX '
                0258       diagTitle = 'Flux of CO2 due to air-sea exch (mol/m2/s)'
                0259       diagUnits = 'mol/m2/sec      '
                0260       diagCode  = 'SM P    L1      '
                0261       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0262      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0263
e0f9a7ba0b Matt*0264       diagName  = 'BLGapco2'
cab777d667 Matt*0265       diagTitle = 'Atmospheric partial pressure of CO2'
                0266       diagUnits = '10^-6 ppm       '
                0267       diagCode  = 'SM      U1      '
                0268       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0269      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0270
                0271 c oxygen
                0272
                0273       diagName  = 'BLGOFLX '
                0274       diagTitle = 'Flux of O2 due to air-sea exch (mol/m2/s)'
e0f9a7ba0b Matt*0275       diagUnits = 'mol/m2/sec      '
cab777d667 Matt*0276       diagCode  = 'SM P    L1      '
                0277       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0278      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0279
b00a067069 Matt*0280       diagName  = 'BLGO2SAT'
                0281       diagTitle = 'surface O2 saturation (percent)'
e0f9a7ba0b Matt*0282       diagUnits = 'percent         '
b00a067069 Matt*0283       diagCode  = 'SM P    L1      '
                0284       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0285      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0286
cab777d667 Matt*0287 c iron
                0288
                0289       diagName  = 'BLGSURFE'
                0290       diagTitle = 'Surface deposition of Fe (mol/m2/s)'
e0f9a7ba0b Matt*0291       diagUnits = 'mol/m2/sec      '
cab777d667 Matt*0292       diagCode  = 'SM P    L1      '
                0293       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0294      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0295
                0296       diagName  = 'BLGFEUP '
                0297       diagTitle = 'Iron uptake by bio'
                0298       diagUnits = 'mol Fe/m3/s     '
                0299       diagCode  = 'SM P    MR      '
                0300       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0301      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0302
c0d1c06c15 Matt*0303       diagName  = 'BLGFESPM'
                0304       diagTitle = 'Iron sinking particle flux'
e337ca9d56 Matt*0305       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0306       diagCode  = 'SM P    MR      '
                0307       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0308      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0309
cab777d667 Matt*0310       diagName  = 'BLGFEREM'
                0311       diagTitle = 'Iron remin of sinking organic matter'
e337ca9d56 Matt*0312       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0313       diagCode  = 'SM P    MR      '
                0314       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0315      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0316
                0317       diagName  = 'BLGFERD '
                0318       diagTitle = 'Iron remin from Diel Vertical Migration'
e337ca9d56 Matt*0319       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0320       diagCode  = 'SM P    MR      '
                0321       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0322      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0323
cab777d667 Matt*0324       diagName  = 'BLGFEREC'
                0325       diagTitle = 'Iron recycling'
e337ca9d56 Matt*0326       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0327       diagCode  = 'SM P    MR      '
                0328       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0329      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0330
cab777d667 Matt*0331       diagName  = 'BLGFEADI'
                0332       diagTitle = 'Iron inorganic adsorption'
e337ca9d56 Matt*0333       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0334       diagCode  = 'SM P    MR      '
                0335       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0336      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0337
cab777d667 Matt*0338       diagName  = 'BLGFEADO'
                0339       diagTitle = 'Iron organic adsorption'
e337ca9d56 Matt*0340       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0341       diagCode  = 'SM P    MR      '
                0342       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0343      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0344
                0345       diagName  = 'BLGFESED'
                0346       diagTitle = 'Iron sediment source'
e337ca9d56 Matt*0347       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0348       diagCode  = 'SM P    MR      '
                0349       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0350      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0351
                0352       diagName  = 'BLGFEBUR'
                0353       diagTitle = 'Iron burial into sediment'
e337ca9d56 Matt*0354       diagUnits = 'mol Fe/m3/s     '
c0d1c06c15 Matt*0355       diagCode  = 'SM P    M1      '
                0356       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0357      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0358
cab777d667 Matt*0359 c nitrogen
c0d1c06c15 Matt*0360
e0f9a7ba0b Matt*0361 #ifndef USE_BLING_V1
cab777d667 Matt*0362       diagName  = 'BLGNUP  '
                0363       diagTitle = 'Nitrogen uptake by bio'
c0d1c06c15 Matt*0364       diagUnits = 'mol N/m3/s      '
                0365       diagCode  = 'SM P    MR      '
                0366       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0367      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0368
cab777d667 Matt*0369       diagName  = 'BLGNFIX '
                0370       diagTitle = 'Nitrogen fixation'
c0d1c06c15 Matt*0371       diagUnits = 'mol N/m3/s      '
                0372       diagCode  = 'SM P    MR      '
                0373       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0374      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0375
cab777d667 Matt*0376       diagName  = 'BLGNSPM '
                0377       diagTitle = 'Nitrogen sinking particle flux'
c0d1c06c15 Matt*0378       diagUnits = 'mol N/m3/s      '
                0379       diagCode  = 'SM P    MR      '
                0380       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0381      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0382
cab777d667 Matt*0383       diagName  = 'BLGNREM '
                0384       diagTitle = 'Nitrogen remin of sinking org matter'
c0d1c06c15 Matt*0385       diagUnits = 'mol N/m3/s      '
                0386       diagCode  = 'SM P    MR      '
                0387       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0388      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0389
cab777d667 Matt*0390       diagName  = 'BLGNDVM '
                0391       diagTitle = 'Nitrogen vertical flux due to DVM'
4ac06494d5 Matt*0392       diagUnits = 'mol N/m3/s      '
                0393       diagCode  = 'SM P    MR      '
                0394       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0395      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0396
cab777d667 Matt*0397       diagName  = 'BLGNRD  '
                0398       diagTitle = 'Nitrogen remin by Diel Vertical Migration'
c0d1c06c15 Matt*0399       diagUnits = 'mol N/m3/s      '
                0400       diagCode  = 'SM P    MR      '
                0401       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0402      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0403
                0404       diagName  = 'BLGNREC '
                0405       diagTitle = 'Nitrogen recycling'
                0406       diagUnits = 'mol N/m3/s      '
                0407       diagCode  = 'SM P    MR      '
                0408       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0409      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0410
cab777d667 Matt*0411       diagName  = 'BLGNDENB'
                0412       diagTitle = 'Benthic denitrification'
c0d1c06c15 Matt*0413       diagUnits = 'mol N/m3/s      '
                0414       diagCode  = 'SM P    MR      '
                0415       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0416      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0417
cab777d667 Matt*0418       diagName  = 'BLGNDENP'
                0419       diagTitle = 'Pelagic denitrification'
c0d1c06c15 Matt*0420       diagUnits = 'mol N/m3/s      '
                0421       diagCode  = 'SM P    MR      '
                0422       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0423      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0424
cab777d667 Matt*0425       diagName  = 'BLGDONP '
                0426       diagTitle = 'Dissolved Organic Nitrogen Production'
c0d1c06c15 Matt*0427       diagUnits = 'mol N/m3/s      '
                0428       diagCode  = 'SM P    MR      '
                0429       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0430      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0431
cab777d667 Matt*0432       diagName  = 'BLGDONR '
                0433       diagTitle = 'Dissolved Organic Nitrogen Remineralization'
                0434       diagUnits = 'mol N/m3/s      '
c0d1c06c15 Matt*0435       diagCode  = 'SM P    MR      '
                0436       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0437      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0438
e0f9a7ba0b Matt*0439       diagName  = 'BLGNONEN'
                0440       diagTitle = 'NO3 no neg'
                0441       diagUnits = '-               '
                0442       diagCode  = 'SM P    MR      '
                0443       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0444      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0445
                0446 #endif
                0447
                0448 c phosphorus
cab777d667 Matt*0449
                0450       diagName  = 'BLGPUP  '
                0451       diagTitle = 'Phosphorus uptake by bio'
c0d1c06c15 Matt*0452       diagUnits = 'mol P/m3/s      '
                0453       diagCode  = 'SM P    MR      '
                0454       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0455      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0456
                0457       diagName  = 'BLGPSPM '
                0458       diagTitle = 'Phosphorus sinking particle flux'
                0459       diagUnits = 'mol P/m3/s      '
                0460       diagCode  = 'SM P    MR      '
                0461       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0462      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0463
cab777d667 Matt*0464       diagName  = 'BLGPREM '
                0465       diagTitle = 'Phosphorus remin of sinking org matter'
c0d1c06c15 Matt*0466       diagUnits = 'mol P/m3/s      '
                0467       diagCode  = 'SM P    MR      '
                0468       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0469      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0470
cab777d667 Matt*0471       diagName  = 'BLGPDVM '
                0472       diagTitle = 'Phosphorus vertical flux due to DVM'
c0d1c06c15 Matt*0473       diagUnits = 'mol P/m3/s      '
                0474       diagCode  = 'SM P    MR      '
                0475       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0476      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0477
cab777d667 Matt*0478       diagName  = 'BLGPRD  '
                0479       diagTitle = 'Phosphorus remin by Diel Vertical Migration'
c0d1c06c15 Matt*0480       diagUnits = 'mol P/m3/s      '
                0481       diagCode  = 'SM P    MR      '
                0482       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0483      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0484
cab777d667 Matt*0485       diagName  = 'BLGPREC '
                0486       diagTitle = 'Phosphorus recycling'
c0d1c06c15 Matt*0487       diagUnits = 'mol P/m3/s      '
                0488       diagCode  = 'SM P    MR      '
                0489       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0490      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0491
cab777d667 Matt*0492       diagName  = 'BLGDOPP '
                0493       diagTitle = 'Dissolved Organic Phosphorus Production'
                0494       diagUnits = 'mol P/m3/s      '
c0d1c06c15 Matt*0495       diagCode  = 'SM P    MR      '
                0496       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0497      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0498
e0f9a7ba0b Matt*0499       diagName  = 'BLGDOPR '
                0500       diagTitle = 'Dissolved Organic Phosphorus Remineralization'
                0501       diagUnits = 'mol P/m3/s      '
                0502       diagCode  = 'SM P    MR      '
                0503       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0504      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0505
                0506 c mld
c0d1c06c15 Matt*0507
e0f9a7ba0b Matt*0508       diagName  = 'BLGMLD  '
                0509       diagTitle = 'Diagnosed mixed layer depth'
                0510       diagUnits = 'm               '
fe1e55ba7d Matt*0511       diagCode  = 'SM P    L1      '
c0d1c06c15 Matt*0512       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0513      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0514
e0f9a7ba0b Matt*0515 c light
                0516
c0d1c06c15 Matt*0517       diagName  = 'BLGIRRIS'
fe1e55ba7d Matt*0518       diagTitle = 'Local irradiance'
e0f9a7ba0b Matt*0519       diagUnits = 'W/m2            '
c0d1c06c15 Matt*0520       diagCode  = 'SM P    MR      '
                0521       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0522      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0523
                0524       diagName  = 'BLGIEFF '
                0525       diagTitle = 'Effective irradiance'
e0f9a7ba0b Matt*0526       diagUnits = 'W/m2            '
c0d1c06c15 Matt*0527       diagCode  = 'SM P    MR      '
                0528       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0529      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0530
                0531       diagName  = 'BLGIMEM '
                0532       diagTitle = 'Irradiance memory'
e0f9a7ba0b Matt*0533       diagUnits = 'W/m2            '
c0d1c06c15 Matt*0534       diagCode  = 'SM P    MR      '
                0535       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0536      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0537
cab777d667 Matt*0538 c bgc
                0539
                0540       diagName  = 'BLGBIOC '
                0541       diagTitle = 'Tendency of DIC due to biology (mol C/m3/s)'
                0542       diagUnits = 'mol C/m3/sec    '
                0543       diagCode  = 'SMRP    MR      '
c0d1c06c15 Matt*0544       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0545      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0546
cab777d667 Matt*0547       diagName  = 'BLGBIOAL'
                0548       diagTitle = 'Tendency of alkalinity due to biology (mol alk/m3/s)'
e0f9a7ba0b Matt*0549       diagUnits = 'mol alk/m3/sec  '
cab777d667 Matt*0550       diagCode  = 'SMRP    MR      '
4ac06494d5 Matt*0551       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0552      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0553
cab777d667 Matt*0554       diagName  = 'BLGBIOO2'
                0555       diagTitle = 'Tendency of oxygen due to biology (mol O/m3/s)'
                0556       diagUnits = 'mol O/m3/sec    '
                0557       diagCode  = 'SMRP    MR      '
c0d1c06c15 Matt*0558       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0559      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0560
cab777d667 Matt*0561       diagName  = 'BLGBION '
                0562       diagTitle = 'Tendency of nitrate due to biology (mol N/m3/s)'
                0563       diagUnits = 'mol N/m3/sec    '
                0564       diagCode  = 'SMRP    MR      '
                0565       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0566      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0567
                0568       diagName  = 'BLGBIOP '
                0569       diagTitle = 'Tendency of phosphate due to biology (mol P/m3/s)'
                0570       diagUnits = 'mol P/m3/sec    '
                0571       diagCode  = 'SMRP    MR      '
                0572       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0573      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0574
                0575       diagName  = 'BLGBIOFE'
                0576       diagTitle = 'Tendency of iron due to biology (mol Fe/m3/s)'
e0f9a7ba0b Matt*0577       diagUnits = 'mol Fe/m3/sec   '
cab777d667 Matt*0578       diagCode  = 'SMRP    MR      '
c0d1c06c15 Matt*0579       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0580      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0581
                0582       diagName  = 'BLGNSED '
                0583       diagTitle = 'Nitrogen flux in/out of sediments'
                0584       diagUnits = 'mol N/m3/s      '
                0585       diagCode  = 'SM P    M1      '
                0586       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0587      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0588
                0589       diagName  = 'BLGPSED '
                0590       diagTitle = 'Phosphorus flux in/out of sediments'
                0591       diagUnits = 'mol P/m3/s      '
                0592       diagCode  = 'SM P    M1      '
                0593       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0594      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0595
                0596       diagName  = 'BLGOSED '
                0597       diagTitle = 'Oxygen flux in/out of sediments'
                0598       diagUnits = 'mol O/m3/s      '
                0599       diagCode  = 'SM P    M1      '
                0600       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0601      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0602
0b0906f953 Matt*0603       diagName  = 'BLGCCdis'
                0604       diagTitle = 'Calcium carbonate dissolution'
                0605       diagUnits = 'mol C/m3/s      '
                0606       diagCode  = 'SM P    MR      '
                0607       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0608      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0609
                0610       diagName  = 'BLGCCpro'
                0611       diagTitle = 'Calcium carbonate production'
                0612       diagUnits = 'mol C/m3/s      '
                0613       diagCode  = 'SM P    MR      '
                0614       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0615      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0616
c0d1c06c15 Matt*0617 #endif /* ALLOW_DIAGNOSTICS */
                0618
                0619       RETURN
                0620       END