1 |
C $Header: $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "BLING_OPTIONS.h" |
5 |
|
6 |
CBOP |
7 |
subroutine BLING_INIT_FIXED( myThid ) |
8 |
|
9 |
C ========================================================== |
10 |
C | subroutine bling_init_fixed |
11 |
C | o Initialize fixed quantities |
12 |
C ========================================================== |
13 |
|
14 |
implicit none |
15 |
|
16 |
C === Global variables === |
17 |
#include "SIZE.h" |
18 |
#include "EEPARAMS.h" |
19 |
#include "PARAMS.h" |
20 |
#include "GRID.h" |
21 |
#include "BLING_VARS.h" |
22 |
|
23 |
C !INPUT PARAMETERS: |
24 |
C myThid :: my Thread Id number |
25 |
INTEGER myThid |
26 |
CEOP |
27 |
|
28 |
#ifdef ALLOW_BLING |
29 |
INTEGER k |
30 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
31 |
INTEGER iUnit |
32 |
|
33 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
34 |
|
35 |
_BEGIN_MASTER(myThid) |
36 |
|
37 |
C set up coefficients for DIC chemistry |
38 |
C define Schmidt no. coefficients for CO2 |
39 |
sca1 = 2073.1 _d 0 |
40 |
sca2 = -125.62 _d 0 |
41 |
sca3 = 3.6276 _d 0 |
42 |
sca4 = -0.043219 _d 0 |
43 |
C define Schmidt no. coefficients for O2 |
44 |
C based on Keeling et al [GBC, 12, 141, (1998)] |
45 |
CAV new Wanninkhof numbers |
46 |
sox1 = 1929.7 _d 0 |
47 |
sox2 = -117.46 _d 0 |
48 |
sox3 = 3.116 _d 0 |
49 |
sox4 = -0.0306 _d 0 |
50 |
|
51 |
C coefficients for determining saturation O2 |
52 |
oA0= 2.00907 _d 0 |
53 |
oA1= 3.22014 _d 0 |
54 |
oA2= 4.05010 _d 0 |
55 |
oA3= 4.94457 _d 0 |
56 |
oA4= -2.56847 _d -1 |
57 |
oA5= 3.88767 _d 0 |
58 |
oB0= -6.24523 _d -3 |
59 |
oB1= -7.37614 _d -3 |
60 |
oB2= -1.03410 _d -2 |
61 |
oB3= -8.17083 _d -3 |
62 |
oC0= -4.88682 _d -7 |
63 |
|
64 |
C Set other constant/flag |
65 |
|
66 |
_END_MASTER(myThid) |
67 |
|
68 |
#ifdef ALLOW_DIAGNOSTICS |
69 |
IF ( useDIAGNOSTICS ) THEN |
70 |
CALL BLING_DIAGNOSTICS_INIT( myThid ) |
71 |
ENDIF |
72 |
#endif /* ALLOW_DIAGNOSTICS */ |
73 |
|
74 |
#endif /* ALLOW_BLING */ |
75 |
|
76 |
RETURN |
77 |
END |