About the hc ,Tcline and hmin

[liangliang]

everyone:
I know 'hc' is citical depth ,and 'Tcline' is surface/bottom layer width,hc=min(hmin,Tcline)
But in my model ,the water is very shallow ,so there are many negative values in the 'h'.
I want to ask some qustions:
If the min 'h' is negative, is the 'depthmin' set to negative value?
And can the 'hc' be set to a negative value?
Thanks

[tony1230]

And can the 'hc' be set to a negative value?

hc can be any positive value , but negative ones!

[kate]

The newer vertical parameters for Vstretching, Vtransform not equal to one are designed for your case and allow hc to be bigger than hmin.

[staalstrom]

Hello

I have more questions conserning this subject.
I want to set hc>hmin. hmin is now 2m and I will now try hc=20.
I use Vtransform=2 and Vstreching=1. I make initial condition by using hc=20
and set Tcline=20 and hc=20 in this file.
In ocean.in I set Tcline=20.
But when I run my application I get an error that say that hc is set to 2m.

Why is hc set to the shallowest depth when Vtransform=2?

brds
Andre Staalstrom

[staalstrom]

Hello again

And one more question:
Is it necessary to recompile the code to use Vtransform=2?
Or is it enough to set it in ocean.in?

I use roms-3.2 by the way.

brds
Staalstrom

[kate]

You shouldn't need to recompile, but you do need to have a reasonably current version of set_depth.F. You might want to compare yours to that in the latest code. Also, there were versions with the new set_depth, but with historic code that wasn't clued in to not reset hc. I'm not sure which routine that was happening in.