source: trunk/tools/density_binning/binning_neutral_and_co/fsbeta.pro @ 163

Last change on this file since 163 was 163, checked in by pinsard, 15 years ago

going on modify unformal header to idldoc 2. header syntax

File size: 1.2 KB
Line 
1;+
2;
3; compute saline expension coefficients
4; used for the computation of the the local Brunt-Vaisala frequency
5;       
6; @param PFT {in} {required}
7; potential temperature in degrees celcius
8;
9; @param PFS {in} {required}
10; salinity anomaly (s-35) in psu
11;
12; @param PFH {in} {required}
13; depth in meters
14;
15; @returns
16; saline expension coefficients
17;
18;   References :
19;   -----------
20;       McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987.
21;
22; @history
23; Original  : 94-07 (G. Madec, M. Imbard)
24; Additions : 97-07 (G. Madec) introduction of statement functions
25; Idl version : 98-12 (M. Levy)
26;
27; @version
28; $Id$
29;
30;-
31FUNCTION fsbeta,  pft,  pfs,  pfh
32      z = $
33       ( ( -0.415613e-09 * pft + 0.555579e-07 ) * pft $
34                               - 0.301985e-05 ) * pft $
35                               + 0.785567e-03 $
36      +( 0.515032e-08 * pfs + 0.788212e-08 * pft - 0.356603e-06 ) * pfs $
37      +(  (   0.121551e-17 * pfh $
38            - 0.602281e-15 * pfs $
39            - 0.175379e-14 * pft + 0.176621e-12 ) * pfh $
40                                 + 0.408195e-10   * pfs $
41         +( - 0.213127e-11 * pft + 0.192867e-09 ) * pft $ 
42                                 - 0.121555e-07 ) * pfh
43      return,  z
44      end
Note: See TracBrowser for help on using the repository browser.