[MITgcm-support] Heat Budget in tutorial (Global Ocean Simulation at 4 deg)

[Krishnakumar Rajagopalan]

Hi All,
 
It seems that the heat globally added at the ocean surface does not match the 
heat absorbed by the ocean.
 
Let Ho and Ht denote the initial  enthalpy and enthalpy at time, t,  of the 
ocean.  Then (Ht – Ho) represents the change in enthalpy , del H, of the ocean 
which is calculated from  summing  over all the cells, mass_cell*Cp*( T(t) - 
T_initial  ),  where T_initial  and T(t) represent the initial temperature and 
temperature at time, t, of the cell.  The quantity, del H,  at any time , must 
 equal the heat supplied to the ocean till that time,  which in the case of this 
tutorial comes from the restoring term { -lamba*(theta-SST) }  and Qnet.   I 
have integrated the diagnostic,  surfForcT (model surface forcing for 
Temperature), over the area of the ocean and with time to calculate the heat, Q 
entering the ocean. (The surfForcT is output every month). In (attached) figure, 
the red line shows del  H as a function of time and the blue line shows Q 
supplied to the ocean  upto that time, t .  Ideally the two lines must 
coincide.  What could be the reason for this apparent discrepency which actually 
increases with time ?
Thanks in advance
Krishnakumar
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20110516/a19f22ce/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dH_and_Q.jpg
Type: image/jpeg
Size: 32761 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20110516/a19f22ce/attachment-0001.jpg>