Question about turbulent kinetic energy

[nmichelet]

Hello everyone,
I would like to know if it is possible to define a turbulent kinetic energy profile as an input into ROMS.
Thanks

[kate]

To what purpose? Would this be part of a vertical mixing scheme? It is certainly possible to modify ROMS, telling it to read new fields, but then what? Is it an initial condition? Something that needs updating from outside? Something that needs to be advected?

[nmichelet]

Thank you Kate for your answer.

The purpose of this is to define in a channel the initial condition for the velocity and the turbulent kinetic energy. I have some data at the enter of my channel so I would like to put those as input.

[kate]

You should look at get_state.F to see how and when it reads TKE already. It is protected within:

Code: Select all

#  if defined GLS_MIXING || defined MY25_MIXING
#endif

and also:

Code: Select all

IF (get_var(idMtke).and.Perfect3D) THEN
END IF

Clearly, for a new run, it is getting TKE some other way (initialized to zero??). Just do some snooping to figure out what to change.