SHDOM single-scattering adjoint

This folder contains SHDOM files related to a high resolution simulation that was performed by Frank Evans in 2013. The cloud field was generated with a large-eddy simulation. SHDOM ran at 320x320xNZxNSH resolution. The NADIR radiance of the cloud field:

The same cloud field is used for adjoint testing in 2D. Presently, the adjoint is working for single-scattering contributions to the misfit function derivative. This component was tested using a 2D slice ($y=5.5\mathrm{km}$).


In [1]:
%pylab inline

# Imports
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm

import pandas as pd


Populating the interactive namespace from numpy and matplotlib

In [2]:
ls


basic_imports.py
exploring_adjoint.ipynb
les0822nh15t13_aer008_w0.466rad.modis
les0822nh15t13_aer008_w0.466.ref
les0822nh15t13_aer008_w0.646a_idir01.jpg
les0822nh15t13_aer008_w0.646rad.modis
les0822nh15t13_aer008_w0.646.ref
les0822nh15t13_aer008_w1.632rad.modis
les0822nh15t13_aer008_w1.632.ref
notebook-check-data.ipynb
radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.adj
radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.arad
radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.log
radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss2.arad
read_modis_sim.pro

In [ ]:


In [3]:
# Utility functions

def print_head(fname, Nlines=10, indent="\t"):
    "Print the head of the file."
    
    # Print a message and then the first N lines
    print("Showing head: {}".format(fname))
    for i, line in zip(range(Nlines), open(fname, 'r')):
        print(indent + line.strip())
        
    print('\n')
    
def read_log(flog, start_line="X      Z   Extinct", max_lines=10**8):
    """
    Read in a log file but save lines only after reaching a line 
    matching start_line.
    """
    
    # Loop over lines record after header and before footer
    ncols = None 
    lines = []
    for i, line in zip(range(max_lines), open(flog, 'r')):
        new = line.strip().split()
        
        # Check if the line is the header
        if ncols is None:
            if line.strip()[:len(start_line)]==start_line:
                ncols = len(new[:9])
                cols = new[:9]
        else:        
            if ncols == len(new[:9]):
                lines.append(new[:9])
            else:
                print(new)
                break
    # Raise an error if the file had no line with the specified columns    
    if ncols is None: 
        raise ValueError("No column line found containing {}".format(start_line))
        
    # Make dictionary of arrays from keys array with given lines
    N = len(lines) - 1
    value_array = np.array(lines[:], dtype='f8')
    
    log_dict = dict(zip(cols, list(value_array.T)))
    return log_dict

Looking at the data files


In [17]:
# The log file
flog = "radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.log"
print_head(flog, Nlines=500)


