Python 2.7.6 (default, Nov 21 2013, 15:55:38) [GCC]

Type "copyright", "credits" or "license" for more information.


IPython 1.0.0 -- An enhanced Interactive Python.

? -> Introduction and overview of IPython's features.

%quickref -> Quick reference.

help -> Python's own help system.

object? -> Details about 'object', use 'object??' for extra details.

%guiref -> A brief reference about the graphical user interface.


In [1]: # ipython qtconsole --pylab=inline &


In [2]: %load ~/src/py4cats-lite/src/py4cats.py


import os, sys

from string import *

from glob import glob


py4cats.py --- catsPath /users/schreier/src/py4cats-lite

sys.path[0]

# INFO: inserting py4cats src dir to sys.path /users/schreier/src/py4cats-lite/src

===> 18 dirs


from ir import amu, c, h, k, pi

from aeiou import awrite, cstack, minmaxmean, loadRecArray, grep_from_header

from pairTypes import Interval

from moreFun import cosdg

from misc import approx, trapez, runningAverage, boxedSpec

from extract import *

from hitran import extract_hitran

from geisa import extract_geisa

from planck import planck

from cgsUnits import unitConversion

from atmos1D import *

from lines import atlas, get_dictOfLineLists, delete_traceGasLines, read_line_file, strengths_and_widths, meanLineWidths

from lbl2xs import lbl2xs

from lbl2od import lbl2od

from lbl2ac import *

from lbl2od import *

from xSection import xsArray, xsRead, xsSave, xsPlot

from oDepth import *

from ac2wf import *

from od2ri import *


INFO --- py4cats: setup done

18 directories in sys.path, the very first should be Py4CAtS' home:

sys.path[0] = /users/schreier/src/py4cats-lite/src


In [4]: rcParams['figure.figsize']= [10.0, 5.0] # default figures are pretty small


In [5]: mls = atmos1D('/users/schreier/src/py4cats-lite/data/atmos/20/mls.xy',zToA=50) # get atmosphere


Atmos1d: got p, T, air

WARNING --- truncateAtmos1d: zToA very small, assuming kilometer units

atmos1D.truncateAtmos1d --- ToA 50.0 km, skip 4 levels at top


In [6]: atmPlot(mls)



In [7]: mlw = atmos1D('/users/schreier/src/py4cats-lite/data/atmos/25/midLatWin.xy',zToA=50) # another atmosphere


Atmos1d: got p, T, air

WARNING --- truncateAtmos1d: zToA very small, assuming kilometer units

atmos1D.truncateAtmos1d --- ToA 50.0 km, skip 7 levels at top


In [8]: atmPlot([mls,mlw],'O3','mb') # ozone profiles vs pressure [mb]



In [9]: # IASI microwindow used for CO retrievals (HItran-GeiSa-exTRACT)

   ...: dll = higstract('/data/hitran/86/lines', (2100.,2150.), molecule='main') # Dictionary of Line Lists

# first line in spectral range at record number 208088 found in 0.883314 sec

31 2100.017000 4.730E-23 1.633E-06.0645.0000 699.63400.760.000000

# last line accepted

31 2149.993000 5.830E-24 1.778E-07.0642.0000 678.88000.760.000000

# first line not accepted

21 2150.055100 2.520E-26 2.519E-09.0710.1005 921.23790.750.000000

main gases only: deleted 5 of 11 linelists

3784 lines of 6 molecules extracted from /data/hitran/86/lines (returning a dictionary)


In [10]: atlas(dll) # visualization of strengths vs positions



In [11]: atlas(dll['CO'],'a') # CO air broadening



In [12]: # explore cross sections

    ...: xs = lbl2xs(dll['CO'])

lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.01325e+06

temperature [K] 296.0 ---> 296.0

cross section CO 66 lines, 1013.25mb, 296.0K, 3329 0.01sec 38.36ns): 7.26446e-22 < xs < 1.71411e-18 in Interval [2101.1,2149.49]



In [13]: xsPlot(xs)



In [14]: xss = lbl2xs(dll['CO'],[1013.25,100,10,'mb']) # remember default units are cgs

    ...: xsPlot(xss)

    ...:

lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.01325e+06

temperature [K] 296.0 ---> 296.0

cross section CO 66 lines, 1013.25mb, 296.0K, 3329 0.01sec 37.70ns): 7.26446e-22 < xs < 1.71411e-18 in Interval [2101.1,2149.49]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 100000

