global_oce_biogeo_bling/code_tap/BLING_OPTIONS.h

4dea327916 aver*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
                0009 C BLING+Nitrogen is the default model. It is a version
                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
                0013 #undef USE_BLING_V1
                0014
                0015 C Options for BLING+Nitrogen code:
                0016 C SiBLING: add a 9th tracer for silica
                0017 #undef USE_SIBLING
                0018 C apply remineralization from diel vertical migration
                0019 #undef USE_BLING_DVM
                0020 C active tracer for total phytoplankton biomass
                0021 #undef ADVECT_PHYTO
                0022 C sub grid scale sediments - NOT IMPLEMENTED YET
                0023 c #undef USE_SGS_SED
                0024
                0025 C Prevents negative values in nutrient fields
                0026 #define BLING_NO_NEG
                0027
                0028 C Use Liebig function instead of geometric mean of the
                0029 C nutrient limitations to calculate maximum phyto growth rate
                0030 #define MIN_NUT_LIM
                0031
                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 #define SIZE_DEP_LIM
4dea327916 aver*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)
                0038 #define ML_MEAN_LIGHT
                0039
                0040 C Assume that phytoplankton are homogenized in the mixed layer
                0041 #define ML_MEAN_PHYTO
                0042
                0043 C Calculate MLD using a threshold criterion. If undefined,
                0044 C MLD is calculated using the second derivative of rho(z)
                0045 #undef BLING_USE_THRESHOLD_MLD
                0046
                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)
4dea327916 aver*0050 #undef USE_QSW
                0051
                0052 C Light absorption scheme from Manizza et al. (2005),
                0053 C with self shading from phytoplankton
                0054 #define PHYTO_SELF_SHADING
                0055
                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)
                0058
                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)
                0061
                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)
                0064
                0065 C Note: atm pCO2 from EXF file is always used if available;
                0066 C otherwise set to constant value in data.bling
                0067
                0068 C Simplify some parts of the code that are problematic when using the adjoint
                0069 #define BLING_ADJOINT_SAFE
                0070
                0071 C For adjoint safe, do not call bling_dvm
                0072 #ifdef BLING_ADJOINT_SAFE
                0073 # undef USE_BLING_DVM
                0074 #endif
                0075
                0076 C ABIOTIC OPTIONS
                0077 C Compile "Solvesaphe" package (Munhoven 2013) for pH/pCO2
                0078 C  can still select Follows et al (2006) solver in data.bling,
9edc0e3a85 aver*0079 C  but will use solvesaphe dissociation coefficient options
4dea327916 aver*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
                0088 C When calculating the fraction of sinking organic matter, use model biomass diagnostics.
                0089 #define NEW_FRAC_EXP
                0090
                0091 #endif /* ALLOW_BLING */
                0092 #endif /* BLING_OPTIONS_H */