Ocean Modeling Discussion

Thanks a lot this is very helpfull

Hello Kate and thanks for your reply,

I have already tried to turn on SOLVE3D and it works well. My simulation run without problem including the nesting. But I want to run this simulation in 2D.
Any suggestion to solve the problem?

Thanks

Thank you for your answer. I have tried to remove the #ifdef SOLVE3D and the following #endif around the definition. It remove the error about ngc but the other errors remain nesting.f90(1576): error #6404: This name does not have a type, and must have an explicit type. [COUPLING] & COUPLING(dg ...

Hello everyone, I am currently trying to run a simulation with 2 grid in nesting in 2D. But I have a compilation error that I don't understand. Could you help me solve that problem please. nesting.f90(215): error #6404: This name does not have a type, and must have an explicit type. [NGC] ngc ...

Hello, I have a question about the stability criterion in ROMS. I have found that there is a CFL criterion which is calculate with 2D time step like : CFL = sqrt(gh)*Delta_t2d*sqrt(Delta_x^-2+Delta_y^-2) And this CFL has to be lower than 1/2. But is there a stability criterion which is calculated ...

Hello everyone,
I am curious about the vertical diffusion terms in k equation. I understand that there is no numerical scheme to solve those terms, is that correct?
If it is not, how can I change the order of this scheme?
Thanks
Nicolas

Thank you for your help Kate it is very useful

Hello everyone, I am currently trying to find where is the vertical boundary condition in the model, in order to modify the condition on u. Can you help me to find where the next equation is in the model Km(du/dz) = tau_x^b # if defined UV_LOGDRAG ! ! Set logarithmic bottom stress. ! DO j=JstrV-1 ...

Thank you Kate for your answer.

The purpose of this is to define in a channel the initial condition for the velocity and the turbulent kinetic energy. I have some data at the enter of my channel so I would like to put those as input.

Hello everyone,
I would like to know if it is possible to define a turbulent kinetic energy profile as an input into ROMS.
Thanks

Hello,
I have a problem with my simulation and I would like to know if someone have seen this before.
I am running a simulation only with tidal forcing and I am looking at the 3D current. The thing is in a direction the velocity look good, but in the other one it is underestimate a lot.

Thanks for your answer Kate but it does not work, I think the two dates are independent. Well to be sure if I put the right date, I want my simulation to begin the 15th of March 1993. So I calculate the number of days between the 1st of January 1900 (my reference time) and this one and put the ...

Hello, I am currently trying to make a simulation with only tidal forcing, but I encounter a problem to calibrate the SSH. It seem to dissipate to much at the end of it and I don't really know what I need to modify. If you can give me an advice it would be really helpfull. Thanks a lot https://www ...

Thanks Kate

Hello everyone, I would like to know if it is possible to deactivate the temperature and salinity calculation. I tried to set it to a constant value but it start to vary again in my simulation, and it make the model blow up. Can you give me an advice to do that, I search but I found nothing about it ...

Thank you Kate I will take a look at that

Hello everyone, I have several questions about coarse2fine, the grid refinement routine. Concerning the bathymetry how does it interpolate from the coarse one to the fine grid? I also have a pretty good resolution on my bathymetry and I would like my grids to interpolate directly from it and not ...

Thanks for your advices, I will work on this and try to find the solution

I have this line in my ocean.in but no file named ocean_itl.nc, what kind of file is this?
I also have tried to remove the comment line but same thing, no change

Hi, I have checked and the -heap-array option is already activated I have looked at the lines 1390 and 87 in the inp_par.f90, but I don't really understand what it means line 1390 : Npts=load_s1d(Nval, Cval, line, label, igrid, Nfiles, ITL) line 87 : CALL read_PhyPar (model, inp, out, Lwrite ...

Deer Kate, I read the entire post but it is not exactly the same problem. I tried the USE_DEBUG and here the information it gives me Model Input Parameters: ROMS/TOMS version 3.5 Wednesday - December 2, 2015 - 3:32:04 PM ...

Hi, I have the same kind of problem, could you explain me how you solve it. I am currently trying to make a two-way nest simulation, so I followed the example of the Lake Jersey AB example. The program don't even start and it gives me the log file I attached. Do I miss something in the example ...

hi, I am looking for a solution for my problem with coarse2fine, I guess it is not a big one but I don't find what to do. Here is the error it gives me F = coarse2fine ('~/roms/Projects/Iroise_sea/agrif_parent.nc','agrif_test.nc',5,20,40,20,40) Interpolating from coarse to fine ... Computing grid ...

Hi, I already use the tide_ramp option, so I don't think this is the problem. I also try do run it with a little time step, it takes a long time, but it still blow up at the same time. I think there is a problem with the Wet_Dry option, because I have some velocity which blows up near the coast ...

Hi,
I have the same problem, did you find how to fix it?

I attached my configuration, if someone can help me with my blowing up problem.
Thanks

Hi,
I still have some trouble with that, the thing is the running is fine, the results are good, and the velocity suddenly blow up
Any idea why it does that?

Hello, My model is blowing and I don't know why. I tried to change the time step but the only thing it does is to delay the blowing up. I need help on this one Thanks 517 0 08:37:00 7.672425E-03 5.298350E+02 5.298427E+02 1.803663E+12 (071,101,01) 2.688424E-02 1.100356E-02 3.424691E+00 1.221882E+00 ...

Thanks for your reply, I solved the problem.
It was a problem with the boundary condition and the cppkeys I used
Thanks for your help

I am using the 3.7 version of ROMS. I change what you said but it still gives me the same results, here is the .h file Can you tell me if something is wrong? /* ** svn $Id: latte_c.h 342 2007-06-06 18:50:48Z zhang ...

Hello everyone, I am currently trying to run a model with only tidal forcing. The problem is that my results for the free surface or u and v are equal to zero everywhere. Can you help me please, here is my .in file and my .h file Thanks .in file ! Application title. TITLE = Iroise Sea Tidal ...