temperature [K] 296.0 ---> 296.0

cross section CO 66 lines, 100mb, 296.0K, 30529 0.01sec 6.40ns): 7.19519e-23 < xs < 1.60165e-17 in Interval [2101.1,2149.49]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 10000

temperature [K] 296.0 ---> 296.0

cross section CO 66 lines, 10mb, 296.0K, 84609 0.02sec 2.71ns): 7.19545e-24 < xs < 5.56761e-17 in Interval [2101.1,2149.49]




In [15]: xss = lbl2xs(dll['CO'],mls[::4]['p'],mls[::4]['T']) # cross sections for some levels

    ...: xsPlot(xss)

    ...:

lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section CO 66 lines, 1013mb, 294.2K, 3329 0.01sec 35.75ns): 7.31572e-22 < xs < 1.71077e-18 in Interval [2101.1,2149.49]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section CO 66 lines, 628mb, 273.2K, 5057 0.01sec 24.98ns): 4.98408e-22 < xs < 2.68058e-18 in Interval [2101.1,2149.49]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section CO 66 lines, 193.4mb, 219.1K, 14017 0.01sec 13.13ns): 2.07237e-22 < xs < 8.01903e-18 in Interval [2101.1,2149.49]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section CO 66 lines, 13.2mb, 233.7K, 85633 0.01sec 2.52ns): 1.29327e-23 < xs < 5.85091e-17 in Interval [2101.1,2149.49]




In [16]: dodList = lbl2od(mls,dll) # now turn to delta/layer optical depths for all molecules

lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section CO2 705 lines, 1013mb, 294.2K, 3393 0.08sec 31.45ns): 8.88655e-26 < xs < 1.54914e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section CO2 705 lines, 902mb, 289.7K, 3777 0.07sec 25.68ns): 7.24251e-26 < xs < 1.54729e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section CO2 705 lines, 802mb, 285.2K, 4225 0.07sec 23.10ns): 5.87397e-26 < xs < 1.54293e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section CO2 705 lines, 710mb, 279.2K, 4673 0.07sec 21.04ns): 4.57159e-26 < xs < 1.5133e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section CO2 705 lines, 628mb, 273.2K, 5185 0.07sec 19.01ns): 3.53722e-26 < xs < 1.47793e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section CO2 705 lines, 554mb, 267.2K, 5825 0.07sec 17.12ns): 2.71642e-26 < xs < 1.43783e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section CO2 705 lines, 398.1mb, 251.4K, 7745 0.10sec 18.18ns): 1.32559e-26 < xs < 1.31071e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section CO2 705 lines, 281mb, 235.3K, 10433 0.12sec 16.41ns): 6.07967e-27 < xs < 1.15254e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section CO2 705 lines, 193.4mb, 219.1K, 14273 0.12sec 12.06ns): 2.61757e-27 < xs < 9.7741e-23 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section CO2 705 lines, 130mb, 215.7K, 20801 0.12sec 8.38ns): 1.59094e-27 < xs < 1.0715e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section CO2 705 lines, 59.5mb, 219.2K, 42945 0.13sec 4.15ns): 8.09794e-28 < xs < 2.4255e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section CO2 705 lines, 27.7mb, 225.1K, 75201 0.14sec 2.55ns): 4.49008e-28 < xs < 5.07979e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section CO2 705 lines, 13.2mb, 233.7K, 103937 0.14sec 1.88ns): 2.7379e-28 < xs < 9.17419e-22 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section CO2 705 lines, 6.52mb, 245.2K, 119937 0.14sec 1.68ns): 1.85255e-28 < xs < 1.40157e-21 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section CO2 705 lines, 3.33mb, 257.5K, 126465 0.14sec 1.61ns): 1.30066e-28 < xs < 1.86777e-21 in Interval [2100.02,2149.99]


