Ocean Modeling Discussion

Just to verify, you've set it up to run in parallel, right? (v.s. serial runs on different nodes).

Thanks for the feedback!

I'm wondering if anyone has experience installing ROMS on a Mac Pro (4-12 core) for parallel processing? Are there any problems that one should be aware of before considering a Mac Pro for this application?

Hi again,

Good point by you, re: the 2D bit. I will be interested to hear what you discover.

Best,
Rachael

How about bathymetry? Have you compared your ROMS model bathymetry to the actual values at the ADCP locations? Has smoothing (to meet Haney/rx0 criterion) altered the bathymetry in this location by an appreciable amount?

Hi Rodrigo, I am not familiar with the use of passive tracers as a dense mass of Lagrangian floats, but my experience with using passive tracers would advise against it. I could be wrong, but you may want to consider floats instead. In my application, an instability in the passive tracer caused the ...

I keep forgetting to post the solution, so--before too much time goes by--here it goes.... The problem that I described above was caused by 0m water depth in the mask region, which resulted in -dZ under land mask. This -dZ was not a problem for the closed boundary or gradient boundary conditions. It ...

I am having a hard time figuring out why my open boundary run doesn't run in ROMS 3.2. When I say "doesn't run", I mean that it blows up at first time step. Here is what works: 1) closed boundaries 2) open boundary on the east (closed N,S,W). East is only open water. 3) open boundaries with non ...

Thanks for the tip (a.k.a. read the comments, dummy!). I can't imagine anything more frustrating from a developer standpoint than providing the information and having it ignored. My apologies. You can commiserate with my husband about how terrible I am about reading directions. At the very least, I ...

I repeated what I did on Monday and the code seems to be running fine! In all likelihood this means that I wasn't strictly repeating myself :wink:; unfortunately, since I don't know what I am doing differently, I don't have any wisdom to impart. Hopefully, no one else will be able to replicate what ...

I am getting the following error when I try to run the UPWELLING test case on a current download of ROMS that I am configuring on my mac: ROMS/Bin/cpp_clean /Users/rmueller/ROMS-3/Projects/upwelling/Build/mod_kinds.f90 cd /Users/rmueller/ROMS-3/Projects/upwelling/Build; /usr/bin/mpif90 -c -frepack ...

Hi Diego, Another potential contribution is the 'angle' parameter. I see a note in a previous post about adding the scalar value of 'angle' to the output, but no mention of the 'angle' array (unless I missed it elsewhere?). I borrowed some code from the seagrid routine to come up with the following ...

Bug found, captured, and made into butterfly.

I believe the smth_bath call in easygrid should be as follows:
h = smth_bath(h,mask_rho, order, rlim, npass);

instead of
h = smth_bath(h,mask_rho);

The later works, it simply applies the default values every time.

Cheers,
Rachael

To prevent any further assault on his name, I just wanted to publicly acknowledge that the "Kevin" mentioned in an earlier post is NOT a nickname for Hernan! I hope my error doesn't inspire new users to do the same. It's just a mistake that I made to publicly embarrass myself for calling Hernan by ...

Good idea to describe the more general case that I was conveniently avoiding. :-) I'm working in the southern hemisphere too and had to convince myself of the signs. My coast is also aligned N-S rather than E-W (hence, u=0). The full expression for my case is a simplified version of: v = eta_0 ...

In the aforementioned reference to "GRID(ng)%yp(Istr-1,j)" values being null, I forgot to clarify that my y_psi(x_psi) values in my non-analytical grid are non-zero . If I'm understanding all this correctly, my grid's y_psi(x_psi) values should populate the yp(xp) variable in mod_grid.F. It doesn't ...

Hi me, The way I fixed this (thanks to Susan H.!) was to 1) added #define ANA_SRFLUX in CPPDEFS.h file 2) changed the "j" to "i" index in the analytical. F file. I made other changes to the indices too, but those were different kinds of bugs than I was reporting here. Unfortunately, it's not smooth ...

Thanks for the commentary Kevin. Just to prevent further confusion in interpreting your representation and my earlier comment, the equation can be written in two ways: U = Uo*exp(-f*y/c)*cos(omega*time), for f>0 where Uo = c*eta_0/H, c = sqrt(g*H) OR Uo = eta_0*sqrt(g/H) I think there may be another ...

Thanks for mentioning this Alexander. The dimensions actually do work out because the equation I posted is incomplete. There is an amplitude scaling in the code that provides the 'm' dimension. I really should have written: U = eta_o * sqrt(g/H) * exp(-f*y/sqrt(g*H))*cos(omega*time), for f>0 In the ...

This may be fixed in ROMS 3.X, but there is a bug in the KELVIN case for ROMS 2.2. The Kelvin wave solution should be of the form U = sqrt(g/H)*exp(-f*y/sqrt(g*H))*cos(omega*time). The KELVIN case in Analytical.F file has the form of: U = sqrt(g/H)*exp(-f*y/sqrt(g/H))*cos(omega*time). I'm still a ...

Hello everyone, Yes--I know--2.2. :? I have plans to upgrade to 3.1 later this summer, I promise. In the meantime, I would like to create an idealistic run that has a simple kelvin wave forcing along the north/south boundary and walls along the east and west. I am trying to do so analytically ...

Roger that, Enrique. Mike Dinniman (CCPO, ODU) helped point out that the version I am using isn't "out of the box" and his improvement included the addition of thermodynamics. At the time of posting, I thought that this may be a good forum discussion point but it seems that my question is not as ...

Well that's not helpful! eh?

take two.

try:
http://staff.acecrc.org.au/~johunter/is ... vities.pdf

Or go to John Hunter's home page and look for "ISOMIP" files.

I'll answer my own question in case someone else is looking for similar information. Check <http> and refer to: Holland and Jenkins (1999) Grosfeld et al. (1997) Jenkins and Holland (2002) Pereira et al. (2002) Blumberg and Mellor (1987) I haven't yet checked the equations to see where specific ...

I am looking for a good reference on the equations used in the ROMS ice shelf coupling, but I am obviously not looking in the right places. Everything that I find seems to refer to the sea ice coupling rather than the ice shelf. What is the paper(s) from which the ROMS ice shelf equations are ...

Thanks for a heads up on the difference with the orthogonality error. Also, the clarification about the paths is helpful. For some reason, on most websites that discuss installing mex and netcdf, they insist that the .dll and .mexw32 need to be in \bin and bin\win32. It worked, so I never really ...

I recently went through a matlab/mex upgrade and seagrid installation. From this experience, Kate, it seems like your error may be more related to the matlab/mex communication than Jamie's code? Could this be the case? I came across a similar error, and it was because of path settings. In case it is ...