Fix up reprojection example

examples/plot_coadd.py

@@ -24,8 +24,8 @@
 ###############################################################################
 # Download STEREO maps
 maps_in_stereo = []
-for i in range(10):
-    d = datetime.now() - timedelta(days=i) - timedelta(days=2)
+for i in range(2):
+    d = datetime.now() - timedelta(days=i) - timedelta(days=1)
     smap = stereo_start_of_day_map(d)
     # Normalise to AIA
     data = (smap.data - 720.0) * 110 / 170
@@ -36,8 +36,8 @@
 ###############################################################################
 # Download AIA maps
 maps_in_aia = []
-for i in range(17):
-    d = datetime.now() - timedelta(days=i) - timedelta(days=2)
+for i in range(2):
+    d = datetime.now() - timedelta(days=i) - timedelta(days=1)
     maps_in_aia.append(aia_start_of_day_map(d, 193 * u.Angstrom))
 
 norm = maps_in_aia[0].plot_settings['norm']

solarsynoptic/reprojection/reprojection.py

@@ -6,6 +6,7 @@
 from astropy.coordinates import SkyCoord
 from astropy.wcs import WCS
 from reproject import reproject_interp
+from sunpy.sun import constants
 from sunpy import log
 from sunpy.map import Map, make_fitswcs_header
 
@@ -34,6 +35,14 @@ def reproject_carrington(smap, shape_out, latitude_projection='CAR',
     carrington_map : sunpy.map.Genericmap
         Reprojected map.
     """
+    rsun_m = constants.radius.to_value(u.m)
+    if 'rsun_ref' not in smap.meta:
+        log.info('Setting rsun_ref in metadata')
+        smap.meta['rsun_ref'] = rsun_m
+    elif smap.meta['rsun_ref'] != rsun_m:
+        log.info('Overwriting rsun_ref with standard photospheric radius')
+        smap.meta['rsun_ref'] = rsun_m
+
     map_out = None
     if cache:
         from solarsynoptic.reprojection.database import DATABASE