pkg/bling/bling_init_fixed.F

c0d1c06c15 Matt*0001 #include "BLING_OPTIONS.h"
                0002
                0003 CBOP
                0004       subroutine BLING_INIT_FIXED( myThid )
                0005
                0006 C     ==========================================================
                0007 C     | subroutine bling_init_fixed
cab777d667 Matt*0008 C     | o Initialize fixed parameters
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 #include "PARAMS.h"
                0017 #include "GRID.h"
                0018 #include "BLING_VARS.h"
                0019
                0020 C     !INPUT PARAMETERS:
                0021 C     myThid       :: my Thread Id number
                0022       INTEGER myThid
                0023 CEOP
                0024
                0025 #ifdef ALLOW_BLING
e0f9a7ba0b Matt*0026 c     CHARACTER*(MAX_LEN_MBUF) msgBuf
c0d1c06c15 Matt*0027
                0028 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0029
                0030       _BEGIN_MASTER(myThid)
                0031
                0032 C set up coefficients for DIC chemistry
                0033 C define Schmidt no. coefficients for CO2
                0034       sca1 = 2073.1 _d 0
                0035       sca2 = -125.62 _d 0
                0036       sca3 =    3.6276 _d 0
                0037       sca4 =   -0.043219 _d 0
                0038 C define Schmidt no. coefficients for O2
                0039 C based on Keeling et al [GBC, 12, 141, (1998)]
                0040 CAV new Wanninkhof numbers
                0041       sox1 = 1929.7 _d 0
                0042       sox2 = -117.46 _d 0
                0043       sox3 =    3.116 _d 0
                0044       sox4 =   -0.0306 _d 0
                0045
                0046 C coefficients for determining saturation O2
                0047       oA0=  2.00907 _d 0
                0048       oA1=  3.22014 _d 0
                0049       oA2=  4.05010 _d 0
                0050       oA3=  4.94457 _d 0
                0051       oA4= -2.56847 _d -1
                0052       oA5=  3.88767 _d 0
                0053       oB0= -6.24523 _d -3
                0054       oB1= -7.37614 _d -3
                0055       oB2= -1.03410 _d -2
                0056       oB3= -8.17083 _d -3
                0057       oC0= -4.88682 _d -7
                0058
e0f9a7ba0b Matt*0059 #ifdef ALLOW_EXF
079948e6a6 Matt*0060       apco2mask = 'c'
                0061       IF ( useSEAICE ) THEN
                0062        apco2mask = ' '
                0063       ENDIF
                0064 #endif
                0065
c0d1c06c15 Matt*0066 C     Set other constant/flag
                0067
                0068       _END_MASTER(myThid)
                0069
                0070 #ifdef ALLOW_DIAGNOSTICS
                0071       IF ( useDIAGNOSTICS ) THEN
                0072         CALL BLING_DIAGNOSTICS_INIT( myThid )
                0073       ENDIF
                0074 #endif /* ALLOW_DIAGNOSTICS */
                0075
                0076 #endif /* ALLOW_BLING */
                0077
                0078       RETURN
                0079       END