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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...

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 ...