Re: error scalar_field: If the input data is 1-dimensional, you must set sfXArray and sfYArray to 1-dimensional arrays of the same length.

Dear Vanúcia,

You read in the first timestep of SIT_NEW which gave you a 91 x 181 (latitude by longitude) array, and then you called dim_avg_n_Wrap with an index of 0, which causes it to take the average of all latitude values for every longitude value , hence giving you a one-dimensional array of length 181.

When you try to create a contour plot of a 1D array of data, you also have to provide 1D X/Y arrays that are the same length, that define the X and Y locations of each point. This is where the error is coming from.

If you are trying to calculate the average temperature for each lat/lon point across all time values, then this is what you want:

ice =f->SIT_NEW ; 120 x 91 x 181
media= dim_avg_n_Wrap(ice,0) ; 91 x 181
 

Note that this reads all of SIT_NEW’s values, and then averages across the time dimension. This will give you something you can plot over a map.

Good luck,

—Mary

On Jun 1, 2014, at 5:00 PM, Vanúcia Schumacher <vanucia-schumacher@hotmail.com> wrote:

> Hi users,
>
> I needed to interpolate the resolution of my data, and now I'm getting this error:
>
> (0) Error: scalar_field: If the input data is 1-dimensional, you must set sfXArray and sfYArray to 1-dimensional arrays of the same length.
> warning:create: Bad HLU id passed to create, ignoring it
>
> My code:
>
> 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"
>
> begin
>
> f=addfile("sit1.nc","r")
>
> ;print(f)
>
> ice =f->SIT_NEW(0,:,:)
>
> media= dim_avg_n_Wrap(ice,0)
>
> wkstype= "png"
> wkstype@wkWidth = 1000
> wkstype@wkHeight = 1000
> wks = gsn_open_wks(wkstype,"ice_2deg")
>
> gsn_define_colormap(wks,"WhiteBlue")
>
> res = True
> res@cnLevelSelectionMode ="ManualLevels"
> res@cnMinLevelValF =0
> res@cnMaxLevelValF = 2
> res@cnLevelSpacingF = 1
>
> res@cnLinesOn = False
> res@cnFillOn = True
>
>
> plot = gsn_csm_contour_map_ce(wks,media,res)
>
> end
>
> print(f)
>
> Variable: f
> Type: file
> filename: sit1
> path: sit1.nc
> file global attributes:
> history : FERRET V6.9 30-May-14
> Conventions : CF-1.0
> dimensions:
> X2DEG = 181
> Y2DEG = 91
> TIME = 120 // unlimited
> bnds = 2
> variables:
> double X2DEG ( X2DEG )
> units : degrees_east
> point_spacing : even
> axis : X
> standard_name : longitude
>
> double Y2DEG ( Y2DEG )
> units : degrees_north
> point_spacing : even
> axis : Y
> standard_name : latitude
>
> double TIME ( TIME )
> bounds : TIME_bnds
> units : days since 1850-01-01 00:00:00
> calendar : 365_day
> axis : T
> long_name : time
> standard_name : time
>
> double TIME_bnds ( TIME, bnds )
>
> double SIT_NEW ( TIME, Y2DEG, X2DEG )
> missing_value : -9.999999999999999e+33
> _FillValue : -9.999999999999999e+33
> long_name : SIT[GX=X2DEG,GY=Y2DEG]
> history : From sit_OImon_CanCM4_decadal1980_r1i1p1_198101-199012
>
>
> ---
> Vanúcia Schumacher
> Mestranda em Meteorologia - UFV
> Meteorologista -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 Sun Jun 1 22:14:46 2014