lbl2xs: CO2 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section CO2 705 lines, 0.951mb, 275.7K, 128833 0.15sec 1.63ns): 5.75268e-29 < xs < 2.50727e-21 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section N2O 677 lines, 1013mb, 294.2K, 3521 0.06sec 27.19ns): 6.52062e-25 < xs < 9.65226e-21 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section N2O 677 lines, 902mb, 289.7K, 3969 0.06sec 24.18ns): 5.41244e-25 < xs < 1.01982e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section N2O 677 lines, 802mb, 285.2K, 4353 0.07sec 22.10ns): 4.48162e-25 < xs < 1.07297e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section N2O 677 lines, 710mb, 279.2K, 4865 0.07sec 20.03ns): 3.60023e-25 < xs < 1.14044e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section N2O 677 lines, 628mb, 273.2K, 5441 0.07sec 17.78ns): 2.8875e-25 < xs < 1.20775e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section N2O 677 lines, 554mb, 267.2K, 6017 0.07sec 16.50ns): 2.30871e-25 < xs < 1.28637e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section N2O 677 lines, 398.1mb, 251.4K, 8001 0.11sec 21.22ns): 1.28091e-25 < xs < 1.53816e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section N2O 677 lines, 281mb, 235.3K, 10753 0.11sec 15.79ns): 6.94859e-26 < xs < 1.85275e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section N2O 677 lines, 193.4mb, 219.1K, 14721 0.11sec 11.49ns): 3.68081e-26 < xs < 2.28489e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section N2O 677 lines, 130mb, 215.7K, 21441 0.12sec 8.01ns): 2.34531e-26 < xs < 3.2371e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section N2O 677 lines, 59.5mb, 219.2K, 44097 0.12sec 4.01ns): 1.13591e-26 < xs < 6.98811e-20 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section N2O 677 lines, 27.7mb, 225.1K, 76545 0.13sec 2.47ns): 5.81618e-27 < xs < 1.36407e-19 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section N2O 677 lines, 13.2mb, 233.7K, 104705 0.13sec 1.88ns): 3.18616e-27 < xs < 2.19468e-19 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section N2O 677 lines, 6.52mb, 245.2K, 120129 0.14sec 1.69ns): 1.89985e-27 < xs < 2.95006e-19 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section N2O 677 lines, 3.33mb, 257.5K, 126273 0.14sec 1.61ns): 1.18983e-27 < xs < 3.49774e-19 in Interval [2100.02,2149.99]


lbl2xs: N2O 44.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section N2O 677 lines, 0.951mb, 275.7K, 128385 0.14sec 1.61ns): 4.60904e-28 < xs < 4.09404e-19 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section H2O 149 lines, 1013mb, 294.2K, 3713 0.02sec 32.97ns): 3.04999e-25 < xs < 2.95441e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section H2O 149 lines, 902mb, 289.7K, 4097 0.02sec 31.50ns): 2.50578e-25 < xs < 3.10208e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section H2O 149 lines, 802mb, 285.2K, 4545 0.02sec 29.58ns): 2.05071e-25 < xs < 3.13642e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section H2O 149 lines, 710mb, 279.2K, 5057 0.02sec 28.05ns): 1.61898e-25 < xs < 3.11693e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section H2O 149 lines, 628mb, 273.2K, 5697 0.02sec 26.41ns): 1.27115e-25 < xs < 2.94125e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section H2O 149 lines, 554mb, 267.2K, 6337 0.02sec 25.63ns): 9.90593e-26 < xs < 2.95906e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section H2O 149 lines, 398.1mb, 251.4K, 8449 0.03sec 20.46ns): 5.01327e-26 < xs < 2.69644e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section H2O 149 lines, 281mb, 235.3K, 11393 0.03sec 15.27ns): 2.36794e-26 < xs < 2.27903e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section H2O 149 lines, 193.4mb, 219.1K, 15617 0.03sec 11.34ns): 1.03491e-26 < xs < 1.75959e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section H2O 149 lines, 130mb, 215.7K, 22337 0.03sec 7.95ns): 6.29142e-27 < xs < 2.13348e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section H2O 149 lines, 59.5mb, 219.2K, 41729 0.03sec 4.87ns): 3.19831e-27 < xs < 3.95134e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section H2O 149 lines, 27.7mb, 225.1K, 62273 0.03sec 3.14ns): 1.76658e-27 < xs < 6.39428e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section H2O 149 lines, 13.2mb, 233.7K, 75713 0.03sec 2.68ns): 1.06636e-27 < xs < 9.62006e-21 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section H2O 149 lines, 6.52mb, 245.2K, 81537 0.03sec 2.49ns): 7.06861e-28 < xs < 1.40301e-20 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section H2O 149 lines, 3.33mb, 257.5K, 83201 0.03sec 2.48ns): 4.82193e-28 < xs < 1.94888e-20 in Interval [2100.02,2149.99]


