get the potential temperature at the surface
get the potential temperature at the surface
Since vertically the potential temperature is on the "RHO" levels, for which the Nth level is not exactly the ocean surface. How to get the potential temperature at the ocean's very surface ? Can I set up the vertical coordinate so that the top level z_r=0 too (like z_w) ?
Short answer - no. You'd break the whole finite volume formulation for ROMS, plus you'd have to change the vertical operations such as diffusivity. If you truly need a surface temperature and you think there is some sort of thin boundary layer there, it will take some special coding. The ice model does this (and I don't really understand it).