scivision / lowtran / 232

pending completion

Build # 232

Build Type

push

travis-ci

Committed by

scivision

Commit Message

<a href="https://github.com/scivision/lowtran/commit/<a class=hub.com/scivision/lowtran/commit/<a class="double-link" href="https://git"><a class=hub.com/scivision/lowtran/commit/<a class="double-link" href="https://git"><a class=hub.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a">92c683e91<a href="https://github.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a">&lt;a href=&quot;https://github.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a">&lt;a href=&quot;https://github.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a&quot;&gt;&amp;quot;&amp;gt;Squashed commit of the following:

commit &amp;lt;/a&amp;gt;&amp;lt;a class=&amp;quot;double-link&amp;quot; href=&amp;quot;https://github.com/scivision/lowtran/commit/&amp;lt;a class=&amp;quot;double-link&amp;quot; href=&amp;quot;https://github.com/scivision/lowtran/commit/30fcb0362d06ceedf46abeee729b9398d5d186b4&quot;&gt;30fcb0362&lt;/a&gt;&quot;&gt;30fcb0362&lt;/a&gt;&lt;a href=&quot;https://github.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a&quot;&gt;
Author: Michael Hirsch, Ph.D &amp;lt;scivision@users.noreply.github.com&amp;gt;
Date:   Mon Aug 13 18:50:14 2018 -0400

    allow user to specify wavelength step (short, long, step)

commit &lt;/a&gt;&lt;a class=&quot;double-link&quot; href=&quot;https://github.com/scivision/lowtran/commit/&lt;a class=&quot;double-link&quot; href=&quot;https://github.com/scivision/lowtran/commit/4f0b26421207acaa2d09fcd61d357a765a124532">4f0b26421</a>">4f0b26421</a><a href="https://github.com/scivision/lowtran/commit/92c683e91528124b1e7e387a73a0d4d4732c628a">
Author: Michael Hirsch, Ph.D &lt;scivision@users.noreply.github.com&gt;
Date:   Mon Aug 13 18:07:30 2018 -0400

    cleanup

commit </a><a class="double-link" href="https://github.com/scivision/lowtran/commit/<a class="double-link" href="https://github.com/scivision/lowtran/commit/c45c10e76206efade9c0e03543a73aa06601e113">c45c10e76">c45c10e76
Author: Michael Hirsch, Ph.D <scivision@users.noreply.github.com>
Date:   Mon Aug 13 18:03:57 2018 -0400

    modularize

Run Details

115 of 115 new or added lines in 6 files covered. (100.0%)

128 of 162 relevant lines covered (79.01%)

0.79 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

57.78

/lowtran/scenarios.py

from pathlib import Path
from pandas import read_csv
from dateutil.parser import parse
import xarray
from typing import Dict, Any
import numpy as np

from .base import loopangle, loopuserdef, golowtran


def scatter(c1: Dict[str, Any]) -> xarray.Dataset:
    # %% low-level Lowtran configuration for this scenario, don't change
    c1.update({
        'itype':  3,  # 3: observer to space
        'iemsct': 2,  # 2: radiance model
    })
# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
    return loopangle(c1)


def irradiance(c1: Dict[str, Any]) -> xarray.Dataset:
    c1.update({
        'itype': 3,   # 3: observer to space
        'iemsct': 3,  # 3: solar irradiance
    })

    return loopangle(c1)


def radiance(c1: Dict[str, Any]) -> xarray.Dataset:

    c1.update({
        'itype':  3,  # 3: observer to space
        'iemsct': 1,  # 1: thermal radiance model  2: radiance model
    })
# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
    return loopangle(c1)


def transmittance(c1: Dict[str, Any]) -> xarray.Dataset:

    c1.update({
        'itype': 3,   # 3: observer to space
        'iemsct': 0,  # 0: transmittance
    })

    return loopangle(c1)


def horizrad(infn: Path, outfn: Path, c1: dict) -> xarray.Dataset:
    """
    read CSV, simulate, write, plot
    """
    if infn is not None:
        infn = Path(infn).expanduser()

        if infn.suffix == '.h5':
            TR = xarray.open_dataset(infn)
            return TR

    c1.update({'model': 0,  # 0: user meterological data
               'itype': 1,  # 1: horizontal path
               'iemsct': 1,  # 1: radiance model
               'im': 1,    # 1: for horizontal path (see Lowtran manual p.42)
               'ird1': 1,  # 1: use card 2C2)
               })
# %% read csv file
    if not infn:  # demo mode
        c1['p'] = [949., 959.]
        c1['t'] = [283.8, 285.]
        c1['wmol'] = [[93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],
                      [93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]]
        c1['time'] = [parse('2017-04-05T12'),
                      parse('2017-04-05T18')]
    else:  # read csv, normal case
        PTdata = read_csv(infn)
        c1['p'] = PTdata['p']
        c1['t'] = PTdata['Ta']
        c1['wmol'] = np.zeros((PTdata.shape[0], 12))
        c1['wmol'][:, 0] = PTdata['RH']
        c1['time'] = [parse(t) for t in PTdata['time']]
# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
    TR = loopuserdef(c1)

    return TR


def horiztrans(c1: Dict[str, Any]) -> xarray.Dataset:

    c1.update({'model': 5,  # 5: subartic winter
               'itype': 1,  # 1: horizontal path
               'iemsct': 0,  # 0: transmittance model
               'im': 0,  # 1: for user-defined atmosphere on horizontal path (see Lowtran manual p.42)
               'ird1': 0,  # 1: use card 2C2
               })

    return golowtran(c1)


def userhoriztrans(c1: Dict[str, Any]) -> xarray.Dataset:

    c1.update({'model': 0,  # 0: user meterological data
               'itype': 1,  # 1: horizontal path
               'iemsct': 0,  # 0: transmittance model
               'im': 1,  # 1: for user-defined atmosphere on horizontal path (see Lowtran manual p.42)
               'ird1': 1,  # 1: use card 2C2
               })

    return golowtran(c1)

1	from pathlib import Path	1✔
2	from pandas import read_csv	1✔
3	from dateutil.parser import parse	1✔
4	import xarray	1✔
5	from typing import Dict, Any	1✔
6	import numpy as np	1✔
7
8	from .base import loopangle, loopuserdef, golowtran	1✔
9
10
11	def scatter(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
12	# %% low-level Lowtran configuration for this scenario, don't change
13	c1.update({	1✔
14	'itype': 3, # 3: observer to space
15	'iemsct': 2, # 2: radiance model
16	})
17	# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
18	return loopangle(c1)	1✔
19
20
21	def irradiance(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
22	c1.update({	1✔
23	'itype': 3, # 3: observer to space
24	'iemsct': 3, # 3: solar irradiance
25	})
26
27	return loopangle(c1)	1✔
28
29
30	def radiance(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
31
32	c1.update({	1✔
33	'itype': 3, # 3: observer to space
34	'iemsct': 1, # 1: thermal radiance model 2: radiance model
35	})
36	# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
37	return loopangle(c1)	1✔
38
39
40	def transmittance(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
41
42	c1.update({	1✔
43	'itype': 3, # 3: observer to space
44	'iemsct': 0, # 0: transmittance
45	})
46
47	return loopangle(c1)	1✔
48
49
50	def horizrad(infn: Path, outfn: Path, c1: dict) -> xarray.Dataset:	1✔
51	"""
52	read CSV, simulate, write, plot
53	"""
54	if infn is not None:	×
55	infn = Path(infn).expanduser()	×
56
57	if infn.suffix == '.h5':	×
58	TR = xarray.open_dataset(infn)	×
59	return TR	×
60
61	c1.update({'model': 0, # 0: user meterological data	×
62	'itype': 1, # 1: horizontal path
63	'iemsct': 1, # 1: radiance model
64	'im': 1, # 1: for horizontal path (see Lowtran manual p.42)
65	'ird1': 1, # 1: use card 2C2)
66	})
67	# %% read csv file
68	if not infn: # demo mode	×
69	c1['p'] = [949., 959.]	×
70	c1['t'] = [283.8, 285.]	×
71	c1['wmol'] = [[93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.],	×
72	[93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]]
73	c1['time'] = [parse('2017-04-05T12'),	×
74	parse('2017-04-05T18')]
75	else: # read csv, normal case
76	PTdata = read_csv(infn)	×
77	c1['p'] = PTdata['p']	×
78	c1['t'] = PTdata['Ta']	×
79	c1['wmol'] = np.zeros((PTdata.shape[0], 12))	×
80	c1['wmol'][:, 0] = PTdata['RH']	×
81	c1['time'] = [parse(t) for t in PTdata['time']]	×
82	# %% TR is 3-D array with axes: time, wavelength, and [transmission,radiance]
83	TR = loopuserdef(c1)	×
84
85	return TR	×
86
87
88	def horiztrans(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
89
90	c1.update({'model': 5, # 5: subartic winter	1✔
91	'itype': 1, # 1: horizontal path
92	'iemsct': 0, # 0: transmittance model
93	'im': 0, # 1: for user-defined atmosphere on horizontal path (see Lowtran manual p.42)
94	'ird1': 0, # 1: use card 2C2
95	})
96
97	return golowtran(c1)	1✔
98
99
100	def userhoriztrans(c1: Dict[str, Any]) -> xarray.Dataset:	1✔
101
102	c1.update({'model': 0, # 0: user meterological data	1✔
103	'itype': 1, # 1: horizontal path
104	'iemsct': 0, # 0: transmittance model
105	'im': 1, # 1: for user-defined atmosphere on horizontal path (see Lowtran manual p.42)
106	'ird1': 1, # 1: use card 2C2
107	})
108
109	return golowtran(c1)	1✔

scivision / lowtran / 232

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous