GenericMappingTools / GenericMappingTools/gmt

grdcontour -C & externals is wrong

Open
#3,701 18 comments 0 reactions 0 assignees View on GitHub
stale
Dominant language
C
Stars
979
Forks
414
Avg merge
17h 26m
Merged PRs (30d)
54

Description

```
julia> C = makecpt(T=(-7,9,2), Vd=1);
makecpt -T-7/9/2
```
Reading grid from file it's (nearly) OK. Nearly because the aliasing is clearly visible
![p_good](https://user-images.githubusercontent.com/537321/88121141-be95cd00-cbbc-11ea-8d3a-4ad584b6abe3.png)

Now do the same but using data in memory.
```
julia> Gp = gmtread("peaks.grd");

julia> grdcontour(Gp, C=C, N=true, savefig="p_bad.png", show=1, Vd=1)
grdcontour -JX12c/0 -Baf -BWSen -C -N -P -K > C:\TMP\GMTjl_tmp.ps
```
![p_bad](https://user-images.githubusercontent.com/537321/88121181-e127e600-cbbc-11ea-9c40-348da03e3866.png)

Finally, another bug is when trying to do the same as a modern one liner
```
gmt grdcontour peaks.grd -JX12c/0 -Baf -BWSen -Cc2.cpt -N -png C2
```
Here it creates two separate figures. One with the colors and another with the contours.

The grid is here but not zipped. Just remove the ``.zip`` (needed to cheat github)
[peaks.grd.zip](https://github.com/GenericMappingTools/gmt/files/4957187/peaks.grd.zip)

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.