-
Notifications
You must be signed in to change notification settings - Fork 0
/
dncoef.h
120 lines (120 loc) · 5.28 KB
/
dncoef.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
c====================== include file "dncoef.h" ==========================
c
c
c normalized temperatures, salinities and coefficients
c generated by program "eqstat" which fits 3rd order polynomials
c to the equation of state for each model level.
c
data to /
& 13.4982640,13.4947868,13.4913025,13.4878112,13.4843128,
& 13.4781048, 8.4728624, 8.4523506, 5.9326036, 3.9075315,
& 2.8688491, 2.8034689, 2.7186610, 2.6141484, 2.9727823/
c
data so /
& -0.0022500,-0.0022500,-0.0022500,-0.0022500,-0.0022500,
& -0.0022500, 0.0001500, 0.0001500,-0.0001000,-0.0002500,
& -0.0002000,-0.0002000,-0.0002000,-0.0002000,-0.0002000/
c
data (c( 1,n),n=1,9)/
& -.2017392E-03,0.7730067E+00,-.4917107E-05,-.2021351E-02,
& 0.3169250E+00,0.3607234E-07,0.3774674E-02,0.3602436E-04,
& 0.1609359E+01/
c
data (c( 2,n),n=1,9)/
& -.2021980E-03,0.7728398E+00,-.4909621E-05,-.2018858E-02,
& 0.3165936E+00,0.3597763E-07,0.3766929E-02,0.3598190E-04,
& 0.1608967E+01/
c
data (c( 3,n),n=1,9)/
& -.2026566E-03,0.7726730E+00,-.4902137E-05,-.2016363E-02,
& 0.3162610E+00,0.3588284E-07,0.3759186E-02,0.3593946E-04,
& 0.1608578E+01/
c
data (c( 4,n),n=1,9)/
& -.2031148E-03,0.7725063E+00,-.4894655E-05,-.2013868E-02,
& 0.3159271E+00,0.3578796E-07,0.3751446E-02,0.3589704E-04,
& 0.1608193E+01/
c
data (c( 5,n),n=1,9)/
& -.2035729E-03,0.7723398E+00,-.4887175E-05,-.2011371E-02,
& 0.3155921E+00,0.3569301E-07,0.3743707E-02,0.3585463E-04,
& 0.1607810E+01/
c
data (c( 6,n),n=1,9)/
& -.2043827E-03,0.7720452E+00,-.4873944E-05,-.2006951E-02,
& 0.3149963E+00,0.3552478E-07,0.3730017E-02,0.3577963E-04,
& 0.1607142E+01/
c
data (c( 7,n),n=1,9)/
& -.1615928E-03,0.7826355E+00,-.5276493E-05,-.2317934E-02,
& 0.4652976E+00,0.4412553E-07,0.4581339E-02,0.4310578E-04,
& 0.1509747E+01/
c
data (c( 8,n),n=1,9)/
& -.1656607E-03,0.7812812E+00,-.5209958E-05,-.2295928E-02,
& 0.4331218E+00,0.4334451E-07,0.4501743E-02,0.4272019E-04,
& 0.1506151E+01/
c
data (c( 9,n),n=1,9)/
& -.1456635E-03,0.7848106E+00,-.5425527E-05,-.2463241E-02,
& 0.2740912E+00,0.4862742E-07,0.4858158E-02,0.4672413E-04,
& 0.1563903E+01/
c
data (c( 10,n),n=1,9)/
& -.1340739E-03,0.7864566E+00,-.5536224E-05,-.2579984E-02,
& -.8651291E-01,0.5236292E-07,0.4997093E-02,0.4959251E-04,
& 0.1605778E+01/
c
data (c( 11,n),n=1,9)/
& -.1376237E-03,0.7846464E+00,-.5447229E-05,-.2594208E-02,
& -.1002929E+01,0.5256849E-07,0.4820744E-02,0.5014700E-04,
& 0.1611559E+01/
c
data (c( 12,n),n=1,9)/
& -.1550969E-03,0.7792967E+00,-.5154945E-05,-.2494927E-02,
& -.4194413E+01,0.4903922E-07,0.4417669E-02,0.4833834E-04,
& 0.1597625E+01/
c
data (c( 13,n),n=1,9)/
& -.1749168E-03,0.7731871E+00,-.4820115E-05,-.2380375E-02,
& -.6548572E+01,0.4484159E-07,0.3961840E-02,0.4625855E-04,
& 0.1584458E+01/
c
data (c( 14,n),n=1,9)/
& -.1962247E-03,0.7665555E+00,-.4456675E-05,-.2254947E-02,
& -.9547533E+01,0.4009664E-07,0.3474258E-02,0.4398563E-04,
& 0.1573560E+01/
c
data (c( 15,n),n=1,9)/
& -.2218872E-03,0.7586723E+00,-.4028944E-05,-.2083144E-02,
& -.9816443E+01,0.3374213E-07,0.2987149E-02,0.4082538E-04,
& 0.1560522E+01/
c
c the above coefficients were calculated using program "eqstat"
c (employing the UNESCO equation of state)
c and are valid for the following depths and T and S ranges
c k depth tmin tmax smin smax
c 1 12.50e2 -2.000 29.000 28.5000 37.0000
c 2 37.50e2 -2.000 29.000 28.5000 37.0000
c 3 62.50e2 -2.000 29.000 28.5000 37.0000
c 4 87.50e2 -2.000 29.000 28.5000 37.0000
c 5 112.50e2 -2.000 29.000 28.5000 37.0000
c 6 156.74e2 -2.000 29.000 28.5000 37.0000
c 7 257.15e2 -2.000 19.000 33.7000 36.6000
c 8 446.09e2 -2.000 19.000 33.7000 36.6000
c 9 748.76e2 -2.000 14.000 34.0000 35.8000
c 10 1181.16e2 -1.000 9.000 34.2000 35.3000
c 11 1748.76e2 -1.000 7.000 34.5000 35.1000
c 12 2446.09e2 -1.000 7.000 34.6000 35.0000
c 13 3257.15e2 -1.000 7.000 34.6000 35.0000
c 14 4156.74e2 -1.000 7.000 34.6000 35.0000
c 15 5112.50e2 0.000 7.000 34.6000 35.0000
c
c the 3rd order polynomial will return density departures [gm/cm**3] as
c a function of potential temperature [deg C] & salinity [model units]
c
c k level reference densities (in sigma units):
c 24.6324 24.7442 24.8560 24.9678 25.0794 25.2769 28.5082 29.3623
c 30.9232 33.0457 35.7885 38.9109 42.4849 46.3783 50.3290
c
data cksumzt / 0.1982364E+07/