pkg/bling/BLING_OPTIONS.h

c0d1c06c15 Matt*0001 #ifndef BLING_OPTIONS_H
                0002 #define BLING_OPTIONS_H
                0003 #include "PACKAGES_CONFIG.h"
                0004 #include "CPP_OPTIONS.h"
                0005
                0006 #ifdef ALLOW_BLING
                0007 C     Package-specific Options & Macros go here
                0008
a85293d087 Mart*0009 C BLING+Nitrogen is the default model. It is a version
a284455135 Matt*0010 C of BLING with 8 tracers and 3 phyto classes.
                0011 C For the original 6-tracer model of Galbraith et al (2010),
                0012 C define USE_BLING_V1 - but note the different order of tracers in data.ptracers
e0f9a7ba0b Matt*0013 #undef USE_BLING_V1
                0014
a284455135 Matt*0015 C Options for BLING+Nitrogen code:
                0016 C SiBLING: add a 9th tracer for silica
e0f9a7ba0b Matt*0017 #undef USE_SIBLING
a284455135 Matt*0018 C apply remineralization from diel vertical migration
e0f9a7ba0b Matt*0019 #undef USE_BLING_DVM
a284455135 Matt*0020 C active tracer for total phytoplankton biomass
52893fa481 Matt*0021 #undef ADVECT_PHYTO
a284455135 Matt*0022 C sub grid scale sediments - NOT IMPLEMENTED YET
e0f9a7ba0b Matt*0023 c #undef USE_SGS_SED
4ac06494d5 Matt*0024
a284455135 Matt*0025 C Prevents negative values in nutrient fields
c0d1c06c15 Matt*0026 #define BLING_NO_NEG
                0027
a284455135 Matt*0028 C Use Liebig function instead of geometric mean of the
                0029 C nutrient limitations to calculate maximum phyto growth rate
e0f9a7ba0b Matt*0030 #define MIN_NUT_LIM
                0031
a284455135 Matt*0032 C Allow different phytoplankton groups to have different growth rates and
9edc0e3a85 aver*0033 C nutrient/light limitations. Parameters implemented have yet to be tuned
                0034 #undef SIZE_DEP_LIM
a284455135 Matt*0035
                0036 C Assume that phytoplankton in the mixed layer experience
                0037 C the average light over the mixed layer (as in original BLING model)
00fa2d4ddd mmaz*0038 #undef ML_MEAN_LIGHT
c0d1c06c15 Matt*0039
a284455135 Matt*0040 C Assume that phytoplankton are homogenized in the mixed layer
e0f9a7ba0b Matt*0041 #define ML_MEAN_PHYTO
c0d1c06c15 Matt*0042
a284455135 Matt*0043 C Calculate MLD using a threshold criterion. If undefined,
                0044 C MLD is calculated using the second derivative of rho(z)
e0f9a7ba0b Matt*0045 #undef BLING_USE_THRESHOLD_MLD
88bed3f0e8 Matt*0046
a284455135 Matt*0047 C Determine PAR from shortwave radiation Qsw;
                0048 C otherwise determined from date and latitude
9edc0e3a85 aver*0049 C (Do not define if not using pkg/exf)
                0050 #undef USE_QSW
4ac06494d5 Matt*0051
a284455135 Matt*0052 C Light absorption scheme from Manizza et al. (2005),
                0053 C with self shading from phytoplankton
00fa2d4ddd mmaz*0054 #undef PHYTO_SELF_SHADING
                0055
a284455135 Matt*0056 C Note: atm pressure from PKG/EXF is always used for air-sea flux calculation
                0057 C if available; otherwise read from file or set to constant value (1 atm)
c0d1c06c15 Matt*0058
a284455135 Matt*0059 C Note: winds from PKG/EXF are always used if available;
                0060 C otherwise read from file or set to constant value (5 m/s)
4ac06494d5 Matt*0061
a284455135 Matt*0062 C Note: ice fraction from PKG/SEAICE or THSICE is always used if available;
                0063 C otherwise read from file or set to constant value (0)
88bed3f0e8 Matt*0064
a284455135 Matt*0065 C Note: atm pCO2 from EXF file is always used if available;
                0066 C otherwise set to constant value in data.bling
c0d1c06c15 Matt*0067
a284455135 Matt*0068 C Simplify some parts of the code that are problematic when using the adjoint
4ac06494d5 Matt*0069 #define BLING_ADJOINT_SAFE
                0070
a284455135 Matt*0071 C For adjoint safe, do not call bling_dvm
4ac06494d5 Matt*0072 #ifdef BLING_ADJOINT_SAFE
f40bb882f5 Jean*0073 # undef USE_BLING_DVM
4ac06494d5 Matt*0074 #endif
c0d1c06c15 Matt*0075
e0f9a7ba0b Matt*0076 C ABIOTIC OPTIONS
a85293d087 Mart*0077 C Compile "Solvesaphe" package (Munhoven 2013) for pH/pCO2
e0f9a7ba0b Matt*0078 C  can still select Follows et al (2006) solver in data.bling,
9edc0e3a85 aver*0079 C  but will use solvesaphe dissociation coefficient options
e0f9a7ba0b Matt*0080 #undef CARBONCHEM_SOLVESAPHE
                0081
                0082 C In S/R CARBON_CHEM convert ak1 and ak2 to the total pH scale
                0083 C  consistent with other coefficients (currently on the seawater scale).
                0084 C NOTE: Has NO effect when CARBONCHEM_SOLVESAPHE is defined (different
                0085 C  coeffs are used).
                0086 #undef CARBONCHEM_TOTALPHSCALE
                0087
a284455135 Matt*0088 C When calculating the fraction of sinking organic matter, use model biomass diagnostics.
                0089 #define NEW_FRAC_EXP
                0090
c0d1c06c15 Matt*0091 #endif /* ALLOW_BLING */
                0092 #endif /* BLING_OPTIONS_H */