Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sat Sep 22, 2018 3:27 am




Post new topic Reply to topic  [ 2 posts ] 

All times are UTC

Author Message
PostPosted: Fri May 28, 2010 3:37 pm 
Offline

Joined: Tue Mar 24, 2009 3:25 pm
Posts: 4
Location: MIT/WHOI Joint Program
Attachment:
DDF8_proposedgrid.pdf [27 KiB]
Downloaded 95 times


Dear ROMS friends,

After Dr. Arango detailed some of the difficulties associated with altering the periodic boundary conditions in ROMS (see post https://www.myroms.org/forum/viewtopic.php?f=3&t=1877), I turned to thinking about another possible domain configuration that might allow me to study boundary layers over sloping bottoms without the interference of boundaries at either end of the slope. Right now, I'm thinking about a grid like the one in the plot I've attached. The idea would be to have closed boundaries, with a test section (on the left) and a narrow channel (on the right) to allow a recirculation of the boundary layer flow, separated by a narrow island. This way, I can have a net flux in the boundary layer in the text section, but allow that flux to recirculate without disrupting the boundary layer dynamics. I've tried running ROMS using this grid, but I'm having some problems:

(1) My island is currently a seamount, with whose top is at h=5m below the surface. I didn't make it break the surface because I wasn't exactly sure what I would need to do with masking if I did. Is this a bad idea? I've never used a coastal mask (I've only used ROMS in idealized domains to study bottom boundary layers), and I couldn't find any information on the online documentation about how it should be applied. So the question is: should I use a coast here?

(2) I would like to force fluid from the test section into the area behind the island through the opening indicated with an arrow. Is there a way I can specify a velocity within the domain? I only know how to specify the velocity on the boundaries, using u3dbc and related functions, but this would be specifying a barotropic velocity across several internal grid points.

(3) For some reason, when I run the model with the attached grid it integrates in a reasonable way for an hour or so (of model time), then blows up. This is almost 1000 timesteps, so is there any reason why ROMS might integrate happily for 1000 timesteps and only then blow up?

Thanks so much for all your help!

Cheers,
Rebecca


Top
 Profile  
Reply with quote  
PostPosted: Fri May 28, 2010 5:56 pm 
Offline

Joined: Mon Aug 16, 2004 8:47 pm
Posts: 81
Location: U.S. Geological Survey, Woods Hole
Hi Rebecca-
(1) If you break the surface (i.e. negative h), mask_rho needs to have a value of 0 instead of 1. If you are using ana_grid, you can specify the masked areas easily enough using ana_mask.h. For this domain, you should probably have the seamount break the surface.

(2) I don't know if you can do this using UV_PSOURCE, basically specifying a point source of momentum in the middle of the domain. Alternatively, could you provide a flow bc on the northeastern side of your domain, forcing the flow south in the deep basin?

(3) The best way to diagnose a blow-up is to save the fields in the restart file often enough, so you can see what is blowing up where. I would suspect you have a problem near the seamount, which may be fixed by masking...


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

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


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