I'm running a ROMS upwelling configuration for north-eastern off Iberia. To check if the model grid and other grid parameters are OK I setup a Run with no forcing, summer stratified configuration and a realistic smothed batymetry. I was exepcting to get a near-steday flow with near zero horizontal velocities.
In fact I did obtain velocities ~0.0001 m/s that seems to be quite large (see att). Can you tell me the expected model velocities error (with no forcing)
model horizontal velocities error
-
- Posts: 39
- Joined: Wed Feb 18, 2004 3:50 pm
- Location: Instituto de Meteorologia
model horizontal velocities error
- Attachments
-
- VH_02_smth_23dias.jpg (34.49 KiB) Viewed 7767 times
Re: model horizontal velocities error
I would expect something like what you are seeing.
Re: model horizontal velocities error
What velocities do you get if you do a simple geostropic balance using your T, S vertical stratification? The rho that you need for it can be calculated using the expressions in ROMS/Nonlinear/rho_eos.F (which you may need to code up in a separte routine in C++, Matlab or Fortran).
- jivica
- Posts: 172
- Joined: Mon May 05, 2003 2:41 pm
- Location: The University of Western Australia, Perth, Australia
- Contact:
Re: model horizontal velocities error
I should say that this is really good to do at the begining ,
just to make sure you have all setup right (!) and can explain "ghost currents" if you have them
What I do is run model similar like you but with analitical (for my region specific) vertical T & S; just to be sure that there are no horizontal gradient introduced in the initial phase. Than after some time (depending on the grid/region), say 2 weeks, you could notice that you have growing velocities at the bottom. This should represent mayor contribution from HPGE (Horizontal Pressure Gradient Error) because of imperfect vertical coordinate. Those regions with higher bottom velocities should correspond with regions characterised with higher values of rx1 -> (in a fact you can compute it and plot it)
If this is not acceptable for you reconsider to smooth bahtymetry more or change vertical structure (less vertical levels, higher theta_s ...)
Cheers
Ivica
just to make sure you have all setup right (!) and can explain "ghost currents" if you have them
What I do is run model similar like you but with analitical (for my region specific) vertical T & S; just to be sure that there are no horizontal gradient introduced in the initial phase. Than after some time (depending on the grid/region), say 2 weeks, you could notice that you have growing velocities at the bottom. This should represent mayor contribution from HPGE (Horizontal Pressure Gradient Error) because of imperfect vertical coordinate. Those regions with higher bottom velocities should correspond with regions characterised with higher values of rx1 -> (in a fact you can compute it and plot it)
If this is not acceptable for you reconsider to smooth bahtymetry more or change vertical structure (less vertical levels, higher theta_s ...)
Cheers
Ivica
-
- Posts: 39
- Joined: Wed Feb 18, 2004 3:50 pm
- Location: Instituto de Meteorologia
Re: model horizontal velocities error
I'm still getting horizontal velocities ~0.001 m/s with no forcing with bathymetry smoothing,the slope parameter maximum value is r = grad(h)/h = 0.2.
What is "no motion" with ROMS? What do you consider acceptable horizontal velocities with no forcing ?
Isabel
What is "no motion" with ROMS? What do you consider acceptable horizontal velocities with no forcing ?
Isabel