[gradsusr] Generating image overlay for Open Street Map / Mercator Projection

[Jonas Kaufmann]

Hello,

I am trying to generate an image with GrADS that should cover the
whole world. The image should be based on an forecast from the GFS
model. I want to use this image as a layer on top of Open Street Map.

I can already create an image without borders, labels etc that I can
put on my map. However it seems that the projection still is wrong,
e.g. the continent borders drawn by GrADS do not match the borders on
the map. I think this might be related to the map projection used, in
OSM this should be Mercator projection which is not available in
GrADS. The complete script can be found at the bottom of the mail.

Is there a way to make the forecast map fit on OSM?


Thank you for your help!


Jonas


-------

'open gfs'

'set gxout shaded'
'set parea 0 11 0 8.5'
'set mpdraw off'
'set frame off'
'set xlab off'
'set ylab off'
'set mproj latlon'
'set grads off'

'set mpdset hires'
'set map 1 1 1'

'set lon -180 180'
'set lat -90 90'

'd mag(ugrd10m, vgrd10m)'

'printim world.png x512 y512'
'quit'