Hi Dennis,
I understand the error, but now returns this:gsn_csm_contour_map_polar: F=
atal: the input data array must be 1D or 2Dfatal:Illegal right-hand side ty=
pe for assignmentfatal:["Execute.c":7743]:Execute: Error occurred at or nea=
r line 51 in file tsm_polar.ncl
Thank you for helps!
Date: Wed, 16 Jul 2014 12:36:28 -0600
Subject: Re: error month_to_seasonN
From: shea@ucar.edu
To: vanucia-schumacher@hotmail.com
CC: ncl-talk@ucar.edu
[1]
The array you are inputting is rank 4 with dimensions (ncl_join,time,la=
t,lon)
[ncl_join | 25] x [time | 12] x [LAT | 224] x [LON | 256]
[2]
The documentation for the month_to_seasonN states the following:
xMon
=
A one-, three-, or four-dimensional array (xMon(time) or
=
xMon(time,lat,lon) or xMon(time,lev,lat,lon) of any
=
numeric type.
xMon must have named dimensions and the time (leftmost)
=
dimension must be divisible by 12. The data are assumed to be monthly
=
mean data and the first record is assumed to be January.
[3]
Hence, your variable does *not* adhere to the documentation. The function=
*requires* 'time' to be the leftmost dimension. In your case. it is not!
====
I do not know why you used
ListSetType(f,"join")
Delete this line and try.
On Wed, Jul 16, 2014 at 12:13 PM, Vanúcia Schumacher <vanucia-schum=
acher@hotmail.com> wrote:
I changed the resolution of the original data to 2 degrees (using linint2)=
, and now need to plot the DJF seasonal average for the data, but retur=
ns this error, and I do not understand why the time has dimensions of jus=
t 12. How do I solve this?
=
Variable: tType: float
=
Total Size: 68812800 bytes 17203200 valuesNumber of Dimensions: =
4Dimensions and sizes: [ncl_join | 25] x [time | 12] x [LAT | 224] x [LON |=
256]
=
Coordinates: LAT: [-89.59821..89.59821] LON: [-180..=
178.5938]Number Of Attributes: 3 _FillValue : 9.96921e+36
=
units : degC long_name : SST
(0) contributed.ncl: month_to_seasonN: dimension must be a multiple of 12=
fils = systemfunc("ls *.nc")
f = addfiles(fils, "r")
ListSetType(f,"join")
=
t =f[:]->bsst(:,:,:,:)
printVarSummary(t)
;media= dim_avg_n_Wrap(t, 0)
x=month_to_seasonN(t,"DJF")
.....
plot = gsn_csm_contour_map_polar(wks,x,res)
---Vanúcia Schumacher
=
Mestranda em Meteorologia - UFVMeteorologista -UFPel
Departamento de Meteorologia Agrícola - DEA
=
Cel: (31) 9978 2522
=
DEA: (31) 3899 1890
=
=
_______________________________________________
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 16 06:43:16 2014
This archive was generated by hypermail 2.1.8 : Fri Aug 01 2014 - 15:10:55 MDT