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-talkReceived on Fri Jun 12 2009 - 10:56:40 MDT
This archive was generated by hypermail 2.2.0 : Mon Jun 15 2009 - 08:53:05 MDT