Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Fri Sep 22, 2017 4:33 am




Post new topic Reply to topic  [ 4 posts ] 

All times are UTC

Author Message
 Post subject: Blowing-up problem
PostPosted: Mon Nov 01, 2010 9:03 am 
Offline

Joined: Thu Oct 04, 2007 4:14 am
Posts: 31
Location: ECNU/NCSU
I'm trying to setup a tidal model with wetting-drying scheme. The minimum bathymetry in my grid is 0.1m, and the minimum depth for wetting and drying was set to 0.5m. Blowing-up occurs at the step 16979. Only tide forcing was added to the model.

DT == 60.0
NDTFAST == 10

Minimum barotropic Courant Number = 1.00127762E-02
Maximum barotropic Courant Number = 4.44498424E-01
Maximum Coriolis Courant Number = 4.82332138E-03

rx0 = 2.000000E-01 (Beckmann and Haidvogel)
rx1 = 1.800000E+00 (Haney)
Attachment:
Fig1.jpg
Fig1.jpg [ 22.52 KiB | Viewed 936 times ]

I checked the restart file and found that the blowing-up was occurred by abnormal v velocity. (See Figure 2, in grid cell x = 185, y = 123),
Attachment:
Fig2.jpg
Fig2.jpg [ 300.34 KiB | Viewed 936 times ]
I also plotted the zeta and v one step before the blowing-up (Figure3).
Attachment:
File comment: (one step before)
Fig3.jpg
Fig3.jpg [ 60.82 KiB | Viewed 936 times ]

I can’t identify what caused the abnormal v velocity.
Can anyone give some suggestions?


Top
 Profile  
Reply with quote  
 Post subject: Re: Blowing-up problem
PostPosted: Tue Nov 02, 2010 1:10 am 
Offline

Joined: Tue Aug 21, 2007 5:44 pm
Posts: 10
Location: North Carolina State University
Haibo,
Is it caused by DT/DTFAST too large?
CFL condition satisfied yet?


Top
 Profile  
Reply with quote  
 Post subject: Re: Blowing-up problem
PostPosted: Tue Nov 02, 2010 5:14 pm 
Offline
User avatar

Joined: Thu Dec 07, 2006 3:14 pm
Posts: 76
Location: USGS
Which version of ROMS are you using?
Did you try a larger NDTFAST? (30 maybe?)
Zafer


Top
 Profile  
Reply with quote  
 Post subject: Re: Blowing-up problem
PostPosted: Tue Nov 02, 2010 6:36 pm 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1011
Location: IMCS, Rutgers University
Yes, we have mentioned several times in this forum about of the time-stepping properties of ROMS. There are several publications about this, see any of the Shchepetkin et al. papers. They discuss this in great detail. User of ROMS need to read the basic core publications about ROMS before starting running your applications. User should be aware that ROMS has an extensive literature with hundreds of papers. There is a lot of great tips and information in the scientific literature that maybe relevant to any your applications. I encourage users to read this literature.

We actually recommend having a value for NDTFAST to be larger than 30. A value of 10 is too small. You need to fine tune DT accordingly to get the appropriate baroclinic time-step. Of course, you need to first set your barotropic time-step according to the fast gravity-wave dynamics which primarily controls the CFL in ROMS.


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

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 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