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

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