gaps in map projections: SHADE, mp_fland.jnl, mp_land.jnl

[karlmsmith]

Reported by @AndrewWittenberg on 2 Apr 2015 19:56 UTC When working with map projections, SHADE leaves a gap at the longitudinal modulo cut. (FILL does not.) Also since mp_fland.jnl does a SHADE, it leaves a gap too. And there's a gap at the same longitude in the mp_land.jnl outline.

    NOAA/PMEL TMAP
    FERRET v6.933 (beta)
    Linux 2.6.32-504.3.3.el6.x86_64 64-bit - 01/27/15
     2-Apr-15 15:43     

yes? use coads_climatology; can reg; set region/l=1
yes? go mp_orthographic 0 90; set grid sst; go mp_aspect
yes? shade/noax/nokey sst*mp_mask, x_page, y_page
yes? go mp_fland; go mp_land 1

I'm attaching the resulting plot.

Is there a way to easily specify the (N+1)x(N+1) bounds for x_page and y_page, needed to display the actual grid cells as in a regular 2d SHADE? This is helpful when the grid is irregular, such that the bounds are not midway between the cell centers, as Ferret assumes when x_page and y_page are NxN.

Migrated-From: http://dunkel.pmel.noaa.gov/trac/ferret/ticket/2257

[karlmsmith]

Comment by @AnsleyManke on 2 Apr 2015 20:49 UTC Setting a region that wraps more than 360 degrees seems to be a workaround. It takes care of the land overlay as well.

yes? use coads_climatology; can reg; set region/l=1/x=-1:361  ! or 19:381
yes? go mp_orthographic 0 90; set grid sst; go mp_aspect
yes? shade/noax/nokey sst*mp_mask, x_page, y_page
yes? go mp_fland; go mp_land 1

[karlmsmith]

Attachment from @AndrewWittenberg on 2 Apr 2015 19:57 UTC The resulting plot.

[karlmsmith]

Attachment from @AnsleyManke on 2 Apr 2015 20:49 UTC the plot with region /x=-1:361