gradient computation

Dear NCL users,

I have a temperature field that is not equally spaced in y, but it is in x.
I'm not sure how to compute its gradient in y direction. I'm wondering if I
can apply this

dphi = center_finite_diff( lat, 1., False, 0 ) * 6.37e6 * acos( -1.
)/180.
dT_dy = center_finite_diff( T( lon|:, lat|: ), dphi, False, 0 )

or I just use

dT_dy = center_finite_diff( T( lon|:, lat|: ), lat * acos(-1.)/180., False,
0 ) / 6.37e6

Thanks.

Best regards,

-- 
Mateus da Silva Teixeira
Registered Linux User #466740

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jul 06 2009 - 16:59:37 MDT