Ocean Modeling Discussion

Thanks jcwarner, these are really helpful.

Hi,

Does anyone know of publications that specifically discuss the terms from the depth-averaged momentum (M2) diagnostics in their analyses? All the references I found so far seem to use only the M3 diagnostics.

Thank you!

Hi hugobastos, thanks for your reply. This plot of maximum points is a single point, or is the maximum in every timestep? I'm thinking that you probably used the maximum per timestep, which indicates that the location of the point can be varying right? Yes. My tip is that after the ~50day, the point...

Updates: I capped the topography at 1200 m instead of at 2000 m, keeping everything else the same, and the HPGE decreased from 10-30 cm/s to ~3 cm/s. Also, the sudden acceleration to ~40 cm/s in the first 5 days is gone. The HPGE now increases more smoothly (attached). After day ~50, weird banded st...

Thanks mastrorocco and hugobastos for bringing this vertical coordinate issue to my attention. It seems like a very good source of trouble. don't know why you are using Tcline = 0 and I am not sure if this is fine (perhaps someone else can comment on that!). I set Tcline to 0 in order to minimize bo...

Thanks for the reply, rduran. That doesn't seem too bad to me, if you were working with single precision (you are probably using double precision, but just to make a point) then that would be zero. It would still be good to know why its not more like 16 orders of magnitude (zero in double precision)...

Hi rduran, thanks for your suggestion. I looked at the ratio of the residue and the leading term, and the ratio of the residue and the smallest term. The ratios between the residue and the leading term show that the residues are everywhere 9-7 orders of magnitude smaller than the leading terms. The ...

Thanks jivica, this procedure seems very useful. I will follow that in another application where I'm seeking a realistic topography. But in this application with idealized topography, I don't understand why I'm getting this much HPGE with this configuration. I'm using safe rx0 and rx1 values (0.075 ...

Got it, saulo. I will fiddle with the grid a little more to decrease rx1. So based on the information I found so far and considering that I have tried pretty much everything, it looks like my slopes are just too steep. At least the largest errors are not located around the shelfbreak, which is the a...

rduran, did you figure out why the momentum balance terms didn't add up? I am having the same problem. I am computing the following sums for the 2D and 3D momentum balances, respectively (and analogously for the xi-direction balance): residue_2D = vbar_hadv+vbar_cor+vbar_prsgrd+vbar_sstr+vbar_bstr+v...

Thanks everyone for your replies. Following Ianerolle's suggestion, I did two test runs changing only the vertical stretching curve as depicted on my previous post. It seems that, for my application, the configuration with Vtransform = 2 and Vstretching = 4 produces less Horizontal Pressure Gradient...

Hi Ianerolle, thank you for your reply. With the parameters you suggested, the vertical stretching curve becomes the red line in the attached figure. The maximum Haney number rx1 increases from 3.7 to 14.1: (ORIGINAL) N, theta_b, theta_s, tcline: 50.0, 0.0, 0.0, 0.0 (MODIFIED) N, theta_b, theta_s, t...

Hello all, I'm having problems with pressure gradient errors over an idealized continental slope application. My initial field is just flat stratification, with an N2 profile derived from a climatological dataset. I have baroclinic (barotropic) velocity errors of 10 cm/s (4 cm/s) after a steady-stat...

Hello, I am having trouble using the LP_smoothing_rx0() and LP_smoothing_rx0_heuristic() functions of the bathy_smoother submodule of pyroms. It seems that I have built and installed bathy_smoother, lp_solve and lpsolve55 (with its python component) successfully, since bathy_smoother imports the lps...

Understood, Arango. It seems that ticket was just before I started learning ROMS. Thanks.

Hi Ivica, thanks for your reply.

Your suggestion seems to be better than the "cycle_length" approach I was using. Now that I have made this change, GET_2DFLD is reading the forcing file only once, and not every time step as was happening before.

Thanks,

Hello, I am new to ROMS, and would like to know whether I am setting the time dimension of the forcing variables correctly. sustr and svstr are given at a single time, and I want them to be applied at every time step. It appears that the routine get_2dfld.F is reading the forcing file properly, but ...