Ocean Modeling Discussion

I had forgotten this page existed. Sandro and I used it for the Bevano River inlet on the east coast of Italy.

If I remember correctly, I did the above fix at one point and it worked fine.

SWAN still needs the grid, bathymetry, winds, and boundary forcing just like it was a stand alone run.

If it is a ROMS/SWAN coupled configuration try looking at the waves and stokes velocities. They may not depend strongly on the fine scale changes in the grid making it not very sensitive to smoothing the grid. If something is happening making the gradient in the waves high it might cause a blow-up ...

The flags also extend to Sediment (flag=4), and I think potenially biology variables but I would have to go look in the code. The best thing to do is to look in the source code and see if there is a flag for a passive tracer. If not, it wasn't that difficult to add the sediment one you might be able ...

It is like salt. Specify the suspended sediment concentration at each time-step, for each vertical layer, in each river. So, it should be a three dimensional matrix, although if only using one river one dimension will be length one. I copied my .cdl file below for reference. netcdf waipaoa_river.nc ...

Search the forum and the wiki for the correct cpp options to specify a point source (river) discharge for your application (you are right you should specify RIVER_SEDIMENT). The CDL file for creating the ROMS forcing file contained in the svn source code has some of the needed variables in it for an ...

I am using it to study sediment transport in Poverty Bay, New Zealand.

Since I can't edit the WIKI...I'm just going to say that I found an issue with what is in the wiki and different versions of (updated) matlab scripts. I had to add another input variable to one of the function calls in the main matlab script for the OSU tidal predictions stuff. Since I can't leave a ...

As preemptive info you may already know.

You will have to change all the line endings in all the files if you are not running it in Windows.

In Linux/Unix this can be done with the recode command.

From a previous post... Posted: Sun Feb 25, 2007 5:55 pm Post subject: -------------------------------------------------------------------------------- Try setting the point source to the middle of the grid, and ramping up the discharge so it is easy to see. Although not realistic, it allows you to ...

I think you can use the recode command in linux to change them all to the correct line endings.

We are having problems along tiles, and one that is a little harder to track down, in our developmental bottom boundary layer gravity flow code. We want to make sure we are mapping u and v points to the rho points and the rho points to the u and v points properly. I have questions regarding the use ...

Yes it is in the mod_ncparam.f90, but it looks like the left hand side is not lined up with the rest of the code, it is off by one space, which is where emacs tabs it to. Is the tab breaking it? Why wasn't I supposed to use a tab? And yes it is in varinfo.dat varinfo.dat 'CsedW' ! Output 'sediment ...

I doesn't make it that far, it doesn't write out anything.

I do get a

before1 0 (from the same file)


and

idCwbl 0 (from inp_par.f90)

So space things out using multiple spaces? What breaks from hitting tab instead? oceanG gives me At line 618 of file inp_par.f90 Traceback: (Innermost first) Called from line 58 of file inp_par.f90 Called from line 80 of file ocean_control.f90 Called from line 67 of file master.f90 Fortran runtime ...

I am trying to add some new variables to the history and averages files, but am running into a snag. Below I will post the name of the file I modified, and what I added, shortened to a single variable only, even though I am adding a total of 5. The error, I believe, is that it is not reading in the ...

Try setting the point source to the middle of the grid, and ramping up the discharge so it is easy to see. Although not realistic, it allows you to make sure it is being set at the exact location you think, and is being discharged through the correct grid cell face (correct direction). Then, simply ...

Here is what we changed to get riverine sediment. added the river sediment to varinfo.dat 'river_mud1' ... ... 'river_sediment, scalar, series' 'river_time' 'idRtrc(idmud)' 'nulvar' 1.0d0 in get_data.F add # ifdef RIVER_SEDIMENT ELSE IF (SOURCES(ng)%Fsrc(is).eq.4.0_r8) THEN SOURCES(ng)%Lsrc(is,itemp ...

I am trying to nest a finer grid over the Po River delta region in the Northern Adriatic Sea. I am using ROMS version 2.1 and it runs fine until I try to add the second grid. Currently I have the two grids and have changed Ngrids in mod_param.F to 2 and set the size of the second grid. Lm = (/ 158 ...

We set river_Vshape to be a 'number of rivers' by 'number of vertical levels' array. for i=1:Nrivers % 20% of transport in each of top 4 levels, % then 10%, 5%, 3%, 1%, 1% below % N=20 - the number of vertical layers River.vshape(i,N )=0.20; River.vshape(i,N-1)=0.20; River.vshape(i,N-2)=0.20; River ...

About not putting all the flow into the top layer and the model numerics: We may have some evidence for not putting it all in the surface layer. When the river discharges are bumped up unexplainable (at least by me) floating point and overflow errors seem to miraculously occur, sometimes. The only ...

Is it our consensus that the best way to put freshwater in is as a surface plume instead of as a well-mixed layer?

I would have thought to put them in as a well-mixed layer with freshwater spread throughout all vertical levels.

Courtney and Aaron