oceanic heat budget diagnostics

[sang-ki_lee@noaa_gov]

Hi,

I am trying to diagnose the heat budget of upper 5 layers.
It is not clear how I can get the heat tendencies due to vertical
diffusion, convection and horizontal diffusion. Has anyone
done this?

Thank you.

Sang-Ki

 

[sang-ki_lee@noaa_gov]

All heat budget terms are computed in CCSM3. You just need to modify
the subroutine "tracer_update(k)" (baroclinic.F) to print them out.

For instance, horizontal diffusion term is computed in line 897:

call hdifft(k,WORKN)
FT = FT + WORKN

Here WORKN(:,:,1) will be heat diffusion term (K/sec). You need to
save it. An easy way is to use one of the tavg routines in tavg.F.
Of course, you need to add the new variables in "gx1v3_tavg_contents"

It will be good to know how NCAR people compute heat budget terms.

Sang-Ki