Skip to content
This repository has been archived by the owner on May 13, 2020. It is now read-only.

overlayFov groundscatter model bugfix #269

Merged
merged 1 commit into from
Aug 17, 2016

Conversation

aburrell
Copy link
Contributor

@aburrell aburrell commented Aug 17, 2016

Fixed bug that caused fields-of-view to not close when the groundscatter model is used.

To test:

import datetime as dt
import davitpy
import matplotlib.pyplot as plt
plt.ion()
f = plt.figure()
ax = f.add_subplot(1,1,1)
import mpl_toolkits.basemap as basemap
rad = 'han'
stime = dt.datetime(2006, 10, 1)
hard = davitpy.pydarn.radar.site(code=rad, dt=stime)
m = davitpy.utils.mapObj(ax=ax, projection="stere", lon_0=hard.geolon, lat_0=hard.geolat, llcrnrlon=0, llcrnrlat=25, urcrnrlat=70, urcrnrlon=135, resolution="l")
m.drawcoastlines(linewidth=0.5, color="0.6")
fov_test = davitpy.pydarn.radar.radFov.fov(site=hard, ngates=75, altitude=350, coords="geo", model="GS", fov_dir="front")
davitpy.pydarn.plotting.overlayFov(m, rad, fovObj=fov_test)

fov_test = davitpy.pydarn.radar.radFov.fov(site=hard, ngates=75, altitude=350, coords="geo", model="IS", fov_dir="back")
davitpy.pydarn.plotting.overlayFov(m, rad, fovObj=fov_test, fovColor="b")

gs_fov

print fov_test

latCenter: (16, 75)
lonCenter: (16, 75)
latFull: (17, 76)
lonFull: (17, 76)
slantRCenter: (16, 75)
slantRFull: (17, 76)
beams: (16,)
gates: (75,)
coords: geo
field of view: back
model: IS

Add field-of-view direction and model to fov class
@Shirling-VT
Copy link
Contributor

Just tested it and get exactly the same results and plot:

print fov_test

latCenter: (16, 75)
lonCenter: (16, 75)
latFull: (17, 76)
lonFull: (17, 76)
slantRCenter: (16, 75)
slantRFull: (17, 76)
beams: (16,)
gates: (75,)
coords: geo
field of view: back
model: IS

On Wed, Aug 17, 2016 at 10:37 AM, Angeline Burrell <notifications@github.com

wrote:

Fixed bug that caused fields-of-view to not close when the groundscatter
model is used.

To test:

import datetime as dt
import davitpy
import matplotlib.pyplot as plt
plt.ion()
f = plt.figure()
ax = f.add_subplot(1,1,1)
import mpl_toolkits.basemap as basemap
rad = 'han'
stime = dt.datetime(2006, 10, 1)
hard = davitpy.pydarn.radar.site(code=rad, dt=stime)
m = davitpy.utils.mapObj(ax=ax, projection="stere", lon_0=hard.geolon,
lat_0=hard.geolat, llcrnrlon=0, llcrnrlat=25, urcrnrlat=70, urcrnrlon=135,
resolution="l")
m.drawcoastlines(linewidth=0.5, color="0.6")
fov_test = davitpy.pydarn.radar.radFov.fov(site=hard, ngates=75,
altitude=350, coords="geo", model="GS", fov_dir="front")
davitpy.pydarn.plotting.overlayFov(m, rad, fovObj=fov_test)

fov_test = davitpy.pydarn.radar.radFov.fov(site=hard, ngates=75,
altitude=350, coords="geo", model="IS", fov_dir="back")
davitpy.pydarn.plotting.overlayFov(m, rad, fovObj=fov_test, fovColor="b")

[image: gs_fov]

https://cloud.githubusercontent.com/assets/7045886/17740425/761fde40-6490-11e6-820d-d60f6e322bcb.png

You can view, comment on, or merge this pull request online at:

#269
Commit Summary

  • bugfix

File Changes

Patch Links:


You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub
#269, or mute the thread
https://github.com/notifications/unsubscribe-auth/AIEZK9QZTqsQGdFmHIYXtLmWQw-VVhSYks5qgxzEgaJpZM4JmhU4
.

Best wishes,
Xueling

@asreimer
Copy link
Contributor

I'll give this a test in a few hours and report back.

@asreimer
Copy link
Contributor

Testing this pull request, I see that #258 is fixed.

Further testing reveals that this pull request doesn't break anything. Nice work @aburrell!

Merging.

@asreimer asreimer merged commit e376806 into vtsuperdarn:develop Aug 17, 2016
@aburrell aburrell deleted the overlayFov_gs_bugfix branch September 29, 2017 20:20
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants