Ocean Modeling Discussion

I'm recently studying the kshedstrom version <https://github.com/kshedstrom/roms/tree/master>. The model works well with the <define UV_VIS4, define TS_DIF4>. However, if I activate Smagorinsky-like diffusion <undef UV_VIS4, undef TS_DIF4, define UV_VIS2, define TS_DIF2, define UV_SMAGORINSKY, defin...

Compared with POM, ROMS is more applicable for coupling with the other models <WRF, CICE, SWAN>, which is more suitbale for my future research on the lake-air interactions. For now, both POM and ROMS can well capture the lake temperature, while their performances in simulating the free-surface eleva...

The red and wet contours in the right panel both denote the free-surface variations <zeta> in the whole simulation area. The following is my model domain and Lake Nam Co depth. model domain and lake depth.png In winter, the currents mainly flow in response to the strong surface winds, which mainly d...

I'm recently studying the kshedstrom version <https://github.com/kshedstrom/roms/tree/master>, and trying to apply it on a freshwater lake. I plot the evolutions of the 5-day averaged depth-averaged currents <ubar, vbar, vector>, free surface elevation <zeta, contour>, and water temperature along a ...

Thanks, Kate! (1) The model works if I reduce the timestep, 10s for the internal mode and 2s for the external mode. (2) I activate the "UV_LOGDRAG" to use the logarithmic bottom stress in <set_vbc.F>. The maximum and minimum value for drag coefficient <Cdb> are 0.5 and 0.000001, just the d...

I'm recently studying the kshedstrom version <https://github.com/kshedstrom/roms/tree/master>, and trying to apply it on a freshwater lake. The model works with the Cp=3985Joules/kg/degC, which denotes the specific heat for seawater in Subrouting <Modules/mod_scalars.F>. However, if I increase Cp to...

Thank for for the suggestion. The model works if I reduce the dt from 300s to 10s.

I'm recently studying the kshedstrom version <https://github.com/kshedstrom/roms/tree/master>, and trying to set the vertical s-coordinate manually. However, if I set the s-coordinate with an intersified surface layer, the model blow up with the 'REASON: KEchar = NaN, PEchar = NaN' report. The follo...

Great thank! The problem has bees solved!

I'm recently studying the kshedstrom version <https://github.com/kshedstrom/roms/tree/master>, and the model blow up at the first timestep. However, the same input 'grid.nc, initial conditions, and atmospheric forcings' can work with the roms-3.7 version. The following is the running log <intelmpi_3...