Ocean Modeling Discussion

Its best to close the boundaries to evaluate the error for sure. Looking back at some of your plots there seems to be a lot more activity near the boundaries than expected. For an application with no forcing and no horizontal variability in T/S you really don't have to worry about boundary reflectio...

The error decreased because I decreased rx1 alone. Actually my case was I had high vertical resolution, a large number of vertical levels and I would had to increase the horizontal resolution more than I could aford to bring rx1 down (so the reverse of what I said, sorry about that, it was a while a...

The rx1 is a bit tricky as it seems application dependent. A value considered high may work on a given setup while a "low" value may not (a setup that depends largely on the topographical slope, which looks real steep in your case and the sharpness of gradients but also on the algorithms y...

So ressurecting this thread... I just got revision 730. It seems Bulk_flux.F still does not take surface currents into account when computing stress. Is this correct? If so, are there other reasons for not including it (other than the reservations discussed above), such as of numerical design in nat...

Hi, You probably disrupted the python indentation when commenting a line of the grid.py, which if I remember is a necessary step on the instructions. Make sure all lines in this function are 4 columns from beguining of def statement. You can use an editor with automatic python indentation as well. M...

Hi, Hsbl is not the mixed layer. It's the boundary layer depth as defined in the KPP parametrization paper by Large et al 94. Basically some function of the surface forcing and density and momentum profile if I'm not mistaken. You can check that by comparing Hsbl with some MLD density difference cri...

From the code in ana_nudgecoef.h if you want to have just a relaxation boundary condition you must activate the option for climatological nudging (grid interior) such as TCLM_NUDGING, even if you don't want a widespread climatological nudging. So if I want to have a fairly strong nudging at the boun...

Did you try and reduce the hyperviscosity coefficient? the rotation of biharmonic tensor induces large vertical fluxes that violate stability constraints (see Marchesiello et al. 2009). Clipping of hyperviscosity is needed in this case ... Thanks for the reply. I'm using for tracer advection the TS...

Dear all, I want to impose a sponge layer near my North and South wall boundaries. I have a resolution of dy ~= 4km, I decided that the layer will have 160km, so its size is 70 eta points. My VISC4 = 6d9 m^4s-1 (not even so sure how large or small this is, it was on some previous namelist) No grid d...

Hi, You can try CCMP: http://podaac.jpl.nasa.gov/DATA_CATALOG/ccmpinfo.html Its not widely used as of yet. Not sure about quality, but its a blend of several satellite products and ECMWF forecasts. My comparisons to buoy are very favorable though. Note this is an analysis product. Also the time cove...

Dear forum, Can I use the bulk algorithm just to compute my momentum stresses and then set all the other fluxes to zero? I assume the COARE 3.0 bulk algo will only work with a specific meteorological input. Say I have all kinds of meteo data from a reanalysis for input to ROMS. But I only care about...