Tue Mar 18 17:20:31 EST 1997 rundis revision 1.2, normal precision, attempted run number 2 RCS/DISOTEST.f,v --> DISOTEST.f revision 1.2 done RCS/DISORT.f,v --> DISORT.f revision 1.2 done RCS/ErrPack.f,v --> ErrPack.f revision 1.2 done RCS/R1MACH.f,v --> R1MACH.f revision 1.2 done RCS/D1MACH.f,v --> D1MACH.f revision 1.2 done RCS/LINPAK.f,v --> LINPAK.f revision 1.2 done code.f: MAIN testdo: ANSI warning on line 1441 of code.f: Lower case letter(s) found getmom: ANSI warning on line 1588 of code.f: Lower case letter(s) found prtfin: ANSI warning on line 1808 of code.f: Lower case letter(s) found BLOCK DATA chekdo: ANSI warning on line 2347 of code.f: Lower case letter(s) found disort: ANSI warning on line 3131 of code.f: Lower case letter(s) found albtrn: ANSI warning on line 3411 of code.f: Lower case letter(s) found altrin: ANSI warning on line 3577 of code.f: Lower case letter(s) found asymtx: ANSI warning on line 4374 of code.f: Lower case letter(s) found chekin: ANSI warning on line 4653 of code.f: Lower case letter(s) found cmpint: ANSI warning on line 4770 of code.f: Lower case letter(s) found fluxes: ANSI warning on line 5023 of code.f: Lower case letter(s) found lepoly: ANSI warning on line 5149 of code.f: Lower case letter(s) found praltr: ANSI warning on line 5194 of code.f: Lower case letter(s) found pravin: ANSI warning on line 5251 of code.f: Lower case letter(s) found prtinp: ANSI warning on line 5423 of code.f: Lower case letter(s) found prtint: ANSI warning on line 5494 of code.f: Lower case letter(s) found qgausn: ANSI warning on line 5666 of code.f: Lower case letter(s) found setdis: ANSI warning on line 5912 of code.f: Lower case letter(s) found setmtx: ANSI warning on line 6204 of code.f: Lower case letter(s) found slftst: ANSI warning on line 6428 of code.f: Lower case letter(s) found soleig: ANSI warning on line 6647 of code.f: Lower case letter(s) found solve0: ANSI warning on line 6955 of code.f: Lower case letter(s) found solve1: ANSI warning on line 7055 of code.f: Lower case letter(s) found spaltr: ANSI warning on line 7164 of code.f: Lower case letter(s) found surfac: ANSI warning on line 7460 of code.f: Lower case letter(s) found terpev: ANSI warning on line 7518 of code.f: Lower case letter(s) found terpso: ANSI warning on line 7656 of code.f: Lower case letter(s) found upbeam: ANSI warning on line 7769 of code.f: Lower case letter(s) found upisot: ANSI warning on line 7868 of code.f: Lower case letter(s) found usrint: ANSI warning on line 8280 of code.f: Lower case letter(s) found zeroal: ANSI warning on line 8438 of code.f: Lower case letter(s) found zeroit: ANSI warning on line 8463 of code.f: Lower case letter(s) found dref: ANSI warning on line 8572 of code.f: Lower case letter(s) found plkavg: ANSI warning on line 8825 of code.f: Lower case letter(s) found ratio: ANSI warning on line 8923 of code.f: Lower case letter(s) found errmsg: ANSI warning on line 8960 of code.f: Lower case letter(s) found wrtbad: ANSI warning on line 8983 of code.f: Lower case letter(s) found wrtdim: ANSI warning on line 9004 of code.f: Lower case letter(s) found tstbad: ANSI warning on line 9021 of code.f: Lower case letter(s) found r1mach: ANSI warning on line 9101 of code.f: Lower case letter(s) found d1mach: ANSI warning on line 9150 of code.f: Lower case letter(s) found sgbco: ANSI warning on line 9528 of code.f: Lower case letter(s) found sgbfa: ANSI warning on line 9675 of code.f: Lower case letter(s) found sgbsl: ANSI warning on line 9847 of code.f: Lower case letter(s) found sgeco: ANSI warning on line 10106 of code.f: Lower case letter(s) found sgefa: ANSI warning on line 10215 of code.f: Lower case letter(s) found sgesl: ANSI warning on line 10358 of code.f: Lower case letter(s) found sasum: ANSI warning on line 10417 of code.f: Lower case letter(s) found saxpy: ANSI warning on line 10503 of code.f: Lower case letter(s) found sdot: ANSI warning on line 10588 of code.f: Lower case letter(s) found sscal: ANSI warning on line 10654 of code.f: Lower case letter(s) found sswap: ANSI warning on line 10755 of code.f: Lower case letter(s) found isamax: ANSI warning on line 10816 of code.f: Lower case letter(s) found **************************************************************************************************** DISORT: Test Case No. 1a: Isotropic Scattering, Ref. VH1, Table 12: b = .03125, a = .20 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .0313 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .0313 .0313 .20000 .00000 .0313 .0313 .20000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 .0000E+00 7.9945E-02 2.5407E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .0313 2.2984E+00 7.9411E-02 -2.8982E-09 1.8653E+01 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 1.1777E-01 2.6417E-02 1.3404E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .031 1.3383E-02 2.6332E-02 1.1590E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 1b: Isotropic Scattering, Ref. VH1, Table 12: b = .03125, a = 1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .0313 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .0313 .0313 1.00000 .00000 .0313 .0313 1.00000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 5.5075E-05 4.2297E-01 3.9609E-05 ( 1.0000) (*********) ( 1.0001) (*********) .0313 2.2984E+00 4.2029E-01 3.1102E-06 2.9495E-05 ( 1.0000) ( 1.0001) (*********) (*********) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 6.2288E-01 1.3976E-01 7.0916E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( .9999) .031 7.0815E-02 1.3934E-01 6.1346E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0001) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ============================================= ==== 3 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 1c: Isotropic Scattering, Ref. VH1, Table 12: b = .03125, a = .99 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .0313 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .0313 .0313 .99000 .00000 .0313 .0313 .99000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 3.1416E+00 9.0656E-02 6.6687E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .0313 .0000E+00 3.0490E+00 1.8926E-07 5.8894E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 1.0000E+00 1.0000E+00 1.0000E+00 1.3318E-01 2.9988E-02 1.5223E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .031 9.8445E-01 9.6936E-01 8.6395E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 1d: Isotropic Scattering, Ref. VH1, Table 12: b = 32.00000, a = .20 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 32.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 32.0000 32.0000 .20000 .00000 32.0000 32.0000 .20000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 .0000E+00 2.5969E-01 2.5777E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 32.0000 .0000E+00 3.3252E-16 6.1753E-27 3.0317E-16 ( 1.0000) (*********) (*********) (*********) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 2.6297E-01 9.0697E-02 5.0285E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 32.000 1.2298E-15 1.3070E-17 6.8884E-18 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 1e: Isotropic Scattering, Ref. VH1, Table 12: b = 32.00000, a = 1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 32.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 32.0000 32.0000 1.00000 .00000 32.0000 32.0000 1.00000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 -2.6226E-06 3.0738E+00 4.6574E-05 ( 1.0000) (*********) ( 1.0000) (*********) 32.0000 .0000E+00 6.7648E-02 -1.5881E-06 1.3967E-07 ( 1.0000) ( .9993) (*********) (*********) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 1.9332E+00 1.0273E+00 7.9716E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 32.000 2.7112E-02 1.8767E-02 1.1630E-02 .0000E+00 .0000E+00 .0000E+00 ( .9993) ( .9993) ( .9993) ( 1.0000) ( 1.0000) ( 1.0000) ============================================= ==== 1 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 1f: Isotropic Scattering, Ref. VH1, Table 12: b = 32.00000, a = .99 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 32.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 32.0000 32.0000 .99000 .00000 32.0000 32.0000 .99000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 3.1416E+00 2.4962E+00 1.1424E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 32.0000 .0000E+00 4.6005E-03 1.0712E-09 7.9364E-05 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 1.0000E+00 1.0000E+00 1.0000E+00 8.7751E-01 8.1514E-01 7.5272E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 32.000 1.8684E-03 1.2649E-03 7.7929E-04 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 2a, Rayleigh Scattering, Ref. SW, Table 1: tau = .20, mu0 = .080442, ss-albedo = .50 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .2000 6 User polar angle cosines : -.98199 -.53826 -.01801 .01801 .53826 .98199 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .08044 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .2000 .2000 .50000 .00000 .2000 .2000 .50000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .100000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 2.5272E-01 .0000E+00 5.3506E-02 1.6657E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .2000 2.1031E-02 4.4179E-02 6.3229E-10 1.8985E-01 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -.982 .0000E+00 ( 1.0000) -.538 .0000E+00 ( 1.0000) -.018 .0000E+00 ( 1.0000) .018 1.6180E-01 ( 1.0000) .538 2.1150E-02 ( 1.0000) .982 7.8671E-03 ( 1.0000) .200 -.982 7.7190E-03 ( 1.0000) -.538 2.0078E-02 ( 1.0000) -.018 2.5768E-02 ( 1.0000) .018 .0000E+00 ( 1.0000) .538 .0000E+00 ( 1.0000) .982 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 2b, Rayleigh Scattering, Ref. SW, Table 1: tau = .20, mu0 = .080442, ss-albedo =1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .2000 6 User polar angle cosines : -.98199 -.53826 -.01801 .01801 .53826 .98199 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .08044 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .2000 .2000 1.00000 .00000 .2000 .2000 1.00000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .100000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 2.5272E-01 -4.0233E-06 1.2556E-01 4.2650E-06 ( 1.0000) (*********) ( 1.0000) (*********) .2000 2.1031E-02 1.0612E-01 -3.4159E-06 6.5426E-07 ( 1.0000) ( 1.0000) (*********) (*********) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -.982 .0000E+00 ( 1.0000) -.538 .0000E+00 ( 1.0000) -.018 .0000E+00 ( 1.0000) .018 3.4768E-01 ( 1.0000) .538 4.8712E-02 ( 1.0000) .982 1.8939E-02 ( 1.0000) .200 -.982 1.8602E-02 ( 1.0000) -.538 4.6406E-02 ( 1.0000) -.018 6.7759E-02 ( 1.0000) .018 .0000E+00 ( 1.0000) .538 .0000E+00 ( 1.0000) .982 .0000E+00 ( 1.0000) ============================================= ==== 4 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 2c, Rayleigh Scattering, Ref. SW, Table 1: tau = 5.00, mu0 = .080442, ss-albedo = .50 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 5.0000 6 User polar angle cosines : -.98199 -.53826 -.01801 .01801 .53826 .98199 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .08044 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 5.0000 5.0000 .50000 .00000 5.0000 5.0000 .50000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .100000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 2.5272E-01 .0000E+00 6.2473E-02 1.6746E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 5.0000 2.5608E-28 2.5168E-04 -1.2547E-12 1.7546E-04 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -.982 .0000E+00 ( 1.0000) -.538 .0000E+00 ( 1.0000) -.018 .0000E+00 ( 1.0000) .018 1.6257E-01 ( 1.0000) .538 2.4579E-02 ( 1.0000) .982 1.0150E-02 ( 1.0000) 5.000 -.982 1.7000E-04 ( 1.0000) -.538 3.9717E-05 ( 1.0000) -.018 1.3247E-05 ( 1.0000) .018 .0000E+00 ( 1.0000) .538 .0000E+00 ( 1.0000) .982 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 2d, Rayleigh Scattering, Ref. SW, Table 1: tau = 5.00, mu0 = .080442, ss-albedo =1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 5.0000 6 User polar angle cosines : -.98199 -.53826 -.01801 .01801 .53826 .98199 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .08044 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 5.0000 5.0000 1.00000 .00000 5.0000 5.0000 1.00000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .100000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 2.5272E-01 2.9802E-07 2.2591E-01 4.4946E-06 ( 1.0000) (*********) ( 1.0000) (*********) 5.0000 2.5608E-28 2.6800E-02 6.3698E-08 5.5175E-08 ( 1.0000) ( 1.0000) (*********) (*********) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -.982 .0000E+00 ( 1.0000) -.538 .0000E+00 ( 1.0000) -.018 .0000E+00 ( 1.0000) .018 3.6401E-01 ( 1.0000) .538 8.2699E-02 ( 1.0000) .982 4.9237E-02 ( 1.0000) 5.000 -.982 1.0595E-02 ( 1.0000) -.538 7.6934E-03 ( 1.0000) -.018 3.7927E-03 ( 1.0000) .018 .0000E+00 ( 1.0000) .538 .0000E+00 ( 1.0000) .982 .0000E+00 ( 1.0000) ============================================= ==== 2 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 3a, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 1.00000, a = .20 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .20000 .01002 .9980 .9980 .19839 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 2.3842E-07 1.6995E-01 2.5922E+01 ( 1.0000) (*********) ( 1.0000) ( 1.0000) 1.0000 1.4263E-04 4.7030E-02 4.5290E-11 7.4364E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 4.8692E-01 5.0880E-02 1.0335E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.000 7.8852E-03 2.2208E-02 2.9000E-03 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 3b, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 1.00000, a = 1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 1.00000 .01002 .9900 .9900 1.00000 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 -2.6226E-06 1.8269E+00 4.7367E-05 ( 1.0000) (*********) ( 1.0000) (*********) 1.0000 1.4263E-04 1.3146E+00 -1.1493E-05 3.4137E-06 ( 1.0000) ( 1.0000) (*********) (*********) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 3.7289E+00 6.3207E-01 1.5121E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.000 2.1400E-01 5.8401E-01 3.7716E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ============================================= ==== 3 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 3c, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 1.00000, a = .99 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .99000 .01002 .9901 .9901 .98990 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 3.1416E+00 5.8312E-01 8.0694E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 .0000E+00 2.4966E+00 -2.5251E-08 4.3524E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 1.0000E+00 1.0000E+00 1.0000E+00 5.6662E-01 2.3221E-01 7.6027E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.000 9.1118E-01 7.4403E-01 4.0361E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 3d, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 8.00000, a = .20 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 8.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 8.0000 8.0000 .20000 .01002 7.9840 7.9840 .19839 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 2.3842E-07 1.7006E-01 2.5922E+01 ( 1.0000) (*********) ( 1.0000) ( 1.0000) 8.0000 5.6701E-35 1.9981E-05 5.8275E-15 1.9194E-05 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 4.8693E-01 5.0919E-02 1.0369E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 8.000 2.0794E-05 9.7696E-07 1.8800E-07 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 3e, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 8.00000, a = 1.00 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 8.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+01 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 8.0000 8.0000 1.00000 .01002 7.9198 7.9198 1.00000 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 2.3842E-06 2.5496E+00 4.8923E-05 ( 1.0000) (*********) ( 1.0000) (*********) 8.0000 5.6701E-35 5.9196E-01 -4.9312E-06 1.2047E-06 ( 1.0000) ( 1.0000) (*********) (*********) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 .0000E+00 .0000E+00 .0000E+00 3.8577E+00 8.6801E-01 3.7979E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 8.000 2.3636E-01 1.6435E-01 9.1499E-02 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ============================================= ==== 2 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 3f, Henyey-Greenstein Scattering, Ref. VH2, Table 35, g = 0.75, b = 8.00000, a = .99 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 8.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .10000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 8.0000 8.0000 .99000 .01002 7.9206 7.9206 .98990 .7500 Layer Phase Function Moments 1 1.000000 .750000 .562500 .421875 .316406 .237305 .177979 .133484 .100113 .075085 .056314 .042235 .031676 .023757 .017818 .013363 .010023 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 3.1416E+00 1.6875E+00 1.0046E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 8.0000 .0000E+00 1.0117E+00 -9.4839E-08 1.7139E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ********* AZIMUTHALLY AVERAGED INTENSITIES ********* Optical Polar Angle Cosines Depth -1.00000 -.50000 -.10000 .10000 .50000 1.00000 .000 1.0000E+00 1.0000E+00 1.0000E+00 7.5732E-01 5.8726E-01 4.3305E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 8.000 4.1633E-01 2.7596E-01 1.5179E-01 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 4a, Haze-L Scattering, Ref. GS, Table 12 **************************************************************************************************** No. streams = 32 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 1.00000 .00026 .9997 .9997 1.00000 .8042 Layer Phase Function Moments 1 1.000000 .804200 .646094 .481851 .359056 .263045 .191995 .140907 .102831 .076154 .056054 .042019 .031295 .023746 .017919 .013730 .010512 .008100 .006302 .004862 .003851 .002963 .002392 .001829 .001506 .001143 .000960 .000722 .000618 .000461 .000402 .000297 .000263 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 -2.0504E-05 1.7324E-01 4.2449E-06 ( 1.0000) (*********) ( 1.0001) (*********) .5000 1.9055E+00 1.1740E+00 1.1111E-01 4.6398E-06 ( 1.0000) ( 1.0000) ( .9999) (*********) 1.0000 1.1557E+00 1.8126E+00 6.2371E-06 4.3662E-06 ( 1.0000) ( 1.0000) (*********) (*********) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -1.000 .0000E+00 ( 1.0000) -.500 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 9.2981E-02 ( .9999) .500 6.6124E-02 ( 1.0000) 1.000 3.5530E-02 ( .9998) .500 -1.000 2.4944E+00 ( 1.0000) -.500 1.1908E-01 ( 1.0000) -.100 1.3515E-01 ( 1.0001) .100 1.2405E-01 ( .9998) .500 4.0279E-02 ( .9998) 1.000 1.7360E-02 ( 1.0001) 1.000 -1.000 3.3467E+00 ( 1.0000) -.500 2.1961E-01 ( 1.0000) -.100 1.5702E-01 ( 1.0001) .100 .0000E+00 ( 1.0000) .500 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) ============================================= ==== 5 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 4b, Haze-L Scattering, Ref. GS, Table 13 **************************************************************************************************** No. streams = 32 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .90000 .00026 .9998 .9998 .89998 .8042 Layer Phase Function Moments 1 1.000000 .804200 .646094 .481851 .359056 .263045 .191995 .140907 .102831 .076154 .056054 .042019 .031295 .023746 .017919 .013730 .010512 .008100 .006302 .004862 .003851 .002963 .002392 .001829 .001506 .001143 .000960 .000722 .000618 .000461 .000402 .000297 .000263 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 -7.1526E-07 1.2367E-01 3.4372E-01 ( 1.0000) (*********) ( 1.0000) ( 1.0000) .5000 1.9055E+00 1.0152E+00 7.8869E-02 3.5239E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.1557E+00 1.5155E+00 3.9171E-08 3.1945E-01 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -1.000 .0000E+00 ( 1.0000) -.500 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 6.5578E-02 ( 1.0000) .500 4.5664E-02 ( 1.0000) 1.000 2.7424E-02 ( 1.0000) .500 -1.000 2.2230E+00 ( 1.0000) -.500 9.6410E-02 ( 1.0000) -.100 9.6292E-02 ( 1.0000) .100 8.4492E-02 ( 1.0000) .500 2.8022E-02 ( 1.0000) 1.000 1.3509E-02 ( 1.0000) 1.000 -1.000 2.9468E+00 ( 1.0000) -.500 1.6751E-01 ( 1.0000) -.100 1.0821E-01 ( 1.0000) .100 .0000E+00 ( 1.0000) .500 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 4c, Haze-L Scattering, Ref. GS, Table 14-16 **************************************************************************************************** No. streams = 32 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 3 User azimuthal angles : .00 90.00 180.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .90000 .00026 .9998 .9998 .89998 .8042 Layer Phase Function Moments 1 1.000000 .804200 .646094 .481851 .359056 .263045 .191995 .140907 .102831 .076154 .056054 .042019 .031295 .023746 .017919 .013730 .010512 .008100 .006302 .004862 .003851 .002963 .002392 .001829 .001506 .001143 .000960 .000722 .000618 .000461 .000402 .000297 .000263 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.5708E+00 4.7684E-07 2.2549E-01 3.8500E-01 ( 1.0000) (*********) ( 1.0000) ( 1.0000) .5000 5.7786E-01 7.0276E-01 1.2385E-01 3.3732E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 2.1258E-01 8.0329E-01 -5.0575E-08 2.1640E-01 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 90.0 180.0 .000 -1.000 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.500 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.100 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) .100 8.7002E-01 8.8611E-02 6.9529E-02 ( 1.0000) ( 1.0000) ( 1.0000) .500 2.2476E-01 5.7692E-02 4.9328E-02 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 2.2832E-02 2.2832E-02 2.2832E-02 ( 1.0000) ( 1.0000) ( 1.0000) .500 -1.000 4.7604E-02 4.7604E-02 4.7604E-02 ( 1.0000) ( 1.0000) ( 1.0000) -.500 3.0026E+00 5.8260E-02 2.5933E-02 ( 1.0000) ( 1.0001) ( .9997) -.100 1.4114E+00 1.0464E-01 6.2496E-02 ( 1.0000) ( 1.0000) ( 1.0000) .100 6.9740E-01 9.1533E-02 5.9019E-02 ( 1.0000) ( 1.0000) ( 1.0000) .500 1.0907E-01 2.9568E-02 2.4459E-02 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 9.3511E-03 9.3511E-03 9.3511E-03 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 -1.000 8.3754E-02 8.3754E-02 8.3754E-02 ( 1.0000) ( 1.0000) ( 1.0000) -.500 2.6879E+00 9.4332E-02 4.0003E-02 ( 1.0000) ( 1.0000) ( 1.0000) -.100 8.7632E-01 8.9596E-02 4.6679E-02 ( 1.0000) ( 1.0000) ( 1.0000) .100 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) .500 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 .0000E+00 .0000E+00 .0000E+00 ( 1.0000) ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 5a, Cloud C.1 Scattering, Ref. GS, Table 19 **************************************************************************************************** No. streams = 48 No. computational layers = 1 3 User optical depths : .0000 32.0000 64.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 64.0000 64.0000 1.00000 .20938 50.5996 50.5996 1.00000 .8480 Layer Phase Function Moments 1 1.000000 .848000 .776600 .652571 .581667 .535182 .496692 .478467 .462294 .447053 .442190 .428522 .424600 .415889 .408655 .403323 .395697 .389314 .384027 .375897 .371488 .363744 .358356 .351894 .345592 .339706 .333453 .327255 .321561 .315051 .309590 .303222 .297615 .291597 .285812 .280056 .274301 .268600 .263000 .257342 .251864 .246313 .240906 .235506 .230146 .224912 .219645 .214547 .209381 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 3.1416E+00 4.7684E-07 2.6614E+00 9.2949E-06 ( 1.0000) (*********) ( .9999) (*********) 32.0000 3.9786E-14 2.2469E+00 1.7671E+00 9.5699E-06 ( 1.0000) ( .9996) ( .9996) (*********) 64.0000 5.0385E-28 4.7960E-01 1.6262E-06 9.6938E-07 ( 1.0000) ( .9995) (*********) (*********) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -1.000 .0000E+00 ( 1.0000) -.500 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 4.5760E-01 ( .9999) .500 7.6884E-01 ( .9999) 1.000 1.0311E+00 ( .9999) 32.000 -1.000 7.5341E-01 ( .9997) -.500 6.9610E-01 ( .9996) -.100 6.5029E-01 ( .9996) .100 6.2738E-01 ( .9996) .500 5.8157E-01 ( .9996) 1.000 5.2431E-01 ( .9996) 64.000 -1.000 1.9513E-01 ( .9995) -.500 1.3192E-01 ( .9995) -.100 7.2030E-02 ( .9995) .100 .0000E+00 ( 1.0000) .500 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) ============================================= ==== 2 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 5b, Cloud C.1 Scattering, Ref. GS, Table 20 **************************************************************************************************** No. streams = 48 No. computational layers = 1 3 User optical depths : 3.2000 12.8000 48.0000 6 User polar angle cosines : -1.00000 -.50000 -.10000 .10000 .50000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 64.0000 64.0000 .90000 .20938 51.9396 51.9396 .87678 .8480 Layer Phase Function Moments 1 1.000000 .848000 .776600 .652571 .581667 .535182 .496692 .478467 .462294 .447053 .442190 .428522 .424600 .415889 .408655 .403323 .395697 .389314 .384027 .375897 .371488 .363744 .358356 .351894 .345592 .339706 .333453 .327255 .321561 .315051 .309590 .303222 .297615 .291597 .285812 .280056 .274301 .268600 .263000 .257342 .251864 .246313 .240906 .235506 .230146 .224912 .219645 .214547 .209381 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) 3.2000 1.2806E-01 1.7477E+00 2.7049E-01 3.1013E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 12.8000 8.6732E-06 2.3398E-01 3.7426E-02 4.5267E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 48.0000 4.4773E-21 6.3835E-05 1.0290E-05 1.2502E-05 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 3.200 -1.000 1.5021E+01 ( 1.0000) -.500 2.1878E-01 ( 1.0000) -.100 1.3652E-01 ( 1.0000) .100 1.1400E-01 ( 1.0000) .500 8.7124E-02 ( 1.0000) 1.000 8.5503E-02 ( 1.0000) 12.800 -1.000 1.8509E-01 ( 1.0000) -.500 4.9273E-02 ( 1.0000) -.100 2.6545E-02 ( 1.0000) .100 2.0216E-02 ( 1.0000) .500 1.2966E-02 ( 1.0000) 1.000 9.5124E-03 ( 1.0000) 48.000 -1.000 3.4129E-05 ( 1.0000) -.500 1.3992E-05 ( 1.0000) -.100 7.4705E-06 ( 1.0000) .100 5.6561E-06 ( 1.0000) .500 3.5825E-06 ( 1.0000) 1.000 2.5786E-06 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6a: No Scattering; Source = Beam ; Bottom Albedo = 0 **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 .0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .0000 .0000 .00000 .00000 .0000 .0000 .00000 .8480 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 .0000E+00 .0000E+00 2.0000E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .0000 1.0000E+02 .0000E+00 .0000E+00 2.0000E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******* WARNING >>>>>> PRTFIN--all intensities zero or negative ******* WARNING >>>>>> PRTFIN--all intensities near underflowing ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6b: No Scattering; Source = Beam ; Bottom Albedo = 0 **************************************************************************************************** No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 .0000E+00 .0000E+00 2.0000E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 .0000E+00 .0000E+00 7.3576E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 .0000E+00 .0000E+00 2.7067E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******* WARNING >>>>>> PRTFIN--all intensities zero or negative ******* WARNING >>>>>> PRTFIN--all intensities near underflowing ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) .500 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) 1.000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6c: No Scattering; Source = Beam ; Bottom Albedo=0.5 Lambert **************************************************************************************************** No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .5000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 .0000E+00 1.4845E+00 2.0201E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 .0000E+00 2.9991E+00 7.7996E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 .0000E+00 6.7668E+00 4.0601E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 9.7788E-05 ( 1.0000) 1.000 7.9239E-01 ( 1.0000) .500 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 1.4513E-02 ( 1.0000) 1.000 1.3064E+00 ( 1.0000) 1.000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 2.1539E+00 ( 1.0000) 1.000 2.1539E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6d: No Scattering; Source = Beam ; Bottom Albedo = Non-Lambert **************************************************************************************************** No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 .00969 .00678 .00475 .00332 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 .0000E+00 6.1910E-01 2.0090E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 .0000E+00 1.3324E+00 7.5761E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 .0000E+00 3.4779E+00 3.6519E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 5.0735E-05 ( 1.0000) 1.000 2.7240E-01 ( 1.0000) .500 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 7.5297E-03 ( 1.0000) 1.000 4.4911E-01 ( 1.0000) 1.000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 1.1175E+00 ( 1.0000) 1.000 7.4045E-01 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6e: No Scattering; Source = Beam , Bottom Emission; Bott Alb = Non-Lambert **************************************************************************************************** No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 .00969 .00678 .00475 .00332 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 300.00 Top temperature = .00 Top emissivity = 1.0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 .000 .000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 .0000E+00 8.2884E+01 3.1055E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 .0000E+00 1.6521E+02 3.1102E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 .0000E+00 3.5989E+02 6.7953E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 3.2050E-03 ( 1.0000) 1.000 4.6484E+01 ( 1.0000) .500 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 4.7567E-01 ( 1.0000) 1.000 7.6639E+01 ( 1.0000) 1.000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 7.0595E+01 ( 1.0000) 1.000 1.2636E+02 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6f: No Scattering; Source = Beam , Bottom+Top Emission; Bott Alb = Non-Lambert **************************************************************************************************** No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E+01 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 .00969 .00678 .00475 .00332 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 300.00 Top temperature = 250.00 Top emissivity = 1.0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 .000 .000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.2150E+02 8.5299E+01 9.5700E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 1.4249E+02 1.7034E+02 5.2925E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 7.0531E+01 3.7284E+02 8.0792E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 1.0234E+02 ( 1.0000) -.100 1.0234E+02 ( 1.0000) .100 3.5888E-03 ( 1.0000) 1.000 4.7512E+01 ( 1.0000) .500 -1.000 6.2070E+01 ( 1.0000) -.100 6.8953E-01 ( 1.0000) .100 5.3263E-01 ( 1.0000) 1.000 7.8334E+01 ( 1.0000) 1.000 -1.000 3.7647E+01 ( 1.0000) -.100 4.6460E-03 ( 1.0000) .100 7.9049E+01 ( 1.0000) 1.000 1.2915E+02 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6g: No Scattering; Source = Beam , Bottom+Top+Internal Emission; Bott Alb = Non-Lambert **************************************************************************************************** ******* WARNING >>>>>> CHEKIN--vertical temperature step may be too large for good accuracy No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E+01 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 .00969 .00678 .00475 .00332 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 300.00 Top temperature = 250.00 Top emissivity = 1.0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .00000 .00000 1.0000 1.0000 .00000 .8480 250.000 300.000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.2150E+02 3.3611E+02 5.8134E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 3.0477E+02 4.1398E+02 1.2862E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 3.6363E+02 4.4305E+02 -1.7112E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 1.0234E+02 ( 1.0000) -.100 1.0234E+02 ( 1.0000) .100 7.8073E+01 ( 1.0000) 1.000 1.1713E+02 ( 1.0000) .500 -1.000 9.7875E+01 ( 1.0000) -.100 1.0105E+02 ( 1.0000) .100 1.1578E+02 ( 1.0000) 1.000 1.3611E+02 ( 1.0000) 1.000 -1.000 1.1006E+02 ( 1.0000) -.100 1.3863E+02 ( 1.0000) .100 1.3337E+02 ( 1.0000) 1.000 1.4286E+02 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 6h: No Scattering; Source = Beam , Bottom+Top+Internal Emission; Bott Alb = Non-Lambert **************************************************************************************************** ******* WARNING >>>>>> CHEKIN--vertical temperature step may be too large for good accuracy No. streams = 16 No. computational layers = 1 3 User optical depths : .0000 1.0000 10.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E+01 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 .00969 .00678 .00475 .00332 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 300.00 Top temperature = 250.00 Top emissivity = 1.0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 10.0000 10.0000 .00000 .00000 10.0000 10.0000 .00000 .8480 250.000 300.000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.2150E+02 2.3735E+02 4.2378E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 2.5545E+02 2.6113E+02 6.1983E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 10.0000 2.0612E-07 4.4344E+02 4.5620E+02 -3.1772E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 90.0 .000 -1.000 1.0234E+02 ( 1.0000) -.100 1.0234E+02 ( 1.0000) .100 7.1262E+01 ( 1.0000) 1.000 7.8074E+01 ( 1.0000) 1.000 -1.000 8.4999E+01 ( 1.0000) -.100 7.7319E+01 ( 1.0000) .100 7.8831E+01 ( 1.0000) 1.000 8.5642E+01 ( 1.0000) 10.000 -1.000 1.3863E+02 ( 1.0000) -.100 1.4544E+02 ( 1.0000) .100 1.4409E+02 ( 1.0000) 1.000 1.4555E+02 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 7a: Absorption + Henyey-Greenstein Scattering, All Sources , Bottom Albedo = 0 **************************************************************************************************** ******* WARNING >>>>>> CHEKIN--vertical temperature step may be too large for good accuracy No. streams = 12 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 2 User azimuthal angles : .00 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+02 Bottom albedo (Lambertian) = .0000 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 320.00 Top temperature = 100.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .50000 .06872 .9656 .9656 .48221 .8000 300.000 200.000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.1983E+02 4.2957E+02 -8.0427E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 3.5410E+02 4.4702E+02 2.5159E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 3.0133E+02 5.9458E+02 7.1596E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 90.0 .000 -1.000 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) -.100 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) .100 1.4164E+02 1.2778E+02 ( 1.0000) ( 1.0000) 1.000 1.4842E+02 1.4842E+02 ( 1.0000) ( 1.0000) .500 -1.000 1.0635E+02 1.0635E+02 ( 1.0000) ( 1.0000) -.100 1.3422E+02 1.0624E+02 ( 1.0000) ( 1.0000) .100 1.0251E+02 9.4135E+01 ( 1.0000) ( 1.0000) 1.000 1.5886E+02 1.5886E+02 ( 1.0000) ( 1.0000) 1.000 -1.000 9.6416E+01 9.6416E+01 ( 1.0000) ( 1.0000) -.100 8.8783E+01 7.4865E+01 ( 1.0000) ( 1.0000) .100 1.8926E+02 1.8926E+02 ( 1.0000) ( 1.0000) 1.000 1.8926E+02 1.8926E+02 ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 7b: Absorption + Henyey-Greenstein Scattering, All Sources , Bottom Albedo = 1 **************************************************************************************************** ******* WARNING >>>>>> CHEKIN--vertical temperature step may be too large for good accuracy No. streams = 12 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 2 User azimuthal angles : .00 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+02 Bottom albedo (Lambertian) = 1.0000 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 320.00 Top temperature = 100.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .50000 .06872 .9656 .9656 .48221 .8000 300.000 200.000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.1983E+02 3.1256E+02 -1.6836E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 3.5055E+02 2.6813E+02 1.0125E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 2.9206E+02 3.0560E+02 4.0933E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 90.0 .000 -1.000 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) -.100 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) .100 1.3584E+02 1.2198E+02 ( 1.0000) ( 1.0000) 1.000 9.5659E+01 9.5659E+01 ( 1.0000) ( 1.0000) .500 -1.000 1.0597E+02 1.0597E+02 ( 1.0000) ( 1.0000) -.100 1.2878E+02 1.0080E+02 ( 1.0000) ( 1.0000) .100 9.0004E+01 8.1625E+01 ( 1.0000) ( 1.0000) 1.000 8.8762E+01 8.8762E+01 ( 1.0000) ( 1.0000) 1.000 -1.000 9.5365E+01 9.5365E+01 ( 1.0000) ( 1.0000) -.100 7.4755E+01 6.0837E+01 ( 1.0000) ( 1.0000) .100 9.7274E+01 9.7274E+01 ( 1.0000) ( 1.0000) 1.000 9.7274E+01 9.7274E+01 ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 7c: Absorption + Henyey-Greenstein Scattering, All Sources , Bottom Albedo = BDR Function **************************************************************************************************** ******* WARNING >>>>>> CHEKIN--vertical temperature step may be too large for good accuracy No. streams = 12 No. computational layers = 1 3 User optical depths : .0000 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 2 User azimuthal angles : .00 90.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+02 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+02 Legendre coeffs of bottom bidirectional reflectivity : 1.00000 .70000 .49000 .34300 .24010 .16807 .11765 .08235 .05765 .04035 .02825 .01977 .01384 Thermal emission in wavenumber interval : .0000 50000.0000 Bottom temperature = 320.00 Top temperature = 100.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .50000 .06872 .9656 .9656 .48221 .8000 300.000 200.000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+02 3.1983E+02 4.0704E+02 -9.8569E+01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 3.6788E+01 3.5325E+02 4.1106E+02 2.1772E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 1.3534E+01 2.9858E+02 5.3050E+02 6.2394E+02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 90.0 .000 -1.000 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) -.100 1.0180E+02 1.0180E+02 ( 1.0000) ( 1.0000) .100 1.4035E+02 1.2634E+02 ( 1.0000) ( 1.0000) 1.000 1.4053E+02 1.4053E+02 ( 1.0000) ( 1.0000) .500 -1.000 1.0626E+02 1.0626E+02 ( 1.0000) ( 1.0000) -.100 1.3294E+02 1.0482E+02 ( 1.0000) ( 1.0000) .100 9.9279E+01 9.0269E+01 ( 1.0000) ( 1.0000) 1.000 1.4862E+02 1.4862E+02 ( 1.0000) ( 1.0000) 1.000 -1.000 9.6139E+01 9.6139E+01 ( 1.0000) ( 1.0000) -.100 8.4588E+01 6.9729E+01 ( 1.0000) ( 1.0000) .100 1.5490E+02 1.3796E+02 ( 1.0000) ( 1.0000) 1.000 1.7627E+02 1.7627E+02 ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 8A: Ref. OS, Table 1, Line 4 (Two Inhomogeneous Layers) **************************************************************************************************** No. streams = 8 No. computational layers = 2 3 User optical depths : .0000 .2500 .5000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 1 User azimuthal angles : 60.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E-01 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .2500 .2500 .50000 .00000 .2500 .2500 .50000 .0000 2 .2500 .5000 .30000 .00000 .2500 .5000 .30000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 2 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 1.0000E+00 9.2963E-02 1.1247E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .2500 .0000E+00 7.2224E-01 2.7895E-02 6.5182E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 .0000E+00 5.1313E-01 2.1986E-10 5.6336E-01 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 .000 -1.000 3.1831E-01 ( 1.0000) -.200 3.1831E-01 ( 1.0000) .200 5.6257E-02 ( 1.0000) 1.000 1.9442E-02 ( 1.0000) .250 -1.000 2.6271E-01 ( 1.0000) -.200 1.3695E-01 ( 1.0000) .200 1.8491E-02 ( 1.0000) 1.000 5.5219E-03 ( 1.0000) .500 -1.000 2.1001E-01 ( 1.0000) -.200 5.6038E-02 ( 1.0000) .200 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 8b: Ref. OS, Table 1, Line 1 (Two Inhomogeneous Layers) **************************************************************************************************** No. streams = 8 No. computational layers = 2 3 User optical depths : .0000 .2500 .5000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 1 User azimuthal angles : 60.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E-01 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .2500 .2500 .80000 .00000 .2500 .2500 .80000 .0000 2 .2500 .5000 .95000 .00000 .2500 .5000 .95000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 2 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 1.0000E+00 2.2514E-01 5.1269E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .2500 .0000E+00 7.9533E-01 1.2635E-01 3.5665E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 .0000E+00 6.5042E-01 -8.5670E-08 5.6810E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 .000 -1.000 3.1831E-01 ( 1.0000) -.200 3.1831E-01 ( 1.0000) .200 1.2369E-01 ( 1.0000) 1.000 4.9558E-02 ( 1.0000) .250 -1.000 2.7750E-01 ( 1.0000) -.200 1.8395E-01 ( 1.0000) .200 8.3570E-02 ( 1.0000) 1.000 2.5057E-02 ( 1.0000) .500 -1.000 2.4073E-01 ( 1.0000) -.200 1.2929E-01 ( 1.0000) .200 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 8c: Ref. OS, Table 1, Line 13 (Two Inhomogeneous Layers) **************************************************************************************************** No. streams = 8 No. computational layers = 2 3 User optical depths : .0000 1.0000 3.0000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 1 User azimuthal angles : 60.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E-01 Bottom albedo (Lambertian) = .0000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .80000 .00000 1.0000 1.0000 .80000 .0000 2 2.0000 3.0000 .95000 .00000 2.0000 3.0000 .95000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 2 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 1.0000E+00 3.7858E-01 5.6509E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 .0000E+00 4.8616E-01 2.4340E-01 2.7670E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 3.0000 .0000E+00 1.5998E-01 1.1863E-09 1.3568E-02 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 .000 -1.000 3.1831E-01 ( 1.0000) -.200 3.1831E-01 ( 1.0000) .200 1.4933E-01 ( 1.0000) 1.000 1.0477E-01 ( 1.0000) 1.000 -1.000 1.8902E-01 ( 1.0000) -.200 9.8816E-02 ( 1.0000) .200 9.6519E-02 ( 1.0000) 1.000 6.5445E-02 ( 1.0000) 3.000 -1.000 6.8476E-02 ( 1.0000) -.200 2.9670E-02 ( 1.0000) .200 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 9a: Ref. DGIS, Tables VI-VII, beta=l=0 (multiple inhomogeneous layers) **************************************************************************************************** No. streams = 8 No. computational layers = 6 5 User optical depths : .0000 1.0500 2.1000 6.0000 21.0000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 1 User azimuthal angles : 60.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E-01 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .65000 .00000 1.0000 1.0000 .65000 .0000 2 2.0000 3.0000 .70000 .00000 2.0000 3.0000 .70000 .0000 3 3.0000 6.0000 .75000 .00000 3.0000 6.0000 .75000 .0000 4 4.0000 10.0000 .80000 .00000 4.0000 10.0000 .80000 .0000 5 5.0000 15.0000 .85000 .00000 5.0000 15.0000 .85000 .0000 6 6.0000 21.0000 .90000 .00000 6.0000 21.0000 .90000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 2 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 3 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 4 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 5 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 6 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 1.0000E+00 2.2797E-01 8.8212E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0500 .0000E+00 3.5515E-01 8.7510E-02 2.3237E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 2.1000 .0000E+00 1.4427E-01 3.6182E-02 9.3344E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 6.0000 .0000E+00 6.7145E-03 2.1929E-03 3.9278E-03 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 21.0000 .0000E+00 6.1697E-07 9.3618E-15 1.0250E-07 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 .000 -1.000 3.1831E-01 ( 1.0000) -.200 3.1831E-01 ( 1.0000) .200 9.9892E-02 ( 1.0000) 1.000 5.9134E-02 ( 1.0000) 1.050 -1.000 1.5351E-01 ( 1.0000) -.200 5.0953E-02 ( 1.0000) .200 3.6701E-02 ( 1.0000) 1.000 2.3190E-02 ( 1.0000) 2.100 -1.000 7.0661E-02 ( 1.0000) -.200 2.0912E-02 ( 1.0000) .200 1.4854E-02 ( 1.0000) 1.000 9.7231E-03 ( 1.0000) 6.000 -1.000 3.7278E-03 ( 1.0000) -.200 1.0882E-03 ( 1.0000) .200 8.8332E-04 ( 1.0000) 1.000 5.9474E-04 ( 1.0000) 21.000 -1.000 2.8766E-07 ( 1.0000) -.200 1.0592E-07 ( 1.0000) .200 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 9b: Ref. DGIS, Tables VI-VII, beta=0,l=8 (multiple inhomogeneous layers) **************************************************************************************************** No. streams = 8 No. computational layers = 6 5 User optical depths : .0000 1.0500 2.1000 6.0000 21.0000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 1 User azimuthal angles : 60.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = .000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 3.183E-01 Bottom albedo (Lambertian) = .0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .65000 .00000 1.0000 1.0000 .65000 .6697 2 2.0000 3.0000 .70000 .00000 2.0000 3.0000 .70000 .6697 3 3.0000 6.0000 .75000 .00000 3.0000 6.0000 .75000 .6697 4 4.0000 10.0000 .80000 .00000 4.0000 10.0000 .80000 .6697 5 5.0000 15.0000 .85000 .00000 5.0000 15.0000 .85000 .6697 6 6.0000 21.0000 .90000 .00000 6.0000 21.0000 .90000 .6697 Layer Phase Function Moments 1 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 2 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 3 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 4 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 5 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 6 1.000000 .669720 .312678 .096296 .024683 .004295 .000516 .000045 .000003 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 .0000E+00 1.0000E+00 1.0008E-01 8.0458E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0500 .0000E+00 4.5236E-01 4.5201E-02 2.5533E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 2.1000 .0000E+00 2.3647E-01 2.4194E-02 1.3098E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 6.0000 .0000E+00 2.7648E-02 4.1602E-03 1.3623E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 21.0000 .0000E+00 7.4185E-05 -4.9076E-12 1.2202E-05 ( 1.0000) ( 1.0000) (*********) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 .000 -1.000 3.1831E-01 ( 1.0000) -.200 3.1831E-01 ( 1.0000) .200 7.3920E-02 ( 1.0000) 1.000 1.3277E-02 ( 1.0000) 1.050 -1.000 1.9661E-01 ( 1.0000) -.200 5.9237E-02 ( 1.0000) .200 3.0023E-02 ( 1.0000) 1.000 7.0557E-03 ( 1.0000) 2.100 -1.000 1.1548E-01 ( 1.0000) -.200 3.0181E-02 ( 1.0000) .200 1.5267E-02 ( 1.0000) 1.000 4.0693E-03 ( 1.0000) 6.000 -1.000 1.4618E-02 ( 1.0000) -.200 3.8559E-03 ( 1.0000) .200 2.3830E-03 ( 1.0000) 1.000 7.7789E-04 ( 1.0000) 21.000 -1.000 3.3774E-05 ( 1.0000) -.200 1.2086E-05 ( 1.0000) .200 .0000E+00 ( 1.0000) 1.000 .0000E+00 ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 9c: Generalization of 9A to include all possible complexity **************************************************************************************************** No. streams = 8 No. computational layers = 6 5 User optical depths : .0000 1.0500 2.1000 6.0000 21.0000 4 User polar angle cosines : -1.00000 -.20000 .20000 1.00000 3 User azimuthal angles : 60.00 120.00 180.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .5000 Thermal emission in wavenumber interval : 999.0000 1000.0000 Bottom temperature = 700.00 Top temperature = 550.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .65000 .00000 1.0000 1.0000 .65000 .1429 600.000 2 2.0000 3.0000 .70000 .00004 1.9999 2.9999 .69999 .2857 610.000 3 3.0000 6.0000 .75000 .00114 2.9974 5.9974 .74979 .4286 620.000 4 4.0000 10.0000 .80000 .01137 3.9636 9.9610 .79816 .5714 630.000 5 5.0000 15.0000 .85000 .06776 4.7120 14.6730 .84083 .7143 640.000 6 6.0000 21.0000 .90000 .29136 4.4267 19.0997 .86446 .8571 650.000 660.000 Layer Phase Function Moments 1 1.000000 .142857 .020408 .002915 .000416 .000059 .000008 .000001 .000000 2 1.000000 .285714 .081633 .023324 .006664 .001904 .000544 .000155 .000044 3 1.000000 .428571 .183673 .078717 .033736 .014458 .006196 .002656 .001138 4 1.000000 .571429 .326531 .186589 .106622 .060927 .034815 .019895 .011368 5 1.000000 .714286 .510204 .364432 .260308 .185934 .132810 .094865 .067760 6 1.000000 .857143 .734694 .629738 .539775 .462664 .396569 .339917 .291357 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.5708E+00 6.0921E+00 4.6840E+00 3.4959E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0001) 1.0500 1.9235E-01 4.9727E+00 4.2437E+00 8.8110E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( .9999) 2.1000 2.3555E-02 4.4661E+00 4.1693E+00 3.5001E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( .9999) 6.0000 9.6513E-06 4.2270E+00 4.3064E+00 1.9447E-02 ( 1.0000) ( .9999) ( .9999) ( 1.0052) 21.0000 9.0313E-19 4.7374E+00 5.1148E+00 7.1503E-02 ( 1.0000) ( 1.0000) ( .9999) ( .9996) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 120.0 180.0 .000 -1.000 1.9392E+00 1.9392E+00 1.9392E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.200 1.9392E+00 1.9392E+00 1.9392E+00 ( 1.0000) ( 1.0000) ( 1.0000) .200 1.6185E+00 1.5789E+00 1.5655E+00 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 1.4387E+00 1.4387E+00 1.4387E+00 ( 1.0000) ( 1.0000) ( 1.0000) 1.050 -1.000 1.6676E+00 1.6676E+00 1.6676E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.200 1.4445E+00 1.4292E+00 1.4244E+00 ( 1.0000) ( 1.0000) ( 1.0000) .200 1.3834E+00 1.3732E+00 1.3703E+00 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 1.3389E+00 1.3389E+00 1.3389E+00 ( 1.0000) ( 1.0000) ( 1.0000) 2.100 -1.000 1.4851E+00 1.4851E+00 1.4851E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.200 1.3501E+00 1.3458E+00 1.3447E+00 ( 1.0000) ( 1.0000) ( 1.0000) .200 1.3308E+00 1.3292E+00 1.3287E+00 ( 1.0000) ( 1.0000) ( 1.0000) 1.000 1.3279E+00 1.3279E+00 1.3279E+00 ( 1.0000) ( 1.0000) ( 1.0000) 6.000 -1.000 1.3451E+00 1.3451E+00 1.3451E+00 ( .9999) ( .9999) ( .9999) -.200 1.3512E+00 1.3512E+00 1.3512E+00 ( .9999) ( .9999) ( .9999) .200 1.3597E+00 1.3597E+00 1.3597E+00 ( .9999) ( .9999) ( .9999) 1.000 1.3791E+00 1.3791E+00 1.3791E+00 ( .9999) ( .9999) ( .9999) 21.000 -1.000 1.4892E+00 1.4892E+00 1.4892E+00 ( 1.0000) ( 1.0000) ( 1.0000) -.200 1.5426E+00 1.5426E+00 1.5426E+00 ( .9999) ( .9999) ( .9999) .200 1.6281E+00 1.6281E+00 1.6281E+00 ( .9999) ( .9999) ( .9999) 1.000 1.6281E+00 1.6281E+00 1.6281E+00 ( .9999) ( .9999) ( .9999) **************************************************************************************************** DISORT: Test Case No. 10a: like 9c, USRANG = True **************************************************************************************************** No. streams = 4 No. computational layers = 6 3 User optical depths : .0000 2.1000 21.0000 4 User polar angle cosines : -.78868 -.21132 .21132 .78868 2 User azimuthal angles : 60.00 120.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .5000 Thermal emission in wavenumber interval : 999.0000 1000.0000 Bottom temperature = 700.00 Top temperature = 550.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .65000 .00042 .9997 .9997 .64991 .1429 600.000 2 2.0000 3.0000 .70000 .00666 1.9907 2.9904 .69859 .2857 610.000 3 3.0000 6.0000 .75000 .03374 2.9241 5.9145 .74351 .4286 620.000 4 4.0000 10.0000 .80000 .10662 3.6588 9.5733 .78135 .5714 630.000 5 5.0000 15.0000 .85000 .26031 3.8937 13.4670 .80738 .7143 640.000 6 6.0000 21.0000 .90000 .53978 3.0852 16.5522 .80552 .8571 650.000 660.000 Layer Phase Function Moments 1 1.000000 .142857 .020408 .002915 .000416 2 1.000000 .285714 .081633 .023324 .006664 3 1.000000 .428571 .183673 .078717 .033736 4 1.000000 .571429 .326531 .186589 .106622 5 1.000000 .714286 .510204 .364432 .260308 6 1.000000 .857143 .734694 .629738 .539775 <----------------------- FLUXES -----------------------> Optical Compu Downward Downward Downward Upward Mean Planck d(Net Flux) Depth Layer Direct Diffuse Total Diffuse Net Intensity Source / d(Op Dep) .0000 1 1.571E+00 6.092E+00 7.663E+00 4.688E+00 2.975E+00 1.985E+00 1.191E+00 3.495E+00 2.1000 2 2.356E-02 4.470E+00 4.493E+00 4.170E+00 3.233E-01 1.366E+00 1.273E+00 3.512E-01 21.0000 6 9.031E-19 4.738E+00 4.738E+00 5.115E+00 -3.769E-01 1.574E+00 1.518E+00 7.131E-02 ********* I N T E N S I T I E S ********* Polar Azimuth angles (degrees) Optical Angle Depth Cosine 60.00 120.00 .0000 -.7887 1.939E+00 1.939E+00 -.2113 1.939E+00 1.939E+00 .2113 1.615E+00 1.576E+00 .7887 1.469E+00 1.458E+00 60.00 120.00 2.1000 -.7887 1.446E+00 1.438E+00 -.2113 1.351E+00 1.347E+00 .2113 1.331E+00 1.329E+00 .7887 1.327E+00 1.326E+00 60.00 120.00 21.0000 -.7887 1.499E+00 1.499E+00 -.2113 1.541E+00 1.541E+00 .2113 1.628E+00 1.628E+00 .7887 1.628E+00 1.628E+00 **************************************************************************************************** DISORT: Test Case No. 10b: like 9C, USRANG = False **************************************************************************************************** No. streams = 4 No. computational layers = 6 3 User optical depths : .0000 2.1000 21.0000 4 User polar angle cosines : -.78868 -.21132 .21132 .78868 2 User azimuthal angles : 60.00 120.00 Boundary condition flag: IBCND = 0 Incident beam with intensity = 3.142E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.000E+00 Bottom albedo (Lambertian) = .5000 Thermal emission in wavenumber interval : 999.0000 1000.0000 Bottom temperature = 700.00 Top temperature = 550.00 Top emissivity = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor Temperature 1 1.0000 1.0000 .65000 .00042 .9997 .9997 .64991 .1429 600.000 2 2.0000 3.0000 .70000 .00666 1.9907 2.9904 .69859 .2857 610.000 3 3.0000 6.0000 .75000 .03374 2.9241 5.9145 .74351 .4286 620.000 4 4.0000 10.0000 .80000 .10662 3.6588 9.5733 .78135 .5714 630.000 5 5.0000 15.0000 .85000 .26031 3.8937 13.4670 .80738 .7143 640.000 6 6.0000 21.0000 .90000 .53978 3.0852 16.5522 .80552 .8571 650.000 660.000 Layer Phase Function Moments 1 1.000000 .142857 .020408 .002915 .000416 2 1.000000 .285714 .081633 .023324 .006664 3 1.000000 .428571 .183673 .078717 .033736 4 1.000000 .571429 .326531 .186589 .106622 5 1.000000 .714286 .510204 .364432 .260308 6 1.000000 .857143 .734694 .629738 .539775 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.5708E+00 6.0921E+00 4.6883E+00 3.4946E+00 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 2.1000 2.3555E-02 4.4699E+00 4.1702E+00 3.5123E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 21.0000 9.0313E-19 4.7383E+00 5.1153E+00 7.1307E-02 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine 60.0 120.0 .000 -.789 1.9392E+00 1.9392E+00 ( 1.0000) ( 1.0000) -.211 1.9392E+00 1.9392E+00 ( 1.0000) ( 1.0000) .211 1.6151E+00 1.5763E+00 ( 1.0000) ( 1.0000) .789 1.4691E+00 1.4577E+00 ( 1.0000) ( 1.0000) 2.100 -.789 1.4458E+00 1.4378E+00 ( 1.0000) ( 1.0000) -.211 1.3512E+00 1.3468E+00 ( 1.0000) ( 1.0000) .211 1.3308E+00 1.3292E+00 ( 1.0000) ( 1.0000) .789 1.3268E+00 1.3265E+00 ( 1.0000) ( 1.0000) 21.000 -.789 1.4994E+00 1.4994E+00 ( 1.0000) ( 1.0000) -.211 1.5412E+00 1.5412E+00 ( 1.0000) ( 1.0000) .211 1.6282E+00 1.6282E+00 ( 1.0000) ( 1.0000) .789 1.6282E+00 1.6282E+00 ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 11a: One Isotropic-Scattering Layer **************************************************************************************************** No. streams = 16 No. computational layers = 1 4 User optical depths : .0000 .0500 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 2 User azimuthal angles : .00 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 1.000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.592E-01 Bottom albedo (Lambertian) = .5000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .90000 .00000 1.0000 1.0000 .90000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <----------------------- FLUXES -----------------------> Optical Compu Downward Downward Downward Upward Mean Planck d(Net Flux) Depth Layer Direct Diffuse Total Diffuse Net Intensity Source / d(Op Dep) .0000 1 5.000E-01 5.000E-01 1.000E+00 5.014E-01 4.986E-01 2.457E-01 .000E+00 3.088E-01 .0500 1 4.524E-01 5.170E-01 9.694E-01 4.861E-01 4.833E-01 2.407E-01 .000E+00 3.025E-01 .5000 1 1.839E-01 5.465E-01 7.304E-01 3.653E-01 3.651E-01 1.788E-01 .000E+00 2.247E-01 1.0000 1 6.767E-02 4.724E-01 5.400E-01 2.700E-01 2.700E-01 1.256E-01 .000E+00 1.578E-01 ********* I N T E N S I T I E S ********* Polar Azimuth angles (degrees) Optical Angle Depth Cosine .00 90.00 .0000 -1.0000 1.592E-01 1.592E-01 -.1000 1.592E-01 1.592E-01 .1000 2.100E-01 2.100E-01 1.0000 1.409E-01 1.409E-01 .00 90.00 .0500 -1.0000 1.621E-01 1.621E-01 -.1000 1.827E-01 1.827E-01 .1000 2.040E-01 2.040E-01 1.0000 1.368E-01 1.368E-01 .00 90.00 .5000 -1.0000 1.708E-01 1.708E-01 -.1000 1.721E-01 1.721E-01 .1000 1.505E-01 1.505E-01 1.0000 1.065E-01 1.065E-01 .00 90.00 1.0000 -1.0000 1.564E-01 1.564E-01 -.1000 1.224E-01 1.224E-01 .1000 8.595E-02 8.595E-02 1.0000 8.595E-02 8.595E-02 **************************************************************************************************** DISORT: Test Case No. 11b: Same as 11a but treated as multiple layers **************************************************************************************************** No. streams = 16 No. computational layers = 3 4 User optical depths : .0000 .0500 .5000 1.0000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 2 User azimuthal angles : .00 90.00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 1.000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = 1.592E-01 Bottom albedo (Lambertian) = .5000 Does not use delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .0500 .0500 .90000 .00000 .0500 .0500 .90000 .0000 2 .4500 .5000 .90000 .00000 .4500 .5000 .90000 .0000 3 .5000 1.0000 .90000 .00000 .5000 1.0000 .90000 .0000 Layer Phase Function Moments 1 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 2 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 3 1.000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 .000000 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 5.0000E-01 5.0000E-01 5.0138E-01 3.0876E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .0500 4.5242E-01 5.1698E-01 4.8608E-01 3.0247E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) .5000 1.8394E-01 5.4647E-01 3.6529E-01 2.2469E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 1.0000 6.7668E-02 4.7236E-01 2.7002E-01 1.5778E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 90.0 .000 -1.000 1.5915E-01 1.5915E-01 ( 1.0000) ( 1.0000) -.100 1.5915E-01 1.5915E-01 ( 1.0000) ( 1.0000) .100 2.1004E-01 2.1004E-01 ( 1.0000) ( 1.0000) 1.000 1.4086E-01 1.4086E-01 ( 1.0000) ( 1.0000) .050 -1.000 1.6208E-01 1.6208E-01 ( 1.0000) ( 1.0000) -.100 1.8270E-01 1.8270E-01 ( 1.0000) ( 1.0000) .100 2.0398E-01 2.0398E-01 ( 1.0000) ( 1.0000) 1.000 1.3684E-01 1.3684E-01 ( 1.0000) ( 1.0000) .500 -1.000 1.7076E-01 1.7076E-01 ( 1.0000) ( 1.0000) -.100 1.7214E-01 1.7214E-01 ( 1.0000) ( 1.0000) .100 1.5049E-01 1.5049E-01 ( 1.0000) ( 1.0000) 1.000 1.0650E-01 1.0650E-01 ( 1.0000) ( 1.0000) 1.000 -1.000 1.5640E-01 1.5640E-01 ( 1.0000) ( 1.0000) -.100 1.2236E-01 1.2236E-01 ( 1.0000) ( 1.0000) .100 8.5949E-02 8.5949E-02 ( 1.0000) ( 1.0000) 1.000 8.5949E-02 8.5949E-02 ( 1.0000) ( 1.0000) **************************************************************************************************** DISORT: Test Case No. 12a: Overhead Beam Striking Absorbing/Scattering Medium **************************************************************************************************** No. streams = 20 No. computational layers = 1 4 User optical depths : .0000 10.0000 19.9000 20.1000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 1.000E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 20.1000 20.1000 .50000 .12158 18.8782 18.8782 .46764 .9000 Layer Phase Function Moments 1 1.000000 .900000 .810000 .729000 .656100 .590490 .531441 .478297 .430467 .387420 .348678 .313810 .282429 .254186 .228768 .205891 .185302 .166772 .150095 .135085 .121577 <----------------------- FLUXES -----------------------> Optical Compu Downward Downward Downward Upward Mean Planck d(Net Flux) Depth Layer Direct Diffuse Total Diffuse Net Intensity Source / d(Op Dep) .0000 1 1.000E+00 -4.768E-07 1.000E+00 6.484E-03 9.935E-01 8.089E-02 .000E+00 5.083E-01 10.0000 1 4.540E-05 3.170E-03 3.216E-03 2.756E-05 3.188E-03 3.063E-04 .000E+00 1.925E-03 19.9000 1 2.278E-09 7.453E-06 7.455E-06 5.476E-06 1.980E-06 1.491E-06 .000E+00 9.368E-06 20.1000 1 1.865E-09 6.628E-06 6.630E-06 6.630E-06 1.364E-12 1.720E-06 .000E+00 1.081E-05 ********* I N T E N S I T I E S ********* Polar Azimuth angles (degrees) Optical Angle Depth Cosine .00 .0000 -1.0000 .000E+00 -.1000 .000E+00 .1000 4.986E-03 1.0000 -1.264E-03 .00 10.0000 -1.0000 7.572E-03 -.1000 4.157E-05 .1000 2.429E-05 1.0000 4.933E-06 .00 19.9000 -1.0000 1.104E-05 -.1000 2.317E-07 .1000 7.381E-07 1.0000 1.906E-06 .00 20.1000 -1.0000 9.706E-06 -.1000 3.131E-07 .1000 2.110E-06 1.0000 2.110E-06 **************************************************************************************************** DISORT: Test Case No. 12b: Same as 12a but uses shortcut for absorption optical depth .GT. 10 **************************************************************************************************** No. streams = 20 No. computational layers = 3 4 User optical depths : .0000 10.0000 19.9000 20.1000 4 User polar angle cosines : -1.00000 -.10000 .10000 1.00000 1 User azimuthal angles : .00 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 1.000E+00 and polar angle cosine = 1.00000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = 1.0000 Uses delta-M method Calculate fluxes and intensities Relative convergence criterion for azimuth series = .00E+00 Sets radiation = 0 below absorption optical depth 10 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 10.0000 10.0000 .50000 .12158 9.3921 9.3921 .46764 .9000 2 9.9000 19.9000 .50000 .12158 9.2982 18.6903 .46764 .9000 3 .2000 20.1000 .50000 .12158 .1878 18.8782 .46764 .9000 Layer Phase Function Moments 1 1.000000 .900000 .810000 .729000 .656100 .590490 .531441 .478297 .430467 .387420 .348678 .313810 .282429 .254186 .228768 .205891 .185302 .166772 .150095 .135085 .121577 2 1.000000 .900000 .810000 .729000 .656100 .590490 .531441 .478297 .430467 .387420 .348678 .313810 .282429 .254186 .228768 .205891 .185302 .166772 .150095 .135085 .121577 3 1.000000 .900000 .810000 .729000 .656100 .590490 .531441 .478297 .430467 .387420 .348678 .313810 .282429 .254186 .228768 .205891 .185302 .166772 .150095 .135085 .121577 <-------------- FLUXES --------------> Optical Downward Downward Upward d(Net Flux) Depth Direct Diffuse Diffuse / d(Op Dep) .0000 1.0000E+00 -4.7684E-07 6.4843E-03 5.0826E-01 ( 1.0000) ( 1.0000) ( 1.0000) ( 1.0000) 10.0000 4.5400E-05 3.1704E-03 2.7556E-05 1.9247E-03 ( 1.0000) ( 1.0000) ( .9998) ( 1.0000) 19.9000 2.2779E-09 7.3673E-06 2.1269E-08 4.4821E-06 ( 1.0000) ( .9885) ( .0039) ( .4785) 20.1000 1.8650E-09 6.5067E-06 -5.0360E-15 3.9147E-06 ( 1.0000) ( .9817) ( .0000) ( .3621) ******** I N T E N S I T I E S ********* Polar Azimuthal Angles (Degrees) Optical Angle Depth Cosine .0 .000 -1.000 .0000E+00 ( 1.0000) -.100 .0000E+00 ( 1.0000) .100 4.9860E-03 ( 1.0000) 1.000 -1.2643E-03 ( 1.0000) 10.000 -1.000 7.5719E-03 ( 1.0000) -.100 4.1565E-05 ( 1.0000) .100 2.4292E-05 ( 1.0000) 1.000 4.9269E-06 ( .9987) 19.900 -1.000 1.1024E-05 ( .9990) -.100 1.0286E-07 ( .4439) .100 4.5348E-08 ( .0614) 1.000 2.1953E-09 ( .0012) 20.100 -1.000 9.6920E-06 ( .9986) -.100 8.8482E-08 ( .2826) .100 .0000E+00 ( .0000) 1.000 .0000E+00 ( .0000) ============================================= ==== 7 SERIOUSLY NON-UNIT RATIOS ==== ============================================= **************************************************************************************************** DISORT: Test Case No. 13a: Albedo and Transmissivity from Shortcut, Single Layer **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User polar angle cosines : -.50000 .50000 No thermal emission Boundary condition flag: IBCND = 1 Isotropic illumination from top and bottom Bottom albedo (Lambertian) = .5000 Uses delta-M method Calculate albedo and transmissivity of medium vs. incident beam angle Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .99000 .02815 .9721 .9721 .98971 .8000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 ******* Flux Albedo and/or Transmissivity of entire medium ******** Beam Zen Ang cos(Beam Zen Ang) Albedo Transmissivity 60.0000 .500000 .54526 8.4500E-01 **************************************************************************************************** DISORT: Test Case No. 13b: Albedo and Transmissivity by Regular Method, Single Layer **************************************************************************************************** No. streams = 16 No. computational layers = 1 2 User optical depths : .0000 1.0000 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .5000 Uses delta-M method Calculate fluxes and azim-averaged intensities only Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 1.0000 1.0000 .99000 .02815 .9721 .9721 .98971 .8000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 <----------------------- FLUXES -----------------------> Optical Compu Downward Downward Downward Upward Mean Planck d(Net Flux) Depth Layer Direct Diffuse Total Diffuse Net Intensity Source / d(Op Dep) .0000 1 1.000E+00 2.384E-07 1.000E+00 5.453E-01 4.547E-01 2.546E-01 .000E+00 3.199E-02 1.0000 1 1.353E-01 7.097E-01 8.450E-01 4.225E-01 4.225E-01 2.190E-01 .000E+00 2.752E-02 **************************************************************************************************** DISORT: Test Case No. 13c: Albedo and Transmissivity from Shortcut, Multiple Layer **************************************************************************************************** No. streams = 16 No. computational layers = 2 No thermal emission Boundary condition flag: IBCND = 1 Isotropic illumination from top and bottom Bottom albedo (Lambertian) = .5000 Uses delta-M method Calculate albedo and transmissivity of medium vs. incident beam angle Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .5000 .5000 .99000 .02815 .4861 .4861 .98971 .8000 2 .5000 1.0000 .50000 .02815 .4930 .9790 .49286 .8000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 2 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 ******* Flux Albedo and/or Transmissivity of entire medium ******** Beam Zen Ang cos(Beam Zen Ang) Albedo Transmissivity 60.0000 .500000 .27620 5.0332E-01 **************************************************************************************************** DISORT: Test Case No. 13d: Albedo and Transmissivity by Regular Method, Multiple Layer **************************************************************************************************** No. streams = 16 No. computational layers = 2 2 User optical depths : .0000 1.0000 No thermal emission Boundary condition flag: IBCND = 0 Incident beam with intensity = 2.000E+00 and polar angle cosine = .50000 and azimuth angle = .00 plus isotropic incident intensity = .000E+00 Bottom albedo (Lambertian) = .5000 Uses delta-M method Calculate fluxes and azim-averaged intensities only Relative convergence criterion for azimuth series = .00E+00 <------------- Delta-M ---------------> Total Single Total Single Optical Optical Scatter Truncated Optical Optical Scatter Asymm Depth Depth Albedo Fraction Depth Depth Albedo Factor 1 .5000 .5000 .99000 .02815 .4861 .4861 .98971 .8000 2 .5000 1.0000 .50000 .02815 .4930 .9790 .49286 .8000 Layer Phase Function Moments 1 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 2 1.000000 .800000 .640000 .512000 .409600 .327680 .262144 .209715 .167772 .134218 .107374 .085899 .068720 .054976 .043980 .035184 .028148 <----------------------- FLUXES -----------------------> Optical Compu Downward Downward Downward Upward Mean Planck d(Net Flux) Depth Layer Direct Diffuse Total Diffuse Net Intensity Source / d(Op Dep) .0000 1 1.000E+00 -1.192E-07 1.000E+00 2.762E-01 7.238E-01 2.133E-01 .000E+00 2.681E-02 1.0000 2 1.353E-01 3.680E-01 5.033E-01 2.517E-01 2.517E-01 1.210E-01 .000E+00 7.604E-01