======================================================== | | | HYDROLIGHT Version 5.1.0 (MAIN CODE) | | Copyright (c) 2008 by Curtis D. Mobley | | | | This copy of HE5 is licensed solely to | | Student Version for Maine 2011 Summer Class | | | | S/N Maine_500_31July2011 | | | | 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: run2 for iop error exercise 2 wed class (07/20/2011 10:47) Run started on 07/20/2011 at 10:47:16.4 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 440.00 2 80.000 0.1730 2 15.00 2 0.010 3 70.000 0.3407 3 30.00 3 5.000 4 60.000 0.4981 4 45.00 4 5.010 5 50.000 0.6403 5 60.00 5 10.000 6 40.000 0.7631 6 75.00 6 10.010 7 30.000 0.8627 7 90.00 7 15.000 8 20.000 0.9361 8 105.00 8 15.010 9 10.000 0.9811 9 120.00 9 20.000 10 0.000 0.9981 10 135.00 10 20.010 11 150.00 11 25.000 12 165.00 12 25.010 13 180.00 13 30.000 14 195.00 14 30.010 15 210.00 15 35.000 16 225.00 16 35.010 17 240.00 17 40.000 18 255.00 18 40.010 19 270.00 19 45.000 20 285.00 20 45.010 21 300.00 21 50.000 22 315.00 22 50.010 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.2400 1/m The scattering coefficient is b = 1.0000 1/m The phase function for component 1 comes from a file titled Petzold "average particle" from L&W Table 3.10; standard quad partition; 7 Sept 98 ***** 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 "HCNRAD" (Harrison and Coombes Normalized RADiances) The cloud fraction is cloud = 0.0 (0 for a clear sky to 1 for a heavy overcast) The sun is placed in quad Q(r,s) = Q( 7, 1) centered at (theta, phi) = (30.000, 0.000) Diffuse and Direct Sky irradiances are computed using the "RADTRANX" (RADTRAN eXtended for 300-1000 nm) model The RADTRAN clear-sky irradiances are adjusted for cloudiness by the model of Kasten and Czeplak Parameter values used in the RADTRANX model: Solar zenith angle = 30.0 degrees Pressure = 29.9 in. mercury Day of year = not set (mean earth-sun distance used) Air mass type = 1.0 Relative Humidity = 80.0 % Precipitable water = 2.5 cm 24-hr wind speed = 5.0 m/sec Current wind speed = 5.0 m/sec Visibility = 15.0 km Total ozone = 300.0 Dobson units Irradiances are computed just above the sea surface Aerosol optical thickness at 550 nm = 0.261 ***** END ATMOSPHERIC MODEL SPECIFICATIONS ***** ***** BEGIN AIR-WATER SURFACE SPECIFICATIONS ***** The air-water surface is for a wind speed of 5.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 depth 50.00 m ***** END BOTTOM BOUNDARY SPECIFICATIONS ***** Initialization of Hydrolight Completed in 0.1 sec. * * * * * Output for wavelength band 1 (439.5 to 440.5 nm; nominal wavelength = 440.0 nm) * * * * * Index of refraction at this wavelength is 1.340 RADTRANX clear-sky spectral irradiances averaged over the band from 439.5 nm to 440.5 nm: Ed(diffuse) = 4.446E-01 Ed(direct) = 8.663E-01 Ed(total) = 1.311E+00 Absorption Coefficients of Individual Components at 440.0 nm iz Opt Depth Geo Depth Comp 1 Total (m) (1/m) (1/m) 1 0.000 0.000 0.2400 0.2400 3 6.200 5.000 0.2400 0.2400 5 12.400 10.000 0.2400 0.2400 7 18.600 15.000 0.2400 0.2400 9 24.800 20.000 0.2400 0.2400 11 31.000 25.000 0.2400 0.2400 13 37.200 30.000 0.2400 0.2400 15 43.400 35.000 0.2400 0.2400 17 49.600 40.000 0.2400 0.2400 19 55.800 45.000 0.2400 0.2400 21 62.000 50.000 0.2400 0.2400 Scattering Coefficients of Individual Components at 440.0 nm iz Opt Depth Geo Depth Comp 1 Total (m) (1/m) (1/m) 1 0.000 0.000 1.0000 1.0000 3 6.200 5.000 1.0000 1.0000 5 12.400 10.000 1.0000 1.0000 7 18.600 15.000 1.0000 1.0000 9 24.800 20.000 1.0000 1.0000 11 31.000 25.000 1.0000 1.0000 13 37.200 30.000 1.0000 1.0000 15 43.400 35.000 1.0000 1.0000 17 49.600 40.000 1.0000 1.0000 19 55.800 45.000 1.0000 1.0000 21 62.000 50.000 1.0000 1.0000 Backscattering Coefficients of Individual Components at 440.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.01833 0.01833 3 6.200 5.000 0.01833 0.01833 5 12.400 10.000 0.01833 0.01833 7 18.600 15.000 0.01833 0.01833 9 24.800 20.000 0.01833 0.01833 11 31.000 25.000 0.01833 0.01833 13 37.200 30.000 0.01833 0.01833 15 43.400 35.000 0.01833 0.01833 17 49.600 40.000 0.01833 0.01833 19 55.800 45.000 0.01833 0.01833 21 62.000 50.000 0.01833 0.01833 Summary of Inherent Optical Properties at 440.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.2400 1.0000 1.2400 0.8065 0.01833 0.01833 3 6.200 5.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 5 12.400 10.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 7 18.600 15.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 9 24.800 20.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 11 31.000 25.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 13 37.200 30.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 15 43.400 35.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 17 49.600 40.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 19 55.800 45.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 21 62.000 50.000 0.2400 1.0000 1.2400 0.8065 0.01833 0.01833 The shape of the asymptotic radiance distribution is given by u theta Linf(theta) 10 0.00 1.0000E+00 9 10.00 9.4078E-01 8 20.00 7.9982E-01 7 30.00 6.1872E-01 6 40.00 4.4288E-01 5 50.00 2.9925E-01 4 60.00 1.9503E-01 3 70.00 1.2526E-01 2 80.00 8.1069E-02 1 87.50 5.8934E-02 1 92.50 4.7341E-02 2 100.00 3.5647E-02 3 110.00 2.5461E-02 4 120.00 1.9048E-02 5 130.00 1.5026E-02 6 140.00 1.2496E-02 7 150.00 1.0919E-02 8 160.00 9.9724E-03 9 170.00 9.4676E-03 10 180.00 9.2918E-03 Asymptotic AOP values are MUu(inf) = 0.3607 MUd(inf) = 0.7202 MU(inf) = 0.6428 R(inf) = 0.0387 Lu(inf)/Ed(inf) = 0.007020 (1/sr) k(inf) = 0.3734 (1/m) absorp(check) = 0.2400 = k(inf)*MU(inf) Irradiances (units of W/m^2 nm), Mean Cosines (Mubars), and Irradiance Reflectance at 440.0 nm iz zeta z(m) Eou Eod Eo Eu Ed MUBARu MUBARd MUBAR R = Eu/Ed in air 1.7528E-01 1.8980E+00 2.0733E+00 6.2572E-02 1.3109E+00 0.3570 0.6907 0.6021 4.7731E-02 1 0.00 0.00 9.5270E-02 1.4922E+00 1.5874E+00 3.8796E-02 1.2871E+00 0.4072 0.8626 0.7864 3.0141E-02 3 6.20 5.00 2.3728E-02 3.1609E-01 3.3982E-01 8.7011E-03 2.3480E-01 0.3667 0.7428 0.6654 3.7058E-02 5 12.40 10.00 3.9721E-03 5.1788E-02 5.5760E-02 1.4377E-03 3.7553E-02 0.3620 0.7251 0.6477 3.8285E-02 7 18.60 15.00 6.2490E-04 8.1111E-03 8.7360E-03 2.2557E-04 5.8512E-03 0.3610 0.7214 0.6440 3.8552E-02 9 24.80 20.00 9.7019E-05 1.2580E-03 1.3550E-03 3.5000E-05 9.0637E-04 0.3608 0.7205 0.6431 3.8616E-02 11 31.00 25.00 1.5016E-05 1.9464E-04 2.0966E-04 5.4161E-06 1.4020E-04 0.3607 0.7203 0.6428 3.8632E-02 13 37.20 30.00 2.3221E-06 3.0098E-05 3.2420E-05 8.3752E-07 2.1677E-05 0.3607 0.7202 0.6428 3.8637E-02 15 43.40 35.00 3.5902E-07 4.6534E-06 5.0125E-06 1.2949E-07 3.3513E-06 0.3607 0.7202 0.6428 3.8639E-02 17 49.60 40.00 5.5506E-08 7.1944E-07 7.7494E-07 2.0020E-08 5.1812E-07 0.3607 0.7202 0.6428 3.8639E-02 19 55.80 45.00 8.5801E-09 1.1121E-07 1.1979E-07 3.0946E-09 8.0089E-08 0.3607 0.7202 0.6428 3.8639E-02 21 62.00 50.00 1.3264E-09 1.7192E-08 1.8518E-08 4.7840E-10 1.2381E-08 0.3607 0.7202 0.6428 3.8639E-02 K-functions (units of 1/meter) at 440.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.010 0.005 0.11014 0.21958 0.21301 0.18321 0.30164 0.30533 0.27285 5.000 5.010 5.005 0.34174 0.35055 0.34994 0.34710 0.36022 0.36073 0.35628 10.000 10.010 10.005 0.36673 0.36846 0.36833 0.36776 0.37044 0.37055 0.36964 15.000 15.010 15.005 0.37183 0.37222 0.37219 0.37206 0.37267 0.37269 0.37249 20.000 20.010 20.005 0.37298 0.37308 0.37307 0.37304 0.37320 0.37320 0.37315 25.000 25.010 25.005 0.37328 0.37330 0.37330 0.37329 0.37334 0.37334 0.37332 30.000 30.010 30.005 0.37336 0.37337 0.37336 0.37336 0.37338 0.37338 0.37337 35.000 35.010 35.005 0.37341 0.37339 0.37339 0.37339 0.37339 0.37339 0.37339 40.000 40.010 40.005 0.37340 0.37340 0.37340 0.37339 0.37339 0.37339 0.37339 45.000 45.010 45.005 0.37343 0.37341 0.37341 0.37340 0.37339 0.37339 0.37339 50.000 50.010 50.005 0.37339 0.37339 0.37339 0.37339 0.37339 0.37339 0.37339 Radiances Just Above the Water Surface (at z = a) (units of W/m^2 sr nm) at 440.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 1.19974E-01 4.54823E-03 7.46607E-03 1.20143E-02 3.46945E-03 0.062231 9 1 10.0 0.0 180.0 7.40291E-02 4.70602E-03 2.46461E-03 7.17063E-03 3.58982E-03 0.033292 9 4 10.0 45.0 135.0 8.13379E-02 4.64652E-03 2.83367E-03 7.48020E-03 3.54444E-03 0.034838 9 7 10.0 90.0 90.0 1.09894E-01 4.58502E-03 5.53922E-03 1.01242E-02 3.49752E-03 0.050405 8 1 20.0 0.0 180.0 5.81488E-02 4.80731E-03 1.41003E-03 6.21734E-03 3.66709E-03 0.024249 8 4 20.0 45.0 135.0 6.47353E-02 4.68964E-03 1.64739E-03 6.33702E-03 3.57732E-03 0.025448 8 7 20.0 90.0 90.0 9.51901E-02 4.62139E-03 3.43930E-03 8.06069E-03 3.52526E-03 0.036131 7 1 30.0 0.0 180.0 5.04181E-02 4.92178E-03 1.21369E-03 6.13546E-03 3.75440E-03 0.024072 7 4 30.0 45.0 135.0 5.58904E-02 4.81980E-03 1.36213E-03 6.18193E-03 3.67661E-03 0.024371 7 7 30.0 90.0 90.0 8.15581E-02 4.74021E-03 2.28829E-03 7.02850E-03 3.61590E-03 0.028057 6 1 40.0 0.0 180.0 4.78564E-02 4.99557E-03 1.34947E-03 6.34504E-03 3.81069E-03 0.028198 6 4 40.0 45.0 135.0 5.23284E-02 4.87769E-03 1.46358E-03 6.34127E-03 3.72078E-03 0.027969 6 7 40.0 90.0 90.0 7.34682E-02 4.88166E-03 2.11628E-03 6.99794E-03 3.72380E-03 0.028805 5 1 50.0 0.0 180.0 5.01758E-02 5.07130E-03 2.13170E-03 7.20300E-03 3.86846E-03 0.042485 5 4 50.0 45.0 135.0 5.34691E-02 4.95332E-03 2.22084E-03 7.17416E-03 3.77846E-03 0.041535 5 7 50.0 90.0 90.0 7.10737E-02 5.06457E-03 2.87111E-03 7.93568E-03 3.86332E-03 0.040396 4 1 60.0 0.0 180.0 5.95595E-02 5.01032E-03 4.96586E-03 9.97618E-03 3.82195E-03 0.083377 4 4 60.0 45.0 135.0 6.09074E-02 4.92489E-03 4.90370E-03 9.82859E-03 3.75678E-03 0.080511 4 7 60.0 90.0 90.0 7.52126E-02 5.19550E-03 5.72641E-03 1.09219E-02 3.96320E-03 0.076136 3 1 70.0 0.0 180.0 8.26392E-02 4.67870E-03 1.37170E-02 1.83957E-02 3.56898E-03 0.165987 3 4 70.0 45.0 135.0 8.03132E-02 4.69880E-03 1.31983E-02 1.78971E-02 3.58431E-03 0.164336 3 7 70.0 90.0 90.0 8.97883E-02 5.03709E-03 1.43526E-02 1.93897E-02 3.84237E-03 0.159849 2 1 80.0 0.0 180.0 1.40405E-01 4.04040E-03 3.40965E-02 3.81369E-02 3.08208E-03 0.242843 2 4 80.0 45.0 135.0 1.30012E-01 4.08536E-03 3.29784E-02 3.70637E-02 3.11637E-03 0.253655 2 7 80.0 90.0 90.0 1.29168E-01 4.42845E-03 3.45373E-02 3.89658E-02 3.37809E-03 0.267382 1 1 87.5 0.0 180.0 2.27009E-01 3.29309E-03 6.01339E-02 6.34269E-02 2.51202E-03 0.264896 1 4 87.5 45.0 135.0 2.04394E-01 3.19145E-03 5.97626E-02 6.29540E-02 2.43448E-03 0.292389 1 7 87.5 90.0 90.0 1.85954E-01 3.85174E-03 6.11631E-02 6.50148E-02 2.93816E-03 0.328916 Selected Radiances (units of W/m^2 sr nm) and Radiance-Irradiance Ratios at 440.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 1.201E-02 1.200E-01 1.452E-01 1.255E-01 1.880E-01 9.165E-03 5.208E+00 4.548E-03 3.469E-03 1 0.00 0.00 8.432E-03 2.233E-01 1.840E-02 2.318E-02 3.454E-02 6.551E-03 4.601E+00 3 6.20 5.00 1.628E-03 1.584E-01 5.971E-03 8.017E-03 1.229E-02 6.933E-03 5.345E+00 5 12.40 10.00 2.629E-04 2.834E-02 1.191E-03 1.440E-03 1.838E-03 7.002E-03 5.468E+00 7 18.60 15.00 4.105E-05 4.447E-03 2.083E-04 2.318E-04 2.632E-04 7.016E-03 5.495E+00 9 24.80 20.00 6.362E-06 6.868E-04 3.424E-05 3.625E-05 3.863E-05 7.020E-03 5.501E+00 11 31.00 25.00 9.842E-07 1.060E-04 5.459E-06 5.623E-06 5.803E-06 7.020E-03 5.503E+00 13 37.20 30.00 1.522E-07 1.639E-05 8.570E-07 8.700E-07 8.838E-07 7.021E-03 5.503E+00 15 43.40 35.00 2.353E-08 2.533E-06 1.335E-07 1.345E-07 1.356E-07 7.021E-03 5.504E+00 17 49.60 40.00 3.638E-09 3.915E-07 2.072E-08 2.080E-08 2.088E-08 7.021E-03 5.504E+00 19 55.80 45.00 5.623E-10 6.052E-08 3.210E-09 3.216E-09 3.222E-09 7.021E-03 5.504E+00 21 62.00 50.00 8.692E-11 9.355E-09 4.966E-10 4.971E-10 4.976E-10 7.021E-03 5.504E+00 Waveband 1 of 1 completed in 6.0 sec. Energy Eo as a function of depth and wavelength Units are W/m^2 [values are integrated over the wavelength bands] depth 440.0 in air 2.073E+00 0.0 1.587E+00 5.0 3.398E-01 10.0 5.576E-02 15.0 8.736E-03 20.0 1.355E-03 25.0 2.097E-04 30.0 3.242E-05 35.0 5.012E-06 40.0 7.749E-07 45.0 1.198E-07 50.0 1.852E-08 Quantum Eo as a function of depth and wavelength Units are (micromol photons)/(m^2 s) [values are integrated over the wavelength bands] depth 440.0 in air 7.625E+00 0.0 5.838E+00 5.0 1.250E+00 10.0 2.051E-01 15.0 3.213E-02 20.0 4.983E-03 25.0 7.710E-04 30.0 1.192E-04 35.0 1.843E-05 40.0 2.850E-06 45.0 4.405E-07 50.0 6.810E-08 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 Run completed on 07/20/2011 at 10:47:22.5 Total (wall clock) run time = 6.1 sec