======================================================== | | | HYDROLIGHT Version 5.1.0 (MAIN CODE) | | Copyright (c) 2009 by Curtis D. Mobley | | | | This copy of HE5 is licensed solely to | | Curtis Mobley, Sequoia Scientific Inc. | | | | S/N XXXXX-XXXXX-XXXXX | | | | THIS PROGRAM IS EXPERIMENTAL AND IS LICENSED "AS IS" | | WITHOUT REPRESENTATION OF WARRANTY OF ANY KIND, | | EITHER EXPRESS OR IMPLIED. THE ENTIRE RISK AS TO | | THE QUALITY AND PERFORMANCE OF HYDROLIGHT IS WITH | | THE USER. HYDROLIGHT IS NOT FAULT TOLERANT. | | | ======================================================== RUN TITLE: HE5.1: User's Guide Example 1 (12/03/2009 15:49) Run started on 12/03/2009 at 15:49:39.9 The radiance distribution is computed on the following (Theta,Phi,Depth,Wavelength) grid: I Theta Mu J Phi K zeta z L Lambda (deg) (deg) (m) (nm) 1 87.500 0.0436 1 0.00 1 0.000 1 550.00 2 80.000 0.1730 2 15.00 2 0.050 3 70.000 0.3407 3 30.00 3 0.100 4 60.000 0.4981 4 45.00 4 0.150 5 50.000 0.6403 5 60.00 5 0.200 6 40.000 0.7631 6 75.00 6 0.250 7 30.000 0.8627 7 90.00 7 0.500 8 20.000 0.9361 8 105.00 8 0.550 9 10.000 0.9811 9 120.00 9 1.000 10 0.000 0.9981 10 135.00 10 1.050 11 150.00 11 2.000 12 165.00 12 2.050 13 180.00 13 3.000 14 195.00 14 3.050 15 210.00 15 5.000 16 225.00 16 5.050 17 240.00 17 10.000 18 255.00 18 10.050 19 270.00 19 15.000 20 285.00 20 15.050 21 300.00 22 315.00 23 330.00 24 345.00 ***** BEGIN IOP MODEL SPECIFICATIONS ***** The IOP routine "abconst" is being used. Constant absorption and scattering coefficients are used: The absorption coefficient is a = 0.1000 1/m The scattering coefficient is b = 0.3000 1/m The phase function for component 1 will be chosen to have a bb/b ratio of 0.0230 ***** END IOP MODEL SPECIFICATIONS ***** ***** BEGIN INELASTIC SCATTER AND INTERNAL SOURCE SPECIFICATIONS ***** Chlorophyll fluorescence is not included in this run CDOM fluorescence is not included in this run Raman scattering is not included in this run ***** END INELASTIC SCATTER AND INTERNAL SOURCE SPECIFICATIONS ***** ***** BEGIN ATMOSPHERIC MODEL SPECIFICATIONS ***** Normalized sky radiances are computed using sky model "COSRAD" (COSine RADiance model) COSRAD parameters: Rsky = ratio of sky to total plane irradiance = 0.300 C = cosine parameter = 0.000 Edtotl = total (sun + sky) Ed in W/(m^2 nm) = 1.000 The sun is placed in quad Q(r,s) = Q( 7, 1) centered at (theta, phi) = (30.000, 0.000) Sky irradiances are computed using sky irradiance model "COSIRRAD" (COSine IRRADiances; for use with the COSRAD model) diffuse (sky) spectral plane irradiance Ed(diffuse) = 3.000E-01 W/(m^2 nm) direct (sun) spectral plane irradiance Ed(direct) = 7.000E-01 W/(m^2 nm) total (sun + sky) spectral plane irradiance Ed(total) = 1.000E+00 W/(m^2 nm) ratio of sky to total plane irradiance, Rsky = 0.300 ***** END ATMOSPHERIC MODEL SPECIFICATIONS ***** ***** BEGIN AIR-WATER SURFACE SPECIFICATIONS ***** The air-water surface is for a wind speed of 0.0 m/s Index of refraction is constant for all wavelengths with the value: 1.340 ***** END AIR-WATER SURFACE SPECIFICATIONS ***** ***** BEGIN BOTTOM BOUNDARY SPECIFICATIONS ***** The bottom boundary is an infinitely deep, homogeneous water body below optical depth 15.00 ***** END BOTTOM BOUNDARY SPECIFICATIONS ***** Initialization of Hydrolight Completed in 0.1 sec. * * * * * Output for wavelength band 1 (549.5 to 550.5 nm; nominal wavelength = 550.0 nm) * * * * * Index of refraction at this wavelength is 1.340 Absorption Coefficients of Individual Components at 550.0 nm iz Opt Depth Geo Depth Comp 1 Total (m) (1/m) (1/m) 1 0.000 0.000 0.1000 0.1000 3 0.100 0.250 0.1000 0.1000 5 0.200 0.500 0.1000 0.1000 7 0.500 1.250 0.1000 0.1000 9 1.000 2.500 0.1000 0.1000 11 2.000 5.000 0.1000 0.1000 13 3.000 7.500 0.1000 0.1000 15 5.000 12.500 0.1000 0.1000 17 10.000 25.000 0.1000 0.1000 19 15.000 37.500 0.1000 0.1000 Scattering Coefficients of Individual Components at 550.0 nm iz Opt Depth Geo Depth Comp 1 Total (m) (1/m) (1/m) 1 0.000 0.000 0.3000 0.3000 3 0.100 0.250 0.3000 0.3000 5 0.200 0.500 0.3000 0.3000 7 0.500 1.250 0.3000 0.3000 9 1.000 2.500 0.3000 0.3000 11 2.000 5.000 0.3000 0.3000 13 3.000 7.500 0.3000 0.3000 15 5.000 12.500 0.3000 0.3000 17 10.000 25.000 0.3000 0.3000 19 15.000 37.500 0.3000 0.3000 Backscattering Coefficients of Individual Components at 550.0 nm (computed from the scattering coefficients and the phase functions) iz Opt Depth Geo Depth Comp 1 Total (m) (1/m) (1/m) 1 0.000 0.000 0.00690 0.00690 3 0.100 0.250 0.00690 0.00690 5 0.200 0.500 0.00690 0.00690 7 0.500 1.250 0.00690 0.00690 9 1.000 2.500 0.00690 0.00690 11 2.000 5.000 0.00690 0.00690 13 3.000 7.500 0.00690 0.00690 15 5.000 12.500 0.00690 0.00690 17 10.000 25.000 0.00690 0.00690 19 15.000 37.500 0.00690 0.00690 Summary of Inherent Optical Properties at 550.0 nm iz Opt Depth Geo Depth total a total b total c albedo total bb total bb/b (m) (1/m) (1/m) (1/m) (1/m) 1 0.000 0.000 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 3 0.100 0.250 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 5 0.200 0.500 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 7 0.500 1.250 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 9 1.000 2.500 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 11 2.000 5.000 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 13 3.000 7.500 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 15 5.000 12.500 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 17 10.000 25.000 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 19 15.000 37.500 0.1000 0.3000 0.4000 0.7500 0.00690 0.02300 The shape of the asymptotic radiance distribution is given by u theta Linf(theta) 10 0.00 1.0000E+00 9 10.00 9.2619E-01 8 20.00 7.5619E-01 7 30.00 5.5095E-01 6 40.00 3.6746E-01 5 50.00 2.3106E-01 4 60.00 1.4119E-01 3 70.00 8.6207E-02 2 80.00 5.3960E-02 1 87.50 3.8704E-02 1 92.50 3.0997E-02 2 100.00 2.3388E-02 3 110.00 1.6924E-02 4 120.00 1.2937E-02 5 130.00 1.0475E-02 6 140.00 8.9457E-03 7 150.00 8.0011E-03 8 160.00 7.4376E-03 9 170.00 7.1390E-03 10 180.00 7.0350E-03 Asymptotic AOP values are MUu(inf) = 0.3708 MUd(inf) = 0.7480 MU(inf) = 0.6826 R(inf) = 0.0308 Lu(inf)/Ed(inf) = 0.006109 (1/sr) k(inf) = 0.1465 (1/m) absorp(check) = 0.1000 = k(inf)*MU(inf) Irradiances (units of W/m^2 nm), Mean Cosines (Mubars), and Irradiance Reflectance at 550.0 nm iz zeta z(m) Eou Eod Eo Eu Ed MUBARu MUBARd MUBAR R = Eu/Ed in air 1.4211E-01 1.4114E+00 1.5535E+00 4.7106E-02 9.9999E-01 0.3315 0.7085 0.6134 4.7106E-02 1 0.00 0.00 6.0096E-02 1.1186E+00 1.1787E+00 2.4716E-02 9.7917E-01 0.4113 0.8754 0.8098 2.5242E-02 3 0.10 0.25 5.9136E-02 1.0921E+00 1.1513E+00 2.4168E-02 9.4949E-01 0.4087 0.8694 0.8037 2.5454E-02 5 0.20 0.50 5.8078E-02 1.0653E+00 1.1234E+00 2.3610E-02 9.2050E-01 0.4065 0.8641 0.7984 2.5649E-02 7 0.50 1.25 5.4611E-02 9.8486E-01 1.0395E+00 2.1918E-02 8.3771E-01 0.4013 0.8506 0.7848 2.6164E-02 9 1.00 2.50 4.8483E-02 8.5632E-01 9.0481E-01 1.9164E-02 7.1353E-01 0.3953 0.8333 0.7674 2.6858E-02 11 2.00 5.00 3.6867E-02 6.3358E-01 6.7044E-01 1.4293E-02 5.1283E-01 0.3877 0.8094 0.7436 2.7871E-02 13 3.00 7.50 2.7253E-02 4.6026E-01 4.8751E-01 1.0439E-02 3.6527E-01 0.3830 0.7936 0.7278 2.8579E-02 15 5.00 12.50 1.4215E-02 2.3514E-01 2.4935E-01 5.3687E-03 1.8207E-01 0.3777 0.7743 0.7087 2.9486E-02 17 10.00 25.00 2.4769E-03 4.0163E-02 4.2640E-02 9.2303E-04 3.0334E-02 0.3727 0.7553 0.6898 3.0429E-02 19 15.00 37.50 4.0622E-04 6.5514E-03 6.9576E-03 1.5083E-04 4.9140E-03 0.3713 0.7501 0.6846 3.0695E-02 K-functions (units of 1/meter) at 550.0 nm (these are accurate local values only when zupper and zlower are closely spaced) zupper zlower z Kou(z) Kod(z) Ko(z) Ku(z) Kd(z) Knet(z) KLu(z) 0.000 0.125 0.062 0.06202 0.09456 0.09290 0.08860 0.12285 0.12373 0.11551 0.250 0.375 0.312 0.07061 0.09876 0.09732 0.09267 0.12382 0.12463 0.11695 0.500 0.625 0.562 0.07651 0.10175 0.10044 0.09588 0.12469 0.12545 0.11821 1.250 1.375 1.312 0.08884 0.10829 0.10727 0.10328 0.12694 0.12758 0.12142 2.500 2.625 2.562 0.10181 0.11572 0.11497 0.11181 0.12993 0.13043 0.12553 5.000 5.125 5.062 0.11649 0.12492 0.12446 0.12236 0.13423 0.13457 0.13121 7.500 7.625 7.562 0.12500 0.13066 0.13034 0.12887 0.13721 0.13745 0.13504 12.500 12.625 12.562 0.13448 0.13742 0.13725 0.13646 0.14101 0.14115 0.13979 25.000 25.125 25.062 0.14321 0.14397 0.14393 0.14372 0.14494 0.14498 0.14461 37.500 37.625 37.562 0.14558 0.14579 0.14578 0.14572 0.14606 0.14607 0.14597 Radiances Just Above the Water Surface (at z = a) (units of W/m^2 sr nm) at 550.0 nm [(Theta,Phi) are the directions of photon travel. Theta is measured from 0 at the zenith. Solar photons travel in the Phi = 180 direction (the sun is located at Phi = 0). (Thus Phi = 45 represents the 135 degree viewing angle for minimizing sun glitter.)] I J Theta Phi Phi-view L(sky) L(water-lv) L(refl-sky) L(total up) Rrs rho (zenith) (Lrefl/Lsky) 10 1 0.0 0.0 0.0 9.54930E-02 3.09395E-03 2.00939E-03 5.10334E-03 3.09396E-03 0.021042 9 1 10.0 0.0 180.0 9.54946E-02 3.13124E-03 2.01644E-03 5.14769E-03 3.13126E-03 0.021116 9 4 10.0 45.0 135.0 9.54922E-02 3.12182E-03 2.01639E-03 5.13821E-03 3.12183E-03 0.021116 9 7 10.0 90.0 90.0 9.54934E-02 3.09921E-03 2.01645E-03 5.11566E-03 3.09923E-03 0.021116 8 1 20.0 0.0 180.0 9.54945E-02 3.18222E-03 2.03822E-03 5.22044E-03 3.18224E-03 0.021344 8 4 20.0 45.0 135.0 9.54921E-02 3.15854E-03 2.03814E-03 5.19668E-03 3.15856E-03 0.021344 8 7 20.0 90.0 90.0 9.54933E-02 3.13594E-03 2.03821E-03 5.17416E-03 3.13596E-03 0.021344 7 1 30.0 0.0 180.0 9.54947E-02 3.27162E-03 2.13209E-03 5.40371E-03 3.27164E-03 0.022327 7 4 30.0 45.0 135.0 9.55107E-02 3.25817E-03 2.13226E-03 5.39042E-03 3.25819E-03 0.022325 7 7 30.0 90.0 90.0 9.53811E-02 3.24339E-03 2.12952E-03 5.37291E-03 3.24341E-03 0.022326 6 1 40.0 0.0 180.0 9.54943E-02 3.32972E-03 2.45013E-03 5.77985E-03 3.32974E-03 0.025657 6 4 40.0 45.0 135.0 9.54919E-02 3.30991E-03 2.45029E-03 5.76019E-03 3.30993E-03 0.025660 6 7 40.0 90.0 90.0 9.54931E-02 3.28686E-03 2.45027E-03 5.73714E-03 3.28688E-03 0.025659 5 1 50.0 0.0 180.0 9.54942E-02 3.42624E-03 3.38837E-03 6.81461E-03 3.42626E-03 0.035482 5 4 50.0 45.0 135.0 9.54918E-02 3.40019E-03 3.39022E-03 6.79041E-03 3.40021E-03 0.035503 5 7 50.0 90.0 90.0 9.54931E-02 3.42239E-03 3.39010E-03 6.81250E-03 3.42241E-03 0.035501 4 1 60.0 0.0 180.0 9.54944E-02 3.46122E-03 6.03601E-03 9.49723E-03 3.46124E-03 0.063208 4 4 60.0 45.0 135.0 9.54920E-02 3.45646E-03 6.03580E-03 9.49225E-03 3.45648E-03 0.063207 4 7 60.0 90.0 90.0 9.54932E-02 3.45596E-03 6.03511E-03 9.49106E-03 3.45598E-03 0.063199 3 1 70.0 0.0 180.0 9.54945E-02 3.33346E-03 1.35093E-02 1.68428E-02 3.33348E-03 0.141467 3 4 70.0 45.0 135.0 9.54921E-02 3.30464E-03 1.34921E-02 1.67967E-02 3.30466E-03 0.141290 3 7 70.0 90.0 90.0 9.54934E-02 3.33722E-03 1.35060E-02 1.68432E-02 3.33724E-03 0.141434 2 1 80.0 0.0 180.0 9.54943E-02 2.50502E-03 3.51022E-02 3.76073E-02 2.50503E-03 0.367585 2 4 80.0 45.0 135.0 9.54919E-02 2.49900E-03 3.50780E-02 3.75770E-02 2.49902E-03 0.367340 2 7 80.0 90.0 90.0 9.54931E-02 2.55476E-03 3.51045E-02 3.76593E-02 2.55478E-03 0.367613 1 1 87.5 0.0 180.0 9.54945E-02 1.13604E-03 7.38274E-02 7.49634E-02 1.13605E-03 0.773106 1 4 87.5 45.0 135.0 9.54921E-02 1.17381E-03 7.38290E-02 7.50029E-02 1.17382E-03 0.773143 1 7 87.5 90.0 90.0 9.54933E-02 1.18757E-03 7.37753E-02 7.49629E-02 1.18758E-03 0.772570 Selected Radiances (units of W/m^2 sr nm) and Radiance-Irradiance Ratios at 550.0 nm [Here (theta,phi) are directions of photon travel relative to nadir (solar photons travel along phi=180).] iz zeta z Lu(z) Ld(z) Lh(z,phi=0) Lh(z,phi=90) Lh(z,phi=180) Lu/Ed Q = Eu/Lu Lw(z) Rrs = Lw/Ed (m) (theta=180) (theta=0) (theta=90) (theta=90) (theta=90) (1/sr) (sr) (theta=180) (1/sr) in air 5.103E-03 9.549E-02 8.523E-02 8.523E-02 8.531E-02 5.103E-03 9.230E+00 3.094E-03 3.094E-03 1 0.00 0.00 5.689E-03 1.672E-01 1.209E-02 1.428E-02 2.031E-02 5.810E-03 4.344E+00 3 0.10 0.25 5.527E-03 1.749E-01 1.235E-02 1.471E-02 2.126E-02 5.821E-03 4.373E+00 5 0.20 0.50 5.367E-03 1.819E-01 1.240E-02 1.485E-02 2.163E-02 5.831E-03 4.399E+00 7 0.50 1.25 4.907E-03 1.984E-01 1.217E-02 1.472E-02 2.172E-02 5.857E-03 4.467E+00 9 1.00 2.50 4.205E-03 2.134E-01 1.128E-02 1.379E-02 2.056E-02 5.894E-03 4.557E+00 11 2.00 5.00 3.050E-03 2.105E-01 9.023E-03 1.115E-02 1.668E-02 5.948E-03 4.685E+00 13 3.00 7.50 2.187E-03 1.842E-01 6.887E-03 8.549E-03 1.271E-02 5.987E-03 4.774E+00 15 5.00 12.50 1.099E-03 1.172E-01 3.768E-03 4.666E-03 6.743E-03 6.037E-03 4.885E+00 17 10.00 25.00 1.847E-04 2.451E-02 7.240E-04 8.616E-04 1.121E-03 6.088E-03 4.998E+00 19 15.00 37.50 2.999E-05 4.197E-03 1.282E-04 1.451E-04 1.717E-04 6.103E-03 5.029E+00 Waveband 1 of 1 completed in 3.7 sec. Energy Eo as a function of depth and wavelength Units are W/m^2 [values are integrated over the wavelength bands] depth 550.0 in air 1.553E+00 0.0 1.179E+00 0.2 1.151E+00 0.5 1.123E+00 1.2 1.039E+00 2.5 9.048E-01 5.0 6.704E-01 7.5 4.875E-01 12.5 2.494E-01 25.0 4.264E-02 37.5 6.958E-03 Quantum Eo as a function of depth and wavelength Units are (micromol photons)/(m^2 s) [values are integrated over the wavelength bands] depth 550.0 in air 7.141E+00 0.0 5.418E+00 0.2 5.292E+00 0.5 5.164E+00 1.2 4.778E+00 2.5 4.159E+00 5.0 3.082E+00 7.5 2.241E+00 12.5 1.146E+00 25.0 1.960E-01 37.5 3.198E-02 PAR was not calculated for this run (PAR is meaningful only for runs covering the near-UV and visible spectrum) PAR range currently set to 400.0 to 700.0 nm CIE chromaticity coordinates and Secchi depth are not computed since the run did not cover at least 400-700 nm. Run completed on 12/03/2009 at 15:49:43.6 Total (wall clock) run time = 3.8 sec