File indexing completed on 2019-08-13 05:10:33 UTC
view on githubraw file Latest commit 9f0da36f on 2019-07-20 15:57:26 UTC
c0d1c06c15 Matt*0001 #include "BLING_OPTIONS.h"
041f4605e9 Jean*0002 #ifdef ALLOW_EXF
0003 # include "EXF_OPTIONS.h"
0004 #endif
c0d1c06c15 Matt*0005
0006
be72e7ae9e Jean*0007 SUBROUTINE BLING_READPARMS( myThid )
c0d1c06c15 Matt*0008
be72e7ae9e Jean*0009
c0d1c06c15 Matt*0010
0011
be72e7ae9e Jean*0012
c0d1c06c15 Matt*0013
e0f9a7ba0b Matt*0014 IMPLICIT NONE
c0d1c06c15 Matt*0015
0016
0017 #include "SIZE.h"
0018 #include "EEPARAMS.h"
0019 #include "PARAMS.h"
e0f9a7ba0b Matt*0020 #ifdef ALLOW_EXF
041f4605e9 Jean*0021 # include "EXF_PARAM.h"
9f0da36f91 Jean*0022 # include "EXF_INTERP_SIZE.h"
0023 # include "EXF_INTERP_PARAM.h"
e0f9a7ba0b Matt*0024 #endif /* ALLOW_EXF */
c0d1c06c15 Matt*0025 #include "BLING_VARS.h"
0026
0027
0028
0029 INTEGER myThid
0030
0031
0032 #ifdef ALLOW_BLING
0033
0034
0035
530382a175 Jean*0036
c0d1c06c15 Matt*0037
0038 CHARACTER*(MAX_LEN_MBUF) msgBuf
530382a175 Jean*0039 INTEGER errCount
c0d1c06c15 Matt*0040 INTEGER iUnit
0041
0042
0043
e0f9a7ba0b Matt*0044
0045
0046
0047
0048
0049
0050
0051
0052
0053
0054
0055
0056
0057
0058
0059
0060
0061
0062
0063
0064
0065
0066
0067
0068
0069
0070
0071
0072
0073
0074
c0d1c06c15 Matt*0075
0076
be72e7ae9e Jean*0077 NAMELIST /ABIOTIC_PARMS/
e0f9a7ba0b Matt*0078 & permil,
0079 & Pa2Atm,
0080 & epsln
0081 #ifdef CARBONCHEM_SOLVESAPHE
0082 & ,selectBTconst,selectFTconst,
0083 & selectHFconst,selectK1K2const,
0084 & selectPHsolver
0085 #endif
c0d1c06c15 Matt*0086
0087
0088
0089
0090
0091 NAMELIST /BIOTIC_PARMS/
0092 & CtoN,
e0f9a7ba0b Matt*0093 & CtoP,
0094 & NtoP,
c0d1c06c15 Matt*0095 & HtoC,
e0f9a7ba0b Matt*0096 & NO3toN,
c0d1c06c15 Matt*0097 & O2toN,
e0f9a7ba0b Matt*0098 & O2toP,
c0d1c06c15 Matt*0099 & CatoN,
e0f9a7ba0b Matt*0100 & CatoP,
c0d1c06c15 Matt*0101 & masstoN,
e0f9a7ba0b Matt*0102 #ifndef USE_BLING_V1
0103 & Pc_0_diaz,
c0d1c06c15 Matt*0104 & alpha_photo,
e0f9a7ba0b Matt*0105 & gamma_DON,
0106 & k_Fe_diaz,
0107 & k_NO3,
0108 & k_PtoN,
0109 & k_FetoN,
0110 & PtoN_min,
0111 & PtoN_max,
0112 & FetoN_min,
0113 & FetoN_max,
0114 & kappa_eppley_diaz,
0115 & phi_dvm,
0116 & sigma_dvm,
0117 #ifdef USE_SIBLING
0118 & k_Si,
0119 & gamma_Si_0,
0120 & kappa_remin_Si,
0121 & wsink_Si,
00fa2d4ddd mmaz*0122 & SitoN_uptake_min,
0123 & SitoN_uptake_max,
0124 & SitoN_uptake_scale,
0125 & SitoN_uptake_exp,
e0f9a7ba0b Matt*0126 & q_SitoN_diss,
0127 #endif
0128 #else
0129 & alpha_max,
0130 & alpha_min,
0131 & gamma_biomass,
0132 & k_FetoP,
0133 & FetoP_max,
0134 & Fe_lim_min,
0135 #endif
0136 & pivotal,
0137 & Pc_0,
0138 & lambda_0,
0139 & resp_frac,
0140 & chl_min,
c0d1c06c15 Matt*0141 & theta_Fe_max_hi,
0142 & theta_Fe_max_lo,
0143 & gamma_irr_mem,
0144 & gamma_DOP,
0145 & gamma_POM,
0146 & k_Fe,
0147 & k_O2,
0148 & k_PO4,
0149 & kFe_eq_lig_max,
0150 & kFe_eq_lig_min,
0151 & kFe_eq_lig_Femin,
be72e7ae9e Jean*0152 & kFe_eq_lig_irr,
c0d1c06c15 Matt*0153 & kFe_org,
0154 & kFe_inorg,
0155 & FetoC_sed,
0156 & remin_min,
0157 & oxic_min,
0158 & ligand,
0159 & kappa_eppley,
0160 & kappa_remin,
0161 & ca_remin_depth,
0162 & phi_DOM,
0163 & phi_sm,
0164 & phi_lg,
0165 & wsink0z,
0166 & wsink0,
0167 & wsinkacc,
0168 & parfrac,
0169 & alpfe,
0170 & k0,
e0f9a7ba0b Matt*0171 & MLmix_max,
0172 & river_conc_po4,
0173 & river_dom_to_nut
c0d1c06c15 Matt*0174
0175
0176
0177
0178
0179 NAMELIST /BLING_FORCING/
0180 & bling_windFile, bling_atmospFile, bling_iceFile,
54b6d6aa91 Matt*0181 & bling_ironFile, bling_silicaFile,
0182 & bling_psmFile, bling_plgFile, bling_PdiazFile,
c0d1c06c15 Matt*0183 & bling_forcingPeriod, bling_forcingCycle,
47de7c1e0e Matt*0184 & bling_Pc_2dFile, bling_Pc_2d_diazFile,
af8f3b9ffc Matt*0185 & bling_alpha_photo2dFile,bling_phi_DOM2dFile,
0186 & bling_k_Fe2dFile, bling_k_Fe_diaz2dFile,
0187 & bling_gamma_POM2dFile, bling_wsink0_2dFile,
e0f9a7ba0b Matt*0188 & bling_phi_sm2dFile,bling_phi_lg2dFile,
0189 & bling_pCO2,
0190 & river_conc_po4, river_dom_to_nut
0191 #ifdef ALLOW_EXF
0192 & ,apco2file, apco2startdate1, apco2startdate2,
041f4605e9 Jean*0193 & apco2RepCycle, apco2period, apco2StartTime,
530382a175 Jean*0194 & exf_inscal_apco2, exf_outscal_apco2, apco2const,
079948e6a6 Matt*0195 & apco2_exfremo_intercept, apco2_exfremo_slope
0196 #ifdef USE_EXF_INTERPOLATION
e0f9a7ba0b Matt*0197 & ,apco2_lon0, apco2_lon_inc, apco2_lat0, apco2_lat_inc,
079948e6a6 Matt*0198 & apco2_nlon, apco2_nlat, apco2_interpMethod
0199 #endif /* USE_EXF_INTERPOLATION */
e0f9a7ba0b Matt*0200 #endif /* ALLOW_EXF */
c0d1c06c15 Matt*0201
0202
e0f9a7ba0b Matt*0203
0204
0205
0206
cab777d667 Matt*0207
e0f9a7ba0b Matt*0208
0209
0210
0211
0212
0213
0214
0215
0216
0217
0218
0219
0220
cab777d667 Matt*0221
0222
e0f9a7ba0b Matt*0223
0224
0225
0226
0227
0228
0229
0230
0231
0232
0233
0234
0235
0236
0237
0238
0239
0240
0241
0242
0243
0244
0245
0246
0247
0248
0249
0250
0251
0252
0253
0254
0255
0256
0257
0258
0259
00fa2d4ddd mmaz*0260
0261
0262
0263
e0f9a7ba0b Matt*0264
0265
0266
0267
cab777d667 Matt*0268
e0f9a7ba0b Matt*0269
0270
0271
0272
0273
cab777d667 Matt*0274
0275
e0f9a7ba0b Matt*0276
0277
0278
0279
0280
0281
0282
0283
0284
0285
0286
0287
0288
c0d1c06c15 Matt*0289 _RL secperday
e0f9a7ba0b Matt*0290 #ifdef ALLOW_EXF
0291 #ifdef USE_EXF_INTERPOLATION
079948e6a6 Matt*0292 INTEGER j
e0f9a7ba0b Matt*0293 #endif /* USE_EXF_INTERPOLATION */
0294 #endif /* ALLOW_EXF */
c0d1c06c15 Matt*0295
0296 _BEGIN_MASTER(myThid)
530382a175 Jean*0297 errCount = 0
c0d1c06c15 Matt*0298
0299
0300
0301
0302 secperday = 86400. _d 0
0303 permil = 1. _d 0 / 1024.5 _d 0
0304 Pa2Atm = 1.01325 _d 5
e0f9a7ba0b Matt*0305 epsln = 1. _d -30
0306 #ifdef CARBONCHEM_SOLVESAPHE
0307 selectBTconst = 1
0308 selectFTconst = 1
0309 selectHFconst = 1
0310 selectK1K2const = 1
0311 selectPHsolver = 0
0312 #endif
0313
c0d1c06c15 Matt*0314 CtoN = 6.75 _d 0
e0f9a7ba0b Matt*0315 CtoP = 106. _d 0
0316 NtoP = 16. _d 0
c0d1c06c15 Matt*0317 HtoC = 48. _d 0 / 106. _d 0
be72e7ae9e Jean*0318 NO3toN = CtoN * (1. _d 0 + 0.25 _d 0 * HtoC)
c0d1c06c15 Matt*0319 & * 0.8 _d 0 + 0.6 _d 0
e0f9a7ba0b Matt*0320 O2toN = CtoN * (1. _d 0 + 0.25 _d 0 * HtoC)
0321 & + 2. _d 0
0322 O2toP = 150. _d 0
c0d1c06c15 Matt*0323 CatoN = CtoN * 0.015 _d 0
e0f9a7ba0b Matt*0324 CatoP = 106.0 _d 0 * 0.015 _d 0
c0d1c06c15 Matt*0325 masstoN = CtoN * 12.001 _d 0
e0f9a7ba0b Matt*0326 #ifndef USE_BLING_V1
c0d1c06c15 Matt*0327 pivotal = 1.9 _d -3 / 1028. _d 0 / CtoN / permil
0328 Pc_0 = 1.7 _d -5
0329 Pc_0_diaz = 0.01 _d -5
be72e7ae9e Jean*0330 alpha_photo = 0.7 _d -5 * 2.77 _d 18 / 6.022 _d 17
c0d1c06c15 Matt*0331 gamma_DON = 0.25 _d 0 / (365.25 _d 0 * secperday)
0332 gamma_DOP = 0.5 _d 0 / (365.25 _d 0 * secperday)
e0f9a7ba0b Matt*0333 PtoN_min = 1. / 28.
0334 PtoN_max = 1. / 9.
0335 FetoN_min = 2. _d -6 * 6.75
0336 FetoN_max = 25. _d -6 * 6.75
cab777d667 Matt*0337 k_Fe = 1.6 _d -10 / permil
c0d1c06c15 Matt*0338 k_Fe_diaz = 7. _d -10 / permil
0339 k_NO3 = 2. _d -6 / permil
cab777d667 Matt*0340 k_PO4 = 1. _d -8 / permil
be72e7ae9e Jean*0341 k_PtoN = 1.5 _d -6 / permil
0342 k_FetoN = 8. _d -10 / permil
e0f9a7ba0b Matt*0343 remin_min = 0.15 _d 0
0344 kappa_eppley_diaz = 0.18 _d 0
0345 sigma_dvm = 40.0 _d 0
0346 phi_dvm = 0.2 _d 0
0347 #ifdef USE_SIBLING
0348 k_Si = 2.0 _d -6 / permil
0349 gamma_Si_0 = 0.05 / secperday
0350 kappa_remin_Si = 0.075
0351 wsink_Si = 100. / secperday
6ffd1aa797 Jean*0352 SitoN_uptake_min = 4
0353 SitoN_uptake_max = 1
00fa2d4ddd mmaz*0354 SitoN_uptake_scale = 0.6
0355 SitoN_uptake_exp = 3.887
e0f9a7ba0b Matt*0356 q_SitoN_diss = 1.0 _d 0
0357 #endif
0358 #else
0359 pivotal = 1.9 _d -3 / 1028. _d 0 / CtoP / permil
0360 Pc_0 = 1 _d -5
0361 alpha_max = 1.6 _d -5 * 2.77 _d 18 / 6.022 _d 17
0362 alpha_min = 0.4 _d -5 * 2.77 _d 18 / 6.022 _d 17
0363 gamma_biomass = 0.5 _d 0 / secperday
0364 gamma_DOP = 0.25 _d 0 / (365.25 _d 0 * secperday)
0365 k_Fe = 8. _d -10 / permil
0366 k_PO4 = 1. _d -7 / permil
0367 k_FetoP = 7. _d -6 * CtoP
0368 FetoP_max = 28. _d -6 * CtoP
0369 Fe_lim_min = 0. _d 0
0370 remin_min = 0.3 _d 0
0371 #endif
0372 lambda_0 = 0.19 _d 0 / secperday
0373 resp_frac = 0. _d 0
0374 chl_min = 1. _d -5
0375 theta_Fe_max_hi = 0.04 _d 0
0376 theta_Fe_max_lo = 0.01 _d 0
0377 gamma_irr_mem = 1. _d 0 / secperday
0378 gamma_POM = 0.12 _d 0 / secperday
0379 k_O2 = 20. _d -6 / permil
c0d1c06c15 Matt*0380 kFe_eq_lig_max = 8.0 _d 10 * permil
0381 kFe_eq_lig_min = 8.0 _d 9 * permil
0382 kFe_eq_lig_Femin = 0.05 _d -9 / permil
0383 kFe_eq_lig_irr = 0.1 _d 0
0384 kFe_org = 0.5 _d 0 / secperday * permil**(0.58)
0385 kFe_inorg = 1. _d 3 / secperday * permil**(0.5)
e0f9a7ba0b Matt*0386 FetoC_sed = 1. _d -4
c0d1c06c15 Matt*0387 oxic_min = 1. _d -6 / permil
0388 Ligand = 1. _d -9 / permil
0389 kappa_eppley = 0.063 _d 0
be72e7ae9e Jean*0390 kappa_remin = -0.032 _d 0
c0d1c06c15 Matt*0391 ca_remin_depth = 1343. _d 0
0392 phi_DOM = 0.1 _d 0
0393 phi_sm = 0.18 _d 0
0394 phi_lg = 1. _d 0
0395 wsink0 = 16. _d 0 / secperday
e0f9a7ba0b Matt*0396 wsink0z = 80. _d 0
c0d1c06c15 Matt*0397 wsinkacc = 0.05 _d 0 / secperday
0398 parfrac = 0.4 _d 0
0399 alpfe = 0.01 _d 0
0400 k0 = 0.04 _d 0
e0f9a7ba0b Matt*0401 MLmix_max = 200.0 _d 0
0402
0403 bling_windFile = ' '
0404 bling_atmospFile = ' '
0405 bling_iceFile = ' '
0406 bling_ironFile = ' '
0407 bling_silicaFile = ' '
0408 bling_psmFile = ' '
0409 bling_plgFile = ' '
0410 bling_pdiazFile = ' '
47de7c1e0e Matt*0411 bling_Pc_2dFile = ' '
0412 bling_Pc_2d_diazFile = ' '
0413 bling_alpha_photo2dFile= ' '
0414 bling_k_Fe2dFile = ' '
0415 bling_k_Fe_diaz2dFile = ' '
af8f3b9ffc Matt*0416 bling_gamma_POM2dFile = ' '
0417 bling_wsink0_2dFile = ' '
0418 bling_phi_DOM2dFile = ' '
0419 bling_phi_sm2dFile = ' '
0420 bling_phi_lg2dFile = ' '
c0d1c06c15 Matt*0421
e0f9a7ba0b Matt*0422 bling_pCO2 = 278. _d -6
0423 river_conc_po4 = 0.007 _d 0
0424 river_dom_to_nut = 0. _d 0
0425
0426 #ifdef ALLOW_EXF
0427 apco2startdate1 = 0
0428 apco2startdate2 = 0
0429 apco2StartTime = UNSET_RL
0430 apco2period = 0.0 _d 0
0431 apco2RepCycle = repeatPeriod
0432 apco2const = 0.0 _d 0
079948e6a6 Matt*0433 apco2_exfremo_intercept = 0.0 _d 0
e0f9a7ba0b Matt*0434 apco2_exfremo_slope = 0.0 _d 0
0435 apco2file = ' '
0436 exf_inscal_apco2 = 1. _d 0
0437 exf_outscal_apco2 = 1. _d 0
079948e6a6 Matt*0438 #ifdef USE_EXF_INTERPOLATION
530382a175 Jean*0439
079948e6a6 Matt*0440
e0f9a7ba0b Matt*0441 apco2_lon0 = inp_lon0
0442 apco2_lat0 = inp_lat0
9f0da36f91 Jean*0443 apco2_nlon = inp_gNx
0444 apco2_nlat = inp_gNy
e0f9a7ba0b Matt*0445 apco2_lon_inc = inp_dLon
079948e6a6 Matt*0446 DO j=1,MAX_LAT_INC
9f0da36f91 Jean*0447 apco2_lat_inc(j) = inp_dLat(j)
079948e6a6 Matt*0448 ENDDO
8125504ce6 Matt*0449 apco2_interpMethod = 1
079948e6a6 Matt*0450 #endif /* USE_EXF_INTERPOLATION */
e0f9a7ba0b Matt*0451 #endif /* ALLOW_EXF */
c0d1c06c15 Matt*0452
be72e7ae9e Jean*0453
e0f9a7ba0b Matt*0454 bling_forcingPeriod = externForcingPeriod
0455 bling_forcingCycle = externForcingCycle
c0d1c06c15 Matt*0456
0457 WRITE(msgBuf,'(A)') ' BLING_READPARMS: opening data.bling'
0458 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
0459 I SQUEEZE_RIGHT, myThid )
0460
0461 CALL OPEN_COPY_DATA_FILE( 'data.bling', 'BLING_READPARMS',
0462 O iUnit, myThid )
0463
0464
0465
0466
0467 READ(UNIT=iUnit,NML=ABIOTIC_PARMS)
0468
0469
0470 READ(UNIT=iUnit,NML=BIOTIC_PARMS)
0471
0472
0473 READ(UNIT=iUnit,NML=BLING_FORCING)
0474
0475 WRITE(msgBuf,'(A)')
0476 & ' BLING_READPARMS: finished reading data.BLING'
0477 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
0478 I SQUEEZE_RIGHT, myThid )
0479
0480
7a77863887 Mart*0481 #ifdef SINGLE_DISK_IO
c0d1c06c15 Matt*0482 CLOSE(iUnit)
7a77863887 Mart*0483 #else
0484 CLOSE(iUnit,STATUS='DELETE')
0485 #endif /* SINGLE_DISK_IO */
c0d1c06c15 Matt*0486
e0f9a7ba0b Matt*0487 iUnit = standardMessageUnit
0488
0489 #ifdef CARBONCHEM_SOLVESAPHE
0490 IF ( selectPHsolver.GT.0 ) THEN
0491 WRITE(msgBuf,'(A)')
0492 & 'Using Munhoven (2013) Solvesaphe for pH/pCO2'
0493 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0494 ELSEIF ( selectPHsolver.EQ.0 ) THEN
0495 WRITE(msgBuf,'(A)')
0496 & 'Using Follows et al. (2006) for pH/pCO2'
0497 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0498 ENDIF
0499 #endif
c0d1c06c15 Matt*0500
e0f9a7ba0b Matt*0501
0502 #ifdef CARBONCHEM_SOLVESAPHE
0503 WRITE(msgBuf,'(A)')
0504 & 'Using Munhoven (2013) Solvesaphe carbon coefficients'
0505 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0506
0507 IF ( selectK1K2const.EQ.1 ) THEN
0508 WRITE(msgBuf,'(A)')
0509 & 'Using Millero (1995)/Mehrbach K1 and K2 coefficients'
0510 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0511 ELSEIF ( selectK1K2const.EQ.2 ) THEN
0512 WRITE(msgBuf,'(A)')
0513 & 'Using Roy et al. (1993) K1 and K2 coefficients'
0514 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0515 ELSEIF ( selectK1K2const.EQ.3 ) THEN
0516 WRITE(msgBuf,'(A)')
0517 & 'Using Millero (1995) "consensus" K1 and K2 coefficients'
0518 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0519 ELSEIF ( selectK1K2const.EQ.4 ) THEN
0520 WRITE(msgBuf,'(A)')
0521 & 'Using Luecker et al. (2000) K1 and K2 coefficients'
0522 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0523 ELSEIF ( selectK1K2const.EQ.5 ) THEN
0524 WRITE(msgBuf,'(A)')
0525 & 'Using Millero et al. (2010) K1 and K2 coefficients'
0526 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0527 ELSEIF ( selectK1K2const.EQ.6 ) THEN
0528 WRITE(msgBuf,'(A)')
0529 & 'Using Waters et al. (2014) K1 and K2 coefficients'
0530 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0531 ENDIF
0532 IF ( selectHFconst.EQ.1 ) THEN
0533 WRITE(msgBuf,'(A)')
0534 & 'Using Dickson and Riley (1979) KF coefficient'
0535 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0536 ELSEIF ( selectHFconst.EQ.2 ) THEN
0537 WRITE(msgBuf,'(A)')
0538 & 'Using Perez and Fraga (1987) KF coefficient'
0539 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0540 ENDIF
0541 IF ( selectBTconst.EQ.1 ) THEN
0542 WRITE(msgBuf,'(A)')
0543 & 'Using Uppstrom (1974) BT estimation from salinity'
0544 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0545 ELSEIF ( selectBTconst.EQ.2 ) THEN
0546 WRITE(msgBuf,'(A)')
0547 & 'Using Lee et al (2010) BT estimation from salinity'
0548 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0549 ENDIF
0550 IF ( selectFTconst.EQ.1 ) THEN
0551 WRITE(msgBuf,'(A)')
0552 & 'Using Riley (1965) FT estimation from salinity'
0553 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0554 ELSEIF ( selectBTconst.EQ.2 ) THEN
0555 WRITE(msgBuf,'(A)')
0556 & 'Using Culkin (1965) FT estimation from salinity'
0557 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0558 ENDIF
0559 #else
0560
0561 WRITE(msgBuf,'(A)')
0562 & 'Using Millero (1995)/Mehrbach K1 and K2 coefficients'
0563 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0564 WRITE(msgBuf,'(A)')
0565 & 'Using Dickson and Riley (1979) KF coefficient'
0566 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0567 WRITE(msgBuf,'(A)')
0568 & 'Using Uppstrom (1974) BT estimation from salinity'
0569 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
0570 WRITE(msgBuf,'(A)')
0571 & 'Using Riley (1965) FT estimation from salinity'
0572 CALL PRINT_MESSAGE(msgBuf,iUnit,SQUEEZE_RIGHT,myThid)
c0d1c06c15 Matt*0573 #endif
0574
530382a175 Jean*0575 IF ( errCount.GE.1 ) THEN
0576 WRITE(msgBuf,'(A,I3,A)')
0577 & 'BLING_READPARMS: detected', errCount,' fatal error(s)'
0578 CALL PRINT_ERROR( msgBuf, myThid )
0579 CALL ALL_PROC_DIE( 0 )
0580 STOP 'ABNORMAL END: S/R BLING_READPARMS'
0581 ENDIF
0582
e0f9a7ba0b Matt*0583 #ifdef CARBONCHEM_SOLVESAPHE
0584 IF ( selectHFconst.GT.2 .OR. selectBTconst.GT.2
0585 & .OR. selectFTconst.GT.2 ) THEN
0586 WRITE(msgBuf,'(A)')
0587 & 'BLING_READPARMS: selectHF, selectBT or selectFT > 2'
0588 CALL PRINT_ERROR( msgBuf, myThid )
0589 STOP 'ABNORMAL END: S/R BLING_READPARMS: solvesaphe_const error'
0590 ELSEIF ( selectK1K2const.GT.6 ) THEN
0591 WRITE(msgBuf,'(A)')
0592 & 'BLING_READPARMS: selectK1K2const > 6 '
0593 CALL PRINT_ERROR( msgBuf, myThid )
0594 STOP 'ABNORMAL END: S/R BLING_READPARMS: solvesaphe_const error'
0595 ELSEIF ( selectPHsolver.GT.3 ) THEN
0596 WRITE(msgBuf,'(A)')
0597 & 'BLING_READPARMS: selectPHsolver > 3 '
0598 CALL PRINT_ERROR( msgBuf, myThid )
0599 STOP 'ABNORMAL END: S/R BLING_READPARMS: select pH solver error'
0600 ENDIF
0601 #endif
0602
c0d1c06c15 Matt*0603 _END_MASTER(myThid)
0604
0605
0606 _BARRIER
0607
0608 #endif /* ALLOW_BLING */
0609
0610 RETURN
0611 END