Ocean Modeling Discussion

I think you should review the straits that communicate with the Sulu Sea. At 1/12 resolution and rx0 max of 0.2 I suspect you are strongly reducing the sectional area of the straits (and hence the tidal transports). The DT you are using suggests to me that you have an extremely smooth bathymetry ...

Jayasankar, I assume that you are also providing the other variables in the IC file (velocities and sea surface elevation). If you are taking your IC from another model they should be mostly balanced unless the geometry of your ROMS application is very different to that of the other model. You ...

Yes, ROMS SSH is the sum of the MDT, the SLA and tides (if you have them defined in ROMS).
Be careful to add the model MDT and not other product such as that provided by AVISO.
The model MDT can be obtained as the mean of a long (say 5 years or more) forward run.

Jvr

Hanq, As Kate mention you would need to tailor ana_nudgcoef.h for your specific needs. However, if you are using a recent version of ROMS you can use a nudge coefficient netcdf input file, which also allows you to define a depth dependence of the relaxation time. The details are here: https://www ...

Hi Yong Chen, We have a prototype model for the area that we have been fine-tuning and is now performing reasonably well. The figure attached shows the temperature profile close to your station 30 and it compares well with your HYCOM plot. Some specific aspects of our configuration relevant to your ...

I recently did a similar evaluation for both ECCO2 (1/4 deg @ 1 day) and SODA analysis (1/4 deg @ 5 days) for the Mid Atlantic Bight (North east coast of US). When compared with observed SST and SSH both reanalysis simulate adequately the seasonal cycle but don’t do so well in the mesoscale. I was r ...

As far as I can remember the appropriate boundary condition for tides is

u - p/c = Flux (= A*cos(w*t) in your case)

with u = ubar, p = g*zeta, c = sqrt(gh)

so you need to define ubar in ana_m2obc accordingly.

Jvr

You could also compute the wind from the stress and then use BULK_FLUXES as it is. If the you have the pseudo-stress Tx = u*W and Ty = v*W where W = sqrt(u^2 + v^2) Then it follows that W = (Tx^2 + Ty^2)^(1/4), and u = Tx/W v = Ty/W If you have stress, then you also need to divide by rho_air*Cd Jvr