pattern_cor() and missing values

Not so much a question but a bug (or at least something that would be useful in the documentation.

The function pattern_cor() can't deal with missing values:

>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

ncl 1> x = new((/2,3,3/),"float")

ncl 2> x(0,:,:) = (/ (/1,2,3/),(/4,5,6/),(/7,8,9/)/)

ncl 6> y=x(0,:,:)

ncl 7> cor1 = pattern_cor(x(0,:,:),y,1.,0) ;works just fine…..

;now try with missing values….

ncl 8> cor2 = pattern_cor(x(1,:,:),y,1.,0)
Segmentation fault: 11

Will

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 26 00:50:29 2013