wrapping for global OPeNDAP data

[Sergey Varlamov]

A Happy New Year!

Does anybody knows how easily handle the global OPeNDAP data sets
like this on
http://usgodae2.fnmoc.navy.mil:80/dods/GDS/fnmoc_glb_ocn/FNMOC_GLB_OCN_uucurr
?
I need to plot the map that cross the initial longitude line.
GrADS fine handle the global GRIB data, but I failed to do same
with this and some other OPeNDAP (netCDF?) ocean data sets.

The mentioned above data set covers global grid for
lon=0,360, step 0.25 (x=1,1441) with first and last longitudes dublicated.
Using
'set lon -10 10'
or
'set lon 350 370'
do not work in this case plotting only lon=0:10 data (first case,
the left part of map keep undefined values) or 350:360 data for second case.

The only solution that I found in the GrADS discussions was to
overlay these two plots, but it requires a lot of additional efforts
to support same colors or arrows scaling for two parts of image etc.
Trying to create single defined variable from two intermediate ones,
like it is shown below, also fails and plots only one zone:

'set lon 350 370'
'define u1=u_current'
'set lon -10 10'
'define u2=u_current'
'define uu=const(u1,0,-u)+const(u2,0,-u)'
'd uu'


Does anybody knows a way to get a 'single' data set wrapped around the
initial longitude like it goes for the GRIB data sets?

I will higly appreciate recomendations of experts.

Sergey Varlamov
sergey.varlamov at meteo.fr