I set theta s=5 and theta b=0.4 , these parameters should work on the spacing of near-surface.
But the value of s_rho shows that the grid is still equidistant.
I wonder how to solve this problems, look forward to your reply.
Parameter Theta s and theta b don't work
Re: Parameter Theta s and theta b don't work
s_rho is not the final dimensional vertical coordinate. s_rho depends on the choice of TRANSFORM and the final z depths depend on the STRETCHING applied to s_rho.
Read https://www.myroms.org/wiki/Vertical_S-coordinate
There are Matlab and Python routines to compute z offline from the Vtransform, Vstretching choices and the parameters theta* etc.
Setting Hout(idpthR) == T in ocean.in will also write the time varying (because of zeta(t)) z coordinates to the history file.
Read https://www.myroms.org/wiki/Vertical_S-coordinate
There are Matlab and Python routines to compute z offline from the Vtransform, Vstretching choices and the parameters theta* etc.
Setting Hout(idpthR) == T in ocean.in will also write the time varying (because of zeta(t)) z coordinates to the history file.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Re: Parameter Theta s and theta b don't work
wilkin wrote:s_rho is not the final dimensional vertical coordinate. s_rho depends on the choice of TRANSFORM and the final z depths depend on the STRETCHING applied to s_rho.
Read https://www.myroms.org/wiki/Vertical_S-coordinate
There are Matlab and Python routines to compute z offline from the Vtransform, Vstretching choices and the parameters theta* etc.
Setting Hout(idpthR) == T in ocean.in will also write the time varying (because of zeta(t)) z coordinates to the history file.
Thank you, your answer really helps a looot!