<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta http-equiv="content-type"
 content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
Sorry for cross-posting but I think the full error message will be more
meaningful. Here the error I get:<br>
<br>
fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR<br>
&nbsp;ERROR&nbsp;&nbsp;&nbsp; 5 IN CTTMTL/ICAEDG - LOGIC ERROR<br>
fatal:ContourPlotDraw: CTMESH - UNCLEARED PRIOR ERROR<br>
fatal:ContourPlotDraw: draw error<br>
warning:WorkstationDeactivate: workstation not active or not opened<br>
<br>
And the code is:<br>
<br>
<blockquote type="cite">&nbsp;load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
&nbsp;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
  <br>
begin<br>
  <br>
&nbsp;inf = addfile("marmara.nc","r")<br>
  <br>
&nbsp;scale_factor = inf-&gt;z@scale_factor<br>
&nbsp;add_offset = inf-&gt;z@add_offset<br>
&nbsp;topo = scale_factor * (inf-&gt;z) + add_offset<br>
  <br>
&nbsp;lons = inf-&gt;x_range<br>
&nbsp;lats = inf-&gt;y_range<br>
&nbsp;latarray = fspan(lats(0),lats(1),108781)<br>
&nbsp;lonarray = fspan(lons(0),lons(1),108781)<br>
  <br>
&nbsp;wks = gsn_open_wks("x11","Black_Sea")<br>
&nbsp;gsn_define_colormap (wks,"gui_default")<br>
&nbsp;<br>
&nbsp;res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>
&nbsp;res@cnFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp; ; turn on color<br>
&nbsp;res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp; ; no contour lines<br>
&nbsp;res@gsnSpreadColors&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp; ; use full color map<br>
&nbsp;res@lbLabelAutoStride&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp; ; every other label<br>
&nbsp;res@cnRasterModeOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>
&nbsp;res@gsnAddCyclic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp; ; regional data<br>
  <br>
&nbsp;res@cnLevelSelectionMode = "AutomaticLevels" ; automatic levels<br>
&nbsp;;res@cnMinLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100.<br>
&nbsp;;res@cnMaxLevelValF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 3500.<br>
&nbsp;;res@cnLevelSpacingF&nbsp;&nbsp;&nbsp;&nbsp; = 100.<br>
&nbsp;<br>
&nbsp;<br>
&nbsp;res@tiMainString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "Black Sea"<br>
&nbsp;res@sfXArray&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = lonarray&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; uses longitude as
plot x-axis<br>
&nbsp;res@sfYArray&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = latarray&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; uses latitude as plot
y-axis<br>
&nbsp;<br>
&nbsp;res@tfDoNDCOverlay&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; do not transform data<br>
&nbsp;res@pmTickMarkDisplayMode = "Always"&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; get a box around the
plot<br>
&nbsp;<br>
&nbsp;plot = gsn_csm_contour(wks,topo,res)<br>
  <br>
end</blockquote>
<br>
<br>
Regards,<br>
<div class="moz-signature">-- <br>
<b>Ekin Akoglu</b>
<br>
<br>
Institute of Marine Sciences
<br>
Middle East Technical University
<br>
P.O. Box 28, 33731
<br>
Erdemli, Mersin
<br>
Turkey
<br>
<br>
Web: <a href="www.ims.metu.edu.tr">www.ims.metu.edu.tr</a>
<br>
Email: <a href="mailto:ekin@ims.metu.edu.tr">ekin@ims.metu.edu.tr</a>
<br>
Phone: +90 324 521 34 34
<br>
GSM: +90 506 554 03 90</div>
<br>
<div class="moz-signature">-- <br>
<b>Ekin Akoglu</b>
<br>
<br>
Institute of Marine Sciences
<br>
Middle East Technical University
<br>
P.O. Box 28, 33731
<br>
Erdemli, Mersin
<br>
Turkey
<br>
<br>
Web: <a href="www.ims.metu.edu.tr">www.ims.metu.edu.tr</a>
<br>
Email: <a href="mailto:ekin@ims.metu.edu.tr">ekin@ims.metu.edu.tr</a>
<br>
Phone: +90 324 521 34 34
<br>
GSM: +90 506 554 03 90</div>
</body>
</html>