Ocean Modeling Discussion

I suspect the current velocities through that pass are just too high. Try reducing your timestep a little and see if that helps.

Most of the work calculating the 2D diagnostics happens in step2d_LF_AM3.h. I believe the forcing terms for the coupling are computed (at least in part) in rhs3d.F.

You may want to try adding the N2 constituent. For the elevation it is at least as large as the S2 and larger than both the K1 and O1, within Long Island Sound anyway. Is the East River resolved in your grid? The tidal forcing through the East River and the western end of LIS is an important ...

2) If can, How to do add the decay rate in the passive tracer ?

I don't believe that ROMS has this capability, but you can modify the code to do it. I once applied a simple exponential decay to some passive tracers. I think I added code to step3d_t.F.

Do you have UV_VIS2 or UV_VIS4 defined? If not, then NDM2d is off by 1 (it was never incremented) You will need to adjust mod_scalars.F. I would try this: # ifdef DIAGNOSTICS_UV ! ! Indices for 2D momentum diagnostic variables. ! Order is important! ! M2pgrd=1 ic=1 # if defined UV_COR M2fcor=ic+1 ic ...

Check the bibliography link on the ROMS homepage for a paper (with pdf link) on boundary conditions by Marchesiello et al. This paper also has references to many other helpful sources (Orlanski, Flather, Chapman, etc.) I recall finding the Palma & Matano papers quite helpful when first learning ...

I did this for pesticide concentrations some time back. They were just passive tracers and I added a bit of code to step3d_t.F to handle the decay. My situation was a little more complicated than Radium would be (increased decay rate in daylight), but it was actually pretty straightforward.

In addition to the boundary condition problem, you should think about the behavior of tides in a semi-enclosed basin, as Kate said. Changing the depth of the basin will change the speed of the wave, hence the transit time through your basin. Think about what that can do to the relative phases of the ...

This sounds familiar, but I think I was seeing it from the other side. Things would blow up after a restart, but then i would try restarting from a different time--sometimes an hour before or after, sometimes only a few timesteps away--and all would be well. I knew I was running really close to the ...

Way back when I submitted the code to output the diagnostic terms, I included a page or two of documentation. Not very wordy, but it does list the different diagnostic variables, etc.

It can be found by clicking Documentation on the ROMS homepage and choosing Diagnostics.

For a 3D application, the bottom stress is computed in set_vbc.F and is not computed in step2d.F at all. From step2d.F: # ifndef SOLVE3D ! !----------------------------------------------------------------------- ! Add in bottom stress ...

I am (finally) porting my application to version 2.0, but I have run into a problem. The array bound check is barfing on get initial: Run-Time Error 455: Array section bounds inconsistent with parent array In Procedure: get_initial Diagnostics Entered From Subroutine get_initial Line 249 Entered ...