lbl2xs: H2O 18.02amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section H2O 149 lines, 0.951mb, 275.7K, 82945 0.03sec 2.49ns): 2.02625e-28 < xs < 2.91311e-20 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section CH4 8 lines, 1013mb, 294.2K, 3713 0.00sec 37.81ns): 9.01524e-29 < xs < 5.22111e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section CH4 8 lines, 902mb, 289.7K, 4097 0.00sec 32.89ns): 7.9474e-29 < xs < 5.73113e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section CH4 8 lines, 802mb, 285.2K, 4545 0.00sec 46.09ns): 6.9893e-29 < xs < 6.25826e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section CH4 8 lines, 710mb, 279.2K, 5057 0.00sec 41.38ns): 6.08851e-29 < xs < 6.86809e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section CH4 8 lines, 628mb, 273.2K, 5633 0.00sec 36.84ns): 5.28916e-29 < xs < 7.48397e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section CH4 8 lines, 554mb, 267.2K, 6273 0.00sec 33.14ns): 4.57373e-29 < xs < 8.14108e-25 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section CH4 8 lines, 398.1mb, 251.4K, 8321 0.00sec 25.64ns): 3.08593e-29 < xs < 1.01343e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section CH4 8 lines, 281mb, 235.3K, 11137 0.00sec 19.38ns): 2.00359e-29 < xs < 1.2581e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section CH4 8 lines, 193.4mb, 219.1K, 15105 0.00sec 14.88ns): 1.23799e-29 < xs < 1.5575e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section CH4 8 lines, 130mb, 215.7K, 21569 0.00sec 10.85ns): 8.10969e-30 < xs < 2.182e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section CH4 8 lines, 59.5mb, 219.2K, 40321 0.00sec 6.86ns): 3.81304e-30 < xs < 4.39838e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section CH4 8 lines, 27.7mb, 225.1K, 59649 0.00sec 5.09ns): 1.85174e-30 < xs < 7.50115e-24 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section CH4 8 lines, 13.2mb, 233.7K, 71937 0.00sec 4.60ns): 9.32416e-31 < xs < 1.05126e-23 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section CH4 8 lines, 6.52mb, 245.2K, 77249 0.00sec 5.06ns): 4.9057e-31 < xs < 1.27456e-23 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section CH4 8 lines, 3.33mb, 257.5K, 78721 0.00sec 4.92ns): 2.64979e-31 < xs < 1.42111e-23 in Interval [2100.02,2149.99]


lbl2xs: CH4 16.04amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section CH4 8 lines, 0.951mb, 275.7K, 78337 0.00sec 4.56ns): 8.07258e-32 < xs < 1.54975e-23 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section CO 66 lines, 1013mb, 294.2K, 3393 0.01sec 30.08ns): 4.46509e-22 < xs < 1.70278e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section CO 66 lines, 902mb, 289.7K, 3777 0.01sec 26.76ns): 4.06914e-22 < xs < 1.90962e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section CO 66 lines, 802mb, 285.2K, 4225 0.01sec 24.40ns): 3.70456e-22 < xs < 2.13014e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section CO 66 lines, 710mb, 279.2K, 4673 0.01sec 23.18ns): 3.38689e-22 < xs < 2.36888e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section CO 66 lines, 628mb, 273.2K, 5249 0.01sec 23.03ns): 3.09633e-22 < xs < 2.65816e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section CO 66 lines, 554mb, 267.2K, 5825 0.01sec 24.11ns): 2.82563e-22 < xs < 2.99928e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section CO 66 lines, 398.1mb, 251.4K, 7809 0.01sec 22.41ns): 2.22888e-22 < xs < 4.11656e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section CO 66 lines, 281mb, 235.3K, 10497 0.01sec 17.57ns): 1.74306e-22 < xs < 5.69968e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section CO 66 lines, 193.4mb, 219.1K, 14465 0.01sec 12.33ns): 1.34051e-22 < xs < 7.98057e-18 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section CO 66 lines, 130mb, 215.7K, 20929 0.01sec 8.62ns): 9.22919e-23 < xs < 1.17084e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section CO 66 lines, 59.5mb, 219.2K, 41665 0.01sec 4.60ns): 4.11987e-23 < xs < 2.38549e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section CO 66 lines, 27.7mb, 225.1K, 68161 0.01sec 2.99ns): 1.84051e-23 < xs < 4.16413e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section CO 66 lines, 13.2mb, 233.7K, 88449 0.01sec 2.47ns): 8.2756e-24 < xs < 5.87528e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section CO 66 lines, 6.52mb, 245.2K, 98625 0.01sec 2.26ns): 3.79526e-24 < xs < 6.91396e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section CO 66 lines, 3.33mb, 257.5K, 102273 0.01sec 2.21ns): 1.79779e-24 < xs < 7.29271e-17 in Interval [2100.02,2149.99]


