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 .
about conservation of roms
Re: about conservation of roms
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?
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?
Re: about conservation of roms
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?
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?
Re: about conservation of roms
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
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
Re: about conservation of roms
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.
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.
Re: about conservation of roms
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.
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.