Ocean Modeling Discussion

Dear Colleagues We would like to draw your attention to Session 058- Mesoscale ocean processes and their representation in earth system models at the 2014 Ocean Sciences meeting in Honolulu HI, USA, from February 23-28 2013. Please consider contributing a presentation. The session description is ...

Hi all, We are looking for a postdoctoral researcher in physical oceanography, please see the following. The Geophysical Institute at the University of Bergen invites applications for a 2-year post-doctoral research position in physical oceanography, under the “Faroe Bank Channel Overflow: Dynamics a ...

> the circulation is not that affected by local extrema (unless out of control) but is greatly affected by diapycnal mixing. I didn't say otherwise. But, when you have a very sharp front, convection zones or overflows, nonmonotnic schemes will also have diapycnal mixing due to local extrama and it ...

Thank you, I'm aware of the Flavien's work. I didn't say that monotonic schemes don't have any problems. All advection schemes have either dispersion or dissipation or both, but that shouldn't mean that we can use a scheme which is not physical and assume that it's fine because deep waters are not ...

Thanks for the quick response. I used RSUP3 for the lock-exchange example which is a standard two layer dam breaking. There was no topography so I evenly spaced the sigma layers to compare with the z-models and RSUP3 had less numerical mixing in that config. In fact, it has a nice convergence when I ...

Hi, I just have one question and one comment: My comment is we just have a new paper coming out about numerical mixing in idealized test cases and global experiments using ROMS, GOLD, MOM and MITgcm. We found that ROMS performs better than the MOM and MITgcm since it uses a high-order upwind biased ...

Hi all, I have a quick question about UV_U3ADV_SPLIT option. I understand that this option seems a delicate issue because of stability. I'm trying to diagnose spurious mixing in a simple problem and this option improves the model a lot. So I was wondering the last situation of this option. Is it ...

Hi, Did anybody have a chance to look at this issue? I also have a very simple setup similar to lock-exchange and I want to use MPDATA as a monotonic scheme, however I'm getting smaller and larger values than minimum and maximum initial temperatures, respectively. I'm using WALL options rather than ...

Hi everybody, In roms 3.0, there is a nice option that you can create your own Project directory and you can create a folder as Functionals and you can put *.h files you need to modify. So, you don't touch original *.h files in ROMS directory. We were wondering can we do the same thing for *.F files ...

Hi everybody, I need to find out the turbulent Prandtl number in KPP mixing scheme, I checked the code, and I found these lines; !For eddy viscosity Kv(i,j,k)=lmd_iwm + lmd_nu0m*nu_sx + lmd_nu0c*nu_sxc !For eddy diffusivity Kt(i,j,k)=lmd_iws + lmd_nu0s*nu_sx + lmd_nu0c*nu_sxc The first terms (i.e ...

Hi, I want to change my theta_b with depth just as "lanerolle" did it in the past, by chance could you send me your code or explain how you did it? Thanks, Mehmet Hi Everybody, I have coded up a new sigma formulation in to ROMS where theta_b is no longer a constant (but theta_s, Tcline, hC are all ...

Thanks for responding, But my problem was more point source tracer concentration rather than initial conditions. Negative concentration at initial condition is probably because of Gibbs oscillations. Since when I used MPDATA there were no negative values. But when I used point sources, even using ...

Hi everybody, I got a problem about passive tracer concentration. I release some dye in my simulation as t(i,j,k,1,itrc) = 100.0_r8 Dye is located as a box volume. It is 100 in a specific location and zero everywhere. I was using 3rd order upwind scheme and I got negative values, so we thought maybe ...

Hi, I am a new user in ROMS, and I try to model an overflow, I need a constant salt forcing at the northern side of the domain, Is there a way to relax some grids of the domain into specified number, just like regular relaxation, but not only at boundaries also in the domain, or anything works for ...