[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