• API
• FAQ
• Tools
• Trends
• Archive
SHARE
TWEET

# Untitled

a guest Jan 14th, 2014 46 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. <?php
2.
3. function salinity(\$measuredConductivity, \$measuredTemperature, \$measuredDepth) {
4.
5. // PHP version by Ben Ostrowsky
6. // Based on Excel formulas by James Douglass
7. //   http://jimbodouglass.blogspot.com/2010/11/conductivity-to-salinity-conversion-for.html
8. // Based on Fofonoff, N.P. and Millard, R.C., Jr. (1983). Algorithms for computation
9. //   of fundamental properties of seawater. http://hdl.handle.net/1912/2470
10. // Added code to compensate for depth.
11.
12.         \$referenceConductivity = 42900.0; // microsiemens per cm^2
13.
14.         \$conductivityRatio = \$measuredConductivity / \$referenceConductivity;
15.
16.         \$waterDensity = 1.03 * 10^3; // kilograms per m^3 (pure water would be 1.00 * 10^3)
17.         \$gravity = 9.8; // meters per s^2
18.         \$airPressure = 10.0; // decibars
19.
20.         \$computedPressure =  (\$waterDensity * \$gravity * \$measuredDepth / 10.0) + \$airPressure;
21.
22.         \$littlert = (
23.                     0.6766097
24.                 +   0.0200564     * \$measuredTemperature
25.                 +   0.0001104259  * (\$measuredTemperature^2)
26.                 + (-6.9698*10^-7) * (\$measuredTemperature^3)
27.                 + ( 1.0031*10^-9) * (\$measuredTemperature^4)
28.         );
29.
30.         \$Rp =
31.                 1 + (
32.                         (
33.                                 \$computedPressure * (
34.                                         2.070*10^-5   * (\$computedPressure^0)
35.                                         + (-6.370*10^-10) * (\$computedPressure^1)
36.                                         +   3.989*10^-15  * (\$computedPressure^2)
37.                                 )
38.
39.                                 /
40.
41.                                 (
42.                                         1
43.                                         + (3.426*10^-2) * (\$measuredTemperature^1)
44.                                         + (4.464*10^-4) * (\$measuredTemperature^2)
45.                                         + ( 4.215*10^-1 + (-3.107*10^-3) * \$measuredTemperature ) * \$conductivityRatio
46.                                 )
47.                         )
48.                 );
49.
50.         \$bigRt =  \$conductivityRatio / (\$littlert * \$Rp);
51.
52.         \$dS =
53.                 (
54.                         ( \$measuredTemperature - 15.0) / (1 + 0.0162*( \$measuredTemperature - 15.0 ) )
55.                 )
56.                 * (
57.                             0.0005
58.                         + (-0.0056) * (\$bigRt^0.5)
59.                         + (-0.0066) * (\$bigRt^1)
60.                         + (-0.0375) * (\$bigRt^1.5)
61.                         + ( 0.0636) * (\$bigRt^2)
62.                         + (-0.0144) * (\$bigRt^2.5)
63.                 );
64.
65.         \$computedSalinity =
66.                     0.008
67.                 + (-0.1692) * (\$bigRt^0.5)
68.                 +  25.3851  * (\$bigRt^1)
69.                 +  14.0941  * (\$bigRt^1.5)
70.                 + (-7.0261) * (\$bigRt^2)
71.                 +   2.7081  * (\$bigRt^2.5)
72.                 + \$dS;
73.
74.         return \$computedSalinity;
75.
76. }
77.
78. ?>
RAW Paste Data
Top