Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Tue Sep 18, 2018 3:39 pm




Post new topic Reply to topic  [ 6 posts ] 

All times are UTC

Author Message
PostPosted: Tue Dec 18, 2012 8:57 am 
Offline

Joined: Fri Nov 02, 2012 3:25 pm
Posts: 19
Location: xmu
Hi everyone,The roms is applied in estuary, and i set some sections in estuary.
but i found the water flux in every section is not the same and some have big difference.
Is there any nonconservation in computering of roms? and is there some options in cppdef.h
needed to defined to ensure the conservation?
thank you .


Top
 Profile  
Reply with quote  
PostPosted: Tue Dec 18, 2012 5:38 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3501
Location: IMS/UAF, USA
Sasha Shchepetkin has written at length about the conservation properties in ROMS. I suggest you read his papers for all the details. There is no cppdef option to change it.

I'm not sure what you are expecting, though. Do you have an estuary with channels between the parts? As long as flow is allowed to go from one box to the next, there's no guarantee about conservation in any one box. Perhaps you could explain your issue more clearly?


Top
 Profile  
Reply with quote  
PostPosted: Wed Dec 19, 2012 2:33 am 
Offline

Joined: Tue Nov 18, 2008 5:06 pm
Posts: 3
Location: xiamen university
thank you kate,
In the estuary ,i set some sections from head to mouth. And there is no branch, so the
tidally averaged water flux in the each section should be equal, at least there is almost
the same. But the water flux in the each section is different and from head to mouth the
value decrease strikingly. At the first i suspect the way to computer the flux may be have
some error, but i tride many methods,the result doesn't change. So i suspect the nonconservation
exit in the roms. can you give some advice ? thank you
.I found xxx_VOLCONS(use if Eastern edge mass conservation enforcement) in cppdefs.h , is this
opions can solve this problem?


Top
 Profile  
Reply with quote  
PostPosted: Wed Dec 19, 2012 2:54 am 
Offline

Joined: Fri Nov 02, 2012 3:25 pm
Posts: 19
Location: xmu
thank you kate,
In the estuary ,i set some sections from head to mouth. And there is no branch, so the
tidally averaged water flux in the each section should be equal, at least there is almost
the same. But the water flux in the each section is different and from head to mouth the
value decrease strikingly. At the first i suspect the way to computer the flux may be have
some error, but i tride many methods,the result doesn't change. So i suspect the nonconservation
exit in the roms. can you give some advice ? thank you
.I found xxx_VOLCONS(use if Eastern edge mass conservation enforcement) in cppdefs.h , is this
opions can solve this problem?

this is my cppdef:
# define SOLVE3D
# define UV_COR
# define UV_ADV
# define UV_C4ADVECTION
# define UV_LOGDRAG
# define DJ_GRADPS
# define TS_MPDATA
# define AVERAGES
# define SPHERICAL
# define MASKING
# define SALINITY
# define NONLIN_EOS
# define SPLINES
# define ANA_STFLUX
# define ANA_SSFLUX
# define ANA_SPFLUX
# define ANA_BSFLUX
# define ANA_BTFLUX
# define ANA_BPFLUX
# define TS_PSOURCE
# define UV_PSOURCE
# define SPONGE
# define VISC_GRID
# define DIFF_GRID
# define UV_VIS2
# define TS_DIF2
# define MIX_S_UV
# define MIX_S_TS
# define UV_SMAGORINSKY
# define MY25_MIXING
# ifdef MY25_MIXING
# define N2S2_HORAVG
# define KANTHA_CLAYSON
# define K_C4ADVECTION
# endif
# define RAMP_TIDES
# define SSH_TIDES
# define ADD_FSOBC
# define UV_TIDES
# define ADD_M2OBC


# define EAST_FSCHAPMAN
# define EAST_M2FLATHER
# define EAST_TCLAMPED
# define EAST_M3CLAMPED

# define SOUTH_FSCHAPMAN
# define SOUTH_M2FLATHER
# define SOUTH_TCLAMPED
# define SOUTH_M3CLAMPED


Top
 Profile  
Reply with quote  
PostPosted: Wed Dec 19, 2012 6:29 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3501
Location: IMS/UAF, USA
Can you show us a picture of your domain?

How are you computing streamflow at the various locations and where are they on the map?

I never use the VOLCONS options because I run with tides and they are inconsistent with VOLCONS.

From your cppdefs:
* I see you have an older version, before the LBC variable.
* I would not use MIX_S_TS, which is asking for cross-isopycnal mixing in the most unphysical way.
* I believe the gang would recommend using MY25 via the GLS_MIXING option instead.


Top
 Profile  
Reply with quote  
PostPosted: Mon Dec 24, 2012 1:33 am 
Offline

Joined: Mon Apr 06, 2009 5:12 pm
Posts: 43
Location: The Hong Kong University of Sci&Tech
That *_VOLCON options do nothing to your problem here, because it only conserves the inward and outward volume transport at the open boundaries.

That non-conserved estuary circulation looks strange, and I never see that kind of problem in my idealized estuary cases with buoyancy forcing only. I suggest you double check your program that calculating the "flux" (I supposed it as "volume transport"), and I notice that you included tides into your model. Please make sure to use the residual component in your calculation.


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 6 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