I found a bug in obs_roms2roms.m. see code below
% Receiver grid.
if (isvec && D.vector_rotation),
if (D.spherical),
Xname = 'lon_rho'; % If requested, interpolate
Yname = 'lat_rho'; % U- and V-points variables
else % to RHO-points instead to
Xname = 'x_rho'; % facilitate the curvilinear
Yname = 'y_rho'; % rotation to true East and
end % North. This rotation is
Mname = 'mask_rho'; % done somewhere else.
if (is3d),
Zname = 'z_r';
end
end
The highlighted part should be R.vector_rotation, not D.vector_rotation. If your outer grid do not need to rotate, while your inner grid needs to rotate, you will see the problem. It took me a whole day to find out this tiny bug... I hope this can be changed in the newer release.
bug in obs_roms2roms.m
-
- Posts: 60
- Joined: Tue Sep 29, 2009 3:50 pm
- Location: School of Environment System Engineering,UWA