Ocean Modeling Discussion

I hope to perform a hot start from a certain point before the model crash, but this brings up several questions.
In my restart file, I noticed there are multiple ocean_time entries, and each ocean_time appears to contain two sets of data. Does this mean I can restart the model directly from a single ...

My model has encountered a new problem—it blew up after running for over nine months.
(273,054,35) 1.360096E-02 1.134846E-02 1.208489E-01 3.810625E+00
493932 2022-10-13 20:10:00.00 2.261205E-02 1.616870E+04 1.616872E+04 1.026477E+16
(274,054,31) 1.138465E-02 1.370948E-02 1.102423E-01 3.910384E+00 ...

I suggest you try setting h = 10 where mask_rho = 0 and #undef WET_DRY. You have very small (10^-4) values of h in the land that might cause trouble.

I notice that your land/sea mask is such that the Malaysia Peninsula connects to the southern boundary. The Singapore Strait is closed.


Thank ...

Sadly, I tried modifying my grid file and the h file today, and also made corresponding changes to other related files. However, the error still hasn't gone away. That said, I believe the cause of the error lies within the grid. This is my grid file—hopefully the issue can be identified.

Thank you both very much for your kind help. I checked the grid file and, as expected, found the presence of zeros. I will proceed to modify the file. Thanks again!

Estimated maximum CFL values U∆t/∆x for u,v,w velocity are reported in the log as Cu,Cv,Cw. At the start of your run,
Cv = NaN so you have v velocity values that are NaN (in your initial conditions) which via continuity equation make w also NaN. You
cannot have NaN values in initial (or boundary ...

Hello everyone,

I’m encountering a confusing error. I initially built an inaccurate model, so I replaced the original grid, ini file, obc file, and tide file. However, I encountered the following error right at the first step of running the model:

Note: The following floating-point exceptions are ...

The following pictures are the ocean currents, salinity and temperature I simulated. The top two lines are my simulated monthly average data, and the bottom two lines are the satellite observation monthly average data provided by ORAS5
scs_SSUV.png
scs_SSS.png
scs_SST.png
As you can see, in ...

It seems that the reason is that I set NTSAVG to 0, after changing it to 1 the model ran smoothly

As a supplement, I did not extend the running time of the old model, but created a new one. The starting times of the two models are different.

Dear friends :) , thank you all very much for your help and answers before. My model ran successfully and gave good results. Now I want to increase the running time of the model from three months to one year.
But my model got off to a bad start, An error occurred before the his file was generated ...

Thanks to MoeinDst for the index and wilkin for your patient explanation, you really helped me a lot. :D

I chose to modify my tide input file to keep only the values for the boundary grid points and set all the data for the internal points to 0. Everything looked much better, except for the tide ...

Hi everyone, after some time, I'm back with a question.
The descriptions of SSH_TIDE and UV_TIDE say that tidal Forcing is added just on the open boundaries, but my tidal forcing file has tidal Data at every grid point, as does the data in \src\Data\ROMS\Forcing\test_head_frc.nc, which really ...

I have a question whether the fact that I used HYCOM's hourly reanalysis data to generate my boundary file, which includes UV, Ubar, and Vbar data, means that my model calculated two tidal forcing

I added tides to my model. The tide data comes from TPXO9_v5, and I used otps2frc(https://github.com/NakamuraTakashi/ROMS_prep/blob/master/OTPS/mfile/otps2frc_wrapper_tpxo9_atlas30.m) to generate the NETCDF file.

With the addition of tides, my currents are significantly skewed. Here are three ...

Dear jcwarner, thank you very much for your helpful explanation. However, this vortex is still bothering me, as it still exists when I set the OBCFAC to 10 or 100, albeit not as noticeably.

Thank you for your suggestion, JH. I will follow your advice and run the model. I sincerely appreciate the ...

I checked the previous experiment and found that after I changed the value of obcfac to 0, my model had an unexplained vortex in eastern Taiwan. In fact, in most of my other experiments, similar vortexes were present, with a strong current moving south along eastern Taiwan, while the northbound ...

Dear jcwarner, thank you very much for your advice. I will go back and try to modify the relevant parameters.
My initial file contains current, all from hycom

Many days passed, but I still couldn't solve the problem, and Kuroshio still didn't show up in my model.
I tried to modify the relevant parameters of obcfac and TNUDG, changing obcfac to 10 and the latter to 3 days, but the result was even worse, there was a strong southward flow field on the ...

After three months of operation, my model had this problem. As shown in the picture, it lacked the Kuroshio flow from Taiwan Island to the Philippines, which caused a large amount of warm current to flow into the southern waters of Indo-China Peninsula. I want to know how to solve this problem ...

I used the data from TPXO9 to generate a tidal forcing file. I defined the following cpp option and ran it for 90 days. There is a big deviation between the tidal level change of the site and the forecast data of TPXO9. Could someone please help me find out what the problem is? Appreciate any ...

Hello. I've never used this module before, so I don't know what kind of variables are saved as output. But if the output file gives information enough, you can apply a simple harmonic analysis. Some good packages in Matlab and python were developed to make harmonic analyses and even predictions ...

Hi everyone, I ran my model and compared the resulting temperature and salinity to the HYCOM data and found that the deviation was increasing over time. I don't know how to solve this problem :( , here are my cpp options, with my script file attached.

Appreciate any suggestions

#define PERFECT ...

Hello, everyone. I have a problem. I use AVERAGES DETIDE to obtain simulated tidal data. Now I generate roms_his and roms_har files, but I don't know how to obtain the amplitude and phase of each component tide . Is there any tool that can help me？
thanks for any suggest.