Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Tue Nov 21, 2017 12:50 am




Post new topic Reply to topic  [ 3 posts ] 

All times are UTC

Author Message
 Post subject: Volume not conserved.
PostPosted: Thu Apr 20, 2017 1:06 am 
Offline

Joined: Tue Nov 10, 2015 3:51 pm
Posts: 7
Location: The University of Western Australia
Hi,
I am running ROMS on a semi-enclosed embayment with result output at each time step (10 seconds), and the result shows that volume is not conserved as I compute the volume fluxes at the four boundaries of a water column, of which the grid size is 555 m * 567 m. The sea surface level of the water column is controlled by the M2 tide with amplitude roughly 2 m, but the total boundary volume flux stays positive all the time, indicating that the volume of the water column is continuously increasing (see the image below).

Image

Is the grid too big to test the volume conservation, or there is some other issue missed? Any suggestion is welcomed and appreciated.

Cheers.
Wencai


Top
 Profile  
Reply with quote  
PostPosted: Thu Apr 20, 2017 4:32 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3251
Location: IMS/UAF, USA
I can't see your image.

I don't usually worry about volume conservation, but I don't generally have major draining or filling over time, just ups and downs. What are your boundary conditions?


Top
 Profile  
Reply with quote  
PostPosted: Fri Apr 21, 2017 1:27 am 
Offline

Joined: Tue Nov 10, 2015 3:51 pm
Posts: 7
Location: The University of Western Australia
Attachment:
File comment: variation of SSL vs. total boundary volume flux of the water column.
zeta_Q.png
zeta_Q.png [ 20.65 KiB | Viewed 360 times ]


Hi kate,

I have uploaded the image to attachment, and the lateral boundary condition type is shown below:

! Set lateral boundary conditions keyword.
! W S E N
! e o a o
! s u s r
! t t t t
! h h
!
! 1 2 3 4

LBC(isFsur) == Cha Clo Cha Cha ! free-surface
LBC(isUbar) == Fla Clo Fla Fla ! 2D U-momentum
LBC(isVbar) == Fla Clo Fla Fla ! 2D V-momentum
LBC(isUvel) == Rad Clo Rad Rad ! 3D U-momentum
LBC(isVvel) == Rad Clo Rad Rad ! 3D V-momentum
LBC(isMtke) == Gra Clo Gra Gra ! mixing TKE

! Set lateral open boundary edge volume conservation switch for nonlinear model and adjoint-based algorithms.
VolCons(west) == F ! western boundary
VolCons(east) == F ! eastern boundary
VolCons(south) == F ! southern boundary
VolCons(north) == F ! northern boundary

The water column where I test the volume conservation is far away from the open boundaries of the study area, so I think it may have little to do with the boundary conditions, isn't it?

Cheers.
Wencai


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

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