ftcurvps freezing up

I am having a problem where my program freezes up
under certain conditions when using ftcurvps.
Below are a few lines from my code.

I have a 20 year period of daily max and min temperatures
(includes 5 leap year days in that period). I want to
create a time series of hourly temperatures. So I assume
the minimum temperatures occur at the 10th hour and
maximums at the 22nd hour (Z time) and apply ftcurvps
to estimate the full hourly time series.

Normally, this works great. But in some rare cases, it
hits certain grid cells and then just freezes up and
the program won't continue. The values have nothing
wrong with them, no strange values, no missing values,
just refuses to perform the task.

What might be wrong?

Michael

nyr=20

xi=new((/365*2*nyr+5*2/),float)

xi(0)=9.
xi(1)=21.
do i=1,365*nyr+5-1
   xi(i*2)=xi(i*2-2)+24.
   xi(i*2+1)=xi(i*2-1)+24.
end do

npts=365*24*nyr+5*24
xo=fspan(0.,int2flt(npts)-1.,npts)

yi=new((/365*2*nyr+5*2/),float)

yi(0::2)=mintpt
yi(1::2)=maxtpt

temp = ftcurvps(xi,yi,npts,1.,xo)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 12 08:24:41 2009