Re: Problem

Adam,

It just stops at :

v=f->a

and not go beyond this point.

matt

________________________________
From: Adam Phillips <asphilli_at_ucar.edu>
To: Mathew Freeman <mattfree0180_at_yahoo.com>
Cc: ncl-talk_at_ucar.edu
Sent: Friday, June 12, 2009 11:56:40 AM
Subject: Re: Problem

Hi Mathew,
Can you tell us which line is causing the segmentation fault? You can use print statements to isolate the line. If the error occurs within the double do loop, let us know the i and j values when the script fails. (Again, you can use print statements for that.)
Adam

Mathew Freeman wrote:
> Hello everyone,
> I am using the following script to convert the daily 0.25 degree data into monthly data but getting segmentation fault error. Can anyone please guide to this case. The script is as follows:
>
>
> begin
> f=addfile("0904.nc","r")
> t=f->time
> v=f->a
> v2=v(latitude|:,longitude|:,time|:)
> delete(v)
> t=fspan(0,dimsizes(t)-1,dimsizes(t))
> utc_date = ut_calendar(t, 0)
>
> year = floattointeger(utc_date(:,0))
> month = floattointeger(utc_date(:,1))
>
> dims=dimsizes(v2)
> v1=new((/12*(2002-1970+1),dims(0),dims(1)/),"float")
> time1=fspan(0,12*(2002-1970+1)-1,12*(2002-1970+1))
> time1_at_units="months since 1970-01-01"
> k=0
> do i=1970,2002
> do j=1,12
>
> indx = ind(year.eq.i.and. month.eq.j)
> printVarSummary(indx)
> print(max(indx))
> print(min(indx))
> v1(k,:,:)=dim_avg(v2(:,:,indx))
> delete(indx)
> ; if (j.eq.1) then
> print(i*100+j)
> ; end if
> k=k+1
> end do
> end do
>
> v1!0="time"
> v1&time=time1
> v1!1="lat" v1&lat=f->latitude
> v1!2="lon"
> v1&lon=f->longitude
>
> system("rm 0904_mon.nc")
> fo=addfile("0904_mon.nc","c")
> fo->a=v1
> end
>
> thankyou
>
> matt
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- --------------------------------------------------------------
Adam Phillips asphilli_at_ucar.edu
National Center for Atmospheric Research tel: (303) 497-1726
ESSL/CGD/CAS fax: (303) 497-1333
P.O. Box 3000
Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jun 12 2009 - 11:09:46 MDT