Ocean Modeling Discussion

Hi Guangyu,

You might find this paper useful: https://doi.org/10.1016/j.joes.2019.12.002

There are two different approaches to force the tides: 1) through the boundaries by downscaling from a tide-resolving model (for which you turn on add_fsobc and add_m2obc) and 2)from a tide model by ...

Giannis,
Another thing you might want to look at are your barotropic boundary conditions. If you are using ADD_FSOBC & ADD_U2OBC as well as SSH_TIDES/UV_TIDES, and if the regional model you are downscaling from is not de-tied (i.e. it has tides) you might be giving ROMS the tidal forcing twice ...

Hello,
The way I understand it is:

dstart = time_start - time_ref
tide_start = tide_base_date - time_start

Where time_ref='19960101', time_start='19960501' and tide_base_date is the base date you used when generating the tide forcing file. Basically if your tide base date pre-dates your ...

Hello,
I am having issues trying to use the lagrangian floats implemented in ROMS.
I compiled ROMS using svn version 655 with #define FLAOTS. My input file
looks like this:


Lfloats == T

! Switch to control the printing of floats positions to standard output
! file, [1:Ngrids].
!
Fprint ...

Ok, I see your point it would be ideal for me to test both ways.
I don't have any experience filtering the tides and from what I have read so far it doesn't seem trivial. Do you have any material that may be helpful, or could you point me in the right direction?
Any suggestion is welcome and ...

Thanks for the reply Kate,
Boundary conditions from NCOM are given every 3 hours and I would say tides are really good (they also assimilate altimetry) but I can't say if better than TPXO, probably not.
So it would probably be better if I provide tide currents from OTPS or adcirc and filter tide ...

Quick question,
I want to set up an operational system simulating full ocean physics in the Caribbean, meaning I want to simulate tides along with wind induced currents and so on. I am providing tidal forcing (tidal ssh and tidal induced currents) at the open boundaries via a netcdf tide forcing ...

Try putting a dot before the multiplication in {sqrt(-1)*v}
So it should look like {sqrt(-1).*v}
I think the reason dimensions do not agree is because you are doing a matrix multiplication of a single value and a matrix instead of an element-wise (or point-wise) multiplication of the matrix with the ...

gloeg002:
I have roughly 1 year of experience with ROMS but I'm pretty sure your problem is with your vertical grid. As you can see the model is blowing up at your first level (111 , 96, 01) , but more importantly your maximum grid stiffness ratios are off the roof! This is specially important if ...

I had a similar problem in a real application. I would suggest as Kate mentioned to provide boundary values (including zeta), the better the temporal update on the values the better solution you may get. What I do in my model is extract initial and boundary values at 3hr intervals from a coarser ...

Safoora,
etopo2 is a 2min Global bathymetry by ngdc, WOAXXXX and COADS are the World Global Atlas and the Comprehensive Ocean Atmospheric Data Set, both global. These two are good for basin scale resolution or lower. If you are working with regional models closer to shore with resolutions higher ...

I have never tried to load two coastlines at the same time in Seagrid, but I do know that coastline files are just a longitude, latitude pair where coastline segments are separated by Nans. So I would just try to separate file 1 and 2 by a pair of Nans so that Seagrid knows these are separated ...

Hi,
Since you have an error in:

% Create S-coordinate system: based on model topography h(i,j),
% fast-time-averaged free-surface field and vertical coordinate
% transformation metrics compute evolving depths of of the three-
% dimensional model grid.
%
hinv=1./h;

and:

% Get the sigma depths ...

Hi Jerry,
I got this same error a couple of weeks ago using ROMS AGRIF version. You error is most probably because you are not giving the correct path for the netcdf libraries. Make sure your netcdf libraries are set correctly in your jobcomp script. From waht you posted my guess is your netcdf ...

Hi Iassad,
Here are some observations on your romstools_param.m that may be the source of your problem:

"grdname=['c:\Litos\ROMS\tools\ROMS_TOOLS\roms_queiroz_grd.nc'];"
I do not know why you specifed this path, but the routine make_grid.m may not find the grid if you save it in another directory ...