Re: Error with dim_area_avg

Please send the printVarSummary of the changed sce variable and your =
complete script,
it is impossible to say what is going wrong without knowing it.


Am 24.07.2014 um 00:05 schrieb Ipshita Majhi <ipmajhi@alaska.edu>:

> Yes I did see the sce printVarSummary. I added latitude and longitude =
to it
> and now I want extract sce data for particular latitude and longitude
>
> I saw that after adding latitude and longitude
> the values are
>
> (87,87),0
>
> Does that make it three dimension
>
>
> On Wed, Jul 23, 2014 at 2:03 PM, Karin Meier-Fleischer =
<meier-fleischer@dkrz.de> wrote:
> Alessandra is right, there is one } to much at the end of the =
following line:
>
> z2 = sce({10:30},{40:60},:,:}) ; select specific region
>
> it should be
>
> z2 = sce({10:30},{40:60},:,:) ; select specific region
>
> BUT your printVarSummary of the variable sce shows that the variable =
has only
> 2 dimensions (rows and cols) so either way the line z2 is incorrect. =

>
> Bye,
> Karin
>
> Am 23.07.2014 um 23:14 schrieb Ipshita Majhi <ipmajhi@alaska.edu>:
>
>> I am also adding the PrintVarSummary in addition to the other error
>> Variable:sce
>> Type:byte
>> Total Size:7744 bytes
>> :7744 bytes
>> Number of dimensions :2
>> Dimesions and sizes: [rows | 88]X [cols | 88]
>> Coordinates :
>> Number of attributes:9
>> time:7
>> long_name:NOAA/NCDC Climate Date Record of snow cover extent
>> standard_name:surface_snow_
>> binary_mask
>> grid_mapping:coord_system
>> coordinates:longitude latitude time snow_cover_threshold
>> valid range: (0,1)
>> flag_values : (0,1)
>> flag_meanings:no_snow snow_covered
>> _FillValue: -9
>>
>> *************
>> Variable:lat
>> Type:double
>> Total Size : 61952 bytes
>> 7744 values
>> Number of Dimensions :2
>> Dimensions and sizes: [rows | 88] X [cols | 88]
>> Coordinates:
>> Number of Attributes:5
>> standard_name: latitude
>> long_name: latitude of cell center in half mesh grid, based on NOAA =
NMC Limited Area Fine Mesh Grid
>> Units : degree_north
>> valid_range : (0,90)
>> _FillValue: -999
>>
>> ************
>> Variable: lon
>> Type : double
>> Total size: 61952 bytes
>> 7744 values
>> Number of Dimensions: 2
>> Dimensions and sizes :[rows | 88] X [ cols | 88]
>> Coordinates:
>> Number of Attributes :5
>> standard_name: longitude
>> long_name: longitude of cell center in half mesh grid, base on NOAA =
NMC Limited Area Fine Mesh grid unites degrees_east
>> valid_range(-180,180)
>> _FillValue:: -999
>>
>>
>> On Wed, Jul 23, 2014 at 1: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


Received on Wed Jul 23 10:10:47 2014