Meet a problem while using contact.m

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
Zrj_tio98
Posts: 2
Joined: Mon Nov 22, 2021 11:05 pm
Location: NUIST

Meet a problem while using contact.m

#1 Unread post by Zrj_tio98 »

Dear all:
I am using Matlab scripts contact.m to set contact points between ROMS nested grids.
I have created my parent grid and used coarse2fine.m to create a child grid successfully.
However, when I try to use contact.m, I face a problem in line 715: FCr.Values = G(dg).angle(:); R.angle = FCr(XrF, YrF);
It tells me that Interpolation is unavailable because Interpolation requires at least two sampling points per grid dimension.
Through debugging, I change the code to FCr.Values = G(dg).angle; R.angle = FCr(XrF, YrF);
And it works. Does anyone ever face a problem like this? Do the contact files I made correct? Thanks!
Attachments
EAST_grid_nest1.nc
fine grid
(5.38 MiB) Downloaded 164 times
EAST_grid_contact.nc
contact file
(2.17 MiB) Downloaded 150 times
EAST_grid.nc
parent grid
(4.92 MiB) Downloaded 173 times

User avatar
wilkin
Posts: 922
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: Meet a problem while using contact.m

#2 Unread post by wilkin »

I get that error recomputing the contacts points file for my DOPPIO and PIONEER application. I had previously computed the file some time ago with an earlier version of contact.m.

So, yes, it looks like a bug has crept in.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

Post Reply