# A weird issue related to subplot

**URL:** <https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460>\
**Category:** GMT.jl Q&A\
**Created:** [December 21, 2021, 4:38pm UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460 "2021-12-21T16:38:23Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![leon6j](https://avatars.discourse-cdn.com/v4/letter/l/a6a055/32.png) [@leon6j](https://forum.generic-mapping-tools.org/u/leon6j)\
**Post date:** [December 21, 2021, 4:38pm UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/1 "2021-12-21T16:38:23Z")

</div>

I’m trying to plot 3 maps as 3 subplots in one figure. Below is my code. Why doesn’t my `region` setting being honored when my land\_clip and my variable contour are plotted using `grdimage`? BTW, each of the 3 maps plotted very well as separate ones.

```
subplot(grid=(3,1), region=(0, 100, 0, 80), panels_size=8, margins="5p", autolabel=true, col_axes=(bott=true,), row_axes=(left=true,), axes="wstr", savefig=F1);

    subplot(:set); 
    coast(region=(180,320,15,80));
    grdimage! ...
    ... 

    subplot(:set); 
    coast(region=(180,320,15,80));
    grdimage! ...
    ... 

    subplot(:set); 
    coast(region=(180,320,15,80));
    grdimage! ...
    ...

```

Here is the plot:  
[test2|196x500](https://forum.generic-mapping-tools.org/uploads/short-url/skvWnnFImoYHKCHHVsgclbDpXJZ.jpeg)

Many thanks!!

---

<div class="post-metadata">

**Author:** ![Joaquim](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/joaquim/32/16_2.png) [@Joaquim](https://forum.generic-mapping-tools.org/u/Joaquim)\
**Post date:** [December 21, 2021, 11:12pm UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/2 "2021-12-21T23:12:49Z")

</div>

`grdimage` after a `coast` obviously hides the coastlines. But why don’t you post something we can try?

---

<div class="post-metadata">

**Author:** ![leon6j](https://avatars.discourse-cdn.com/v4/letter/l/a6a055/32.png) [@leon6j](https://forum.generic-mapping-tools.org/u/leon6j)\
**Post date:** [December 21, 2021, 11:48pm UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/3 "2021-12-21T23:48:59Z")

</div>

Many thanks, @Joaquim.

Below is a shorter version of my code:

```
using GMT;

Land_clip = grdclip("@earth_relief_05m", low=[0 NaN]); # Remove ocean areas

subplot(grid=(3,1), panels_size=8, margins="5p",
autolabel=true, col_axes=(bott=true,), row_axes=(left=true,), axes="wstr",
savefig="Variable.jpg");

subplot(:set);
coast!(region=(180,320,15,80));
grdimage!(Land_clip, color=:geo, nan_alpha=true, shade=false,
xaxis=(annot=15, grid=30, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
yaxis=(annot=15, grid=15, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
par=(MAP_GRID_PEN_PRIMARY="0.005p",), coast=true);

subplot(:set);
coast!(region=(180,320,15,80));
grdimage!(Land_clip, color=:geo, nan_alpha=true, shade=false,
xaxis=(annot=15, grid=30, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
yaxis=(annot=15, grid=15, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
par=(MAP_GRID_PEN_PRIMARY="0.005p",), coast=true);

subplot(:set);
coast!(region=(180,320,15,80));
grdimage!(Land_clip, color=:geo, nan_alpha=true, shade=false,
xaxis=(annot=15, grid=30, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
yaxis=(annot=15, grid=15, FONT_ANNOT_PRIMARY="12p,Helvetica,black",),
par=(MAP_GRID_PEN_PRIMARY="0.005p",), coast=true);

subplot(:show)
```

---

<div class="post-metadata">

**Author:** ![Joaquim](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/joaquim/32/16_2.png) [@Joaquim](https://forum.generic-mapping-tools.org/u/Joaquim)\
**Post date:** [December 22, 2021, 1:07am UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/4 "2021-12-22T01:07:29Z")

</div>

> [@leon6j](#):
>
> ```auto
> coast!(region=(180,320,15,80));
> grdimage!(Land_clip, color=:geo, nan_alpha=true, shade=false,
> 
> ```

This is wrong. I warned a couple of times. `grdimage` **overlaps** the coastlines.

> `grdimage` after a `coast` obviously hides the coastlines.

---

<div class="post-metadata">

**Author:** ![Joaquim](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/joaquim/32/16_2.png) [@Joaquim](https://forum.generic-mapping-tools.org/u/Joaquim)\
**Post date:** [December 22, 2021, 1:14am UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/5 "2021-12-22T01:14:50Z")

</div>

> [@leon6j](#):
>
> `coast!(region=(180,320,15,80));`

> [@leon6j](#):
>
> `par=(MAP_GRID_PEN_PRIMARY="0.005p",), coast=true);`

This is twice the coastlines

---

<div class="post-metadata">

**Author:** ![Joaquim](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/joaquim/32/16_2.png) [@Joaquim](https://forum.generic-mapping-tools.org/u/Joaquim)\
**Post date:** [December 22, 2021, 2:07am UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/6 "2021-12-22T02:07:24Z")

</div>

But the true problem is this bug (the region is not what was asked)

```auto
grdimage(Land_clip, region=(-180,-40,15,80), color=:geo, xaxis=(annot=15, grid=30), yaxis=(annot=15, grid=15), Vd=1, show=1)
 grdimage -R-180/-40/15/80 -JX14c/0 -Bpya15g15 -Bpxa15g30 -n+a -Cgeo -P -K > C:\TEMP\GMTjl_tmp.ps

```

No idea where it comes from. I suggest that you crop the grid at the same time you clip it to avoid this. i.e:

```auto
Land_clip = grdclip("@earth_relief_05m", low=[0 NaN], region=(-180,-40,15,80));

```

 ![GMTjl_tmp-fs8](https://canada1.discourse-cdn.com/flex047/uploads/gmt/original/2X/c/ce8ec207267a78e8ad31fbf4007b909929f93919.png)

---

<div class="post-metadata">

**Author:** ![leon6j](https://avatars.discourse-cdn.com/v4/letter/l/a6a055/32.png) [@leon6j](https://forum.generic-mapping-tools.org/u/leon6j)\
**Post date:** [December 22, 2021, 2:19am UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/7 "2021-12-22T02:19:13Z")

</div>

> [@Joaquim](#):
>
> `Land_clip = grdclip("@earth_relief_05m", low=[0 NaN], region=(-180,-40,15,80));`

Thank you so much! This solution fixes the issue.

---

<div class="post-metadata">

**Author:** ![Joaquim](https://yyz1.discourse-cdn.com/flex047/user_avatar/forum.generic-mapping-tools.org/joaquim/32/16_2.png) [@Joaquim](https://forum.generic-mapping-tools.org/u/Joaquim)\
**Post date:** [December 22, 2021, 2:45am UTC](https://forum.generic-mapping-tools.org/t/a-weird-issue-related-to-subplot/2460/8 "2021-12-22T02:45:54Z")

</div>

Note, this problem does dot occur if one use a projection other than a linear one.