search

timduly4 / pyglow

Upper atmosphere climatological models in Python
MIT License
104 stars 57 forks source link

Not found point.TEC property #153

Open leeliangchao opened 1 year ago

leeliangchao commented 1 year ago

When I want to use IRI 2016 to generate TEC, I found that there is no option for TEC in point.[property], how to use IRI to generate TEC, thank you

# Cyclic calculation of TEC
for i, lat in enumerate(latitudes):
    for j, lon in enumerate(longitudes):
        point = Point(date, lat, lon, 0)
        point.run_iri(version=2016)
        tec = point.??????
        tec_data[i, j] = tec
bharding512 commented 1 year ago

pyglow does not currently calculate TEC. You'd have to calculate your own integral by instantiating a number of points along the line of sight and running IRI at each point.

It could be useful to add TEC to pyglow, but I'm not sure how to reconcile it with the fact that a "point" (1-D) philosophically cannot have a TEC (which is a LoS integral). I'm open to ideas though.