Problem with gsn_csm_contour_map

From: Andy Barrett <apbarret_at_nyahnyahspammersnyahnyah>
Date: Thu, 09 Nov 2006 11:21:42 -0700

Hi

I'm trying to plot seasonal means of evaporation from ERA-40. I can
make plots when allow gsn_csm_contour_map to select contour levels
automatically but I get the following error message if I try to define
contour levels using cnLevels.

fatal:ContourPlotDraw: CPCLAM - CONTRADICTORY AREA-IDENTIFIER INFORMATION
fatal:ContourPlotDraw: draw error
warning:WorkstationDeactivate: workstation not active or not opened

I can plot other fields such as precipitation but run into the same
error message if I try other evaporation fields. So I understand that
this is a data problem but I am not sure how to fix it.

Ultimately, I want to plot fields for the north polar region but I get
the same error message when I try to do this.
 
The printVarSummary result and max/min for the field and a snippit of
the code I'm using are below.

Thanks

Andy Barrett
NSIDC-CIRES
University of Colorado at Boulder

--- The printVarSummary for the field is

Variable: t1
Type: float
Total Size: 204800 bytes
            51200 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 160] x [lon | 320]
Coordinates:
            lat: [-89.14152..89.14152]
            lon: [ 0..358.875]
Number Of Attributes: 12
  season : DJF
  average_op_ncl : dim_avg over dimension: time
  sum_op_ncl : dim_sum over dimension: initial_time0
  center : European Center for Medium-Range Weather Forecasts - Reading
  long_name : Evaporation
  units : m
  level_indicator : 1
  gds_grid_type : 4
  parameter_table_version : 128
  parameter_number : 182
  forecast_time : 6
  _FillValue : -999

(0) Evap: max/min: -319.673 / 3.39301

---- The code snippit for making the plot is...

;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","era40_evap") ; open a ps file
  gsn_define_colormap(wks,"gui_default")

  plot = new(2,graphic)

  res = True

  res_at_mpProjection = "Mercator" ; choose projection
  res_at_mpGridAndLimbOn = True ; turn on lat/lon lines
  res_at_mpPerimOn = False ; turn off box around plot
  res_at_mpGridLatSpacingF = 20. ; spacing for lat lines
  res_at_mpGridLonSpacingF = 30. ; spacing for lon lines
  res_at_mpFillOn = False

  res_at_cnFillOn = True ; color plot desired
  res_at_cnLineLabelsOn = False ; turn off contour lines
  res_at_txFontHeightF = 0.015

  res_at_vpXF = 0.1 ; make plot bigger
  res_at_vpYF = 0.9
  res_at_vpWidthF = 0.8
  res_at_vpHeightF = 0.8

  res_at_gsnSpreadColors = True ; use full range of color map
  res_at_cnLevels = (/-200,-175,-150,-125,-100,-75,-50,-25,0,50/)
  res_at_lbLabelFontHeightF = 0.015 ; label bar font height
  res_at_lbLabelBarOn = True

  res_at_tiMainFontHeightF = .018 ; font height

  res_at_tiMainString = "Winter" ; add a title
  contour = gsn_csm_contour_map(wks,t1,res) ; create the plot

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 09 2006 - 11:21:42 MST

This archive was generated by hypermail 2.2.0 : Thu Nov 09 2006 - 11:46:51 MST