[ncl-talk] warning:ScalarFieldSetValues: 2d coordinate array sfYArray has an incorrect dimension size: defaulting sfYArray
Amit Kumar
amitkumar.777 at hotmail.com
Thu Jul 7 01:19:51 MDT 2016
Hi All,
I did following two modifications and this error was removed.
res at gsnAddCyclic = False
plot(i) = gsn_csm_contour_map(wks,var,res)
But now I'm getting 1st plots data in all the plots. It is displaying only data of 1st nc file in all the plots. It is not reading rest files.
Please help me out.
Regards
Amit Kumar,
________________________________
From: Amit Kumar <amitkumar.777 at hotmail.com>
Sent: Thursday, July 7, 2016 10:52 AM
To: ncl-talk at ucar.edu
Subject: warning:ScalarFieldSetValues: 2d coordinate array sfYArray has an incorrect dimension size: defaulting sfYArray
Hi All,
I'm new to NCL, I'm trying to make a panel plot of HEM for a day. so i concatenated all the files into one file using cdo.
amit at amit-HP-xw6400-Workstation:~/Desktop/images/output/UK$ cdo infov out2.nc
-1 : Date Time Level Gridsize Miss : Minimum Mean Maximum : Parameter name
1 : 2016-06-30 10:00:00 0 3124 0 : 0.0000 0.18826 28.871 : HEM
2 : 2016-06-30 10:30:00 0 3124 0 : 0.0000 0.40692 41.558 : HEM
3 : 2016-06-30 11:00:00 0 3124 0 : 0.0000 0.27961 39.316 : HEM
4 : 2016-06-30 11:30:00 0 3124 0 : 0.0000 0.39080 54.272 : HEM
5 : 2016-06-30 12:00:00 0 3124 0 : 0.0000 0.24352 21.233 : HEM
6 : 2016-06-30 12:30:00 0 3124 0 : 0.0000 0.34891 26.116 : HEM
7 : 2016-06-30 13:00:00 0 3124 0 : 0.0000 0.31744 20.495 : HEM
8 : 2016-06-30 13:30:00 0 3124 0 : 0.0000 0.26712 15.045 : HEM
9 : 2016-06-30 14:00:00 0 3124 0 : 0.0000 0.14747 16.468 : HEM
10 : 2016-06-30 14:30:00 0 3124 0 : 0.0000 0.21267 46.071 : HEM
11 : 2016-06-30 15:00:00 0 3124 0 : 0.0000 0.46982 65.314 : HEM
12 : 2016-06-30 15:30:00 0 3124 0 : 0.0000 0.57988 103.28 : HEM
13 : 2016-06-30 16:00:00 0 3124 0 : 0.0000 0.34465 232.33 : HEM
14 : 2016-06-30 16:30:00 0 3124 0 : 0.0000 0.65428 84.749 : HEM
15 : 2016-06-30 17:00:00 0 3124 0 : 0.0000 1.4974 132.02 : HEM
16 : 2016-06-30 17:30:00 0 3124 0 : 0.0000 1.3078 121.53 : HEM
17 : 2016-06-30 18:00:00 0 3124 0 : 0.0000 1.5186 116.90 : HEM
18 : 2016-06-30 18:30:00 0 3124 0 : 0.0000 2.8212 139.77 : HEM
19 : 2016-06-30 19:00:00 0 3124 0 : 0.0000 4.9052 283.62 : HEM
20 : 2016-06-30 19:30:00 0 3124 0 : 0.0000 3.9615 161.71 : HEM
21 : 2016-06-30 20:00:00 0 3124 0 : 0.0000 3.8913 144.44 : HEM
22 : 2016-06-30 20:30:00 0 3124 0 : 0.0000 7.5982 232.91 : HEM
23 : 2016-06-30 21:00:00 0 3124 0 : 0.0000 2.8257 94.449 : HEM
24 : 2016-06-30 21:30:00 0 3124 0 : 0.0000 2.3140 122.16 : HEM
25 : 2016-06-30 22:00:00 0 3124 0 : 0.0000 2.4510 107.22 : HEM
26 : 2016-06-30 22:30:00 0 3124 0 : 0.0000 3.0760 114.25 : HEM
27 : 2016-06-30 23:00:00 0 3124 0 : 0.0000 2.5792 92.298 : HEM
28 : 2016-06-30 23:30:00 0 3124 0 : 0.0000 3.4650 120.34 : HEM
cdo infon: Processed 87472 values from 1 variable over 28 timesteps ( 0.00s )
when i'm trying to plot this in a panel i'm getting this error, warning:ScalarFieldSetValues: 2d coordinate array sfYArray has an incorrect dimension size: defaulting sfYArray. and the plots are blank.
my ncl script is as follows:
load "/usr/share/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/share/ncarg/nclscripts/csm/contributed.ncl"
begin
a = addfile("/home/amit/Desktop/images/output/UK/out2.nc","r")
s = "/media/amit/E/India_shapefile/Districtbound(Satmet).shp"
var = doubletofloat(a->HEM(0,:,:))
lat = a->latitude
lon = a->longitude
wks = gsn_open_wks("png","HEM_30JUN2016") ; send graphics to PNG file
res = True ; plot mods desired
res at gsnDraw = False
res at gsnFrame = False
res at gsnAddCyclic = True
res at gsnStringFontHeightF = 0.015
res at cnFillOn = True ; turn on color fill
res at cnLinesOn = False
res at gsnMaximize = True ; maximize plot in frame
res at mpOutlineOn = False ; Use outlines from shapefile
res at mpFillOn = False
res at cnLevelSelectionMode = "ExplicitLevels"
res at cnLevels = (/0.1,2.5,15.5,64.5,115.5,204.5,400.5,600.5,800.5/) ; ; 13 contour values
res at cnFillPalette = (/"Snow","LightGoldenRodYellow","Yellow"\ ; 13 contour colors
,"Orange","Green","Green4" \
,"Cyan","Blue2","Blue3","Blue4"/)
res at tiMainString = "HEM_30JUN2016"
res at tiMainFontHeightF = 0.016
res at lbLabelBarOn = False
res at lbOrientation = "Vertical"
res at sfXArray = lon
res at sfYArray = lat
res at mpMinLonF = 78.0
res at mpMaxLonF = 81.0
res at mpMinLatF = 29.0
res at mpMaxLatF = 31.0
res at pmTickMarkDisplayMode = "Always"
res at tiXAxisString = "Longitude"
res at tiYAxisString = "Lattitude"
res at tiXAxisFontHeightF = 0.02
res at tiYAxisFontHeightF = 0.02
;create empty graphic plot array
nplot = 28
plot = new(nplot,"graphic")
;create the plots
do i=0,nplot-1
if (i.ne.0) then
res at gsnLeftString = ""
res at gsnRightString = ""
end if
plot(i) = gsn_csm_contour_map_ce(wks,var,res)
end do
;panel res
pres = True
pres at gsnPanelTop = 0.96
pres at gsnPanelBottom = 0.012
pres at gsnPanelLabelBar = True
pres at pmLabelBarOrthogonalPosF = -0.02
pres at txstring = "Panel: 3 X 3 + Common label and title"
pres at txFontHeightF = 0.020
pres at txPosXF = 0.5
pres at txPosYF = 0.85
pres at txJust = "CenterCenter"
gsn_panel(wks,plot,(/7,4/),pres)
end
variable summary is as follows:
Variable: var
Type: float
Total Size: 12496 bytes
3124 values
Number of Dimensions: 2
Dimensions and sizes: [44] x [71]
Coordinates:
Number Of Attributes: 1
_FillValue : -9e+33
Please help me out of this
Regards
Amit Kumar,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160707/2ff0d24e/attachment.html
More information about the ncl-talk
mailing list