Last change
on this file since 68 was
2,
checked in by post_it, 17 years ago
|
Initial import from ~/POST_IT/
|
File size:
1.1 KB
|
Rev | Line | |
---|
[2] | 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.