Showing head: radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.log
	
	Name of this SHDOMADJ run:
	radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1
	Input particle properties filename (or NONE for Rayleigh only)
	les0822nh15t13y89_ocaer1.part
	Number of input scattering tables
	2
	Name of each scattering table
	water_w0.646.scat
	aero_oc_w0.646.scat
	Index or type of each scattering table
	1  2
	Polarization flag (T or F) for each scattering table
	T T
	Wavelength for Rayleigh molecular scattering (micron)
	0.6460000000000000
	Surface pressure for Rayleigh molecular scattering (mb or hPa)
	1010.000000000000
	Number of extra Z levels
	0
	Surface properties file name (or NONE)
	NONE
	Correlated k-distribution file name (or NONE)
	NONE
	Number of Stokes parameters (1, 3, or 4; 1=unpolarized)
	1
	Number of discrete ordinates in mu and phi
	16           32
	Independent pixel approximation (0,1,2,3)
	0
	Delta-M method (T,F)
	T
	Thermal, solar, or both source (T, S, B)
	S
	Solar flux and direction (F, mu, phi[deg])
	1.000000000000000       -0.9221000000000000         339.4000000000000
	Isotropic sky radiance
	0.000000000000000
	Ground albedo
	5.0000000000000003E-002
	Wavelength (micron)
	0.6460000000000000
	Cell splitting and spherical harmonic series accuracy
	2.9999999999999999E-002   3.0000000000000001E-003
	Acceleration flag (T or F), solution accuracy, and max iterations
	T   3.0000000000000001E-005            0
	Radiance grid output Z level, number in X and Y
	15.00000000000000               320            1
	Number of radiance directions output
	3
	Radiance mu, phi[deg] pairs
	0.8660300000000000         180.0000000000000
	1.000000000000000         0.000000000000000
	0.8660300000000000         0.000000000000000
	Radiance output file name
	radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.arad
	Observed radiance input file name (or NONE)
	radobs_les0822nh15t13y89_ocaer1_w0.646_ns1z3.arad
	Observed radiance uncertainty (sigma)
	2.9999999999999999E-002
	Adjoint output file name
	radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.adj
	Approximate maximum memory to use (MB for 4 byte reals)
	2000.000000000000
	Ratio of total grid points to base grid points (>=1.0)
	5.000000000000000
	Ratio of average number of spherical harmonic terms to total possible (<=1.0)
	0.9500000000000000
	Ratio of number of grid cells to grid points (1.0 to 2.0)
	1.500000000000000
	
	Rayleigh total atmosphere optical depth:    5.0269465934799501E-002
	Starting solution iterations for SHDOM run:
	radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1
	rms of I radiance residuals:   0.43786
	Finite differencing IX=            1
	Finite differencing IX=           41
	Finite differencing IX=           81
	Finite differencing IX=          121
	Finite differencing IX=          161
	Finite differencing IX=          201
	Finite differencing IX=          241
	Finite differencing IX=          281
	X      Z   Extinct FDextinct ADextinct  FDssalb   ADssab    FDsscatI  ADsscatI  FDsscatQ  ADsscatQ  FDsscatU  ADsscatU
	0.000  0.000   0.070  -0.00001  -0.00001  0.000001  0.000001  0.000025  0.000025
	0.000  0.010   0.070  -0.00003  -0.00003  0.000005  0.000005  0.000081  0.000081
	0.000  0.032   0.070  -0.00004  -0.00004  0.000007  0.000007  0.000111  0.000111
	0.000  0.057   0.069  -0.00005  -0.00005  0.000008  0.000008  0.000113  0.000113
	0.000  0.088   0.069  -0.00007  -0.00007  0.000011  0.000011  0.000180  0.000180
	0.000  0.125   0.068  -0.00010  -0.00010  0.000015  0.000015  0.000301  0.000301
	0.000  0.165   0.067  -0.00011  -0.00011  0.000021  0.000021  0.000466  0.000466
	0.000  0.205   0.067  -0.00018  -0.00018  0.000030  0.000030  0.000762  0.000762
	0.000  0.245   0.066  -0.00026  -0.00026  0.000043  0.000043  0.001200  0.001200
	0.000  0.285   0.065  -0.00022  -0.00022  0.000059  0.000059  0.001716  0.001716
	0.000  0.325   0.064  -0.00014  -0.00014  0.000074  0.000074  0.002134  0.002134
	0.000  0.365   0.064  -0.00014  -0.00014  0.000086  0.000086  0.002457  0.002457
	0.000  0.405   0.063  -0.00009  -0.00009  0.000099  0.000099  0.002855  0.002855
	0.000  0.445   0.062  -0.00002  -0.00002  0.000115  0.000115  0.003374  0.003374
	0.000  0.485   0.062  -0.00005  -0.00005  0.000125  0.000125  0.003715  0.003715
	0.000  0.525   0.061  -0.00005  -0.00005  0.000139  0.000139  0.004198  0.004198
	0.000  0.565   0.060   0.00001   0.00001  0.000158  0.000158  0.004894  0.004894
	0.000  0.605   0.059   0.00003   0.00003  0.000174  0.000174  0.005481  0.005481
	0.000  0.645   0.059   0.00064   0.00064  0.000233  0.000233  0.007711  0.007711
	0.000  0.685   0.058   0.00115   0.00115  0.000297  0.000297  0.010090  0.010090
	0.000  0.725   0.057   0.00160   0.00160  0.000348  0.000348  0.011998  0.011998
	0.000  0.765   0.057   0.00313   0.00313  0.000449  0.000449  0.015763  0.015763
	0.000  0.805   0.056   0.00472   0.00472  0.000549  0.000549  0.019438  0.019438
	0.000  0.845   0.055   0.00475   0.00475  0.000551  0.000551  0.019450  0.019450
	0.000  0.885   0.055   0.00427   0.00427  0.000513  0.000513  0.017979  0.017979
	0.000  0.925   0.054   0.00357   0.00357  0.000465  0.000465  0.016162  0.016162
	0.000  0.965   0.053   0.00267   0.00267  0.000404  0.000404  0.013996  0.013996
	0.000  1.005   0.052   0.00179   0.00179  0.000320  0.000320  0.010992  0.010992
	0.000  1.045   0.052   0.00091   0.00091  0.000236  0.000236  0.008067  0.008067
	0.000  1.085   0.051   0.00037   0.00037  0.000173  0.000173  0.005907  0.005907
	0.000  1.125   0.050   0.00011   0.00011  0.000132  0.000132  0.004476  0.004476
	0.000  1.165   0.050  -0.00014  -0.00014  0.000099  0.000099  0.003337  0.003337
	0.000  1.205   0.049  -0.00042  -0.00042  0.000067  0.000067  0.002221  0.002221
	0.000  1.245   0.048  -0.00058  -0.00058  0.000042  0.000042  0.001355  0.001355
	0.000  1.285   0.048  -0.00061  -0.00061  0.000026  0.000026  0.000821  0.000821
	0.000  1.325   0.047  -0.00065  -0.00065  0.000015  0.000015  0.000444  0.000444
	0.000  1.365   0.046  -0.00061  -0.00061  0.000009  0.000009  0.000248  0.000248
	0.000  1.405   0.046  -0.00051  -0.00051  0.000007  0.000007  0.000172  0.000172
	0.000  1.445   0.045  -0.00044  -0.00044  0.000005  0.000005  0.000107  0.000107
	0.000  1.485   0.044  -0.00042  -0.00042  0.000004  0.000004  0.000074  0.000074
	0.000  1.525   6.331  -0.00026  -0.00026  0.001453  0.001453  0.011993  0.011993
	0.000  1.565  20.798  -0.00023  -0.00023  0.006042  0.006042  0.064990  0.064990
	0.000  1.605  30.529  -0.00026  -0.00026  0.017429  0.017429  0.325657  0.325657
	0.000  1.645  32.122  -0.00043  -0.00043  0.062453  0.062453  1.594797  1.594797
	0.000  1.686  27.786  -0.00053  -0.00053  0.184585  0.184585  5.457382  5.457382
	0.000  1.722   7.103   0.00142   0.00142  0.093610  0.093610  2.662625  2.662625
	0.000  1.753   0.040   0.00501   0.00501  0.000664  0.000664  0.013813  0.013813
	0.000  1.778   0.039   0.00449   0.00449  0.000547  0.000547  0.011317  0.011317
	0.000  1.799   0.039   0.00415   0.00415  0.000474  0.000474  0.009757  0.009757
	0.000  1.819   0.039   0.00418   0.00418  0.000451  0.000451  0.009196  0.009196
	0.000  1.839   0.038   0.00436   0.00436  0.000446  0.000446  0.008988  0.008988
	0.000  1.859   0.038   0.00452   0.00452  0.000437  0.000437  0.008731  0.008731
	0.000  1.879   0.038   0.00464   0.00464  0.000429  0.000429  0.008491  0.008491
	0.000  1.899   0.037   0.00466   0.00466  0.000421  0.000421  0.008242  0.008242
	0.000  1.919   0.037   0.00470   0.00470  0.000414  0.000414  0.007998  0.007998
	0.000  1.939   0.037   0.00475   0.00475  0.000407  0.000407  0.007748  0.007748
	0.000  1.959   0.036   0.00481   0.00481  0.000401  0.000401  0.007507  0.007507
	0.000  1.979   0.036   0.00490   0.00490  0.000394  0.000394  0.007272  0.007272
	0.000  1.999   0.035   0.00498   0.00498  0.000386  0.000386  0.007001  0.007001
	0.000  2.019   0.035   0.00502   0.00502  0.000379  0.000379  0.006763  0.006763
	0.000  2.039   0.035   0.00533   0.00533  0.000386  0.000386  0.006680  0.006680
	0.000  2.060   0.035   0.00619   0.00619  0.000433  0.000433  0.007268  0.007268
	0.000  2.085   0.035   0.00753   0.00753  0.000511  0.000511  0.008267  0.008267
	0.000  2.116   7.548   0.00269   0.00269  0.096058  0.096058  2.421359  2.421359
	0.000  2.153   2.166   0.00807   0.00807  0.049806  0.049806  0.959019  0.959019
	0.000  2.193   0.034   0.01468   0.01468  0.001077  0.001077  0.014527  0.014527
	0.000  2.233   2.880   0.00694   0.00694  0.075139  0.075139  1.556807  1.556807
	0.000  2.273   0.033   0.01516   0.01516  0.001266  0.001266  0.017414  0.017414
	0.000  2.313   0.033   0.01479   0.01479  0.001310  0.001310  0.019117  0.019117
	0.000  2.353   0.032   0.01384   0.01384  0.001302  0.001302  0.019089  0.019089
	0.000  2.393   0.032   0.01369   0.01369  0.001277  0.001277  0.018550  0.018550
	0.000  2.433   0.032   0.01435   0.01435  0.001254  0.001254  0.018026  0.018026
	0.000  2.473   0.031   0.01513   0.01513  0.001234  0.001234  0.017461  0.017461
	0.000  2.513   0.031   0.01568   0.01568  0.001210  0.001210  0.016718  0.016718
	0.000  2.553   0.031   0.01625   0.01625  0.001189  0.001189  0.016090  0.016090
	0.000  2.593   0.030   0.01738   0.01738  0.001190  0.001190  0.015953  0.015953
	0.000  2.633   0.030   0.01854   0.01854  0.001192  0.001192  0.015816  0.015816
	0.000  2.673   0.029   0.01931   0.01931  0.001189  0.001189  0.015641  0.015641
	0.000  2.713   0.029   0.01914   0.01914  0.001180  0.001180  0.015349  0.015349
	0.000  2.753   0.029   0.01916   0.01916  0.001172  0.001172  0.015056  0.015056
	0.000  2.793   0.028   0.29280   0.29280  0.018793  0.018793  0.245885  0.245885
	0.000  4.000   0.018   1.11947   1.11947  0.042883  0.042883  0.554364  0.554364
	0.000  6.000   0.011   2.04891   2.04891  0.041939  0.041939  0.443700  0.443700
	0.000  8.000   0.005   9.28452   9.28452  0.072127  0.072127  0.479430  0.479430
	0.000 15.000   0.001  16.02494  16.02494  0.021535  0.021535  0.070929  0.070929
	2.500  0.000   0.070  -0.00145  -0.00145  0.000268  0.000268  0.003729  0.003729
	2.500  0.010   0.070  -0.00480  -0.00480  0.000895  0.000895  0.012602  0.012602
	2.500  0.032   0.070  -0.00724  -0.00724  0.001382  0.001382  0.019992  0.019992
	2.500  0.057   0.069  -0.00845  -0.00845  0.001655  0.001655  0.024715  0.024715
	2.500  0.088   0.069  -0.01010  -0.01010  0.001992  0.001992  0.030523  0.030523
	2.500  0.125   0.068  -0.01161  -0.01161  0.002322  0.002322  0.036787  0.036787
	2.500  0.165   0.067  -0.01199  -0.01199  0.002454  0.002454  0.040255  0.040255
	2.500  0.205   0.067  -0.01188  -0.01188  0.002414  0.002414  0.039955  0.039955
	2.500  0.245   0.066  -0.01183  -0.01183  0.002381  0.002381  0.039781  0.039781
	2.500  0.285   0.065  -0.01174  -0.01174  0.002357  0.002357  0.039567  0.039567
	2.500  0.325   0.064  -0.01133  -0.01133  0.002344  0.002344  0.038704  0.038704
	2.500  0.365   0.064  -0.01090  -0.01090  0.002332  0.002332  0.037861  0.037861
	2.500  0.405   0.063  -0.01045  -0.01045  0.002317  0.002317  0.037077  0.037077
	2.500  0.445   0.062  -0.00967  -0.00967  0.002292  0.002292  0.035787  0.035787
	2.500  0.485   0.062  -0.00868  -0.00868  0.002267  0.002267  0.034549  0.034549
	2.500  0.525   0.061  -0.00776  -0.00776  0.002213  0.002213  0.033607  0.033607
	2.500  0.565   0.060  -0.00782  -0.00782  0.002077  0.002077  0.031572  0.031572
	2.500  0.605   0.059  -0.00809  -0.00809  0.001942  0.001942  0.029537  0.029537
	2.500  0.645   0.059  -0.00857  -0.00857  0.001836  0.001836  0.028561  0.028561
	2.500  0.685   0.058  -0.00880  -0.00880  0.001749  0.001749  0.027877  0.027877
	2.500  0.725   0.057  -0.00839  -0.00839  0.001715  0.001715  0.027803  0.027803
	2.500  0.765   0.057  -0.00852  -0.00852  0.001695  0.001695  0.028407  0.028407
	2.500  0.805   0.056  -0.00901  -0.00901  0.001675  0.001675  0.028945  0.028945
	2.500  0.845   0.055  -0.00969  -0.00969  0.001661  0.001661  0.029028  0.029028
	2.500  0.885   0.055  -0.00959  -0.00959  0.001682  0.001682  0.029704  0.029704
	2.500  0.925   0.054  -0.00895  -0.00895  0.001712  0.001712  0.030536  0.030536
	2.500  0.965   0.053  -0.00823  -0.00823  0.001741  0.001741  0.030912  0.030912
	2.500  1.005   0.052  -0.00785  -0.00785  0.001767  0.001767  0.031214  0.031214
	2.500  1.045   0.052  -0.00735  -0.00735  0.001798  0.001798  0.031646  0.031646
	2.500  1.085   0.051  -0.00632  -0.00632  0.001830  0.001830  0.032198  0.032198
	2.500  1.125   0.050  -0.00486  -0.00486  0.001861  0.001861  0.032757  0.032757
	2.500  1.165   0.050  -0.00326  -0.00326  0.001867  0.001867  0.032481  0.032481
	2.500  1.205   0.049  -0.00274  -0.00274  0.001801  0.001801  0.030449  0.030449
	2.500  1.245   0.048  -0.00285  -0.00285  0.001733  0.001733  0.028443  0.028443
	2.500  1.285   0.048  -0.00260  -0.00260  0.001708  0.001708  0.027840  0.027840
	2.500  1.325   0.047  -0.00223  -0.00223  0.001684  0.001684  0.027312  0.027312
	2.500  1.365   0.046  -0.00179  -0.00179  0.001666  0.001666  0.026969  0.026969
	2.500  1.405   0.046  -0.00094  -0.00094  0.001689  0.001689  0.027989  0.027989
	2.500  1.445   0.045  -0.00019  -0.00019  0.001714  0.001714  0.029058  0.029058
	2.500  1.485   0.044   0.00018   0.00018  0.001723  0.001723  0.029486  0.029486
	2.500  1.525   0.044   0.00056   0.00056  0.001732  0.001732  0.029901  0.029901
	2.500  1.565   0.043   0.00095   0.00095  0.001740  0.001740  0.030272  0.030272
	2.500  1.605   0.042   0.00123   0.00123  0.001730  0.001730  0.030011  0.030011
	2.500  1.645   0.041   0.00157   0.00157  0.001728  0.001728  0.029825  0.029825
	2.500  1.686   0.041   0.00185   0.00185  0.001644  0.001644  0.028232  0.028232
	2.500  1.722   0.040   0.00187   0.00187  0.001419  0.001419  0.024256  0.024256
	2.500  1.753   0.040   0.00175   0.00175  0.001173  0.001173  0.019972  0.019972
	2.500  1.778   0.039   0.00161   0.00161  0.000974  0.000974  0.016531  0.016531
	2.500  1.799   0.039   0.00153   0.00153  0.000853  0.000853  0.014433  0.014433
	2.500  1.819   0.039   0.00158   0.00158  0.000825  0.000825  0.013913  0.013913
	2.500  1.839   0.038   0.00169   0.00169  0.000825  0.000825  0.013874  0.013874
	2.500  1.859   0.038   0.00180   0.00180  0.000821  0.000821  0.013754  0.013754
	2.500  1.879   0.038   0.00190   0.00190  0.000817  0.000817  0.013656  0.013656
	2.500  1.899   0.037   0.00201   0.00201  0.000814  0.000814  0.013546  0.013546
	2.500  1.919   0.037   0.00212   0.00212  0.000793  0.000793  0.012881  0.012881
	2.500  1.939   0.037   0.00218   0.00218  0.000769  0.000769  0.012117  0.012117
	2.500  1.959   0.036   0.00225   0.00225  0.000746  0.000746  0.011382  0.011382
	2.500  1.979   0.036   0.00231   0.00231  0.000722  0.000722  0.010657  0.010657
	2.500  1.999   0.035   0.00238   0.00238  0.000697  0.000697  0.009852  0.009852
	2.500  2.019   0.035   0.00232   0.00232  0.000677  0.000677  0.009303  0.009303
	2.500  2.039   0.035   0.00218   0.00218  0.000674  0.000674  0.008840  0.008840
	2.500  2.060   0.035   0.00225   0.00225  0.000741  0.000741  0.009247  0.009247
	2.500  2.085   0.035   0.00233   0.00233  0.000850  0.000850  0.009931  0.009931
	2.500  2.116   0.034   0.00224   0.00224  0.000969  0.000969  0.010268  0.010268
	2.500  2.153   0.034   0.00258   0.00258  0.001113  0.001113  0.011859  0.011859
	2.500  2.193   0.034   0.00273   0.00273  0.001162  0.001162  0.012505  0.012505
	2.500  2.233   0.033   0.00325   0.00325  0.001161  0.001161  0.012710  0.012710
	2.500  2.273   0.033   0.00563   0.00563  0.001226  0.001226  0.014852  0.014852
	2.500  2.313   0.033   0.00827   0.00827  0.001293  0.001293  0.017061  0.017061
	2.500  2.353   0.032   0.01042   0.01042  0.001313  0.001313  0.017736  0.017736
	2.500  2.393   0.032   0.01180   0.01180  0.001304  0.001304  0.017533  0.017533
	2.500  2.433   0.032   0.01198   0.01198  0.001294  0.001294  0.017332  0.017332
	2.500  2.473   0.031   0.01193   0.01193  0.001279  0.001279  0.016940  0.016940
	2.500  2.513   0.031   0.01184   0.01184  0.001259  0.001259  0.016410  0.016410
	2.500  2.553   0.031   0.01167   0.01167  0.001241  0.001241  0.015976  0.015976
	2.500  2.593   0.030   0.01154   0.01154  0.001258  0.001258  0.016588  0.016588
	2.500  2.633   0.030   0.01140   0.01140  0.001277  0.001277  0.017248  0.017248
	2.500  2.673   0.029   0.01114   0.01114  0.001277  0.001277  0.017302  0.017302
	2.500  2.713   0.029   0.01084   0.01084  0.001248  0.001248  0.016365  0.016365
	2.500  2.753   0.029   0.01083   0.01083  0.001218  0.001218  0.015437  0.015437
	2.500  2.793   0.028   0.17463   0.17463  0.022760  0.022760  0.337792  0.337792
	2.500  4.000   0.018   0.87586   0.87586  0.048460  0.048460  0.612745  0.612745
	2.500  6.000   0.011   1.74001   1.74001  0.045131  0.045131  0.466419  0.466419
	2.500  8.000   0.005   7.38290   7.38290  0.063415  0.063415  0.381025  0.381025
	2.500 15.000   0.001  14.79615  14.79615  0.020729  0.020729  0.076545  0.076545
	5.000  0.000   0.070  -0.00091  -0.00091  0.000214  0.000214  0.004860  0.004860
	5.000  0.010   0.070  -0.00300  -0.00300  0.000708  0.000708  0.016098  0.016098
	5.000  0.032   0.070  -0.00454  -0.00454  0.001076  0.001076  0.024421  0.024421
	5.000  0.057   0.069  -0.00530  -0.00530  0.001264  0.001264  0.028624  0.028624
	5.000  0.088   0.069  -0.00627  -0.00627  0.001507  0.001507  0.034011  0.034011
	5.000  0.125   0.068  -0.00714  -0.00714  0.001736  0.001736  0.038983  0.038983
	5.000  0.165   0.067  -0.00735  -0.00735  0.001811  0.001811  0.040456  0.040456
	5.000  0.205   0.067  -0.00721  -0.00721  0.001804  0.001804  0.040018  0.040018
	5.000  0.245   0.066  -0.00711  -0.00711  0.001810  0.001810  0.039806  0.039806
	5.000  0.285   0.065  -0.00702  -0.00702  0.001813  0.001813  0.039634  0.039634
	5.000  0.325   0.064  -0.00693  -0.00693  0.001811  0.001811  0.039401  0.039401
	5.000  0.365   0.064  -0.00686  -0.00686  0.001810  0.001810  0.039208  0.039208
	5.000  0.405   0.063  -0.00678  -0.00678  0.001805  0.001805  0.038990  0.038990
	5.000  0.445   0.062  -0.00672  -0.00672  0.001798  0.001798  0.038767  0.038767
	5.000  0.485   0.062  -0.00666  -0.00666  0.001790  0.001790  0.038546  0.038546
	5.000  0.525   0.061  -0.00656  -0.00656  0.001785  0.001785  0.038302  0.038302
	5.000  0.565   0.060  -0.00647  -0.00647  0.001780  0.001780  0.038062  0.038062
	5.000  0.605   0.059  -0.00637  -0.00637  0.001775  0.001775  0.037817  0.037817
	5.000  0.645   0.059  -0.00626  -0.00626  0.001771  0.001771  0.037566  0.037566
	5.000  0.685   0.058  -0.00614  -0.00614  0.001767  0.001767  0.037318  0.037318
	5.000  0.725   0.057  -0.00601  -0.00601  0.001761  0.001761  0.037027  0.037027
	5.000  0.765   0.057  -0.00588  -0.00588  0.001758  0.001758  0.036781  0.036781
	5.000  0.805   0.056  -0.00573  -0.00573  0.001753  0.001753  0.036489  0.036489
	5.000  0.845   0.055  -0.00560  -0.00560  0.001750  0.001750  0.036240  0.036240
	5.000  0.885   0.055  -0.00543  -0.00543  0.001750  0.001750  0.035950  0.035950
	5.000  0.925   0.054  -0.00526  -0.00526  0.001749  0.001749  0.035671  0.035671
	5.000  0.965   0.053  -0.00503  -0.00503  0.001759  0.001759  0.035379  0.035379
	5.000  1.005   0.052  -0.00478  -0.00478  0.001772  0.001772  0.035075  0.035075
	5.000  1.045   0.052  -0.00453  -0.00453  0.001784  0.001784  0.034780  0.034780
	5.000  1.085   0.051  -0.00408  -0.00408  0.001828  0.001828  0.034466  0.034466
	5.000  1.125   0.050  -0.00362  -0.00362  0.001873  0.001873  0.034154  0.034154
	5.000  1.165   0.050  -0.00314  -0.00314  0.001917  0.001917  0.033845  0.033845
	5.000  1.205   0.049  -0.00234  -0.00234  0.002009  0.002009  0.033515  0.033515
	5.000  1.245   0.048  -0.00156  -0.00156  0.002094  0.002094  0.033194  0.033194
	5.000  1.285   0.048  -0.00105  -0.00105  0.002129  0.002129  0.032852  0.032852
	5.000  1.325   0.047  -0.00046  -0.00046  0.002173  0.002173  0.032500  0.032500
	5.000  1.365   0.046   0.00002   0.00002  0.002201  0.002201  0.032180  0.032180
	5.000  1.405   0.046   0.00027   0.00027  0.002186  0.002186  0.031833  0.031833
	5.000  1.445   0.045   0.00061   0.00061  0.002171  0.002171  0.031481  0.031481
	5.000  1.485   0.044   0.00118   0.00118  0.002146  0.002146  0.031126  0.031126
	5.000  1.525   0.044   0.00203   0.00203  0.002094  0.002094  0.030766  0.030766
	5.000  1.565   0.043   0.00289   0.00289  0.002043  0.002043  0.030395  0.030395
	5.000  1.605   0.042   0.00364   0.00364  0.001985  0.001985  0.030012  0.030012
	5.000  1.645   0.041   0.00461   0.00461  0.001915  0.001915  0.029825  0.029825
	5.000  1.686   0.041   0.00533   0.00533  0.001765  0.001765  0.028232  0.028232
	5.000  1.722   0.040   0.00507   0.00507  0.001497  0.001497  0.024256  0.024256
	5.000  1.753   0.040   0.00453   0.00453  0.001220  0.001220  0.019972  0.019972
	5.000  1.778   0.039   0.00398   0.00398  0.001000  0.001000  0.016531  0.016531
	5.000  1.799   0.039   0.00368   0.00368  0.000871  0.000871  0.014435  0.014435
	5.000  1.819   0.039   0.00376   0.00376  0.000841  0.000841  0.013927  0.013927
	5.000  1.839   0.038   0.00397   0.00397  0.000840  0.000840  0.013899  0.013899
	5.000  1.859   0.038   0.00418   0.00418  0.000835  0.000835  0.013790  0.013790
	5.000  1.879   0.038   0.00437   0.00437  0.000831  0.000831  0.013702  0.013702
	5.000  1.899   0.037   0.00457   0.00457  0.000826  0.000826  0.013603  0.013603
	5.000  1.919   0.037   0.00475   0.00475  0.000822  0.000822  0.013492  0.013492
	5.000  1.939   0.037   0.00486   0.00486  0.000819  0.000819  0.013403  0.013403
	5.000  1.959   0.036   0.00499   0.00499  0.000815  0.000815  0.013303  0.013303
	5.000  1.979   0.036   0.00511   0.00511  0.000811  0.000811  0.013202  0.013202
	5.000  1.999   0.035   0.00525   0.00525  0.000807  0.000807  0.013088  0.013088
	5.000  2.019   0.035   0.00531   0.00531  0.000802  0.000802  0.012980  0.012980
	5.000  2.039   0.035   0.00556   0.00556  0.000824  0.000824  0.013325  0.013325
	5.000  2.060   0.035   0.00642   0.00642  0.000934  0.000934  0.015088  0.015088
	5.000  2.085   0.035   0.00783   0.00783  0.001116  0.001116  0.017983  0.017983
	5.000  2.116   0.034   0.00961   0.00961  0.001337  0.001337  0.021506  0.021506
	5.000  2.153   0.034   0.01144   0.01144  0.001539  0.001539  0.024578  0.024578
	5.000  2.193   0.034   0.01240   0.01240  0.001608  0.001608  0.025434  0.025434
	5.000  2.233   0.033   0.01276   0.01276  0.001604  0.001604  0.025084  0.025084
	5.000  2.273   0.033   0.01348   0.01348  0.001664  0.001664  0.024889  0.024889
	5.000  2.313   0.033   0.01423   0.01423  0.001726  0.001726  0.024691  0.024691
	5.000  2.353   0.032   0.01460   0.01460  0.001753  0.001753  0.024495  0.024495
	5.000  2.393   0.032   0.01471   0.01471  0.001776  0.001776  0.024294  0.024294
	5.000  2.433   0.032   0.01484   0.01484  0.001797  0.001797  0.024091  0.024091
	5.000  2.473   0.031   0.01511   0.01511  0.001761  0.001761  0.023749  0.023749
	5.000  2.513   0.031   0.01541   0.01541  0.001702  0.001702  0.023323  0.023323
	5.000  2.553   0.031   0.01568   0.01568  0.001651  0.001651  0.022915  0.022915
	5.000  2.593   0.030   0.01575   0.01575  0.001593  0.001593  0.021984  0.021984
	5.000  2.633   0.030   0.01581   0.01581  0.001532  0.001532  0.021004  0.021004
	5.000  2.673   0.029   0.01587   0.01587  0.001484  0.001484  0.020122  0.020122
	5.000  2.713   0.029   0.01570   0.01570  0.001411  0.001411  0.018469  0.018469
	5.000  2.753   0.029   0.01554   0.01554  0.001339  0.001339  0.016838  0.016838
	5.000  2.793   0.028   0.20665   0.20665  0.019244  0.019244  0.200518  0.200518
	5.000  4.000   0.018   1.02063   1.02063  0.047958  0.047958  0.436663  0.436663
	5.000  6.000   0.011   1.74305   1.74305  0.040229  0.040229  0.342797  0.342797
	5.000  8.000   0.005   8.32753   8.32753  0.068575  0.068575  0.457894  0.457894
	5.000 15.000   0.001  15.53560  15.53560  0.021386  0.021386  0.090833  0.090833
	7.500  0.000   0.070   0.00000   0.00000  0.000000  0.000000  0.000008  0.000008
	7.500  0.010   0.070   0.00000   0.00000  0.000001  0.000001  0.000023  0.000023
	7.500  0.032   0.070   0.00000   0.00000  0.000001  0.000001  0.000031  0.000031
	7.500  0.057   0.069   0.00000   0.00000  0.000001  0.000001  0.000031  0.000031
	7.500  0.088   0.069   0.00000   0.00000  0.000001  0.000001  0.000028  0.000028
	7.500  0.125   0.068   0.00000   0.00000  0.000001  0.000001  0.000021  0.000021
	7.500  0.165   0.067   0.00000   0.00000  0.000000  0.000000  0.000013  0.000013
	7.500  0.205   0.067   0.00000   0.00000  0.000000  0.000000  0.000008  0.000008
	7.500  0.245   0.066   0.00000   0.00000  0.000000  0.000000  0.000004  0.000004
	7.500  0.285   0.065   0.00000   0.00000  0.000000  0.000000  0.000002  0.000002
	7.500  0.325   0.064   0.00000   0.00000  0.000000  0.000000  0.000001  0.000001
	7.500  0.365   0.064   0.00000   0.00000  0.000000  0.000000  0.000001  0.000001
	7.500  0.405   0.063   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.445   0.062   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.485   0.062   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.525   0.061   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.565   0.060   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.605   0.059   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.645   0.059   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.685   0.058   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.725   0.057   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.765   0.057   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.805   0.056   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.845   0.055   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.885   0.055   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.925   0.054   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  0.965   0.053   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.005   0.052   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.045   0.052   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.085   0.051   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.125   0.050   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.165   0.050   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.205   0.049   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.245   0.048   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.285   0.048   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.325   0.047   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.365   0.046   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.405   0.046   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.445   0.045   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.485   0.044   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.525   0.044   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.565   0.043   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.605   0.042   0.00000   0.00000  0.000000  0.000000  0.000000  0.000000
	7.500  1.645  13.922   0.00000   0.00000  0.000020  0.000020  0.000000  0.000000
	7.500  1.686  13.627   0.00000   0.00000  0.000030  0.000030  0.000000  0.000000
	7.500  1.722   0.040  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.753   0.040  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.778   0.039  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.799   0.039  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.819   0.039  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.839   0.038  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.859   0.038  -0.00001  -0.00001  0.000000  0.000000  0.000000  0.000000
	7.500  1.879   0.038  -0.00001  -0.00001  0.000000  0.000000  0.000001  0.000001
	7.500  1.899   0.037  -0.00001  -0.00001  0.000000  0.000000  0.000001  0.000001
	7.500  1.919   0.037  -0.00001  -0.00001  0.000000  0.000000  0.000001  0.000001
	7.500  1.939   0.037   0.00000   0.00000  0.000000  0.000000  0.000001  0.000001
	7.500  1.959   3.754   0.00000   0.00000  0.000008  0.000008  0.000087  0.000087
	7.500  1.979  11.113   0.00000   0.00000  0.000025  0.000025  0.000290  0.000290
	7.500  1.999  16.797   0.00000   0.00000  0.000037  0.000037  0.000558  0.000558
	7.500  2.019  21.602   0.00000   0.00000  0.000047  0.000047  0.000878  0.000878
	7.500  2.039  26.094   0.00000   0.00000  0.000059  0.000059  0.001293  0.001293
	7.500  2.060  31.458   0.00000   0.00000  0.000091  0.000091  0.002265  0.002265
	7.500  2.085  37.802   0.00000   0.00000  0.000174  0.000174  0.004797  0.004797
	7.500  2.116  41.296   0.00000   0.00000  0.000307  0.000307  0.011118  0.011118
	7.500  2.153  35.998  -0.00001  -0.00001  0.000622  0.000622  0.022547  0.022547
	7.500  2.193  19.576  -0.00001  -0.00001  0.000442  0.000442  0.013923  0.013923
	7.500  2.233   5.480  -0.00001  -0.00001  0.000055  0.000055  0.001344  0.001344
	7.500  2.273   8.643   0.00000   0.00000  0.000062  0.000062  0.001546  0.001546
	7.500  2.313  27.426   0.00000   0.00000  0.000187  0.000187  0.005583  0.005583
	7.500  2.353  43.808   0.00000   0.00000  0.000826  0.000826  0.026383  0.026383
	7.500  2.393  52.030  -0.00003  -0.00003  0.005827  0.005827  0.173528  0.173528
	7.500  2.433  52.951  -0.00016  -0.00016  0.038880  0.038880  1.093014  1.093014
	7.500  2.473  41.103  -0.00048  -0.00048  0.170771  0.170771  4.562814  4.562814
	7.500  2.513  18.746   0.00067   0.00067  0.252853  0.252853  6.610350  6.610350
	7.500  2.553   0.031   0.01843   0.01843  0.001133  0.001133  0.015217  0.015217
	7.500  2.593   0.030   0.01866   0.01866  0.001130  0.001130  0.014640  0.014640
	7.500  2.633   0.030   0.01948   0.01948  0.001168  0.001168  0.015858  0.015858
	7.500  2.673   0.029   0.01897   0.01897  0.001161  0.001161  0.015789  0.015789
	7.500  2.713   0.029   0.01848   0.01848  0.001130  0.001130  0.015201  0.015201
	7.500  2.753   0.029   0.01959   0.01959  0.001099  0.001099  0.014617  0.014617
	7.500  2.793   0.028   0.29912   0.29912  0.017022  0.017022  0.222949  0.222949
	7.500  4.000   0.018   1.15987   1.15987  0.040118  0.040118  0.480068  0.480068
	7.500  6.000   0.011   2.23957   2.23957  0.045016  0.045016  0.487460  0.487460
	7.500  8.000   0.005   9.23782   9.23782  0.070440  0.070440  0.523122  0.523122
	7.500 15.000   0.001  14.77270  14.77270  0.019614  0.019614  0.092468  0.092468
	10.000  0.000   0.070  -0.00052  -0.00052  0.000078  0.000078  0.000000  0.000000
	10.000  0.010   0.070  -0.00172  -0.00172  0.000259  0.000259  0.000000  0.000000
	10.000  0.032   0.070  -0.00254  -0.00254  0.000397  0.000397  0.000001  0.000001
	10.000  0.057   0.069  -0.00294  -0.00294  0.000467  0.000467  0.000001  0.000001
	10.000  0.088   0.069  -0.00345  -0.00345  0.000557  0.000557  0.000003  0.000003
	10.000  0.125   0.068  -0.00387  -0.00387  0.000642  0.000642  0.000005  0.000005
	10.000  0.165   0.067  -0.00396  -0.00396  0.000666  0.000666  0.000008  0.000008
	10.000  0.205   0.067  -0.00387  -0.00387  0.000643  0.000643  0.000029  0.000029
	10.000  0.245   0.066  -0.00380  -0.00380  0.000616  0.000616  0.000056  0.000056
	10.000  0.285   0.065  -0.00370  -0.00370  0.000604  0.000604  0.000175  0.000175
	10.000  0.325   0.064  -0.00360  -0.00360  0.000619  0.000619  0.001139  0.001139
	10.000  0.365   0.064  -0.00391  -0.00391  0.000640  0.000640  0.002205  0.002205
	10.000  0.405   0.063  -0.00414  -0.00414  0.000704  0.000704  0.004269  0.004269
	10.000  0.445   0.062  -0.00344  -0.00344  0.000888  0.000888  0.009978  0.009978
	10.000  0.485   0.062  -0.00235  -0.00235  0.001097  0.001097  0.016402  0.016402
	10.000  0.525   0.061  -0.00199  -0.00199  0.001199  0.001199  0.018729  0.018729
	10.000  0.565   0.060  -0.00170  -0.00170  0.001290  0.001290  0.020697  0.020697
	10.000  0.605   0.059  -0.00223  -0.00223  0.001326  0.001326  0.021723  0.021723
	10.000  0.645   0.059  -0.00299  -0.00299  0.001330  0.001330  0.022275  0.022275
	10.000  0.685   0.058  -0.00321  -0.00321  0.001308  0.001308  0.022412  0.022412
	10.000  0.725   0.057  -0.00374  -0.00374  0.001264  0.001264  0.022227  0.022227
	10.000  0.765   0.057  -0.00413  -0.00413  0.001273  0.001273  0.023936  0.023936
	10.000  0.805   0.056  -0.00472  -0.00472  0.001280  0.001280  0.025410  0.025410
	10.000  0.845   0.055  -0.00495  -0.00495  0.001297  0.001297  0.026089  0.026089
	10.000  0.885   0.055  -0.00497  -0.00497  0.001329  0.001329  0.026960  0.026960
	10.000  0.925   0.054  -0.00496  -0.00496  0.001361  0.001361  0.027796  0.027796
	10.000  0.965   0.053  -0.00499  -0.00499  0.001367  0.001367  0.027737  0.027737
	10.000  1.005   0.052  -0.00543  -0.00543  0.001354  0.001354  0.027273  0.027273
	10.000  1.045   0.052  -0.00563  -0.00563  0.001351  0.001351  0.026963  0.026963
	10.000  1.085   0.051  -0.00601  -0.00601  0.001316  0.001316  0.025480  0.025480
	10.000  1.125   0.050  -0.00661  -0.00661  0.001279  0.001279  0.023949  0.023949
	10.000  1.165   0.050  -0.00690  -0.00690  0.001295  0.001295  0.023920  0.023920
	10.000  1.205   0.049  -0.00673  -0.00673  0.001343  0.001343  0.024680  0.024680
	10.000  1.245   0.048  -0.00650  -0.00650  0.001392  0.001392  0.025471  0.025471
	10.000  1.285   0.048  -0.00568  -0.00568  0.001452  0.001452  0.026579  0.026579
	10.000  1.325   0.047  -0.00479  -0.00479  0.001481  0.001481  0.027170  0.027170
	10.000  1.365   0.046  -0.00440  -0.00440  0.001481  0.001481  0.027256  0.027256
	10.000  1.405   0.046  -0.00371  -0.00371  0.001508  0.001508  0.028268  0.028268
	10.000  1.445   0.045  -0.00295  -0.00295  0.001537  0.001537  0.029329  0.029329
	10.000  1.485   0.044  -0.00234  -0.00234  0.001547  0.001547  0.029723  0.029723
	10.000  1.525   0.044  -0.00175  -0.00175  0.001555  0.001555  0.030026  0.030026
	10.000  1.565   0.043  -0.00115  -0.00115  0.001561  0.001561  0.030291  0.030291
	10.000  1.605   0.042  -0.00095  -0.00095  0.001552  0.001552  0.030012  0.030012
	10.000  1.645   0.041  -0.00117  -0.00117  0.001550  0.001550  0.029825  0.029825
	10.000  1.686   0.041  -0.00133  -0.00133  0.001475  0.001475  0.028232  0.028232
	10.000  1.722   0.040  -0.00133  -0.00133  0.001273  0.001273  0.024256  0.024256
	10.000  1.753   0.040  -0.00123  -0.00123  0.001052  0.001052  0.019972  0.019972
	10.000  1.778   0.039  -0.00111  -0.00111  0.000874  0.000874  0.016531  0.016531
	10.000  1.799   0.039  -0.00104  -0.00104  0.000765  0.000765  0.014435  0.014435
	10.000  1.819   0.039  -0.00107  -0.00107  0.000740  0.000740  0.013927  0.013927
	10.000  1.839   0.038  -0.00113  -0.00113  0.000741  0.000741  0.013899  0.013899
	10.000  1.859   0.038  -0.00118  -0.00118  0.000737  0.000737  0.013790  0.013790
	10.000  1.879   0.038  -0.00124  -0.00124  0.000734  0.000734  0.013702  0.013702
	10.000  1.899   0.037  -0.00129  -0.00129  0.000731  0.000731  0.013603  0.013603
	10.000  1.919   0.037  -0.00124  -0.00124  0.000731  0.000731  0.013492  0.013492
	10.000  1.939   0.037  -0.00117  -0.00117  0.000732  0.000732  0.013403  0.013403
	10.000  1.959   0.036  -0.00109  -0.00109  0.000733  0.000733  0.013303  0.013303
	10.000  1.979   0.036  -0.00101  -0.00101  0.000733  0.000733  0.013202  0.013202
	10.000  1.999   0.035  -0.00091  -0.00091  0.000734  0.000734  0.013088  0.013088
	10.000  2.019   0.035  -0.00089  -0.00089  0.000731  0.000731  0.012980  0.012980
	10.000  2.039   0.035  -0.00091  -0.00091  0.000752  0.000752  0.013325  0.013325
	10.000  2.060   0.035  -0.00103  -0.00103  0.000855  0.000855  0.015088  0.015088
	10.000  2.085   0.035  -0.00123  -0.00123  0.001023  0.001023  0.017983  0.017983
	10.000  2.116   0.034  -0.00147  -0.00147  0.001229  0.001229  0.021506  0.021506
	10.000  2.153   0.034  -0.00174  -0.00174  0.001409  0.001409  0.024578  0.024578
	10.000  2.193   0.034  -0.00186  -0.00186  0.001462  0.001462  0.025434  0.025434
	10.000  2.233   0.033  -0.00188  -0.00188  0.001447  0.001447  0.025084  0.025084
	10.000  2.273   0.033  -0.00185  -0.00185  0.001438  0.001438  0.024889  0.024889
	10.000  2.313   0.033  -0.00181  -0.00181  0.001430  0.001430  0.024691  0.024691
	10.000  2.353   0.032  -0.00172  -0.00172  0.001425  0.001425  0.024495  0.024495
	10.000  2.393   0.032  -0.00150  -0.00150  0.001426  0.001426  0.024294  0.024294
	10.000  2.433   0.032  -0.00126  -0.00126  0.001426  0.001426  0.024091  0.024091
	10.000  2.473   0.031  -0.00107  -0.00107  0.001422  0.001422  0.023884  0.023884
	10.000  2.513   0.031  -0.00088  -0.00088  0.001419  0.001419  0.023678  0.023678
	10.000  2.553   0.031  -0.00070  -0.00070  0.001414  0.001414  0.023468  0.023468
	10.000  2.593   0.030  -0.00066  -0.00066  0.001400  0.001400  0.023256  0.023256



