Re: FW: lspoly didn't wok when the number of data points equal to the number of coefficients

Forgot NCL driver ... attached

On 01/03/2014 01:25 PM, Dennis Shea wrote:
> Attached is the SLATEC fortran code and a sample test script
> that calls the fortran subroutine. You will have to change the
> NCL driver script.
>
> %> WRAPIT slatec_DPOLFT.f
>
> Please note: This is a double precision version of the
> fortran code. The data passed to it must be dpoble
> precision. Hence, if x and y are type float
>
> xd = xf*1d0 ; new array in double precisiom
> yd = yf*1d0
>
> LSPOLY::POLFT (M, xd, yd, wgt ...)
>
> ===
>
> Please read the Mini-Language Manual for calling fortran code
>
> http://www.ncl.ucar.edu/Document/Manuals/
>
> Good luck
>
[SNIP]

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk