SOS Calculating REOF

Hi everyone:

I have some big problems to ask you for helps when i calculated the REOF. Tow questions are that:

(1) One question:When it runs at the eof_rot = eofunc_varimax_Wrap(eof,1), there is error happened, with following that "[main] ncl 6100 _cygtls::handle_exceptions:error while dumping state (probably corrupted stack)". I don't know how to solve this problem.

(2)As
mentioned a sample at the website
"http://www.ncl.ucar.edu/Document/Functions/Contributed/eofunc_varimax_reorder.shtml",
if the REOF is got by "eof_rot = eofunc_varimax_Wrap", and then how to get eof_tc_rot. Is it like this eof_ts_rot=eofunc_ts_Wrap (x, eof_rot, optETS)?

 Please help me, it cost me lots of time, i feel headache right now, even crazy.
Thanks so much ahead.

Simbol Chen

The segment of code calculating the REOF is followed by:

; -----------------------------------------------------------------
; Caculating the EOF
; -----------------------------------------------------------------
  x = wsic(lat|:,lon|:,time|:)
  eof = eofunc_Wrap(x, nyrs, optEOF)
  eof_ts = eofunc_ts_Wrap (x, eof, optETS)
  
; -----------------------------------------------------------------
; Caculating the REOF
; Reorder the results returned by eof_varimax into descending order by percent variance explained.
; -----------------------------------------------------------------
  print("Caculating REOF")
  eof_rot = eofunc_varimax_Wrap(eof,1)
  eof_varimax_reorder ( eof_rot ) ; reorder the information
                                                                                              
_________________________________________________________________
Ô¼»á˵²»ÇåµØ·½£¿À´ÊÔÊÔ΢ÈíµØͼ×îÐÂmsn»¥¶¯¹¦ÄÜ£¡
http://ditu.live.com/?form=TL&swm=1

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 16 22:53:08 2010