Masking Ocean areas in rotated Cartesian coordinates?

Q&A
1

I have an oblique rectangular domain, where I would like to plot a 2D data on top of topography.

I am playing with the domain in a Cartesian coordinate system:
Topography:

image
image1253×260 18.3 KB

Data:
image
image932×238 14.9 KB

What I am doing is finding out all Ocean areas in the topography, and try to mask them out in plotting the data, as seen in the bottom in the figure below. My approach is creating a grd file with NaNs for Ocean areas, then blending it with the data.grd:
gmt grdblend temp.grd coast_mask_extlarge.grd -Cu -I0.02/0.02 -Gtemp2grd
which, however, totally messed up my colormap.
image
With coast command, we can select this Ocean area (the bottom red part) and fill in some color, but how to do it in this Cartesian plot?

Thank you!

2

I have a workaround now, perhaps this is the correct way.
The idea is to create another grid file with NaNs and some certain number representing ocean. Then plotting this grid with -Q flag to visualize crust part, while giving some color to ocean.

gmt makecpt -Clightblue
gmt grdimage mask.grd -C -Q

I haven’t tried to use GMT to create such mask.grid from topography, instead I created a binary using Python, then convert it to mask.grid.