Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sun Aug 25, 2019 2:38 am




Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Author Message
PostPosted: Mon Jul 22, 2019 3:36 pm 
Offline

Joined: Mon Dec 31, 2018 8:42 pm
Posts: 5
Location: China Institute of Water Resources and Hydropower
contact.m works well for creating the nesting grid in the Cartesian coordinate system. However, when I made a nesting grid in the spherical coordinate system, I found there may be one bug in contact.m. Apologies if this is correct as it is.

lines 721-723
FSr.Values = G(dg).angle(:); R.angle = FCr(XrF, YrF);
FSr.Values = G(dg).f(:); R.f = FCr(XrF, YrF);
FSr.Values = G(dg).h(:); R.h = FCr(XrF, YrF);


should be
FSr.Values = G(dg).angle(:); R.angle = FSr(XrF, YrF);
FSr.Values = G(dg).f(:); R.f = FSr(XrF, YrF);
FSr.Values = G(dg).h(:); R.h = FSr(XrF, YrF);

However, one new problem occurs after I made these changes.

Subscripted assignment between dissimilar structures.
Error in contact (line 329)
S.refined(cr) = R;


Then I checked the function [S,G] = contact(Gnames, Cname, Lmask, MaskInterp, Lplot)
I found output S includs following two elements
S.refined(cr).xi_rho(:,:) - Receiver grid curvilinear
S.refined(cr).eta_rho(:,:) (XI,ETA) coordinates

However, R (R = refine_coordinates(cr, dg, rg, G, S, MaskInterp)) includes eight elements
R.xi_rho
R.eta_rho
R.xi_psi
R.eta_psi
R.xi_u
R.eta_u
R.xi_v
R.eta_v

Thanks a lot if you can look at them.


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 3 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group