Пример #1
0
def main():
    p = ArgumentParser(description='Lowtran 7 interface')
    p.add_argument('-z',
                   '--obsalt',
                   help='altitude of bother observers on horizontal path [km]',
                   type=float,
                   default=0.3)
    p.add_argument('-r',
                   '--range_km',
                   help='range between observers on horizontal path [km]',
                   type=float,
                   default=1.0)
    p.add_argument(
        '-a',
        '--zenang',
        help='zenith angle [deg]  can be single value or list of values',
        type=float,
        default=0.)
    p.add_argument('-s',
                   '--short',
                   help='shortest wavelength nm ',
                   type=float,
                   default=200)
    p.add_argument('-l',
                   '--long',
                   help='longest wavelength nm ',
                   type=float,
                   default=30000)
    p.add_argument('-step',
                   help='wavelength step size cm^-1',
                   type=float,
                   default=20)
    P = p.parse_args()

    c1 = {
        'zmdl': P.obsalt,
        'h1': P.obsalt,
        'range_km': P.range_km,
        'wlshort': P.short,
        'wllong': P.long,
        'wlstep': P.step,
    }

    atmos = {
        'p': 949.,
        't': 283.8,
        'wmol': [93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]
    }

    c1.update(atmos)

    TR = lowtran.userhoriztrans(c1).squeeze()

    plothoriz(TR, c1)

    show()
Пример #2
0
def main():
    p = ArgumentParser(description='Lowtran 7 interface')
    p.add_argument('-z',
                   '--obsalt',
                   help='altitude of bother observers on horizontal path [km]',
                   type=float,
                   default=0.3)
    p.add_argument('-r',
                   '--range_km',
                   help='range between observers on horizontal path [km]',
                   type=float,
                   default=1.0)
    p.add_argument(
        '-a',
        '--zenang',
        help='zenith angle [deg]  can be single value or list of values',
        type=float,
        default=0.)
    p.add_argument('-w',
                   '--wavelen',
                   help='wavelength range nm (start,stop)',
                   type=float,
                   nargs=2,
                   default=(200, 30000))
    P = p.parse_args()

    c1 = {
        'zmdl': P.obsalt,
        'h1': P.obsalt,
        'range_km': P.range_km,
        'wlnmlim': P.wavelen,
    }

    atmos = {
        'p': 949.,
        't': 283.8,
        'wmol': [93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]
    }

    c1.update(atmos)

    TR = lowtran.userhoriztrans(c1).squeeze()

    plothoriz(TR, c1)

    show()
Пример #3
0
from lowtran.plots import plothoriz

if __name__=='__main__':

    from argparse import ArgumentParser
    p = ArgumentParser(description='Lowtran 7 interface')
    p.add_argument('-z','--obsalt',help='altitude of observer [km]',type=float,default=0.05)
    p.add_argument('-a','--zenang',help='zenith angle [deg]  can be single value or list of values',type=float,default=0.)
    p.add_argument('-w','--wavelen',help='wavelength range nm (start,stop)',type=float,nargs=2,default=(200,30000))
    p=p.parse_args()

    c1={'model':0, # 0: user meterological data
        'itype':1, # 1: horizontal path
        'iemsct':0, # 0: transmittance model
        'im': 1, # 1: for horizontal path (see Lowtran manual p.42)
        'ird1': 1, # 1: use card 2C2
        'range_km':p.obsalt,
        'zmdl':p.obsalt,
        'h1': p.obsalt,
        'wlnmlim': p.wavelen,
        'p':949.,
        't':283.8,
        'wmol':[93.96,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.]
        }

    TR = lowtran.golowtran(c1)

    plothoriz(TR, c1, False)

    show()
from lowtran.plots import plothoriz

if __name__=='__main__':

    from argparse import ArgumentParser
    p = ArgumentParser(description='Lowtran 7 interface')
    p.add_argument('-z','--obsalt',help='altitude of observer [km]',type=float,default=0.05)
    p.add_argument('-a','--zenang',help='zenith angle [deg]  can be single value or list of values',type=float,default=0.)
    p.add_argument('-w','--wavelen',help='wavelength range nm (start,stop)',type=float,nargs=2,default=(200,30000))
    p=p.parse_args()

    c1={'model':0, # 0: user meterological data
        'itype':1, # 1: horizontal path
        'iemsct':0, # 0: transmittance model
        'im': 1, # 1: for horizontal path (see Lowtran manual p.42)
        'ird1': 1, # 1: use card 2C2
        'range_km':p.obsalt,
        'zmdl':p.obsalt,
        'p':949.,
        't':283.8,
        'wmol':[93.96,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.]
        }

    TR = lowtran.golowtran(p.obsalt,p.zenang,p.wavelen,c1)

 
    plothoriz(TR,p.zenang,c1,False)


    show()
        type=float,
        default=0.)
    p.add_argument('-w',
                   '--wavelen',
                   help='wavelength range nm (start,stop)',
                   type=float,
                   nargs=2,
                   default=(200, 30000))
    p = p.parse_args()

    c1 = {
        'zmdl': p.obsalt,
        'h1': p.obsalt,
        'range_km': p.range_km,
        'wlnmlim': p.wavelen,
    }

    atmos = {
        'p': 949.,
        't': 283.8,
        'wmol': [93.96, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.]
    }

    c1.update(atmos)

    TR = lowtran.userhoriztrans(c1).squeeze()

    plothoriz(TR, c1)

    show()