!-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/
let mmmm = if mask lt 0 then 1
let xpts = x[gx=mmmm]+0*y[gy=mmmm]
let ypts = y[gy=mmmm]+0*y[gx=mmmm]
plot/vs/ov/SYMBOL=dot/nolab/nokey/size=0.15/opacity=90 mmmm*xpts*mp_mask, mmmm*ypts*mp_mask, x_page, y_page
!-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/
2. Regarding viewports, the following script might be helpful.
!-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/-/
set mem/size=1000
set window/quality=high
let mean =SST[l=@ave]
set view view_2
go mp_stereographic_south 180 0; GO mp_grid mean ; SET GRID mean
fill/NOAXIS/lev=(-inf)(-25,25,2.5)(inf)/palette=blue_orange/nolab/set mean*mp_mask, x_page, y_page
[ferret_users] stippling on the polar stereographic plots