In [5]:
# Make a log-dictionary
log_dict = read_log(flog)


['FDsfcalb', 'ADsfcalb:', '1.373958E+03', '1.373958E+03']

In [6]:
nz = np.unique(log_dict["Z"]).size
nx = np.unique(log_dict["X"]).size

x = log_dict['X'].reshape(nx, nz)
z = log_dict["Z"].reshape(nx, nz)
ext = log_dict['Extinct'].reshape(nx, nz)

In [7]:
f0 = plt.figure(0, (15,8), facecolor='white')
ax0 = f0.add_axes([.1, .1, .8, .8])
contour0 = ax0.contourf(x, z, ext, 100, cmap=cm.gray)
ax0.set_ybound((0,5))
f0.colorbar(contour0, ax=ax0)


Out[7]:
<matplotlib.colorbar.Colorbar instance at 0x979c4560>

In [18]:
# The the adjoint results
fadj = "radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.adj"
print_head(fadj, Nlines=100)


Showing head: radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.adj
	! SHDOMADJ optical property adjoint output
	!   X       Y       Z     ADEXTINCT   ADALBEDO   ADSINGSCAT(1..NSTSING,1..NRADDIR)
	0.0000  0.0000  0.0000 -7.220E-06  1.384E-06  2.540E-05  2.610E-05  8.162E-06
	0.0000  0.0000  0.0095 -2.529E-05  4.514E-06  8.070E-05  8.761E-05  2.636E-05
	0.0000  0.0000  0.0315 -4.138E-05  6.660E-06  1.110E-04  1.388E-04  3.794E-05
	0.0000  0.0000  0.0574 -5.251E-05  7.601E-06  1.131E-04  1.745E-04  4.170E-05
	0.0000  0.0000  0.0878 -7.434E-05  1.058E-05  1.803E-04  2.251E-04  4.762E-05
	0.0000  0.0000  0.1245 -9.725E-05  1.517E-05  3.012E-04  2.890E-04  5.036E-05
	0.0000  0.0000  0.1650 -1.076E-04  2.087E-05  4.663E-04  3.559E-04  4.724E-05
	0.0000  0.0000  0.2050 -1.844E-04  3.026E-05  7.619E-04  4.354E-04  4.552E-05
	0.0000  0.0000  0.2450 -2.578E-04  4.331E-05  1.200E-03  5.195E-04  3.654E-05
	0.0000  0.0000  0.2850 -2.246E-04  5.936E-05  1.716E-03  6.427E-04  3.027E-05
	0.0000  0.0000  0.3250 -1.412E-04  7.413E-05  2.134E-03  8.084E-04  3.024E-05
	0.0000  0.0000  0.3650 -1.356E-04  8.567E-05  2.457E-03  9.419E-04  2.570E-05
	0.0000  0.0000  0.4050 -8.608E-05  9.922E-05  2.855E-03  1.072E-03  2.404E-05
	0.0000  0.0000  0.4450 -2.021E-05  1.154E-04  3.374E-03  1.179E-03  2.132E-05
	0.0000  0.0000  0.4850 -4.900E-05  1.254E-04  3.715E-03  1.217E-03  1.692E-05
	0.0000  0.0000  0.5250 -4.602E-05  1.389E-04  4.198E-03  1.239E-03  1.550E-05
	0.0000  0.0000  0.5650  8.894E-06  1.581E-04  4.894E-03  1.258E-03  1.395E-05
	0.0000  0.0000  0.6050  3.360E-05  1.740E-04  5.481E-03  1.257E-03  1.246E-05
	0.0000  0.0000  0.6450  6.393E-04  2.332E-04  7.711E-03  1.248E-03  1.077E-05
	0.0000  0.0000  0.6850  1.151E-03  2.967E-04  1.009E-02  1.239E-03  8.987E-06
	0.0000  0.0000  0.7250  1.600E-03  3.482E-04  1.200E-02  1.229E-03  7.861E-06
	0.0000  0.0000  0.7650  3.131E-03  4.495E-04  1.576E-02  1.221E-03  6.620E-06
	0.0000  0.0000  0.8050  4.719E-03  5.492E-04  1.944E-02  1.211E-03  5.392E-06
	0.0000  0.0000  0.8450  4.754E-03  5.510E-04  1.945E-02  1.202E-03  4.632E-06
	0.0000  0.0000  0.8850  4.267E-03  5.130E-04  1.798E-02  1.188E-03  3.606E-06
	0.0000  0.0000  0.9250  3.566E-03  4.646E-04  1.616E-02  1.140E-03  2.605E-06
	0.0000  0.0000  0.9650  2.669E-03  4.045E-04  1.400E-02  1.017E-03  1.825E-06
	0.0000  0.0000  1.0050  1.788E-03  3.196E-04  1.099E-02  8.265E-04  1.018E-06
	0.0000  0.0000  1.0450  9.144E-04  2.357E-04  8.067E-03  6.200E-04  5.329E-07
	0.0000  0.0000  1.0850  3.657E-04  1.734E-04  5.907E-03  4.573E-04  2.676E-07
	0.0000  0.0000  1.1250  1.133E-04  1.320E-04  4.476E-03  3.492E-04  1.058E-07
	0.0000  0.0000  1.1650 -1.448E-04  9.918E-05  3.337E-03  2.762E-04  5.673E-08
	0.0000  0.0000  1.2050 -4.173E-04  6.713E-05  2.221E-03  2.181E-04  3.659E-08
	0.0000  0.0000  1.2450 -5.764E-04  4.184E-05  1.355E-03  1.617E-04  2.135E-08
	0.0000  0.0000  1.2850 -6.080E-04  2.605E-05  8.213E-04  1.206E-04  1.112E-08
	0.0000  0.0000  1.3250 -6.461E-04  1.485E-05  4.440E-04  9.330E-05  3.615E-09
	0.0000  0.0000  1.3650 -6.059E-04  8.969E-06  2.482E-04  7.702E-05  0.000E+00
	0.0000  0.0000  1.4050 -5.095E-04  6.656E-06  1.724E-04  6.902E-05  0.000E+00
	0.0000  0.0000  1.4450 -4.415E-04  4.650E-06  1.073E-04  6.204E-05  0.000E+00
	0.0000  0.0000  1.4850 -4.212E-04  3.592E-06  7.449E-05  5.658E-05  0.000E+00
	0.0000  0.0000  1.5250 -2.561E-04  1.453E-03  1.199E-02  9.328E-03  0.000E+00
	0.0000  0.0000  1.5650 -2.330E-04  6.042E-03  6.499E-02  5.776E-02  1.270E-07
	0.0000  0.0000  1.6050 -2.632E-04  1.743E-02  3.257E-01  2.791E-01  1.147E-05
	0.0000  0.0000  1.6450 -4.341E-04  6.245E-02  1.595E+00  1.244E+00  9.209E-05
	0.0000  0.0000  1.6855 -5.292E-04  1.846E-01  5.457E+00  3.932E+00  4.222E-04
	0.0000  0.0000  1.7222  1.418E-03  9.361E-02  2.663E+00  1.810E+00  2.968E-04
	0.0000  0.0000  1.7526  5.010E-03  6.642E-04  1.381E-02  9.481E-03  2.308E-06
	0.0000  0.0000  1.7780  4.485E-03  5.474E-04  1.132E-02  7.807E-03  2.373E-06
	0.0000  0.0000  1.7992  4.151E-03  4.744E-04  9.757E-03  6.762E-03  2.920E-06
	0.0000  0.0000  1.8190  4.177E-03  4.515E-04  9.196E-03  6.487E-03  6.316E-06
	0.0000  0.0000  1.8390  4.364E-03  4.455E-04  8.988E-03  6.444E-03  9.434E-06
	0.0000  0.0000  1.8590  4.521E-03  4.372E-04  8.731E-03  6.364E-03  1.244E-05
	0.0000  0.0000  1.8790  4.644E-03  4.294E-04  8.491E-03  6.294E-03  1.539E-05
	0.0000  0.0000  1.8990  4.663E-03  4.213E-04  8.242E-03  6.220E-03  1.840E-05
	0.0000  0.0000  1.9190  4.704E-03  4.143E-04  7.998E-03  6.141E-03  6.618E-05
	0.0000  0.0000  1.9390  4.751E-03  4.074E-04  7.748E-03  6.073E-03  1.238E-04
	0.0000  0.0000  1.9590  4.812E-03  4.007E-04  7.507E-03  6.003E-03  1.775E-04
	0.0000  0.0000  1.9790  4.897E-03  3.941E-04  7.272E-03  5.935E-03  2.300E-04
	0.0000  0.0000  1.9990  4.977E-03  3.863E-04  7.001E-03  5.849E-03  2.871E-04
	0.0000  0.0000  2.0190  5.019E-03  3.791E-04  6.763E-03  5.749E-03  3.885E-04
	0.0000  0.0000  2.0388  5.329E-03  3.862E-04  6.680E-03  5.839E-03  6.661E-04
	0.0000  0.0000  2.0600  6.193E-03  4.332E-04  7.268E-03  6.531E-03  1.034E-03
	0.0000  0.0000  2.0854  7.533E-03  5.115E-04  8.267E-03  7.679E-03  1.624E-03
	0.0000  0.0000  2.1158  2.692E-03  9.606E-02  2.421E+00  2.335E+00  6.271E-01
	0.0000  0.0000  2.1525  8.069E-03  4.981E-02  9.590E-01  9.298E-01  3.443E-01
	0.0000  0.0000  2.1930  1.468E-02  1.077E-03  1.453E-02  1.541E-02  7.689E-03
	0.0000  0.0000  2.2330  6.943E-03  7.514E-02  1.557E+00  1.436E+00  8.072E-01
	0.0000  0.0000  2.2730  1.516E-02  1.266E-03  1.741E-02  1.729E-02  9.232E-03
	0.0000  0.0000  2.3130  1.479E-02  1.310E-03  1.912E-02  1.720E-02  8.744E-03
	0.0000  0.0000  2.3530  1.384E-02  1.302E-03  1.909E-02  1.706E-02  8.447E-03
	0.0000  0.0000  2.3930  1.369E-02  1.277E-03  1.855E-02  1.692E-02  8.151E-03
	0.0000  0.0000  2.4330  1.435E-02  1.254E-03  1.803E-02  1.677E-02  7.865E-03
	0.0000  0.0000  2.4730  1.513E-02  1.234E-03  1.746E-02  1.662E-02  7.818E-03
	0.0000  0.0000  2.5130  1.568E-02  1.210E-03  1.672E-02  1.648E-02  7.832E-03
	0.0000  0.0000  2.5530  1.625E-02  1.189E-03  1.609E-02  1.633E-02  7.839E-03
	0.0000  0.0000  2.5930  1.738E-02  1.190E-03  1.595E-02  1.618E-02  8.099E-03
	0.0000  0.0000  2.6330  1.854E-02  1.192E-03  1.582E-02  1.603E-02  8.381E-03
	0.0000  0.0000  2.6730  1.931E-02  1.189E-03  1.564E-02  1.587E-02  8.499E-03
	0.0000  0.0000  2.7130  1.914E-02  1.180E-03  1.535E-02  1.572E-02  8.548E-03
	0.0000  0.0000  2.7530  1.916E-02  1.172E-03  1.506E-02  1.557E-02  8.592E-03
	0.0000  0.0000  2.7930  2.928E-01  1.879E-02  2.459E-01  2.428E-01  1.375E-01
	0.0000  0.0000  4.0000  1.119E+00  4.288E-02  5.544E-01  4.239E-01  2.061E-01
	0.0000  0.0000  6.0000  2.049E+00  4.194E-02  4.437E-01  3.395E-01  1.858E-01
	0.0000  0.0000  8.0000  9.285E+00  7.213E-02  4.794E-01  3.637E-01  3.140E-01
	0.0000  0.0000 15.0000  1.602E+01  2.154E-02  7.093E-02  6.365E-02  7.774E-02
	0.0625  0.0000  0.0000 -1.450E-05  2.022E-06  3.898E-05  3.226E-05  1.621E-05
	0.0625  0.0000  0.0095 -5.198E-05  7.255E-06  1.444E-04  1.123E-04  5.601E-05
	0.0625  0.0000  0.0315 -8.232E-05  1.328E-05  2.817E-04  1.922E-04  9.432E-05
	0.0625  0.0000  0.0574 -1.021E-04  1.907E-05  4.315E-04  2.553E-04  1.227E-04
	0.0625  0.0000  0.0878 -1.636E-04  2.955E-05  7.370E-04  3.484E-04  1.527E-04
	0.0625  0.0000  0.1245 -1.980E-04  4.840E-05  1.328E-03  4.875E-04  1.839E-04
	0.0625  0.0000  0.1650 -1.232E-04  7.422E-05  2.187E-03  6.426E-04  1.985E-04
	0.0625  0.0000  0.2050 -1.096E-04  8.572E-05  2.540E-03  7.453E-04  1.953E-04
	0.0625  0.0000  0.2450 -1.235E-04  9.611E-05  2.860E-03  8.490E-04  1.743E-04
	0.0625  0.0000  0.2850 -1.411E-04  1.055E-04  3.153E-03  9.359E-04  1.533E-04
	0.0625  0.0000  0.3250 -1.860E-04  1.145E-04  3.493E-03  9.664E-04  1.174E-04
	0.0625  0.0000  0.3650 -2.563E-04  1.228E-04  3.815E-03  9.854E-04  7.617E-05
	0.0625  0.0000  0.4050 -2.974E-04  1.325E-04  4.176E-03  1.002E-03  5.797E-05



