[ncl-talk] Issue with shapefile_mask_data

[Amit Kumar]

Hi All,

I want to do shape file masking but I'm getting
shapefile_mask_data: Error: not a valid rectilinear, curvilinear, or unstructured grid

My variable description by ncl_filedump is as follows:

Variable: f
Type: file
   dimensions:
      time = 1
      lon = 2805
      lat = 2816
   variables:
      double time ( time )
         units :    hours since 0001-01-01 00:00:00.0
         calendar :    gregorian

      double latitude ( lat, lon )
         units :    degrees_north

      double longitude ( lat, lon )
         units :    degrees_east

      double HEM_MNT ( time, lat, lon )
         coordinates :    latitude longitude
         fillvalue :    -9.e+33
         missing_value :    -9e+33
         _FillValue :    -9e+33


and my script section for masking is as follows:

  var                 = doubletofloat(a->HEM_MNT)
  lat                 =  a->latitude
  lon                 =  a->longitude
  opt             = True
  opt at return_mask = True
  shp_mask = shapefile_mask_data(var(0,:,:),filename,opt)  -------> error in this line

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160809/b075a652/attachment.html