diff --git a/Gallery/MapProjections/NCL_sat_1.py b/Gallery/MapProjections/NCL_sat_1.py index 0bd2ba001..167ecd01a 100644 --- a/Gallery/MapProjections/NCL_sat_1.py +++ b/Gallery/MapProjections/NCL_sat_1.py @@ -48,8 +48,11 @@ # Set figure size fig = plt.figure(figsize=(8, 8)) -# Set global axes with an orthographic projection -proj = ccrs.Orthographic(central_longitude=270, central_latitude=45) +# Set global axes with a nearside perspective projection (equivalent to NCL's +# satellite projection) +proj = ccrs.NearsidePerspective(central_longitude=270.0, + central_latitude=45.0, + satellite_height=12742000) ax = plt.axes(projection=proj) ax.set_global() diff --git a/Gallery/MapProjections/NCL_sat_2.py b/Gallery/MapProjections/NCL_sat_2.py index 2dfb7914f..3c3fba5a1 100644 --- a/Gallery/MapProjections/NCL_sat_2.py +++ b/Gallery/MapProjections/NCL_sat_2.py @@ -50,8 +50,11 @@ # Set figure size fig = plt.figure(figsize=(8, 8)) -# Set global axes with an orthographic projection -proj = ccrs.Orthographic(central_longitude=270, central_latitude=45) +# Set global axes with a nearside perspective projection (equivalent to NCL's +# satellite projection) +proj = ccrs.NearsidePerspective(central_longitude=270.0, + central_latitude=45.0, + satellite_height=12742000) ax = plt.axes(projection=proj) ax.set_global() diff --git a/Gallery/MapProjections/NCL_sat_3.py b/Gallery/MapProjections/NCL_sat_3.py index f8e67376b..b185f007d 100644 --- a/Gallery/MapProjections/NCL_sat_3.py +++ b/Gallery/MapProjections/NCL_sat_3.py @@ -41,10 +41,10 @@ plt.figure(figsize=(8, 8)) -# Create an axis with an orthographic projection -ax = plt.axes(projection=ccrs.Orthographic(central_longitude=-35, - central_latitude=60), - anchor='C') +# Create an axis with an orthographic projection (equivalent to NCL's satellite +# projection where mpSatelliteDistF <= 1.0) +ax = plt.axes( + projection=ccrs.Orthographic(central_longitude=-35, central_latitude=60)) # Set extent of map ax.set_extent((-80, -10, 30, 80), crs=ccrs.PlateCarree())