In [9]:
# Read the adjoint file

cols = ['X', 'Y', 'Z', 'ADEXTINCT', "ADALBEDO"]
value_list = []
for line in open(fadj):
    if line[0] =='!':
        pass
    else:
        new = line.strip().split()
        value_list.append(new[:len(cols)])
        
values = np.array(value_list, dtype="f8").T
adjoint_dict = dict(zip(cols, values))

In [10]:
# Get values to plot
nx = np.unique(adjoint_dict['X']).size
nz = np.unique(adjoint_dict['Z']).size
x = adjoint_dict['X'].reshape(nx, nz)
z = adjoint_dict['Z'].reshape(nx, nz)
adj_ext = adjoint_dict['ADEXTINCT'].reshape(nx, nz)
adj_alb = adjoint_dict['ADALBEDO'].reshape(nx, nz)

print("maxes: {}, {}.".format(abs(adj_ext).max(),abs(adj_alb).max()))


maxes: 16.91, 0.501.

In [11]:
#plot the adjoint of the extinction
f1 = plt.figure(1, (15, 8), facecolor='white')
ax1 = f1.add_axes([.1, .1, .8, .8])
contour1 =ax1.contourf(x, z, adj_alb, levels=np.linspace(0, .6,50), cmap=cm.gray)
ax1.set_ybound(0,5)
f1.colorbar(contour1, ax=ax1)

