Perhaps there is an easy answer to this, but I seem to be struggling to find it on the forum. Is there a way to specify the vertical profile of a tracer as a function of depth, rather than by layer? Currently, I'm specifying temperature on specific layers in TRA_PRO_Z.F using something like this:
IF(IOPT.EQ.0) THEN
IF(ITRAC.EQ.1 .OR. .NOT.SEDI) THEN
IF(IPLAN>=0.5*NPLAN) THEN
TRA_PROF_Z = 2.D0
ELSE
TRA_PROF_Z = 1.D0
ENDIF
ENDIF...
It works just fine to generate a two layer stratification where the thermocline essentially follows the terrain. However, I'd rather have a more realistic condition where the thermocline is at a constant depth. I'd like to keep the sigma grid if possible. Any thoughts? Thanks in advance!