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