Thanks for your replay;
Some times I look part of the code and the solution looks simple, but I'm sure the Kernel is a lot more complex than it

. That's ok.. I will try another solution

.
Now let me explain my idea:
I'm modeling an inverse estuary (Spencer Gulf, South Australia), where we have a very strong density gradient during the summer and winter. Basically, we have very dense water in the head of the gulf and lighter water in the ocean and it produce strong currents.
I used summer temperature and salinity distribution with TS_FIXED to get the associated sea level and currents. It works fine. I'm happy with the results and the input that it gave to my work. After that I was trying to use some tracers to explore some ideas. But that's fine, I don't want to create problems.
The diagnostic run was just the first step. Now I will run the model with the proper, wind stress, heat and water flux and I hope to get the summer scenario. After that I will be able to explore my ideas with the passive tracer option.
Again, many thanks for your explanation.
Carlos Teixeira