Getting Error in lat, lon

Dear NCL users,
                 I want to get a regional contour plot [lon 60:100,lat
-5:28] of my 2D variable, using binary (direct access) file.
While running my script am getting error regarding my lat, lon. My script
 is as follows, and I have already ftp my binary file (2012_RH.dat) to ftp
ftp.cgd.ucar.edu.I also want to put the red shades only in my color scheme.
**************************************************
 ;Plot relative humidity
 ; drawing a color filled contours using color map

;******************************************
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

;******************************************************

  nrec=0
  dims=(/561,3601/)
  rh=fbindirread("H:/2012_RH.dat",0,dims,"float")
  lat_g = latGlobeF(561, "lat", "latitude", "degrees_north")
  lon_g=lonGlobeF(3601, "lon","longitude","degrees_east")
  lon_g=(/lon_g-180. /)
  lon_g&lon=lon_g
  rlon=lon_g({60.1:100})
  rlat=lat_g({-5.1:28})
   ;******************************************
   wks=gsn_open_wks("pdf","color")

   res =True
   res@tiMainString = "RH: 2012"
   res@cnFillOn= True
   plot= gsn_csm_contour_map(wks,rh,res)

Thanks,
Ipsita

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Dec 2 07:23:21 2012