Re: Error with dim_area_avg

Hi Ipshita,
I think you might need to take a step back and really look at the data you
are trying to read in. This should always be the first thing you do when
constructing a program.

- From your printVarSummary output your sce array is of type byte. In
general you will almost always want to use the byte2flt function to convert
the byte values to floats.
http://www.ncl.ucar.edu/Document/Functions/Contributed/byte2flt.shtml

- Again from the printVarSummary your data array is 2 dimensions, but you
are specifying more than 2 dimensions when referring to sce.

- Your data is also on a curvilinear grid, in that it has 2 dimensional
latitudes and longitudes. Different syntax is needed when dealing with
rectilinear data vs. curvilinear data. You cannot use coordinate
subscripting on curvilinear data as you are doing here:
z2 = sce({10:30},{40:60},:) ; select specific region

You should try using region_ind to subset your data:
http://www.ncl.ucar.edu/Document/Functions/Contributed/region_ind.shtml
Note the warning on the page that not all curvilinear grids work with
region_ind. You will want to make sure that the subsetting works but say
drawing a quick plot of the resulting array.

2 final things:
1)If you haven't done so already please read through the NCL mini-language
manual.
http://www.ncl.ucar.edu/Document/Manuals/

2) Instead of immediately replying to every response that comes in via
ncl-talk trying to help you, try to work through some things on your own.
This is the best way to learn. Use the NCL website: Look at the examples,
look up functions, etc. You can search past ncl-talk postings by using the
search engine on the NCL website.

If after applying what I say above you are still having issues with the
script please do what Karin just requested: Send your complete script (as
an attachment) along with the exact error messages that you are getting to
ncl-talk so everyone can better assist you.

Good luck,
Adam



On Wed, Jul 23, 2014 at 3:28 PM, Ipshita Majhi <ipmajhi@alaska.edu> wrote:

> I don't know the thing is that I had a 2D file called snow cover extent
> (sce), and then I added lat and lon to it
> then the cooridinates of sce shows
> latitude longitude time and snow cover extent
>
> I made the changes you requested still the errors persist!
>
>
>
> On Wed, Jul 23, 2014 at 1:20 PM, Alessandra Giannini <
> alesall@iri.columbia.edu> wrote:
>
>>
>> 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
>
>


--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Division, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>

Received on Wed Jul 23 10:22:33 2014