ax1.set_title("Single-scattering adjoint derivative w.r.t. albedo", fontsize="x-large", fontweight="bold")


Out[11]:
<matplotlib.text.Text at 0x4f0a910>

In [19]:
# The first adjoint radiance field
farad1 = "radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.arad"
print_head(farad1, Nlines=100)


Showing head: radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss1.arad
	! Polarized Spherical Harmonic Discrete Ordinate Radiative Transfer Radiance Output
	!  L= 15  M= 15  NLM=  256   NMU= 16  NPHI= 32  NANG=  354   NSH=   1757885
	!  NSTOKES=1   NX= 320   NY=   1   NZ=  85    NPTS=   27285   NCELLS=   26880
	!  PARTICLE_FILE=les0822nh15t13y89_ocaer1.part
	!  CORRELATED_K-DIST_FILE=NONE   NUM_G= 1
	!  SOURCE_TYPE=SOLAR               DELTA-M METHOD
	!  GRID_TYPE=EVEN (X,Y)    PROPERTY-FILE (Z)    INDEPENDENT_PIXEL=2
	!  SURFACE_TYPE=FIXED LAMBERTIAN         HORIZ_BOUNDARY_COND=0
	!  GROUND_ALBEDO=0.0500000  SKY_RAD= 0.00000E+00
	!  SOLAR_FLUX= 0.100000E+01   SOLAR_MU=-0.9221000   SOLAR_AZ= 339.400
	!  UNITS=WATTS/(M^2 MICRON STER)    WAVELENGTH=      0.65
	!  SPLITTING_ACCURACY= 0.300E-01   SPHERICAL_HARMONIC_ACCURACY= 0.300E-02
	!  SOLUTION_ACCURACY= 0.300E-04
	!  MAXIMUM_ITERATIONS=   0   NUMBER_ITERATIONS=   0
	!  RADIANCE AT Z= 15.000
	320    1    3 ! NXOUT, NYOUT, NRADDIR
	!  ID      mu       phi
	!   1  0.8660300  180.00
	!   2  1.0000000    0.00
	!   3  0.8660300    0.00
	! ID  IY  IX    RADIANCE (STOKES I, Q, U, V)
	1   1   1   0.134613E-01
	1   1   2   0.115650E-01
	1   1   3   0.125559E-01
	1   1   4   0.125286E-01
	1   1   5   0.121647E-01
	1   1   6   0.117333E-01
	1   1   7   0.124280E-01
	1   1   8   0.115808E-01
	1   1   9   0.988294E-02
	1   1  10   0.123670E-01
	1   1  11   0.118903E-01
	1   1  12   0.117026E-01
	1   1  13   0.114063E-01
	1   1  14   0.110016E-01
	1   1  15   0.109395E-01
	1   1  16   0.111800E-01
	1   1  17   0.119092E-01
	1   1  18   0.124901E-01
	1   1  19   0.131205E-01
	1   1  20   0.126538E-01
	1   1  21   0.138799E-01
	1   1  22   0.130093E-01
	1   1  23   0.133729E-01
	1   1  24   0.134214E-01
	1   1  25   0.127005E-01
	1   1  26   0.161253E-01
	1   1  27   0.191886E-01
	1   1  28   0.176198E-01
	1   1  29   0.179799E-01
	1   1  30   0.172136E-01
	1   1  31   0.174579E-01
	1   1  32   0.181361E-01
	1   1  33   0.185459E-01
	1   1  34   0.184355E-01
	1   1  35   0.184394E-01
	1   1  36   0.180923E-01
	1   1  37   0.178597E-01
	1   1  38   0.178598E-01
	1   1  39   0.178598E-01
	1   1  40   0.178598E-01
	1   1  41   0.178598E-01
	1   1  42   0.178598E-01
	1   1  43   0.178598E-01
	1   1  44   0.178598E-01
	1   1  45   0.178598E-01
	1   1  46   0.178598E-01
	1   1  47   0.178598E-01
	1   1  48   0.178598E-01
	1   1  49   0.178598E-01
	1   1  50   0.178598E-01
	1   1  51   0.178598E-01
	1   1  52   0.178598E-01
	1   1  53   0.178598E-01
	1   1  54   0.178598E-01
	1   1  55   0.178598E-01
	1   1  56   0.178598E-01
	1   1  57   0.178598E-01
	1   1  58   0.178598E-01
	1   1  59   0.178598E-01
	1   1  60   0.178598E-01
	1   1  61   0.178598E-01
	1   1  62   0.178598E-01
	1   1  63   0.178598E-01
	1   1  64   0.178598E-01
	1   1  65   0.178598E-01
	1   1  66   0.178598E-01
	1   1  67   0.178598E-01
	1   1  68   0.174253E-01
	1   1  69   0.151741E-01
	1   1  70   0.132923E-01
	1   1  71   0.148204E-01
	1   1  72   0.169102E-01
	1   1  73   0.174304E-01
	1   1  74   0.174491E-01
	1   1  75   0.174725E-01
	1   1  76   0.184029E-01
	1   1  77   0.181126E-01
	1   1  78   0.179639E-01
	1   1  79   0.178554E-01