lbl2xs: CO 28.01amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section CO 66 lines, 0.951mb, 275.7K, 103041 0.01sec 2.19ns): 4.62603e-25 < xs < 7.36288e-17 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 1.013e+06

temperature [K] 296.0 ---> 294.2

cross section O3 2037 lines, 1013mb, 294.2K, 3841 0.20sec 25.11ns): 5.1466e-23 < xs < 7.71705e-20 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 902000

temperature [K] 296.0 ---> 289.7

cross section O3 2037 lines, 902mb, 289.7K, 4225 0.20sec 22.86ns): 4.66247e-23 < xs < 8.23995e-20 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 802000

temperature [K] 296.0 ---> 285.2

cross section O3 2037 lines, 802mb, 285.2K, 4737 0.20sec 20.52ns): 4.21112e-23 < xs < 8.79562e-20 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 710000

temperature [K] 296.0 ---> 279.2

cross section O3 2037 lines, 710mb, 279.2K, 5249 0.20sec 18.63ns): 3.78781e-23 < xs < 9.40404e-20 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 628000

temperature [K] 296.0 ---> 273.2

cross section O3 2037 lines, 628mb, 273.2K, 5825 0.20sec 16.71ns): 3.39576e-23 < xs < 1.0051e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 554000

temperature [K] 296.0 ---> 267.2

cross section O3 2037 lines, 554mb, 267.2K, 6529 0.20sec 14.97ns): 3.03481e-23 < xs < 1.06935e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 398100

temperature [K] 296.0 ---> 251.45

cross section O3 2037 lines, 398.1mb, 251.4K, 8641 0.35sec 19.63ns): 2.23962e-23 < xs < 1.24598e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 281000

temperature [K] 296.0 ---> 235.3

cross section O3 2037 lines, 281mb, 235.3K, 11649 0.35sec 14.70ns): 1.61319e-23 < xs < 1.55451e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 193400

temperature [K] 296.0 ---> 219.05

cross section O3 2037 lines, 193.4mb, 219.1K, 15937 0.35sec 10.75ns): 1.12592e-23 < xs < 2.15396e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 130000

temperature [K] 296.0 ---> 215.7

cross section O3 2037 lines, 130mb, 215.7K, 23169 0.36sec 7.54ns): 7.59502e-24 < xs < 2.92476e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 59500

temperature [K] 296.0 ---> 219.2

cross section O3 2037 lines, 59.5mb, 219.2K, 47489 0.37sec 3.78ns): 3.47423e-24 < xs < 5.00737e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 27700

temperature [K] 296.0 ---> 225.1

cross section O3 2037 lines, 27.7mb, 225.1K, 81793 0.39sec 2.36ns): 1.61224e-24 < xs < 9.12019e-19 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 13200

temperature [K] 296.0 ---> 233.7

cross section O3 2037 lines, 13.2mb, 233.7K, 111041 0.40sec 1.78ns): 7.63634e-25 < xs < 1.34766e-18 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 6520

temperature [K] 296.0 ---> 245.2

cross section O3 2037 lines, 6.52mb, 245.2K, 126785 0.41sec 1.59ns): 3.73352e-25 < xs < 1.64179e-18 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 3330

temperature [K] 296.0 ---> 257.5

cross section O3 2037 lines, 3.33mb, 257.5K, 132993 0.42sec 1.54ns): 1.88214e-25 < xs < 1.75408e-18 in Interval [2100.02,2149.99]


lbl2xs: O3 48.00amu pressure [g/cm/s**2] 1.01325e+06 ---> 951

temperature [K] 296.0 ---> 275.7

cross section O3 2037 lines, 0.951mb, 275.7K, 135041 0.42sec 1.54ns): 5.25138e-26 < xs < 1.7456e-18 in Interval [2100.02,2149.99]



INFO --- xs2ac: atmosphere with 16 levels with 8 molecules H2O CO2 O3 N2O CO CH4 N2 O2 but 6 of 8 gases without cross sections!



INFO --- xs2ac: x limits identical for all xs 2100.017 2149.993


Absorption coefficient for 8 molecules, 16 levels, and 135041 wavenumber points:

0 0.00km <ac> = 1.08e-05 1.77e-07 < ac < 0.00138 @ 2136.139315

