======================================================== | | | 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: run1 for iop error exercise 2 wed class (07/20/2011 10:42) Run started on 07/20/2011 at 10:42:34.0 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.3000 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.3000 0.3000 3 6.500 5.000 0.3000 0.3000 5 13.000 10.000 0.3000 0.3000 7 19.500 15.000 0.3000 0.3000 9 26.000 20.000 0.3000 0.3000 11 32.500 25.000 0.3000 0.3000 13 39.000 30.000 0.3000 0.3000 15 45.500 35.000 0.3000 0.3000 17 52.000 40.000 0.3000 0.3000 19 58.500 45.000 0.3000 0.3000 21 65.000 50.000 0.3000 0.3000 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.500 5.000 1.0000 1.0000 5 13.000 10.000 1.0000 1.0000 7 19.500 15.000 1.0000 1.0000 9 26.000 20.000 1.0000 1.0000 11 32.500 25.000 1.0000 1.0000 13 39.000 30.000 1.0000 1.0000 15 45.500 35.000 1.0000 1.0000 17 52.000 40.000 1.0000 1.0000 19 58.500 45.000 1.0000 1.0000 21 65.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.500 5.000 0.01833 0.01833 5 13.000 10.000 0.01833 0.01833 7 19.500 15.000 0.01833 0.01833 9 26.000 20.000 0.01833 0.01833 11 32.500 25.000 0.01833 0.01833 13 39.000 30.000 0.01833 0.01833 15 45.500 35.000 0.01833 0.01833 17 52.000 40.000 0.01833 0.01833 19 58.500 45.000 0.01833 0.01833 21 65.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.3000 1.0000 1.3000 0.7692 0.01833 0.01833 3 6.500 5.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 5 13.000 10.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 7 19.500 15.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 9 26.000 20.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 11 32.500 25.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 13 39.000 30.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 15 45.500 35.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 17 52.000 40.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 19 58.500 45.000 0.3000 1.0000 1.3000 0.7692 0.01833 0.01833 21 65.000 50.000 0.3000 1.0000 1.3000 0.7692 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.3227E-01 8 20.00 7.7407E-01 7 30.00 5.7781E-01 6 40.00 3.9603E-01 5 50.00 2.5542E-01 4 60.00 1.5906E-01 3 70.00 9.8002E-02 2 80.00 6.1227E-02 1 87.50 4.3545E-02 1 92.50 3.4567E-02 2 100.00 2.5666E-02 3 110.00 1.8089E-02 4 120.00 1.3435E-02 5 130.00 1.0574E-02 6 140.00 8.8041E-03 7 150.00 7.7167E-03 8 160.00 7.0733E-03 9 170.00 6.7339E-03 10 180.00 6.6166E-03 Asymptotic AOP values are MUu(inf) = 0.3577 MUd(inf) = 0.7393 MU(inf) = 0.6756 R(inf) = 0.0299 Lu(inf)/Ed(inf) = 0.005449 (1/sr) k(inf) = 0.4441 (1/m) absorp(check) = 0.3000 = 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.6837E-01 1.8980E+00 2.0664E+00 5.9045E-02 1.3109E+00 0.3507 0.6907 0.6058 4.5041E-02 1 0.00 0.00 7.5455E-02 1.4759E+00 1.5514E+00 3.0326E-02 1.2822E+00 0.4019 0.8687 0.8069 2.3652E-02 3 6.50 5.00 1.3037E-02 2.1624E-01 2.2927E-01 4.7267E-03 1.6372E-01 0.3626 0.7571 0.6935 2.8871E-02 5 13.00 10.00 1.5136E-03 2.4632E-02 2.6145E-02 5.4276E-04 1.8293E-02 0.3586 0.7427 0.6789 2.9670E-02 7 19.50 15.00 1.6642E-04 2.6989E-03 2.8653E-03 5.9558E-05 1.9972E-03 0.3579 0.7400 0.6762 2.9821E-02 9 26.00 20.00 1.8114E-05 2.9357E-04 3.1168E-04 6.4801E-06 2.1709E-04 0.3577 0.7395 0.6757 2.9851E-02 11 32.50 25.00 1.9677E-06 3.1885E-05 3.3853E-05 7.0387E-07 2.3575E-05 0.3577 0.7394 0.6756 2.9857E-02 13 39.00 30.00 2.1366E-07 3.4620E-06 3.6757E-06 7.6426E-08 2.5596E-06 0.3577 0.7393 0.6756 2.9859E-02 15 45.50 35.00 2.3197E-08 3.7587E-07 3.9907E-07 8.2976E-09 2.7789E-07 0.3577 0.7393 0.6756 2.9859E-02 17 52.00 40.00 2.5183E-09 4.0804E-08 4.3323E-08 9.0078E-10 3.0168E-08 0.3577 0.7393 0.6756 2.9859E-02 19 58.50 45.00 2.7340E-10 4.4300E-09 4.7034E-09 9.7795E-11 3.2752E-09 0.3577 0.7393 0.6756 2.9859E-02 21 65.00 50.00 2.9683E-11 4.8096E-10 5.1065E-10 1.0618E-11 3.5559E-10 0.3577 0.7393 0.6756 2.9859E-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.15198 0.28148 0.27518 0.24029 0.36884 0.37195 0.33896 5.000 5.010 5.005 0.41612 0.42422 0.42376 0.42081 0.43229 0.43263 0.42886 10.000 10.010 10.005 0.43893 0.44035 0.44027 0.43975 0.44182 0.44189 0.44120 15.000 15.010 15.005 0.44307 0.44334 0.44332 0.44322 0.44362 0.44363 0.44350 20.000 20.010 20.005 0.44386 0.44392 0.44392 0.44389 0.44397 0.44398 0.44395 25.000 25.010 25.005 0.44403 0.44404 0.44404 0.44403 0.44405 0.44405 0.44405 30.000 30.010 30.005 0.44406 0.44406 0.44406 0.44406 0.44407 0.44407 0.44407 35.000 35.010 35.005 0.44414 0.44408 0.44408 0.44408 0.44408 0.44408 0.44408 40.000 40.010 40.005 0.44407 0.44408 0.44408 0.44408 0.44408 0.44408 0.44408 45.000 45.010 45.005 0.44409 0.44408 0.44408 0.44408 0.44408 0.44408 0.44408 50.000 50.010 50.005 0.44409 0.44409 0.44409 0.44408 0.44408 0.44408 0.44408 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 3.53561E-03 7.46607E-03 1.10017E-02 2.69701E-03 0.062231 9 1 10.0 0.0 180.0 7.40291E-02 3.68104E-03 2.46461E-03 6.14565E-03 2.80795E-03 0.033292 9 4 10.0 45.0 135.0 8.13379E-02 3.62510E-03 2.83367E-03 6.45877E-03 2.76528E-03 0.034838 9 7 10.0 90.0 90.0 1.09894E-01 3.55939E-03 5.53922E-03 9.09861E-03 2.71515E-03 0.050405 8 1 20.0 0.0 180.0 5.81488E-02 3.76883E-03 1.41003E-03 5.17886E-03 2.87492E-03 0.024249 8 4 20.0 45.0 135.0 6.47353E-02 3.66060E-03 1.64739E-03 5.30798E-03 2.79236E-03 0.025448 8 7 20.0 90.0 90.0 9.51901E-02 3.58069E-03 3.43930E-03 7.02000E-03 2.73140E-03 0.036131 7 1 30.0 0.0 180.0 5.04181E-02 3.86235E-03 1.21369E-03 5.07604E-03 2.94626E-03 0.024072 7 4 30.0 45.0 135.0 5.58904E-02 3.76033E-03 1.36213E-03 5.12246E-03 2.86843E-03 0.024371 7 7 30.0 90.0 90.0 8.15581E-02 3.66544E-03 2.28829E-03 5.95373E-03 2.79605E-03 0.028057 6 1 40.0 0.0 180.0 4.78564E-02 3.91600E-03 1.34947E-03 5.26547E-03 2.98718E-03 0.028198 6 4 40.0 45.0 135.0 5.23284E-02 3.79671E-03 1.46358E-03 5.26029E-03 2.89619E-03 0.027969 6 7 40.0 90.0 90.0 7.34682E-02 3.76727E-03 2.11628E-03 5.88354E-03 2.87373E-03 0.028805 5 1 50.0 0.0 180.0 5.01758E-02 3.96637E-03 2.13170E-03 6.09807E-03 3.02561E-03 0.042485 5 4 50.0 45.0 135.0 5.34691E-02 3.84770E-03 2.22084E-03 6.06854E-03 2.93508E-03 0.041535 5 7 50.0 90.0 90.0 7.10737E-02 3.90717E-03 2.87111E-03 6.77829E-03 2.98045E-03 0.040396 4 1 60.0 0.0 180.0 5.95595E-02 3.90956E-03 4.96586E-03 8.87542E-03 2.98227E-03 0.083377 4 4 60.0 45.0 135.0 6.09074E-02 3.81927E-03 4.90370E-03 8.72297E-03 2.91340E-03 0.080511 4 7 60.0 90.0 90.0 7.52126E-02 4.00919E-03 5.72641E-03 9.73560E-03 3.05827E-03 0.076136 3 1 70.0 0.0 180.0 8.26392E-02 3.64448E-03 1.37170E-02 1.73615E-02 2.78006E-03 0.165987 3 4 70.0 45.0 135.0 8.03132E-02 3.64279E-03 1.31983E-02 1.68411E-02 2.77877E-03 0.164336 3 7 70.0 90.0 90.0 8.97883E-02 3.89256E-03 1.43526E-02 1.82451E-02 2.96930E-03 0.159849 2 1 80.0 0.0 180.0 1.40405E-01 3.14460E-03 3.40965E-02 3.72411E-02 2.39875E-03 0.242843 2 4 80.0 45.0 135.0 1.30012E-01 3.16861E-03 3.29784E-02 3.61470E-02 2.41706E-03 0.253655 2 7 80.0 90.0 90.0 1.29168E-01 3.42805E-03 3.45373E-02 3.79654E-02 2.61497E-03 0.267382 1 1 87.5 0.0 180.0 2.27009E-01 2.56427E-03 6.01339E-02 6.26981E-02 1.95606E-03 0.264896 1 4 87.5 45.0 135.0 2.04394E-01 2.47916E-03 5.97626E-02 6.22417E-02 1.89114E-03 0.292389 1 7 87.5 90.0 90.0 1.85954E-01 2.98867E-03 6.11631E-02 6.41518E-02 2.27980E-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.100E-02 1.200E-01 1.449E-01 1.251E-01 1.875E-01 8.392E-03 5.367E+00 3.536E-03 2.697E-03 1 0.00 0.00 6.555E-03 2.233E-01 1.478E-02 1.871E-02 2.834E-02 5.112E-03 4.626E+00 3 6.50 5.00 8.832E-04 1.152E-01 3.379E-03 4.555E-03 7.032E-03 5.395E-03 5.352E+00 5 13.00 10.00 9.950E-05 1.477E-02 4.677E-04 5.653E-04 7.224E-04 5.439E-03 5.455E+00 7 19.50 15.00 1.088E-05 1.643E-03 5.714E-05 6.353E-05 7.205E-05 5.447E-03 5.474E+00 9 26.00 20.00 1.183E-06 1.789E-04 6.582E-06 6.961E-06 7.410E-06 5.449E-03 5.478E+00 11 32.50 25.00 1.285E-07 1.942E-05 7.362E-07 7.576E-07 7.813E-07 5.449E-03 5.479E+00 13 39.00 30.00 1.395E-08 2.108E-06 8.112E-08 8.231E-08 8.357E-08 5.449E-03 5.479E+00 15 45.50 35.00 1.514E-09 2.289E-07 8.873E-09 8.938E-09 9.005E-09 5.449E-03 5.479E+00 17 52.00 40.00 1.644E-10 2.485E-08 9.668E-10 9.704E-10 9.740E-10 5.449E-03 5.479E+00 19 58.50 45.00 1.785E-11 2.697E-09 1.052E-10 1.054E-10 1.055E-10 5.449E-03 5.479E+00 21 65.00 50.00 1.938E-12 2.929E-10 1.143E-11 1.144E-11 1.145E-11 5.449E-03 5.479E+00 Waveband 1 of 1 completed in 6.3 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.066E+00 0.0 1.551E+00 5.0 2.293E-01 10.0 2.615E-02 15.0 2.865E-03 20.0 3.117E-04 25.0 3.385E-05 30.0 3.676E-06 35.0 3.991E-07 40.0 4.332E-08 45.0 4.703E-09 50.0 5.106E-10 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.599E+00 0.0 5.705E+00 5.0 8.432E-01 10.0 9.615E-02 15.0 1.054E-02 20.0 1.146E-03 25.0 1.245E-04 30.0 1.352E-05 35.0 1.468E-06 40.0 1.593E-07 45.0 1.730E-08 50.0 1.878E-09 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:42:40.5 Total (wall clock) run time = 6.4 sec