Ocean Modeling Discussion

Thanks for your reply. I looked at the fields of temp, u, v, ssh etc (not AKv yet) and the pattern is similar to what was seen in the TKE: 5 simulations with three slightly different results, and two being reproducible. I didn't see strong differences at the boundaries, as well as near the rivers. R...

Hi all, A free run (Exp1 - red line) was generated in Sep 2021. In March 2022, I ran Exp2 (blue line) turning on passive tracers (all needed cppdefs), and adding dye_01 to ini, bnd and river files. The experiments generated different results (fig attached - time series of total KE in cm^2/s^2). I th...

Thanks for your message.

In which field of NL_fDataPenalty should I get the final NL model data penalty? > ncdump -v NL_fDataPenalty ocean_mod_0006.nc NL_fDataPenalty = 785.246249424019, 591.14144905359, 0, 0, 0, 0, 188.059436356603, 6.0453640138258, 0, 0, 0, 0 ; Whereas, the total cost function shows: > ncdump -v Jact ocea...

Hello,

Do in situ u and v components of velocity data need to be rotated to grid angle before assimilation?

Thanks

Rafael

Thanks Joao and John for your replies. The problem was in obs_Zgrid. I've changed the depth of temp from 0.02 to 0.5 m in the scripts I'm using and I got better values for obs_Zgrid (below). Before: obs_Zgrid = -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -0, -...

I'm also having SST gridded satellite data getting rejected. About 195 out of 204 SST data points are being rejected. And only 4 out of 204 SSH data points were rejected. I'm setting the SST (temp) depth to 2 cm and the minimum error to 0.2oC. Please let me know if you have any ideas. Fragments of t...

Hi, I'm running simulations for the northeast region of New Zealand using JRA55-do (https://climate.mri-jma.go.jp/~htsujino/jra55do.html), which has corrected atmospheric variables for ocean forcing. The dataset does not have upward shortwave and longwave fluxes. Therefore, I'm using the option LONG...

Thanks for your reply John. I checked both your routines: roms_{get,write}_era5_NCARds633_frcfile.m. They are quite nice. Though, we have annual files of JRA55-do already download on the server. Then I decided to modify d_ecmwf2roms.m to read and write JRA55-do variables for bulk_fluxes option. The ...

Hello everyone,

Is there a script in matlab or python to write ROMS forcing files using JRA55-do dataset?

Thanks in advance,
Rafael

Thanks Kate, Actually, at the beginning I used residual (or monthly averages) of vbar_northward (at rho points)*dx*(h+zeta residual). Then I summed the residual transport through the section, positive transport plus negative transport to obtain the residual transport. I was wondering if the eulerian...

Dear users, I am running a seasonal ROMS simulation with tides and wind for the Todos os Santos bay, located northeastern Brazil. At the beginning, I started calculating the residual transport by multiplying the monthly (average of 29 tidal cycles) residual velocity and the area. The area is equal t...