Last change
on this file since 122 was
2,
checked in by post_it, 17 years ago
|
Initial import from ~/POST_IT/
|
File size:
1.1 KB
|
Line | |
---|
1 | FUNCTION cpsw, T, S, p |
---|
2 | ; |
---|
3 | ; Specific heat of sea water (J/KG C) |
---|
4 | ; |
---|
5 | CP1 = 0. |
---|
6 | CP2 = 0. |
---|
7 | SR=SQRT(ABS(S)) |
---|
8 | ; SPECIFIC HEAT CP0 FOR P=0 (MILLERO ET AL. 1973) |
---|
9 | A = (-1.38E-3*T+0.10727)*T-7.644 |
---|
10 | B = (5.35E-5*T-4.08E-3)*T+0.177 |
---|
11 | C = (((2.093236E-5*T-2.654387E-3)*T+0.1412855)*T-3.720283)*T+4217.4 |
---|
12 | CP0 = (B*SR + A) * S + C |
---|
13 | ; CP1 PRESSURE AND TEMPERATURE TERMS FOR S = 0 |
---|
14 | A = (((1.7168E-8*T+2.0357E-6)*T-3.13885E-4)*T+1.45747E-2)*T-0.49592 |
---|
15 | B = (((2.2956E-11*T-4.0027E-9)*T+2.87533E-7)*T-1.08645E-5)*T+2.4931E-4 |
---|
16 | C = ((6.136E-13*T-6.5637E-11)*T+2.6380E-9)*T-5.422E-8 |
---|
17 | CP1 = ((C*P+B)*P+A)*P |
---|
18 | ; CP2 PRESSURE AND TEMPERATURE TERMS FOR S > 0 |
---|
19 | A = (((-2.9179E-10*T+2.5941E-8)*T+9.802E-7)*T-1.28315E-4)*T+4.9247E-3 |
---|
20 | B = (3.122E-8*T-1.517E-6)*T-1.2331E-4 |
---|
21 | A = (A+B*SR)*S |
---|
22 | B = ((1.8448E-11*T-2.3905E-9)*T+1.17054E-7)*T-2.9558E-6 |
---|
23 | B = (B+9.971E-8*SR)*S |
---|
24 | C = (3.513E-13*T-1.7682E-11)*T+5.540E-10 |
---|
25 | C = (C-1.4300E-12*T*SR)*S |
---|
26 | CP2 = ((C*P+B)*P+A)*P |
---|
27 | ; SPECIFIC HEAT RETURN |
---|
28 | cp = CP0 + CP1 + CP2 |
---|
29 | return, cp |
---|
30 | END |
---|
Note: See
TracBrowser
for help on using the repository browser.