main: run brdf routines side by side main: wmic,wavenumber (2pi/wmic) 0.7600 8.267E+00 main: nreal,nim 1.330E+00 -1.300E-07 main: windspeed (m/sec) 8.0000 main: mre,mim 1.330D+00 -1.300D-07 main:(for nn=nn1) suntheta,sensorview 1.710D+01 1.353D+01 main: for Breon geometry, azimuth of sun and sensor main: phisun 180.0000 main: phiview 0.0000 main: windaz 45.0000 main: windazdiff (windaz+180-phisun) 45.0000 main: the final (1,1) reflectance term main: rspgiss from coxmunk_vfunction_plus.f90 main: rspnew from coxmunk_new.f90 main: rspbreon from breon.f90 main: rspguerin from guerin.f90 main: nn,lats,zensun,rspgiss,rspnew,rspbreon,rspguerin 1 25.14480 17.10090 0.12066 0.12066 0.13212 0.12829 2 22.78770 17.63180 0.12107 0.12107 0.13236 0.12865 3 16.75890 19.97350 0.12346 0.12346 0.13429 0.13094 4 12.48190 22.48920 0.12644 0.12644 0.13665 0.13377 5 7.35300 26.08590 0.13172 0.13172 0.14097 0.13885 6 2.17110 30.13240 0.13976 0.13976 0.14790 0.14670 7 -2.32180 33.88180 0.14974 0.14974 0.15667 0.15649 8 -7.14550 38.06720 0.16508 0.16508 0.17042 0.17164 9 -12.95030 43.23340 0.19332 0.19332 0.19609 0.19964 10 -16.81730 46.75710 0.22161 0.22161 0.22192 0.22773 11 -22.72300 52.18220 0.28921 0.28921 0.28356 0.29476 12 -27.92050 57.02640 0.39088 0.39088 0.37532 0.39512 13 -31.29330 60.16130 0.49105 0.49105 0.46450 0.49341 14 -37.79220 66.25930 0.82646 0.82646 0.75682 0.81957 15 -43.03410 71.19920 1.34283 1.34283 1.19309 1.31565 16 -46.26030 74.24110 1.87224 1.87224 1.63581 1.82466 main: the M matrix (1,1) term main: m11divnew is m11/dmodmt from coxmunk_new.f90 main: m11breon is from breon.f90 (uses hafermanmatrix.f90) main: m11guerin is from guerin.f90 main: nn,lats,zensun,m11divnew,m11breon,m11guerin 1 25.14480 17.10090 0.02012 0.02012 0.02012 2 22.78770 17.63180 0.02013 0.02013 0.02013 3 16.75890 19.97350 0.02017 0.02017 0.02017 4 12.48190 22.48920 0.02025 0.02025 0.02025 5 7.35300 26.08590 0.02041 0.02041 0.02041 6 2.17110 30.13240 0.02071 0.02071 0.02071 7 -2.32180 33.88180 0.02116 0.02116 0.02116 8 -7.14550 38.06720 0.02194 0.02194 0.02194 9 -12.95030 43.23340 0.02350 0.02350 0.02350 10 -16.81730 46.75710 0.02512 0.02512 0.02512 11 -22.72300 52.18220 0.02895 0.02895 0.02895 12 -27.92050 57.02640 0.03437 0.03437 0.03437 13 -31.29330 60.16130 0.03933 0.03933 0.03933 14 -37.79220 66.25930 0.05370 0.05370 0.05370 15 -43.03410 71.19920 0.07171 0.07171 0.07171 16 -46.26030 74.24110 0.08684 0.08684 0.08684 main: the M matrix (1,2) term main: m12divnew is m11/dmodmt from coxmunk_new.f90 main: m12breon is from breon.f90 (uses hafermanmatrix.f90) main: m12guerin is from guerin.f90 main: nn,lats,zensun,m12divnew,m12breon,m12guerin 1 25.14480 17.10090 -0.00223 -0.00223 -0.00223 2 22.78770 17.63180 -0.00237 -0.00237 -0.00237 3 16.75890 19.97350 -0.00308 -0.00308 -0.00308 4 12.48190 22.48920 -0.00395 -0.00395 -0.00395 5 7.35300 26.08590 -0.00541 -0.00541 -0.00541 6 2.17110 30.13240 -0.00740 -0.00740 -0.00740 7 -2.32180 33.88180 -0.00960 -0.00960 -0.00960 8 -7.14550 38.06720 -0.01252 -0.01252 -0.01252 9 -12.95030 43.23340 -0.01691 -0.01691 -0.01691 10 -16.81730 46.75710 -0.02047 -0.02047 -0.02047 11 -22.72300 52.18220 -0.02701 -0.02701 -0.02701 12 -27.92050 57.02640 -0.03407 -0.03407 -0.03407 13 -31.29330 60.16130 -0.03932 -0.03932 -0.03932 14 -37.79220 66.25930 -0.05113 -0.05113 -0.05113 15 -43.03410 71.19920 -0.06209 -0.06209 -0.06209 16 -46.26030 74.24110 -0.06940 -0.06940 -0.06940 main: the final sun glint reflectance matrix (1,1) term main: sunglint11giss is from gisscoxmunk_vfunction_plus.f90 main: sunglint11new is from coxmunk_new.f90 main: sunglint11breon is from breon.f90 (uses hafermanmatrix.f90) main: suunglint11guerin is from guerin.f90 main: nn,lats,zensun, singlint11new,sunglint11breon,sunglint11guerin 1 25.14480 17.10090 0.12066 0.12066 0.13212 0.12829 2 22.78770 17.63180 0.12107 0.12107 0.13236 0.12865 3 16.75890 19.97350 0.12346 0.12346 0.13429 0.13094 4 12.48190 22.48920 0.12644 0.12644 0.13665 0.13377 5 7.35300 26.08590 0.13172 0.13172 0.14097 0.13885 6 2.17110 30.13240 0.13976 0.13976 0.14790 0.14670 7 -2.32180 33.88180 0.14974 0.14974 0.15667 0.15649 8 -7.14550 38.06720 0.16508 0.16508 0.17042 0.17164 9 -12.95030 43.23340 0.19332 0.19332 0.19609 0.19964 10 -16.81730 46.75710 0.22161 0.22161 0.22192 0.22773 11 -22.72300 52.18220 0.28921 0.28921 0.28356 0.29476 12 -27.92050 57.02640 0.39088 0.39088 0.37532 0.39512 13 -31.29330 60.16130 0.49105 0.49105 0.46450 0.49341 14 -37.79220 66.25930 0.82646 0.82646 0.75682 0.81957 15 -43.03410 71.19920 1.34283 1.34283 1.19307 1.31565 16 -46.26030 74.24110 1.87224 1.87224 1.63570 1.82466 main: the final sun glint reflectance matrix (1,2) term main: sunglint12giss is from gisscoxmunk_vfunction_plus.f90 main: sunglint12new is from coxmunk_new.f90 main: sunglint12breon is from breon.f90 (uses hafermanmatrix.f90) main: suunglint12guerin is from guerin.f90 main: nn,lats,zensun, singlint12new,sunglint12breon,sunglint12guerin 1 25.14480 17.10090 -0.01335 -0.01335 -0.01462 -0.01419 2 22.78770 17.63180 -0.01424 -0.01424 -0.01557 -0.01513 3 16.75890 19.97350 -0.01883 -0.01883 -0.02048 -0.01997 4 12.48190 22.48920 -0.02465 -0.02465 -0.02664 -0.02608 5 7.35300 26.08590 -0.03489 -0.03489 -0.03734 -0.03678 6 2.17110 30.13240 -0.04991 -0.04991 -0.05281 -0.05238 7 -2.32180 33.88180 -0.06791 -0.06791 -0.07105 -0.07097 8 -7.14550 38.06720 -0.09420 -0.09420 -0.09725 -0.09794 9 -12.95030 43.23340 -0.13907 -0.13907 -0.14107 -0.14363 10 -16.81730 46.75710 -0.18056 -0.18056 -0.18082 -0.18555 11 -22.72300 52.18220 -0.26984 -0.26984 -0.26457 -0.27502 12 -27.92050 57.02640 -0.38747 -0.38747 -0.37204 -0.39167 13 -31.29330 60.16130 -0.49094 -0.49094 -0.46439 -0.49330 14 -37.79220 66.25930 -0.78700 -0.78700 -0.72068 -0.78043 15 -43.03410 71.19920 -1.16278 -1.16278 -1.03310 -1.13924 16 -46.26030 74.24110 -1.49617 -1.49617 -1.30714 -1.45815 main: nn,lats,zensun,shadowgiss,shadowbreon,shadowguerin 1 25.1448 17.1009 1.0000 1.0000 1.0000 2 22.7877 17.6318 1.0000 1.0000 1.0000 3 16.7589 19.9735 1.0000 1.0000 1.0000 4 12.4819 22.4892 1.0000 1.0000 1.0000 5 7.3530 26.0859 1.0000 1.0000 1.0000 6 2.1711 30.1324 1.0000 1.0000 1.0000 7 -2.3218 33.8818 1.0000 1.0000 1.0000 8 -7.1455 38.0672 1.0000 1.0000 1.0000 9 -12.9503 43.2334 1.0000 1.0000 1.0000 10 -16.8173 46.7571 1.0000 1.0000 1.0000 11 -22.7230 52.1822 1.0000 1.0000 1.0000 12 -27.9205 57.0264 1.0000 1.0000 1.0000 13 -31.2933 60.1613 1.0000 1.0000 1.0000 14 -37.7922 66.2593 0.9999 1.0000 1.0000 15 -43.0341 71.1992 0.9984 1.0000 1.0000 16 -46.2603 74.2411 0.9943 0.9999 0.9999 main: fresnelbreon is the fresnel value main: gcbreon is the breon Gram-Charlier series sum main: comparebreon and comparenew should be similar main: nn,lats,zensun,gcbreon,comparebreon,comparenew 1 25.1448 17.1009 0.0201 1.0804 6.0785 5.9859 2 22.7877 17.6318 0.0201 1.0788 6.0956 6.0031 3 16.7589 19.9735 0.0201 1.0736 6.1970 6.1044 4 12.4819 22.4892 0.0202 1.0671 6.3170 6.2243 5 7.3530 26.0859 0.0204 1.0572 6.5171 6.4247 6 2.1711 30.1324 0.0207 1.0459 6.7985 6.7063 7 -2.3218 33.8818 0.0211 1.0348 7.1125 7.0208 8 -7.1455 38.0672 0.0219 1.0219 7.5405 7.4498 9 -12.9503 43.2334 0.0235 1.0053 8.2081 8.1195 10 -16.8173 46.7571 0.0251 0.9934 8.7730 8.6868 11 -22.7230 52.1822 0.0289 0.9742 9.8768 9.7967 12 -27.9205 57.0264 0.0343 0.9557 11.1700 11.1001 13 -31.2933 60.1613 0.0393 0.9428 12.2054 12.1462 14 -37.7922 66.2593 0.0536 0.9154 14.8730 14.8511 15 -43.0341 71.1992 0.0716 0.8901 17.8922 17.9328 16 -46.2603 74.2411 0.0868 0.8740 20.4895 20.5931 main: exvalbreon and exvalnew are exp(-xx)/sigma2 values main: theta86new is the eqn 86 theta from getmtmatrix main: betadegbreon is the surface slope in degrees main: cs4breon is the cos(suface slope)**4 term main: term4new is RDZ4/(FACT1*FACT1) value main: term5new is dmodmt value main: nn,lats,zensun,exvalbreon,exvalnew,theta86new,betadegbreon,cs4breon main: term4new,term5new 1 25.1448 17.1009 45.1884 44.5004 15.3147 1.7862 0.9981 0.9981 0.0674 2 22.7877 17.6318 45.1081 44.4241 15.7776 1.8542 0.9979 0.9979 0.0750 3 16.7589 19.9735 44.8284 44.1580 17.8990 2.0745 0.9974 0.9974 0.1171 4 12.4819 22.4892 44.4545 43.8024 20.1508 2.3384 0.9967 0.9967 0.1750 5 7.3530 26.0859 43.8246 43.2029 23.3561 2.7298 0.9955 0.9955 0.2808 6 2.1711 30.1324 43.0337 42.4498 26.9732 3.1592 0.9939 0.9939 0.4272 7 -2.3218 33.8818 42.1850 41.6411 30.3120 3.5698 0.9923 0.9923 0.5766 8 -7.1455 38.0672 41.1283 40.6335 34.0351 4.0321 0.9901 0.9901 0.7405 9 -12.9503 43.2334 39.6597 39.2315 38.6197 4.6137 0.9871 0.9871 0.9048 10 -16.8173 46.7571 38.5451 38.1663 41.7350 5.0221 0.9847 0.9847 0.9743 11 -22.7230 52.1822 36.6522 36.3550 46.5115 5.6707 0.9806 0.9806 0.9944 12 -27.9205 57.0264 34.7224 34.5049 50.7348 6.2917 0.9761 0.9761 0.9225 13 -31.2933 60.1613 33.3149 33.1534 53.4343 6.7270 0.9727 0.9727 0.8387 14 -37.7922 66.2593 30.1877 30.1432 58.5977 7.6616 0.9648 0.9648 0.6259 15 -43.0341 71.1992 27.0767 27.1382 62.6287 8.5705 0.9561 0.9561 0.4446 16 -46.2603 74.2411 24.9761 25.1024 65.0568 9.1843 0.9497 0.9497 0.3421 main: compare refl(1,2) / refl(1,1) ratios main: ratiom12m11giss is from gisscoxmunk_vfuncrtion_plus.f90 main: ratio of R1(2)/R1(1) main: ratiom12m11breon is from breon.f90 ((which uses main: the hafermanmatrix.f90 routine) main: ratiom12m11mt is from Mishchenko and Travis MT paper main: nn,lats,zensun,ratiom12m11giss,ratiom12m11breon,ratiom12m11mt 1 25.1448 17.1009 -0.1106 -0.1106 -0.1106 2 22.7877 17.6318 -0.1176 -0.1176 -0.1176 3 16.7589 19.9735 -0.1525 -0.1525 -0.1525 4 12.4819 22.4892 -0.1949 -0.1949 -0.1949 5 7.3530 26.0859 -0.2649 -0.2649 -0.2649 6 2.1711 30.1324 -0.3571 -0.3571 -0.3571 7 -2.3218 33.8818 -0.4535 -0.4535 -0.4535 8 -7.1455 38.0672 -0.5706 -0.5706 -0.5706 9 -12.9503 43.2334 -0.7194 -0.7194 -0.7194 10 -16.8173 46.7571 -0.8148 -0.8148 -0.8148 11 -22.7230 52.1822 -0.9330 -0.9330 -0.9330 12 -27.9205 57.0264 -0.9913 -0.9913 -0.9913 13 -31.2933 60.1613 -0.9998 -0.9998 -0.9998 14 -37.7922 66.2593 -0.9523 -0.9523 -0.9523 15 -43.0341 71.1992 -0.8659 -0.8659 -0.8659 16 -46.2603 74.2411 -0.7991 -0.7991 -0.7991 main: factor62mult is MT equation 62 scalar part * dmodmt main: rspdivbreon is rsp / fresnel from breon.f90 main: ratio62=factor62multnew(nn)/rspdivbreon(nn) main: for perfect agreement ratio62 should be 1.00 main: nn,lats,zensun,factor62multnew,rspdivbreon 1 25.1448 17.1009 5.9976 6.5801 0.9115 2 22.7877 17.6318 6.0157 6.5898 0.9129 3 16.7589 19.9735 6.1204 6.6704 0.9175 4 12.4819 22.4892 6.2451 6.7633 0.9234 5 7.3530 26.0859 6.4539 6.9212 0.9325 6 2.1711 30.1324 6.7472 7.1541 0.9431 7 -2.3218 33.8818 7.0755 7.4174 0.9539 8 -7.1455 38.0672 7.5240 7.7826 0.9668 9 -12.9503 43.2334 8.2256 8.3594 0.9840 10 -16.8173 46.7571 8.8215 8.8500 0.9968 11 -22.7230 52.1822 9.9908 9.8125 1.0182 12 -27.9205 57.0264 11.3716 10.9362 1.0398 13 -31.2933 60.1613 12.4866 11.8291 1.0556 14 -37.7922 66.2593 15.3935 14.1125 1.0908 15 -43.0341 71.1992 18.7567 16.6566 1.1261 16 -46.2603 74.2411 21.6838 18.8549 1.1500 main: recalculate from scalar and matrix terms main: the fianl sun glint reflection matrix (1,1) terms main: from coxmunk_new, breon.f90, guerin.f90 main: windspeedbreon,windaz 8.0000 45.0000 main: a1=factor62new,b1=m11new,c1=a1*b1 main: a2=factor62multnew,b2=m11divnew,c2=a2*b2 main: a3=rspdivbreon,b3=m11breon,c3=a3*b3 main: c4=rspguerin main: nn,lats,zensun, c1,c2,c3,c4 1 25.1448 17.1009 0.1207 0.1207 0.1324 0.1283 2 22.7877 17.6318 0.1211 0.1211 0.1326 0.1287 3 16.7589 19.9735 0.1235 0.1235 0.1346 0.1309 4 12.4819 22.4892 0.1264 0.1264 0.1369 0.1338 5 7.3530 26.0859 0.1317 0.1317 0.1413 0.1388 6 2.1711 30.1324 0.1398 0.1398 0.1482 0.1467 7 -2.3218 33.8818 0.1497 0.1497 0.1570 0.1565 8 -7.1455 38.0672 0.1651 0.1651 0.1708 0.1716 9 -12.9503 43.2334 0.1933 0.1933 0.1965 0.1996 10 -16.8173 46.7571 0.2216 0.2216 0.2223 0.2277 11 -22.7230 52.1822 0.2892 0.2892 0.2841 0.2948 12 -27.9205 57.0264 0.3909 0.3909 0.3759 0.3951 13 -31.2933 60.1613 0.4910 0.4910 0.4652 0.4934 14 -37.7922 66.2593 0.8266 0.8266 0.7578 0.8196 15 -43.0341 71.1992 1.3450 1.3450 1.1944 1.3156 16 -46.2603 74.2411 1.8830 1.8830 1.6374 1.8247