the depth of rho points for ROMS

[peigen]

Hi all,

I calculated the depth of rho point for ROMS using "z=zlevs(h,zeta,theta_s,theta_b,hc,N,vtype);"
But the result is not reasonable (shown as attached figure) that the numbers do not keep increasing from top (20th row) to bottom (1st row).
Could you please let me know why that happen?
Thank you!


Best regards,
Peigen

[arango]

The function zlevs is not part of our official distribution code and we are not responsible for third party functions coded by others.

In the ROMS Matlab repository we offer: set_depth.m. This is the official Matlab version of such a function for computing depths at any C-grid location, if you use the version of the code that we distribute.

[m.hadfield]

What Hernan said, but do remember that

• ROMS levels are numbered from the bottom to the top; and
• The ROMS s coordinate is -1 at z=-h and 0 at z-zeta.

I periodically confuse myself by forgetting that s is negative!