In [13]:
# The second adjoint radiance field
farad2 = "radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss2.arad"
print_head(farad2)


Showing head: radobs_les0822nh15t13y89_ocaer1_w0.646_ns1ss2.arad
	! Polarized Spherical Harmonic Discrete Ordinate Radiative Transfer Radiance Output
	!  L= 15  M= 15  NLM=  256   NMU= 16  NPHI= 32  NANG=  354   NSH=   2422485
	!  NSTOKES=1   NX= 320   NY=   1   NZ=  85    NPTS=   30052   NCELLS=   31290
	!  PARTICLE_FILE=les0822nh15t13y89_ocaer1.part
	!  CORRELATED_K-DIST_FILE=NONE   NUM_G= 1
	!  SOURCE_TYPE=SOLAR               DELTA-M METHOD
	!  GRID_TYPE=EVEN (X,Y)    PROPERTY-FILE (Z)    INDEPENDENT_PIXEL=2
	!  SURFACE_TYPE=FIXED LAMBERTIAN         HORIZ_BOUNDARY_COND=0
	!  GROUND_ALBEDO=0.0500000  SKY_RAD= 0.00000E+00
	!  SOLAR_FLUX= 0.100000E+01   SOLAR_MU=-0.9221000   SOLAR_AZ= 339.400



