Ocean Modeling Discussion

The reason may be the Coriolis parameter (f(i,j-1)+f(i,j)) in tl_balance.F90 and ad_balance.F90 becomes to or near to zero in the equator, which may cause unstable results. code in tl_balance.F90, which using (f(i,j-1)+f(i,j)) : Compute balanced, surface U-momentum from baroclinic and barotropic ...

Did you smooth your bathymetry? Thanks LiuZHQ! I smoothed the bathymetry with rfac=0.35(maximum slope of topography allowed in smoothing), and I have a test with the bathymetry to be a constant 800m, it is no useful. I changed the initial state from levitas climate data to hycom exp08, but still no ...

Thanks Kate! I have change the time step from 720s to 90s, and the Maximum grid stiffness ratios: rx0 = 3.069334E-01 (Beckmann and Haidvogel) rx1 = 9.416266E+00 (Haney). The SST and speed at the beginging is 屏幕快照 2012-08-24 下午10.58.48.png at time record 2 (the total time = 2*90s)is: 屏幕快照 2012-08-24 下 ...

I have a succesful modelling of an area using a grid of 116x156x32, but when I changes the grid size of the same area to 0378x0358x032, It blows up just at the beginging, Could any one help me? the following is the output the model ...

Hi leonjld, In my test, I get the error: >> seagrid2roms('seagrid.mat', 'test.grd'); ## SeaGrid Source File : seagrid.mat ## ROMS Destination File: test.grd ## Defining Global Attributes... ??? Undefined function or method 'endef' for input arguments of type 'double'. Error in ==> ncvar.subsref at ...