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
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 [ ]:
Content source: wgm2111/wgm-shdom-files-AGU2015
Similar notebooks: