/[MITgcm]/MITgcm_contrib/bling/pkg/bling_init_fixed.F
ViewVC logotype

Annotation of /MITgcm_contrib/bling/pkg/bling_init_fixed.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (hide annotations) (download)
Fri May 23 17:33:43 2014 UTC (11 years, 2 months ago) by mmazloff
Branch: MAIN
Adding package BLING

1 mmazloff 1.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     sox1 = 1638.0 _d 0
46     sox2 = -81.83 _d 0
47     sox3 = 1.483 _d 0
48     sox4 = -0.008004 _d 0
49    
50     C coefficients for determining saturation O2
51     oA0= 2.00907 _d 0
52     oA1= 3.22014 _d 0
53     oA2= 4.05010 _d 0
54     oA3= 4.94457 _d 0
55     oA4= -2.56847 _d -1
56     oA5= 3.88767 _d 0
57     oB0= -6.24523 _d -3
58     oB1= -7.37614 _d -3
59     oB2= -1.03410 _d -2
60     oB3= -8.17083 _d -3
61     oC0= -4.88682 _d -7
62    
63     C Set other constant/flag
64    
65     _END_MASTER(myThid)
66    
67     #ifdef ALLOW_DIAGNOSTICS
68     IF ( useDIAGNOSTICS ) THEN
69     CALL BLING_DIAGNOSTICS_INIT( myThid )
70     ENDIF
71     #endif /* ALLOW_DIAGNOSTICS */
72    
73     #endif /* ALLOW_BLING */
74    
75     RETURN
76     END

  ViewVC Help
Powered by ViewVC 1.1.22