1 1.00km <ac> = 6.7e-06 1.08e-07 < ac < 0.000962 @ 2136.144496

2 2.00km <ac> = 4e-06 6.61e-08 < ac < 0.000617 @ 2136.145977

3 3.00km <ac> = 2.09e-06 3.78e-08 < ac < 0.000343 @ 2136.144867

4 4.00km <ac> = 1.15e-06 2.28e-08 < ac < 0.000187 @ 2136.147827

5 5.00km <ac> = 6.17e-07 1.18e-08 < ac < 9.88e-05 @ 2136.142276

6 7.50km <ac> = 2.17e-07 3.12e-09 < ac < 2.57e-05 @ 2136.144126

7 10.00km <ac> = 1.02e-07 7.65e-10 < ac < 4.93e-06 @ 2115.629689

8 12.50km <ac> = 6.58e-08 2.55e-10 < ac < 3.74e-06 @ 2119.680617

9 15.00km <ac> = 5.21e-08 1.03e-10 < ac < 2.22e-06 @ 2119.682097

10 20.00km <ac> = 6.86e-08 3.45e-11 < ac < 1.97e-06 @ 2123.086490

11 25.00km <ac> = 7.21e-08 1.32e-11 < ac < 3.87e-06 @ 2123.086860

12 30.00km <ac> = 4.76e-08 3.41e-12 < ac < 3.86e-06 @ 2123.086860

13 35.00km <ac> = 2.81e-08 8.49e-13 < ac < 2.8e-06 @ 2123.086860

14 40.00km <ac> = 1.15e-08 1.7e-13 < ac < 1.23e-06 @ 2123.086860

15 50.00km <ac> = 1.19e-09 5.49e-15 < ac < 1.22e-07 @ 2123.086860


Optical depth array with shape (135041, 15)

0.0142581 <= xy[:,00] <= 115.605 mean 0.87425 norm 1782.84

0.0086901 <= xy[:,01] <= 78.5223 mean 0.53469 norm 1123.83

0.00526501 <= xy[:,02] <= 47.6434 mean 0.30413 norm 649.556

0.00303478 <= xy[:,03] <= 26.4214 mean 0.16174 norm 342.335

0.00177531 <= xy[:,04] <= 14.1232 mean 0.088241 norm 178.259

0.00186857 <= xy[:,05] <= 15.3974 mean 0.10424 norm 183.743

0.000489132 <= xy[:,06] <= 3.80657 mean 0.039857 norm 53.8481

0.000127527 <= xy[:,07] <= 1.0789 mean 0.020932 norm 27.1892

4.47596e-05 <= xy[:,08] <= 0.742381 mean 0.014739 norm 16.1836

3.45108e-05 <= xy[:,09] <= 0.887937 mean 0.030161 norm 24.685

1.22676e-05 <= xy[:,10] <= 1.46058 mean 0.035166 norm 35.4468

4.14934e-06 <= xy[:,11] <= 1.93304 mean 0.029928 norm 37.5762

1.0656e-06 <= xy[:,12] <= 1.66462 mean 0.018932 norm 28.1874

2.54791e-07 <= xy[:,13] <= 1.00711 mean 0.0099106 norm 16.0767

8.77812e-08 <= xy[:,14] <= 0.675009 mean 0.0063629 norm 10.6116



In [18]: odPlot(dodList) # quicklook plot



In [19]: vGrid, radUp = dod2ri(dodList) # downwelling radiation seen by an uplooking observer

    ...: radDown = dod2ri(vGrid, dod, mls['T'], 180.) # upwelling radiation seen by a nadir view

    ...:


In [20]: vGrid, radNadir = dod2ri(dodList, 180., mls['T'][0]) # forgot surface emission, assume T_surface=T_BoA=mls['T'][0]=mls[0]['T']

Adding Planck BB background for 294.2K ===> 3.21273 <= B <= 3.82278


In [21]: plot(vGrid,radUp,'r',label='upLooking');

    ...: plot(vGrid,radDown,'g',label='down without surface');

    ...: plot(vGrid,radNadir,'b',label='down with surface');

    ...: legend();

    ...: xlabel('wavenumber $[cm^{-1}]$');

    ...: ylabel('radiance $[erg/s/(cm^2.sr.cm^{-1})]$');

    ...: title('midlatitude summer')

    ...:

Out[21]: <matplotlib.text.Text at 0xbe9d8d0>


In [22]: