Main menu

Navigation

oceanic heat budget diagnostics

2 posts / 0 new
Last post
sang-ki.lee@...
oceanic heat budget diagnostics

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

Dr. Sang-Ki Lee
NOAA/AOML
4301 Rickenbacker Causeway Miami,
FL 33149 USA
Tel) 305-361-4521
Fax) 305-361-4412
E-mail) eval(unescape('%64%6f%63%75%6d%65%6e%74%2e%77%72%69%74%65%28%27%3c%61%20%68%72%65%66%3d%22%6d%61%69%6c%74%6f%3a%73%61%6e%67%2d%6b%69%2e%6c%65%65%40%6e%6f%61%61%2e%67%6f%76%22%3e%73%61%6e%67%2d%6b%69%2e%6c%65%65%40%6e%6f%61%61%2e%67%6f%76%3c%2f%61%3e%27%29%3b'))

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

Dr. Sang-Ki Lee
NOAA/AOML
4301 Rickenbacker Causeway Miami,
FL 33149 USA
Tel) 305-361-4521
Fax) 305-361-4412
E-mail) eval(unescape('%64%6f%63%75%6d%65%6e%74%2e%77%72%69%74%65%28%27%3c%61%20%68%72%65%66%3d%22%6d%61%69%6c%74%6f%3a%73%61%6e%67%2d%6b%69%2e%6c%65%65%40%6e%6f%61%61%2e%67%6f%76%22%3e%73%61%6e%67%2d%6b%69%2e%6c%65%65%40%6e%6f%61%61%2e%67%6f%76%3c%2f%61%3e%27%29%3b'))

Log in or register to post comments

Who's new

  • robert.grumbine@...
  • lwheele@...
  • junya.wu@...
  • afroberts@...
  • richard.allard@...