[ncl-install] Week average with years with varying weeks
Ipshita Majhi
ipmajhi at alaska.edu
Wed Nov 12 09:36:56 MST 2014
Thank you
On Wed, Nov 12, 2014 at 9:23 AM, Mary Haley <haley at ucar.edu> wrote:
> This forum is for people having trouble installing NCL.
>
> Please repost your question to ncl-talk at ucar.edu.
>
> --Mary
>
>
> On Tue, Nov 11, 2014 at 10:07 PM, Ipshita Majhi <ipmajhi at alaska.edu>
> wrote:
>
>> Hi,
>>
>> I wrote a code to extract each year and then take a long term weekly
>> average i.e for 52 weeks. There are some years which have 53 years . I get
>> a matrix of 46 by 52 but it is not writing it out in a matrix.
>> Error: Subscript out of range, error in subscript #0
>> Then when I average it, it gives me only one value instead of 52 and if
>> I use
>> dim_avg_n then it gives other sets of concern.
>>
>> Error : Invalid dimension sizes
>>
>> ;*******************************************
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>> ;*******************************************
>>
>> ;Reading in the data
>>
>> a=addfile("~/Documents/Snow_cover_IMS/nhsce_v01r01_19661004_20140602.nc
>> ","r")
>>
>> sce=a->snow_cover_extent
>> time_sce=a->time
>>
>>
>> ;Now converting the time from georgian to ut_calender
>> ;Doing an area average for sce
>>
>> sce_area=wgt_areaave_Wrap(sce,1.0,1.0, 0)
>>
>> ut_sce=cd_calendar(time_sce,0)
>>
>>
>> ;Extracting just the year information from the time series
>>
>> xyears=ut_sce(:,0); Years from the SCE dataset
>>
>> ;This part of code is to create a matrix with each years weekly data
>> ;The years go from 1967 to 2013 i.e 47 years with 52 weeks for each year
>>
>> sce_week=new((/46,52/),"float")
>> years=ispan(1967,2013,1)
>>
>>
>> do i=1,46
>>
>>
>> yind=years(i)
>>
>>
>> ;Finding the index matching the year
>>
>> time_ind=ind(xyears.eq.2012)
>>
>> ; Extracting that years snow cover extent
>>
>> sce_year=sce_area(time_ind)
>> sce_52=sce_year(0:51)
>>
>> ;Now checking for years with 53 weeks
>> z=dimsizes(sce_year)
>>
>> sce_week(i,0:51)=sce_52(i)
>>
>> opt = True
>> opt at title = "sce_week_matrix.txt"
>>
>> ;write_matrix(sce_week, "52f1.10", opt)
>>
>> end do
>>
>> sce_week_avg=dim_avg_n_Wrap( sce_week, 2 )
>>
>> Thank YOu
>> Best Regards
>> Ipshita
>>
>> _______________________________________________
>> ncl-install mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-install
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-install/attachments/20141112/36d3c345/attachment-0001.html
More information about the ncl-install
mailing list