[ncl-talk] Error in computing heat sources (Q1) and sinks (Q2)

[Lyndz]

Dear Sir Dennis,

Thank you for this!

I figured out the mistake in my code.

[1] The unit od dTdt is actually K/day not K/s.

[2] This line:
p2 = conform(H,p,1)
p2 = p2*100
dHdp = center_finite_diff_n(H,p2,False,0,*1*)  ;dim should be 1

[3] Pressure should be in Pascals for the integration

After integration, the output will be in W/m2 (Q1) or [kg/kg](W/m2)
Attached is the working code.


Sincerely,
Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190913/275bc875/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: heatsource.ncl
Type: application/octet-stream
Size: 5366 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190913/275bc875/attachment.obj>