[ncl-talk] CSV file plotting issue.

[Ajay Bankar]

Dear NCL Users,
                       I'm trying to make spatial plot of rainfall from
station data available in csv file. I wrote the following script but data
is not masking with shapefile. Also I'm getting the warning as:
*warning:tofloat: A bad value was passed to (string) tofloat, input strings
must contain numeric digits, replacing with missing value.*
I have attached plot and sample data file for reference.

Thanks for any help!

load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/scratch/rakesh/WRF-AJAY/WRFV3/RUN_14AUG/NCL_PROGRAM/shp.ncl"

begin
 ;-- shapefile name
 shp_filename =
"/scratch/Administrative_Boundary_Headquarters/District_Boundary.shp"
 infile = "/scratch/rakesh/WRF-AJAY/DATASETS/STATION/station_data_final.csv"
 lines = asciiread(infile, -1, "string")
;print(lines)
;delim = ","

rain = tofloat( str_get_field(lines(1:), 7, ","))
rain at lat1d  = tofloat( str_get_field(lines(1:), 5, ","))
rain at lon1d  = tofloat( str_get_field(lines(1:), 6, ","))
;print(rain)

;-- define sub-region, here for Karnataka
     minlat =  11
     maxlat =  19
     minlon =  73.5
     maxlon =  79.1

;-- open workstation
     wks_type          = "pdf"
     wks  = gsn_open_wks(wks_type,"1_station_plot")
 ;-- resource settings
     res                        =  True
     res at gsnDraw                =  False                       ;-- don't
draw plot yet
     res at gsnFrame               =  False                       ;-- don't
advance frame yet
     res at gsnMaximize            =  True                        ;-- maximize
plot in frame

  ;  res at tfDoNDCOverlay        = True

    cmap     := read_colormap_file("BlAqGrYeOrReVi200")
    cmap(0,:) = (/0,0,0,0/)    ; make first color fully transparent
    res at cnFillOn               =  True
    res at cnLinesOn              =  False
    res at cnLineLabelsOn         =  False
    res at cnFillPalette          = cmap
;"BlueYellowRed"
    res at cnLevelSelectionMode   = "ManualLevels"               ;-- manually
set contour levels
    res at cnMinLevelValF         = 5                            ; 260.0
                ;-- contour min lev
    res at cnMaxLevelValF         = 95                           ; 280.0
                ;-- contour max lev
    res at cnLevelSpacingF        = 10                           ;   0.5
                ;-- contour spacing
    res at gsnAddCyclic           = False

    res at mpGridLineColor        =  0                          ;"grey40"
    res at mpFillOn               =  False
    res at mpOutlineOn            =  False
    res at mpGeophysicalLineColor = "black"
    res at mpLimitMode            = "LatLon"
    res at mpMinLatF              = minlat
    res at mpMaxLatF              = maxlat
    res at mpMinLonF              = minlon
    res at mpMaxLonF              = maxlon
    ;res at mpAreaMaskingOn        = True

   res at lbBoxMinorExtentF      =  0.1                         ;-- decrease
height of labelbar boxes
   res at pmLabelBarOrthogonalPosF = 0.08                      ;-- move
labelbar to the left side of plot
   res at lbBoxMinorExtentF      =  0.2                         ;-- decrease
height of labelbar boxes

   res at pmTickMarkDisplayMode  = "Always"
   res at tmXTOn                 =  False
   res at tiMainString           = "STATION DATA PLOT"

 ;-- create plot
    plot = gsn_csm_contour_map(wks,rain,res)
    var_mask = shapefile_mask_data(rain,shp_filename,True)

 ;-- create contours of masked data
    plot_mask = gsn_csm_contour_map(wks,var_mask,res)

 ;-- add plot_mask to plot
     lnres                    =  True
     lnres at gsLineColor        = "black"
     lnres at gsLineThicknessF   =  1.0

     plot at lines = gsn_add_shapefile_polylines(wks, plot_mask, shp_filename,
lnres)

  ;-- draw the plot and advance the frame
     draw(plot_mask)
     frame(wks)
 end
-- 

*Thanks & Regards,*

*Ajay Bankar.*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191207/830f7fcb/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: station_plot.pdf
Type: application/pdf
Size: 1037167 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191207/830f7fcb/attachment-0001.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sample.csv
Type: text/csv
Size: 174753 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191207/830f7fcb/attachment-0001.csv>