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 0.0000 main: windazdiff (windaz+180-phisun) 0.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.13678 0.12803 2 22.78770 17.63180 0.12107 0.12107 0.13725 0.12840 3 16.75890 19.97350 0.12346 0.12346 0.13996 0.13069 4 12.48190 22.48920 0.12644 0.12644 0.14334 0.13355 5 7.35300 26.08590 0.13172 0.13172 0.14935 0.13868 6 2.17110 30.13240 0.13976 0.13976 0.15850 0.14664 7 -2.32180 33.88180 0.14974 0.14974 0.16984 0.15659 8 -7.14550 38.06720 0.16508 0.16508 0.18730 0.17202 9 -12.95030 43.23340 0.19332 0.19332 0.21943 0.20057 10 -16.81730 46.75710 0.22161 0.22161 0.25166 0.22927 11 -22.72300 52.18220 0.28921 0.28921 0.32869 0.29794 12 -27.92050 57.02640 0.39088 0.39088 0.44468 0.40124 13 -31.29330 60.16130 0.49105 0.49105 0.55910 0.50297 14 -37.79220 66.25930 0.82646 0.82646 0.94327 0.84372 15 -43.03410 71.19920 1.34283 1.34283 1.53940 1.37082 16 -46.26030 74.24110 1.87224 1.87224 2.16062 1.91970 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.13678 0.12803 2 22.78770 17.63180 0.12107 0.12107 0.13725 0.12840 3 16.75890 19.97350 0.12346 0.12346 0.13996 0.13069 4 12.48190 22.48920 0.12644 0.12644 0.14334 0.13355 5 7.35300 26.08590 0.13172 0.13172 0.14935 0.13868 6 2.17110 30.13240 0.13976 0.13976 0.15850 0.14664 7 -2.32180 33.88180 0.14974 0.14974 0.16984 0.15659 8 -7.14550 38.06720 0.16508 0.16508 0.18730 0.17202 9 -12.95030 43.23340 0.19332 0.19332 0.21943 0.20057 10 -16.81730 46.75710 0.22161 0.22161 0.25166 0.22927 11 -22.72300 52.18220 0.28921 0.28921 0.32869 0.29794 12 -27.92050 57.02640 0.39088 0.39088 0.44468 0.40124 13 -31.29330 60.16130 0.49105 0.49105 0.55910 0.50297 14 -37.79220 66.25930 0.82646 0.82646 0.94327 0.84372 15 -43.03410 71.19920 1.34283 1.34283 1.53937 1.37082 16 -46.26030 74.24110 1.87224 1.87224 2.16047 1.91970 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.01513 -0.01416 2 22.78770 17.63180 -0.01424 -0.01424 -0.01614 -0.01510 3 16.75890 19.97350 -0.01883 -0.01883 -0.02135 -0.01993 4 12.48190 22.48920 -0.02465 -0.02465 -0.02794 -0.02603 5 7.35300 26.08590 -0.03489 -0.03489 -0.03956 -0.03673 6 2.17110 30.13240 -0.04991 -0.04991 -0.05660 -0.05236 7 -2.32180 33.88180 -0.06791 -0.06791 -0.07703 -0.07102 8 -7.14550 38.06720 -0.09420 -0.09420 -0.10688 -0.09816 9 -12.95030 43.23340 -0.13907 -0.13907 -0.15786 -0.14429 10 -16.81730 46.75710 -0.18056 -0.18056 -0.20504 -0.18680 11 -22.72300 52.18220 -0.26984 -0.26984 -0.30668 -0.27798 12 -27.92050 57.02640 -0.38747 -0.38747 -0.44080 -0.39774 13 -31.29330 60.16130 -0.49094 -0.49094 -0.55897 -0.50286 14 -37.79220 66.25930 -0.78700 -0.78700 -0.89823 -0.80344 15 -43.03410 71.19920 -1.16278 -1.16278 -1.33297 -1.18702 16 -46.26030 74.24110 -1.49617 -1.49617 -1.72651 -1.53410 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.1136 6.1053 5.9859 2 22.7877 17.6318 0.0201 1.1133 6.1246 6.0031 3 16.7589 19.9735 0.0201 1.1123 6.2340 6.1044 4 12.4819 22.4892 0.0202 1.1109 6.3649 6.2243 5 7.3530 26.0859 0.0204 1.1085 6.5846 6.4247 6 2.1711 30.1324 0.0207 1.1055 6.8930 6.7063 7 -2.3218 33.8818 0.0211 1.1022 7.2390 7.0208 8 -7.1455 38.0672 0.0219 1.0981 7.7122 7.4498 9 -12.9503 43.2334 0.0235 1.0922 8.4540 8.1195 10 -16.8173 46.7571 0.0251 1.0877 9.0855 8.6868 11 -22.7230 52.1822 0.0289 1.0799 10.3282 9.7967 12 -27.9205 57.0264 0.0343 1.0716 11.8027 11.1001 13 -31.2933 60.1613 0.0393 1.0654 12.9998 12.1462 14 -37.7922 66.2593 0.0536 1.0511 16.1442 14.8511 15 -43.0341 71.1992 0.0716 1.0361 19.8321 17.9328 16 -46.2603 74.2411 0.0868 1.0254 23.0666 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.3880 44.5004 15.3147 1.7862 0.9981 0.9981 0.0674 2 22.7877 17.6318 45.3229 44.4241 15.7776 1.8542 0.9979 0.9979 0.0750 3 16.7589 19.9735 45.0958 44.1580 17.8990 2.0745 0.9974 0.9974 0.1171 4 12.4819 22.4892 44.7918 43.8024 20.1508 2.3384 0.9967 0.9967 0.1750 5 7.3530 26.0859 44.2785 43.2029 23.3561 2.7298 0.9955 0.9955 0.2808 6 2.1711 30.1324 43.6320 42.4498 26.9732 3.1592 0.9939 0.9939 0.4272 7 -2.3218 33.8818 42.9357 41.6411 30.3120 3.5698 0.9923 0.9923 0.5766 8 -7.1455 38.0672 42.0650 40.6335 34.0351 4.0321 0.9901 0.9901 0.7405 9 -12.9503 43.2334 40.8477 39.2315 38.6197 4.6137 0.9871 0.9871 0.9048 10 -16.8173 46.7571 39.9179 38.1663 41.7350 5.0221 0.9847 0.9847 0.9743 11 -22.7230 52.1822 38.3272 36.3550 46.5115 5.6707 0.9806 0.9806 0.9944 12 -27.9205 57.0264 36.6890 34.5049 50.7348 6.2917 0.9761 0.9761 0.9225 13 -31.2933 60.1613 35.4831 33.1534 53.4343 6.7270 0.9727 0.9727 0.8387 14 -37.7922 66.2593 32.7678 30.1432 58.5977 7.6616 0.9648 0.9648 0.6259 15 -43.0341 71.1992 30.0124 27.1382 62.6287 8.5705 0.9561 0.9561 0.4446 16 -46.2603 74.2411 28.1175 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.8123 0.8804 2 22.7877 17.6318 6.0157 6.8330 0.8804 3 16.7589 19.9735 6.1204 6.9522 0.8804 4 12.4819 22.4892 6.2451 7.0943 0.8803 5 7.3530 26.0859 6.4539 7.3324 0.8802 6 2.1711 30.1324 6.7472 7.6667 0.8801 7 -2.3218 33.8818 7.0755 8.0413 0.8799 8 -7.1455 38.0672 7.5240 8.5531 0.8797 9 -12.9503 43.2334 8.2256 9.3545 0.8793 10 -16.8173 46.7571 8.8215 10.0357 0.8790 11 -22.7230 52.1822 9.9908 11.3740 0.8784 12 -27.9205 57.0264 11.3716 12.9572 0.8776 13 -31.2933 60.1613 12.4866 14.2382 0.8770 14 -37.7922 66.2593 15.3935 17.5892 0.8752 15 -43.0341 71.1992 18.7567 21.4915 0.8727 16 -46.2603 74.2411 21.6838 24.9041 0.8707 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 0.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.1371 0.1280 2 22.7877 17.6318 0.1211 0.1211 0.1375 0.1284 3 16.7589 19.9735 0.1235 0.1235 0.1402 0.1307 4 12.4819 22.4892 0.1264 0.1264 0.1436 0.1335 5 7.3530 26.0859 0.1317 0.1317 0.1496 0.1387 6 2.1711 30.1324 0.1398 0.1398 0.1588 0.1466 7 -2.3218 33.8818 0.1497 0.1497 0.1702 0.1566 8 -7.1455 38.0672 0.1651 0.1651 0.1877 0.1720 9 -12.9503 43.2334 0.1933 0.1933 0.2198 0.2006 10 -16.8173 46.7571 0.2216 0.2216 0.2521 0.2293 11 -22.7230 52.1822 0.2892 0.2892 0.3293 0.2979 12 -27.9205 57.0264 0.3909 0.3909 0.4454 0.4012 13 -31.2933 60.1613 0.4910 0.4910 0.5599 0.5030 14 -37.7922 66.2593 0.8266 0.8266 0.9445 0.8437 15 -43.0341 71.1992 1.3450 1.3450 1.5411 1.3708 16 -46.2603 74.2411 1.8830 1.8830 2.1627 1.9197