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

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

conversion from IDL help header syntax to IDLDOC 2. header syntax

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