Ocean Modeling Discussion

Hi, I wanted to look at the vertical structure of stresses, and so I computed rho*AKv*du/dz (or dv/dz). The vertical gradient is computed using the finite difference method, not splines. I compared the values at the surface and bottom boundaries to sustr/bustr, which are output into history files ...

How about just calculate the barotropic PGF from the standard zeta output? I think it is better than changing the ROMS codes. you can get baroclinic PGF by subtracting barotropic PGF from total PGF. Thank you for your reply. Actually I did the opposite recently. I computed the baroclinic pressure ...

Hi, Is there any way to separate the barotropic and baroclinic pressure gradient in the 2D diagnostic output? There is only one term (ubar_prsgrd/vbar_prsgrd) in the diagnostic output. I learned that this term contains both the barotropic and baroclinic parts. I found where the baroclinic pressure ...

Hi Dr. Wilkins, yes your suggestion works! This time I computed 'ubar_sustr+ubar_bustr+ubar_prsgrd+ubar_cor+ubar_hadv+ubar_hvisc-ubar_accel' and the result is zero. Thank you!

Hi, I'm trying to do some momentum balance analysis using ROMS. I output the 2D momentum terms into the diagnostic output files, but I'm not quite sure how the momentum balance is completed. I chose a point in the model, and added up all the terms (ubar_sustr+ubar_bustr+ubar_prsgrd+ubar_cor+ubar ...

Hi, I'm wondering if anyone has a script for computing the 2-D momentum equation terms for ROMS. I know this can be realized by turning on the diagnostic output option, but my model run really takes time and I don't want to repeat it to just produce the diagnostic files. So a well-written code for ...

Thanks arango, and then I can try with the volcon.

Hi rduran,

I did see some discussions on VOLCON before. But i remember that somebody mentioned using VOLCON would increase the computational time by several times. This could be an alternative if there are no other solutions. Thank you for your hint!

Hi Simon, Thank you for your suggestion for using the radiation boundary for M2. Actually when I increased the integration time, the surface elevation began to stabilize in the 2nd year, although the values are still high compared to observations (but the variability matches the obs well).So I guess ...

Hi all, I'm running a model for the Gulf of Mexico. Now I'm trying to set up the open boundary conditions. The model has two open boundaries, the southern boundary opened at Yucatan and the eastern boundary opened at Florida strait. Since now I'm mainly focusing on the circulation of northwestern ...

Hi Kate, Arango and Shchepet, Thanks for your replies! Previously the model was working fine, and after I changed the values of some parameters in the ocean.in file, such as NTIMES, NHIS, and DSTART, that error occurred. I'm using gfortran for the compilation and openmpi to do the parallel computing ...

Hi, I tried to submit a job and found the following error in the .err file. At line 2033 of file inp_par.f90 Fortran runtime error: Bad real number in item 1 of list input I looked at line 2033 of the inp_par.f90 file and the statements are as follows: 2031 string=Vstring(is:ie) 2032 LenS=LEN_TRIM ...

Thanks, mathiu, and I'm going to doublecheck it.

The river inflow in ROMS is quite complicated. In the netcdf file of rivers you have the following arrays: double river_Xposition(river) ; river_Xposition:long_name = "river ETA-position at RHO-points" ; river_Xposition:units = "nondimensional" ; river_Xposition:field = "river_Xposition, scalar ...