Ocean Modeling Discussion

Dear community,

Could anyone direct me to a python script/package/project to write netcdf boundary files?
I have a bunch of depth-averaged velocity and sea surface elevation timeseries with which I would like to force a 3D ROMS simulation.

Thanks in advance,

-T-

Dear Community, Here is the thing. I would like to generate eddy/vortex-like structures which would travel in a simplified one-direction open-boundary flow. Firstly, does anyone have any idea if it is even possible to generate/impose such structures? Secondly, if it is indeed possible, do you have ...

Thank you all,

Here is an early christmas present, a python translation of John's tpxo8_to_roms package.
I've started testing it...so far, so good.

Enjoy

-T-

Dear all, I realize this thread is rather old but I noticed this odd line in TPXO2ROMS_v4pt1.m: line 253: ncwrite(fnOut,'tide_Vphase',vamp) Should not it be "ncwrite(fnOut,'tide_Vphase',vpha)" rather? BTW, if some could direct me to a python package permitting similar functionality (i.e. creating ...

Dear community, I would like borrow your brains for a moment. Although I am guessing that this is out of ROMS' capability scope, I was wondering if there would be a way, given a small-enough time step and a large-enough number of iterations, to impose a turbulence spectrum (i.e. Von-Karman like ...

Greetings everyone, I was wondering if anyone has ever used CoastGrid to generate grid files (i.e. *.nc) !? I find CoastGrid interface pretty slick and straight forward and easy to use. it is fairly simple to generate the *.dat grid files yet I miss one major thing: how do you convert those files in ...

Unfortunately, it is a bit more complex than that. During my PhD, I have developed a methodology for representing hydrodynamic effects induced by individual tidal current turbine within ROMS (i.e. http://mhk.pnnl.gov/wiki/images/1/1b/Roc_et_al._2013.pdf). Now I am investigating the hydrodynamics ...

Hi ROMS community, I would like to create something like a 3D mask which I could use as an input file in order to turn on and off my custom momentum and turbulence source terms. What is the best approach to achieve such thing? FYI, I have successfully implemented these custom source terms yet I use ...

Yeah I spotted this one as well however the problem doesn't come from that...but thank you anyway. I have actually cracked my problem though. In my input file (i.e. *.in) the open-boundary conditions for temperature and salinity were set to "Clamped" which imply to provide some values, analytic or ...

Hi everybody, I am trying to run a very simple case (a simple rectangular channel with a constant inflow). My simulation is compiling without any trouble but when I am running the executable this message appears: "READ_PHYPAR - could not find input file: ocean_frc.nc" I guess this is due to a ...

Thank you very much Kate !

-T-

Hi Guys, I face troubles restarting a tidal flow simulation from an *_his.nc file. The model restart perfectly using the *_his.nc but the tidal forcing doesn't follow and damp to zero after the first steps. Notice that the tidal forcing is specified through the ana_fobc file as follow: fac=0.75_r8 ...

Hi everyone,

Is there a flag or an option allowing to choose between x-y and/or lon-lat coordinates for the output files, particularly the history files ?

Cheers,


-T-

Hi folks,

I was wondering if there is any usual rule thumb (or related references) to parameterize a priori the set of background vertical mixing coefficients ( AKT_BAK, AKV_BAK, AKK_BAK and AKP_BAK) ?

Cheers,

-T-

Thank you very much for the clarification.

-T-

Hello everyone,

Quick question, does anyone know what are the units of gls in the output files ? m3/s2 as defined in varinfo.dat or the units of the length scale corresponding to the chosen two-equation turbulence model (e.g.: s-1 for k-omega, m2/s3 for k-epsilon, m3/s2 for k-kl,...) ?

Cheers


-T-

That's what I thought as well !

Thank you so much for the help

Absolutely spot-on, I'd like to add a source term on the RHS of the previously mention equation number 13, function of the squared production of turbulent kinetic energy.
So how do I proceed ?

Cheers

Hi everyone, I would like to add a turbulence dissipation rate production term into the GLS turbulence closure and more precisely in the GLS equation. Physically speaking it would be a quadratic function of the production of turbulence kinetic enegy. So far I am struggling with numerical and ...

Spot on John !

Many thanks to you both.


Thomas

Hi Kate,

Sorry for the lateness of my answer.

I have just checked my analytical.f90 and yes indeed, my ana_initial module is coming from here.
By the way no trace of ana_NLMinitial_tile in analytical.f90. I have got instead ana_initial_tile.

Thanks a lot for your help.

Hi every one, I have just upgraded my ROMS code through SVN and rerun a case which was working perfectly before the upgrading and unfortunately this following error message now pops up during the compilation: analytical.f90(336): error #6645: The name of the module procedure conflicts with a name in ...

Hi Howard,

I'm struggling with the similar problem right now.
Have you managed to sort this out ?

Cheers

Hi everyone, It's probably a dummy question but can anyone explain me the physical meaning of "wrk" in rhs3d.F in the part reserved for the bodyforce expression. ! ! Apply bottom stress as a bodyforce: determine the thickness (m) ! of the bottom layer; then add in bottom stress as a bodyfoce. ! DO j ...

Thank you Kate. Indeed I am not using periodic boundary conditions but I am using DISTRIBUTE (running ROMS in parallel on 2 processors) and one can find in ANA_FLUX.h: #if defined EW_PERIODIC || defined NS_PERIODIC || defined DISTRIBUTE USE exchange_3d_mod, ONLY : exchange_r3d_tile #endif and ...

Hi everybody, I'm trying to run in parallel (mpi on 2 intel's cores duo) a 3 D simulation of a simple East-West channel including areas of different roughness height (Z0). To set up these areas I am using ANA_BMFLUX. Unfortunately this message pops up each time I try to compile the code: analytical ...

Hi everybody, I'd like to add a mid depth local pressure jump, defined as a function of the velocity, in my 3D simulation. As far as I've understood the best way to do this is use the "bodyforce" option and express directly my pressure jump expression into the code. My questions are: which files ...

Many thanks for the piece of advice, now it's compiling and running correctly. But a new error has popped up, from now I can use the "make clean" command neither recompiling the code. It can be read: " makefile:229: INCLUDING FILE Build/make_macros.mk WHICH CONTAINS APPLICATION-DEPENDENT MAKE ...

Hi everyone, I'm trying to compile ROMS on 4 processors with the command "make -j 4". My exploitation system is Fedora 9 on a 64 bits architecture and I'm using intel fortran compiler and intel MPI library. In the makefile I've only switched on the "USE_MPI" and "USE_LARGE" options and define the ...

My mistake ! My supervisor pointed out the problem in my simulation. I was calculating west_ubar as function of another unknown(zeta_west-zeta_east). Thus non linear behaviours pop up. By using a fixed (zeta_west-zeta_east) and using the proper time step, the simulation reaches a steady solution ...

Hi, I'm currently facing the same problems. I'm trying to run a barotropic simulation of a simple east-west channel with a homogeneous bathymetry. In this case zeta_west is gently increasing up to a threshold of 0.4m and zeta_east remains the same (thanks to ANA_FSOBC). The boundary conditions are ...

Yes it does !

Thanks a lot.

Hi everybody, I am trying to run a very simple case (a simple channel with a rdrg2 local variation). My simulation is compiling without any trouble but when I am running the executable this message appears: "READ_PHYPAR - could not find input file: ocean_ini.nc" I guess this is due to a consistency ...

Hi everybody, I am a new user of ROMS and I bring with me whole bunch of questions. To get started, I am trying to build up a simple case composed of a simple channel where an additional shear stress is imposed on a particular area in order to simulate a power extraction. This additional shear ...