Re: ncl fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR

Hi Ekin,

A 601x181 grid [108781 pts] is not particularly large.
A commonly encountered 0.5 degree grid is 360x720. Some
of the high resolution ocean model grids are 2400x3600
and NCL can plot.

I am not sure what a "custom grid" is. Please send [jpeg/png/gif]
or point us toward the Matlab generated plot. If the data are at
"nodes", say from a finite element grid then NCL can plot

http://www.ncl.ucar.edu/Applications/
Click: Non-uniform grids/random data
[In particular. Triangular meshes]

Other info
    http://www.ncl.ucar.edu/Document/Graphics/contour_grids.shtml

I don't know ehere the error message was generated.
I speculate memory was 'clobbered'

On 6/8/10 10:53 AM, Ekin Akoglu wrote:
> Dear Dennis,
>
> I know it is a huge data to plot. I have already looked at the data
> before posting and I am aware this way it cannot be plotted. I plotted
> it in MATLAB with a custom grid without problem. However, I could not
> figure it out how to do a similar thing (custom grid via interpolating
> the data) through the documentation on NCL web site. But I suggest that
> NCL may try to be a little bit more "logical" while giving errors so
> that people can understand. If it has been so, I wouldn't have posted
> this as a problem.
>
> Thank you again for your help.
>
> Regards,
>
> Ekin
>
> On 06/08/2010 05:59 PM, Dennis Shea wrote:
>> The most important thing in data processing is to look at your data.
>>
>> %> ncdump -v x_range,y_range,spacing,dimension marmara.nc
>>
>> xysize = 108781
>> [SNIP]
>> x_range = 26, 36 ;
>>
>> y_range = 40, 43 ;
>>
>> spacing = 0.0166666666666667, 0.0166666666666667 ;
>>
>> dimension = 601, 181 ;
>>
>> ===
>> the way you did things you were trying to plot a
>> 108781*108781 array.
>>
>> The array is (601,181) or (181,601) not sure.
>>
>> The attached script works. You may have to change things.
>> Again, please look at your data files.
>>
>> Good luck
>> D
>>
>>
>>
>>
>>>> On 06/07/2010 05:37 AM, Ekin Akoglu wrote:
>>>>> Sorry for cross-posting but I think the full error message will be
>>>>> more
>>>>> meaningful. Here the error I get:
>>>>>
>>>>> fatal:ContourPlotDraw: CTTMTL/ICAEDG - LOGIC ERROR
>>>>> ERROR 5 IN CTTMTL/ICAEDG - LOGIC ERROR
>>>>> fatal:ContourPlotDraw: CTMESH - UNCLEARED PRIOR ERROR
>>>>> fatal:ContourPlotDraw: draw error
>>>>> warning:WorkstationDeactivate: workstation not active or not opened
>>>>>
>>>>> And the code is:
>>>>>
>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>>>
>>>>>> begin
>>>>>>
>>>>>> inf = addfile("marmara.nc","r")
>>>>>>
>>>>>> scale_factor = inf->z@scale_factor
>>>>>> add_offset = inf->z@add_offset
>>>>>> topo = scale_factor * (inf->z) + add_offset
>>>>>>
>>>>>> lons = inf->x_range
>>>>>> lats = inf->y_range
>>>>>> latarray = fspan(lats(0),lats(1),108781)
>>>>>> lonarray = fspan(lons(0),lons(1),108781)
>>>>>>
>>>>>> wks = gsn_open_wks("x11","Black_Sea")
>>>>>> gsn_define_colormap (wks,"gui_default")
>>>>>>
>>>>>> res = True
>>>>>> res@cnFillOn = True ; turn on color
>>>>>> res@cnLinesOn = False ; no contour lines
>>>>>> res@gsnSpreadColors = True ; use full color map
>>>>>> res@lbLabelAutoStride = True ; every other label
>>>>>> res@cnRasterModeOn = True
>>>>>> res@gsnAddCyclic = False ; regional data
>>>>>>
>>>>>> res@cnLevelSelectionMode = "AutomaticLevels" ; automatic levels
>>>>>> ;res@cnMinLevelValF = 100.
>>>>>> ;res@cnMaxLevelValF = 3500.
>>>>>> ;res@cnLevelSpacingF = 100.
>>>>>>
>>>>>>
>>>>>> res@tiMainString = "Black Sea"
>>>>>> res@sfXArray = lonarray ; uses longitude as plot x-axis
>>>>>> res@sfYArray = latarray ; uses latitude as plot y-axis
>>>>>>
>>>>>> res@tfDoNDCOverlay = True ; do not transform data
>>>>>> res@pmTickMarkDisplayMode = "Always" ; get a box around the plot
>>>>>>
>>>>>> plot = gsn_csm_contour(wks,topo,res)
>>>>>>
>>>>>> end
>>>>>
>>>>>
>>>>> Regards,
>>>>> --
>>>>> *Ekin Akoglu*
>>>>>
>>>>> Institute of Marine Sciences
>>>>> Middle East Technical University
>>>>> P.O. Box 28, 33731
>>>>> Erdemli, Mersin
>>>>> Turkey
>>>>>
>>>>> Web: www.ims.metu.edu.tr
>>>>> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
>>>>> Phone: +90 324 521 34 34
>>>>> GSM: +90 506 554 03 90
>>>>>
>>>>> --
>>>>> *Ekin Akoglu*
>>>>>
>>>>> Institute of Marine Sciences
>>>>> Middle East Technical University
>>>>> P.O. Box 28, 33731
>>>>> Erdemli, Mersin
>>>>> Turkey
>>>>>
>>>>> Web: www.ims.metu.edu.tr
>>>>> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
>>>>> Phone: +90 324 521 34 34
>>>>> GSM: +90 506 554 03 90
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>>
>>> --
>>> *Ekin Akoglu*
>>>
>>> Institute of Marine Sciences
>>> Middle East Technical University
>>> P.O. Box 28, 33731
>>> Erdemli, Mersin
>>> Turkey
>>>
>>> Web: www.ims.metu.edu.tr
>>> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
>>> Phone: +90 324 521 34 34
>>> GSM: +90 506 554 03 90
>
>
> --
> *Ekin Akoglu*
>
> Institute of Marine Sciences
> Middle East Technical University
> P.O. Box 28, 33731
> Erdemli, Mersin
> Turkey
>
> Web: www.ims.metu.edu.tr
> Email: ekin@ims.metu.edu.tr <mailto:ekin@ims.metu.edu.tr>
> Phone: +90 324 521 34 34
> GSM: +90 506 554 03 90
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 8 11:57:37 2010