In [14]:
# A sample modis file
fmodis_blue = "les0822nh15t13_aer008_w0.646rad.modis"
print_head(fmodis_blue)


Showing head: les0822nh15t13_aer008_w0.646rad.modis
	40  40  23
	! Simulated MODIS reflectances from SHDOM radiance output
	! IA   Rad(1:Ncross)
	1 1.0000  0.00   0.00 20.000 20.000 ! idir, mu, theta, phi_shdom, size_cross, size_along
	1 0.0765 0.1141 0.2162 0.1455 0.1407 0.0961 0.0401 0.0496 0.0690 0.0696 0.0722 0.0709 0.0778 0.0677 0.0694 0.0692 0.0725 0.1055 0.0934 0.0724 0.0701 0.0680 0.0680 0.0679 0.0680 0.0682 0.0684 0.0691 0.0747 0.0888 0.1230 0.0937 0.1201 0.0896 0.0944 0.0857 0.0592 0.0672 0.0700 0.0722
	2 0.0702 0.0723 0.0972 0.0969 0.2869 0.2755 0.0567 0.0559 0.0687 0.0636 0.0761 0.0977 0.0996 0.1314 0.0961 0.0598 0.0686 0.0689 0.0840 0.0941 0.0648 0.0692 0.0710 0.0670 0.0679 0.0678 0.0682 0.0687 0.0697 0.0782 0.1025 0.0723 0.0575 0.0718 0.0841 0.0720 0.0803 0.0733 0.0694 0.0699
	3 0.0697 0.0700 0.0712 0.0723 0.1212 0.1625 0.0773 0.0977 0.0573 0.0622 0.0979 0.0866 0.1896 0.2539 0.1146 0.1198 0.1418 0.0673 0.0754 0.1100 0.0931 0.0671 0.0703 0.0683 0.0682 0.0679 0.0679 0.0683 0.0690 0.0737 0.0956 0.0804 0.0713 0.0749 0.0772 0.0540 0.0561 0.0685 0.0757 0.0685
	4 0.0694 0.0696 0.0701 0.0710 0.0913 0.1209 0.3438 0.3862 0.1215 0.0807 0.2028 0.1568 0.0961 0.1207 0.1067 0.0881 0.0820 0.0575 0.0992 0.0793 0.0712 0.0679 0.0682 0.0682 0.0680 0.0680 0.0680 0.0681 0.0685 0.0693 0.0719 0.0790 0.0823 0.1356 0.1346 0.0854 0.0563 0.0640 0.0692 0.0697
	5 0.0604 0.0688 0.0695 0.0700 0.0711 0.0791 0.2566 0.4028 0.1625 0.1041 0.2124 0.2489 0.0943 0.0469 0.0659 0.0868 0.1006 0.1008 0.0864 0.0592 0.0674 0.0679 0.0683 0.0682 0.0680 0.0680 0.0681 0.0680 0.0681 0.0688 0.0703 0.0903 0.1479 0.1072 0.1253 0.0785 0.0636 0.0553 0.0680 0.0623
	6 0.0458 0.0667 0.0691 0.0695 0.0700 0.0715 0.0777 0.1253 0.0972 0.0941 0.2480 0.1939 0.0607 0.0509 0.0740 0.0917 0.0764 0.0891 0.0833 0.0703 0.0682 0.0671 0.0684 0.0683 0.0681 0.0681 0.0682 0.0684 0.0680 0.0683 0.0696 0.0764 0.1183 0.1811 0.1659 0.1587 0.0729 0.0854 0.1028 0.0425



