Spherical harmonic routines problem

From: Bao, Qing (baoqing AT mail.iap.ac.cn)
Date: Sat May 28 2005 - 08:52:37 MDT


Hi, all!

Does anyone use NCL functions rhomb_trunC or tri_trunc, maybe there are some problems with the functions. The following is my NCL code to deal with array tmp(108,128):

begin
  temp = fspan(19,1000,128*108)
  tmp = onedtond(temp,(/108,128/))
    ab = shagC(tmp)
    abc = rhomb_trunC (ab,42)
; print(ab)
    bb = ab - abc
    bb = fabs(bb)
    print(sum(bb) )
end

The arrays 'ab' and 'abc' are same, the result 'bb' is 0.

My NCL version is 4.2.0.a032, Linux

Best regards!

Qing

**************************************
* Bao, Qing *
* Ph.D. Candidate *
* LASG *
* Institute of Atmospheric Physics *
* Chinese Academy of Sciences *
* P.O.BOX 9804 *
* Beijing, 100029 *
* P. R. China *
* E-mail: baoqing@mail.iap.ac.cn *
**************************************

_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk



This archive was generated by hypermail 2b29 : Tue May 31 2005 - 08:13:59 MDT