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
This archive was generated by hypermail 2.2.0 : Tue Jul 07 2009 - 11:13:18 MDT