In [15]:
# Reflectance from the MODIS case, but at full 320 x 320 resolution 
fhres_blue = "les0822nh15t13_aer008_w0.646.ref"
print_head(fhres_blue)


Showing head: les0822nh15t13_aer008_w0.646.ref
	0   0 0.0499 0.1004 0.2231 0.0945 0.0426 0.1586 0.0464
	1   0 0.0495 0.1004 0.2272 0.0976 0.0393 0.2537 0.0463
	2   0 0.0493 0.1004 0.2193 0.0935 0.0390 0.2863 0.0463
	3   0 0.0493 0.1004 0.2151 0.0913 0.0391 0.2772 0.0463
	4   0 0.0498 0.1004 0.2003 0.0937 0.0394 0.2357 0.0462
	5   0 0.0533 0.1003 0.1858 0.0979 0.0399 0.1731 0.0462
	6   0 0.0583 0.1003 0.1592 0.1063 0.0407 0.1175 0.0461
	7   0 0.0602 0.1002 0.1362 0.1148 0.0413 0.0771 0.0460
	8   0 0.0602 0.1002 0.1257 0.1197 0.0413 0.0479 0.0459
	9   0 0.0602 0.1001 0.1255 0.1194 0.0413 0.0425 0.0459



In [ ]:


In [ ]:


In [ ]: