getting data plotted on a map

Hi,

I'm trying to get some WRF output plotted onto a map, but am unable to
do so. The error messages I get are
(0) is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must
either be the same dimension sizes as the data, or one element larger in
both directions. Your data will most likely not be overlaid on the map
correctly.
(0) check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units'
attribute equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
north'

However, I believe that I have addressed each of the possible error
messages. I have put in printVarSummary commands for the array of
interest, plus the lat2d and lon 2d arrays, and it appears to me like it
all looks fine. Would you possibly be able to help me find what I'm
missing? If so, the CODE is directly below, while the output from the 3
printVarSummary lines is immediately below that.

Thanks!

Joe Grim

CODE
;*************************************************
; NCL Graphics: raster_1.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
;************************************************
; open netCDF file
;************************************************
  a = addfile("/d2/CFDRC/output/prelim_precip_output.nc","r")
  hourly_mean_precip = a->hourly_mean_precip
  lat2d = a->lat2d
  lon2d = a->lon2d
;************************************************
; create plots
;************************************************
printVarSummary(hourly_mean_precip)
printVarSummary(lat2d)
printVarSummary(lon2d)
do ii=0,23
  char_ii=sprinti("%0.2i",ii)
  wks = gsn_open_wks("ps","hourly_mean_precip_"+char_ii+"Z") ;
open a ps file
  res = True ; plot mods desired

  gsn_define_colormap(wks,"WhBlGrYeRe") ; Choose colormap.

  res_at_gsnMaximize = True ; Maximize plot in frame.

  res_at_gsnStringFontHeightF = 0.013
  res_at_gsnRightStringOrthogonalPosF = 0.02
  res_at_gsnLeftStringOrthogonalPosF = 0.02

  res_at_cnFillOn = True ; Color plot desired.
  res_at_cnLinesOn = False ; Turn off contour lines .
  res_at_cnLineLabelsOn = False ; Turn off contour labels.

  res_at_mpProjection = "LambertConformal"
  res_at_mpLambertParallel1F = 33.0 ; two parallels
  res_at_mpLambertParallel2F = 45.0
  res_at_mpLambertMeridianF = avg(lon2d) ; central meridian
  res_at_mpLimitMode = "LatLon"

  res_at_mpMinLatF = min(lat2d) ; map area
  res_at_mpMaxLatF = max(lat2d) ; lat2ds
  res_at_mpMinLonF = min(lon2d) ; and
  res_at_mpMaxLonF = max(lon2d) ; lon2ds

  res_at_gsnSpreadColors = True ; Use full range of colormap.
  res_at_gsnSpreadColorStart = 0 ; Start color in colormap.
  res_at_gsnSpreadColorEnd = 190 ; End color in colormap.
  res_at_cnInfoLabelOn = False ; Turn off info label.

  res_at_mpGeophysicalLineColor = "red" ; color of continental outlines
  res_at_mpUSStateLineColor = "red"
  res_at_mpGridLineDashPattern = 2 ; lat/lon lines as dashed
  res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"

  res_at_pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks.

  res_at_cnFillMode = "RasterFill"
  res_at_cnLevelSelectionMode = "ManualLevels"
  res_at_cnMinLevelValF = 0.1 ; min contour level
  res_at_cnMaxLevelValF = max(hourly_mean_precip) ; max
contour level
  res_at_cnLevelSpacingF = 1.0 ; contour spacing

  res_at_pmLabelBarWidthF = 0.6
  res_at_lbLabelFontHeightF = 0.013
  res_at_lbLabelStride = 2
  res_at_lbBoxLinesOn = False

  ;plot =
gsn_csm_contour(wks,hourly_mean_precip(ii,:,:),res)
  map = gsn_csm_contour_map(wks,hourly_mean_precip(ii,:,:),res)

end do

end

OUTPUT from printVarSummary lines
Variable: hourly_mean_precip
Type: float
Total Size: 960000 bytes
            240000 values
Number of Dimensions: 3
Dimensions and sizes: [ncl5 | 24] x [ncl6 | 100] x [ncl7 | 100]
Coordinates:
Number Of Attributes: 5
  long_name : Mean precipitation rate for time of day
  units : mm/hour
  lon2d : <ARRAY of 10000 elements>
  lat2d : <ARRAY of 10000 elements>
  _FillValue : -999

Variable: lat2d
Type: float
Total Size: 40000 bytes
            10000 values
Number of Dimensions: 2
Dimensions and sizes: [g3_x_0 | 100] x [g3_y_1 | 100]
Coordinates:
Number Of Attributes: 11
  La1 : 26.856
  Lo1 : -82.448
  Lov : -80.6
  Dx : 3000
  Dy : 3000
  Latin1 : 60
  Latin2 : 30
  units : degrees_north
  GridType : Lambert Conformal Secant or Tangent, Conical or bipolar
  long_name : latitude
  corners : ( 26.856, 26.86411, 29.78162, 29.77304 )

Variable: lon2d
Type: float
Total Size: 40000 bytes
            10000 values
Number of Dimensions: 2
Dimensions and sizes: [g3_x_0 | 100] x [g3_y_1 | 100]
Coordinates:
Number Of Attributes: 11
  La1 : 26.856
  Lo1 : -82.448
  Lov : -80.6
  Dx : 3000
  Dy : 3000
  Latin1 : 60
  Latin2 : 30
  units : degrees_east
  GridType : Lambert Conformal Secant or Tangent, Conical or bipolar
  long_name : longitude
  corners : ( -82.448, -79.20031, -79.14114, -82.52611 )
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 01 2009 - 16:16:58 MDT