[gradsusr] Question on computing Kinetic Energy

[Saiprasanth Bhalachandran]

Hello,
My name is Sai and I am trying to compute the Kinetic energy of the upper-level divergent winds in a tropical cyclone. However, I am unsure if I am doing the integration correctly. Can one of you please correct me if I am wrong and tell me how to go about this?

As you know, the formula for Kinetic energy is :
[cid:4A2F25CE-A70F-4E04-94B8-90D9C61AB884 at itap.purdue.edu]

My domain is a 4 degree x 4 degree box and here is my grads script. I am particularly not confident about the integrations. Can one of you please give it a look at tell me?
Thank you - Sai

"reinit"
"open basic_std.ctl"
"set grid    off     5     1"
"set xlopts 1 12 0.2"
"set ylopts 1 15 0.25"
"set clopts 1 5 0.15"
"set font 5"
"set display color white"
'clear'
"set lat -2 2"
"set lon -2 2"

time=1
while(time<61)
"set t "time
"set t 24"
"set lev 300"
"den = (lev*100)/(287*(1+0.608*qq)*tt)"
"speed = mag(uu,vv)"

*Note: The area is 4degrees by 4degrees. 111.325 * 4 * 1000 is the dx and dy

"KElat = gint(0.5 *den*(speed*speed)*(445300),lat=-2,lat=2)"
"KE    = gint(KElat*(445300),lon=-2,lon=2)"


time=time+1
endwhile

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20171214/e78d3a54/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2017-12-14 at 9.34.31 AM.png
Type: image/png
Size: 12271 bytes
Desc: Screen Shot 2017-12-14 at 9.34.31 AM.png
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20171214/e78d3a54/attachment-0001.png