Update from 5.2.0 to 5.2.1 breaks function readAsciiTable()

From: Marco Steinacher <steinacher_at_nyahnyahspammersnyahnyah>
Date: Wed Mar 30 2011 - 01:18:42 MDT

Hi

The update from NCL 5.2.0 to 5.2.1 applies the following change to the
function readAsciiTable() [ncarg/nclscripts/csm/contributed.ncl line 7893]:

- head = readAsciiHead( fName, opt(0))
- nh = dimsizes(head) ; # header records
+ if(opt(0).gt.0) then
+ head = readAsciiHead( fName, opt(0))
+ nh = dimsizes(head) ; # header records
+ else
+ nh = 0 ; no header to read
+ end if

According to the documentation opt(0) can be a number or a string to
specify the end of the header. In 5.2.1 the function does not work
correctly anymore if opt(0) is a string.

To fix this I changed the if-statement in the code above to:

   if((nopt.eq.1 .and. typeof(opt).eq."string") .or. opt(0).gt.0) then

With that change my scripts work again.

Cheers,
Marco

-- 
***************************************
Marco Steinacher
Climate and Environmental Physics
Physics Institute, University of Bern
Sidlerstrasse 5, CH-3012 Bern
Phone ++41 (0)31 631 34 02
Fax   ++41 (0)31 631 87 42
steinacher@climate.unibe.ch
http://www.climate.unibe.ch/~steinach
OpenPGP Key ID: 0x62937F7F
***************************************
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Mar 30 01:26:51 2011

This archive was generated by hypermail 2.1.8 : Tue Apr 05 2011 - 09:01:21 MDT