Ocean Modeling Discussion

Actually, I have a similar problem : my velocities, temperature (my ocean will boil soon !!!) and salinity behave abnormal near the boundaries (although the boundaries seem OK !?!) You wrote : The velocity seem to be go bad on the boundaries (from the bottom to the surface), I suspect my tidal ...

Have you checked the value of your CFL criteria? Minimum barotropic Courant Number = 1.17951147E-02 Maximum barotropic Courant Number = 4.07452685E-01 in the print you can read as your model is running (and that I put in a 'log_file'). If your Courant Number is too high, that means your grid can be ...

Thank you for your advice. :) Actually, I understood I had to use the advection scheme if I wanted the more efficient solution for configuration with steep bathymetry, following some advices from K. Hedstrom. 2016. Technical Manual for a Coupled Sea-Ice/Ocean Circulation Model (Version 4). U.S. Dept ...

Hi Kate, thank you for your answer. I've tried to activate LIMIT_BSTRESS as you suggested. Unfortunately, that does not work either. (It blows up with the same abnormal velocity values). I then tried # define LIMIT_VVISC #define RI_SPLINES I tried to change my drag scheme (with #define UV_LOGDRAG ...

Hi everyone, I keep working on my configuration and I'm facing for several days a new blowing-up event after four time step (one time step is 180 s). I look deeper into my output files (netcdf and log error) and I had some clues. But I'm afraid that I have not the keys to understand well. Actually ...

Hi lalvarez! I'm working with the HYCOM files too (but the GLBv0.08 analysis). I tried to download one file from the GLBu0.08 to test it but I got this message : Requested time range 1999-04-18T00:00:00Z - 1999-06-20T00:00:00Z does not intersect actual time range 2016-04-18T00:00:00Z - 2018-06-20T00 ...

Thank you for your advices. Actually, I had some problems with my atmospheric input files because of interpolation online. In fact, my input files were written with the latitude decreasing instead of increasing. So when ROMS was scanning (in subroutine hindices in interpolate.F) my latitude values ...

Hello modellers, I've made some modifications to ROMS take into account the "hour" units as time axis. Now, here is what I've added. In mod_scalars.F : I added line 705 real(r8), parameter :: day2sec = 86400.0_r8 real(r8), parameter :: hour2day = 1.0 / 24.0_r8 ! LINE 705 real(r8), parameter ...

Thank you for your advice! :) None of these files seems to be weird (No Nan, and values are coherent). But I noticed something concerning my wind input. On the netCDF in the initial file, all values seem OK. But as I make my ROMS model write the values of wind it reads in netCDF (adding WRITE ...

That does not matter the fact that the K from (Akt) is in capital letter in one line and not in the other?
(Anyway, I tested with Akt at the place of AKt, and that does not solve my problem... )
I keep on searching...

Thank you Kate for your answer. I tried to smooth (even exageratedly) but, as you expected, that does not solve anything. I tried to compute an other domain with a smoother bathymetry as joining the coast but that does not change anything neither. With a colleague, we have noticed that my ocean_his ...

Hi modellers, I'm trying to model the ocean dynamics in a “small” domain along the coast of West South Africa in the Benguela Current. I chose it a bit randomly as a first test domain as I want to study other domains afterwards. At the fist time step (t=1), ROMS blows up while indicating that : F ...

Silly mistake! I was right on that point!
Thank you very much
It works ... further now!

Hello modellers, I'm trying to run ROMS using the BULK FLUX Formulae to simulate my atmosphere. For that, I have build the netcdf files using the pyroms tools and inspiring from the MERRA example (providing with pyroms package). atm_forc_GFS_2018-01-01_Pair.nc atm_forc_GFS_2018-01-01_Tair.nc atm ...

Thanks a lot! Indeed, I did not notice the modification in the index as downloading the HYCOM files. Now, my ssh interpolation works correctly. And at the same time, that explains the error I received after the ssh interpolation, as computing the interpolation of the following variables. Traceback ...

Hello everyone, I try to interpolate some HYCOM output fields (temperature, salinity, ssh and wind velocity) to a ROMS grid thanks the pyroms tools. To achieve that, I try to reproduce the HYCOM_Arctic example from the examples directory. Yet, I succeed in interpolating the ssh from HYCOM GLBa0.08 ...

Thanks a lot!

Hello everyone, I'm trying to use the pyroms tools to create my own files for initializing, defining boundaries and forcing. As the files I want to use are from HYCOM outputs (from GLBv0.08), I first want to train on the Arctic_HYCOM example. To do it, I need a definition of the grid ARCTIC2 (that I ...