Re: Error with dim_area_avg

From: Alessandra Giannini <alesall_at_nyahnyahspammersnyahnyah>
Date: Wed Jul 23 2014 - 15:20:37 MDT


Hi Ipshita,

I see an error in your very last script [but not in the prior version]. Could it be so simple as the following:

the last parenthesis in your definition of z2 should be a ")", not a "}" ??

z2=sce({10:30},{40:60},:,:}

should actually be

z2=sce({10:30},{40:60},:,:)

Ah, also in the previous script, I think you have an extra "}" at the end:

> z2 = sce({10:30},{40:60},:,:})

there shouldn't be a "}" before the last ")"

i.e. it should read

> z2 = sce({10:30},{40:60},:,:)


best, alessandra








On Jul 23, 2014, at 5:13 PM, Ipshita Majhi <ipmajhi@alaska.edu> wrote:

> Here are the error messages
>
> fatal: error in statement
> fatal: syntax error : line 3 in file
> area_avg_ims_sce.ncl
>
> *
> ^
> fatal : error in statement
> fatal: syntax error :line 22 in file
> area_avg_ims_sce.ncl before or near }
>
> z2=sce({10:30},{40:60},:,:}
>
> fatal:sytax error : line 32 in file area_avg_ims_sce.ncl before or near ; create
> timeXlong array
> z2ma=dim_avg_wrap(z2(time
> fatal:syntax error: line -1 before or near
> ^
> fatal: error in statement
>
>
> On Wed, Jul 23, 2014 at 12:48 PM, Karin Meier-Fleischer <meier-fleischer@dkrz.de> wrote:
> Please send the complete error message.
>
>
> Am 23.07.2014 um 22:32 schrieb Ipshita Majhi <ipmajhi@alaska.edu>:
>
>> Why is this not working? It says the synatax at Z2 is not correct
>>
>> ************************************************
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> ;************************************************
>> begin
>> ;************************************************
>> ; read in netCDF file
>> ;************************************************
>> a =addfile("~/Documents/Snow_cover_IMS/nhsce_v01r01_19661004_20140602.nc","r")
>>
>> sce=a->snow_cover_extent(0,:,:)
>> lat=a->latitude
>> lon=a->longitude
>>
>>
>> ;Attaching latitude and longitude to sce
>> sce@lat2d=lat
>> sce@lon2d=lon
>>
>> z2 = sce({10:30},{40:60},:,:}) ; select specific region
>> delete(sce)
>> printVarSummary(z2)
>>
>> z2Area = wgt_areaave_Wrap(z2,1.0,1.0,0) ; create area-avg. timeseries
>>
>> printVarSummary(z2Area)
>>
>> z2_za = dim_avg_Wrap(z2) ; create time x lat array
>>
>> z2_ma = dim_avg_Wrap(z2(time|:,lon|:,lat|:) ; create time x lon array
>>
>> printVarSummary(z2_za)
>> printVarSummary(z2_ma)
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

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

Received on Wed Jul 23 09:20:47 2014

This archive was generated by hypermail 2.1.8 : Fri Aug 01 2014 - 15:10:55 MDT